Cited by 1,775Single-cell and spatial transcriptomicsBreast Cancer Treatment StudiesCancer Immunotherapy and Biomarkers

A single-cell and spatially resolved atlas of human breast cancers

Open access
Wu, Sunny Z., Al-Eryani, Ghamdan, Roden, Daniel Lee, Junankar, Simon, Harvey, Kate, Andersson, Alma, Thennavan, Aatish, Wang, Chenfei +32 more·2021·Nature Genetics

Abstract

Breast cancers are complex cellular ecosystems where heterotypic interactions play central roles in disease progression and response to therapy. However, our knowledge of their cellular composition and organization is limited. Here we present a single-cell and spatially resolved transcriptomics analysis of human breast cancers. We developed a single-cell method of intrinsic subtype classification (SCSubtype) to reveal recurrent neoplastic cell heterogeneity. Immunophenotyping using cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) provides high-resolution immune profiles, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcome. Mesenchymal cells displayed diverse functions and cell-surface protein expression through differentiation within three major lineages. Stromal-immune niches were spatially organized in tumors, offering insights into antitumor immune regulation. Using single-cell signatures, we deconvoluted large breast cancer cohorts to stratify them into nine clusters, termed 'ecotypes', with unique cellular compositions and clinical outcomes. This study provides a comprehensive transcriptional atlas of the cellular architecture of breast cancer.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41588-021-00911-1

Cites

11 of 11 references verified0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • The proliferation score, defined based on S-phase genes and G2/M cell-cycle genes, was used to identify proliferative cells, given that senescent cells exhibit G1-arrest.

    “The proliferation score, defined based on S-phase genes and G2/M cell-cycle genes [24], was used to identify proliferative cells, given that senescent cells exhibit G1-arrest.”

    From Ensemble of Time-Evolving SASP Gene Sets Identifies IGFBP7 and CDKN1A as a Potential Marker Pair for Senescent Fibroblast Subpopulations Across Tissues (2026) · cites this paper as [24]

    The Methods section, “Calculating Proliferation and differentiation scores” (page unavailable), confirms that the paper calculated a proliferation score. It defines the score as the average normalized expression of 11 specified genes, not explicitly as a score based on S-phase and G2/M gene sets. The citation is attached to the score definition; the citing authors’ subsequent use of the score and G1-arrest rationale are not attributed to this source. (LLM, conf=0.85)

    Experimental — accuracy not yet measured

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Research works printed but not linked: 59

  1. [1]Kim HK Discordance of the PAM50 Intrinsic Subtypes Compared with Immunohistochemistry-Based Surrogate in Breast Cancer Patients: Potential Implication of Genomic Alterations of Discordance. Cancer Res Treat 51, 737– 747 (2019). 30189722 PMC6473265
  2. [2]Picornell AC Breast cancer PAM50 signature: correlation and concordance between RNA-Seq and digital multiplexed gene expression technologies in a triple negative breast cancer series. BMC Genomics 20, 452 (2019). 31159741 PMC6547580
  3. [3]Parker JS Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol 27, 1160– 7 (2009). 19204204 PMC2667820
  4. [4]Perou CM Molecular portraits of human breast tumours. Nature 406, 747– 52 (2000). 10963602
  5. [5]Sorlie T Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A 98, 10869– 74 (2001). 11553815 PMC58566
  6. [6]Sorlie T Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A 100, 8418– 23 (2003). 12829800 PMC166244
  7. [7]Su S CD10(+)GPR77(+) Cancer-Associated Fibroblasts Promote Cancer Formation and Chemoresistance by Sustaining Cancer Stemness. Cell 172, 841– 856 e16 (2018). 29395328
  8. [8]Cazet AS Targeting stromal remodeling and cancer stem cell plasticity overcomes chemoresistance in triple negative breast cancer. Nature Communications 9, 2897 (2018). PMC6057940 30042390
  9. [9]Dushyanthen S Relevance of tumor-infiltrating lymphocytes in breast cancer. BMC Med 13, 202 (2015). 26300242 PMC4547422
  10. [10]Cassetta L Human Tumor-Associated Macrophage and Monocyte Transcriptional Landscapes Reveal Cancer-Specific Reprogramming, Biomarkers, and Therapeutic Targets. Cancer Cell 35, 588– 602 e10 (2019). 30930117 PMC6472943
  11. [11]Katzenelenbogen Y Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer. Cell 182, 872– 885 e19 (2020). 32783915
  12. [12]Medler TR Complement C5a Fosters Squamous Carcinogenesis and Limits T Cell Response to Chemotherapy. Cancer Cell 34, 561– 578 e6 (2018). 30300579 PMC6246036
  13. [13]Nakamura K & Smyth MJ TREM2 marks tumor-associated macrophages. Signal Transduct Target Ther 5, 233 (2020). 33037186 PMC7547002
  14. [15]Wu SZ Stromal cell diversity associated with immune evasion in human triple-negative breast cancer. EMBO J, e104063 (2020). 32790115 PMC7527929
  15. [18]Biffi G IL1-Induced JAK/STAT Signaling Is Antagonized by TGFbeta to Shape CAF Heterogeneity in Pancreatic Ductal Adenocarcinoma. Cancer Discov (2018). PMC6368881 30366930
  16. [19]Elyada E Cross-Species Single-Cell Analysis of Pancreatic Ductal Adenocarcinoma Reveals Antigen-Presenting Cancer-Associated Fibroblasts. Cancer Discov 9, 1102– 1123 (2019). 31197017 PMC6727976
  17. [20]Puram SV Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer. Cell 171, 1611– 1624 e24 (2017). 29198524 PMC5878932
  18. [21]Lambrechts D Phenotype molding of stromal cells in the lung tumor microenvironment. Nat Med 24, 1277– 1289 (2018). 29988129
  19. [22]Azizi E Single-Cell Map of Diverse Immune Phenotypes in the Breast Tumor Microenvironment. Cell 174, 1293– 1308 e36 (2018). 29961579 PMC6348010
  20. [23]Savas P Single-cell profiling of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis. Nat Med 24, 986– 993 (2018). 29942092
  21. [24]Kim C Chemoresistance Evolution in Triple-Negative Breast Cancer Delineated by Single-Cell Sequencing. Cell 173, 879– 893 e13 (2018). 29681456 PMC6132060
  22. [25]Ali HR Imaging mass cytometry and multiplatform genomics define the phenogenomic landscape of breast cancer. Nature Cancer 1, 163– 175 (2020). 35122013
  23. [26]Wagner J A Single-Cell Atlas of the Tumor and Immune Ecosystem of Human Breast Cancer. Cell 177, 1330– 1345 e18 (2019). 30982598 PMC6526772
  24. [29]Pliner HA, Shendure J & Trapnell C Supervised classification enables rapid annotation of cell atlases. Nature Methods 16, 983– 986 (2019). 31501545 PMC6791524
  25. [30]Neftel C An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma. Cell 178, 835– 849 e21 (2019). 31327527 PMC6703186
  26. [31]Tirosh I Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189– 96 (2016). 27124452 PMC4944528
  27. [32]Cancer Genome Atlas N Comprehensive molecular portraits of human breast tumours. Nature 490, 61– 70 (2012). 23000897 PMC3465532
  28. [33]Prat A Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res 12, R68 (2010). 20813035 PMC3096954
  29. [34]Nielsen TO A comparison of PAM50 intrinsic subtyping with immunohistochemistry and clinical prognostic factors in tamoxifen-treated estrogen receptor-positive breast cancer. Clin Cancer Res 16, 5222– 32 (2010). 20837693 PMC2970720
  30. [35]Stoeckius M Simultaneous epitope and transcriptome measurement in single cells. Nat Methods 14, 865– 868 (2017). 28759029 PMC5669064
  31. [36]Stuart T Comprehensive Integration of Single-Cell Data. Cell 177, 1888– 1902 e21 (2019). 31178118 PMC6687398
  32. [37]Glajcar A, Szpor J, Hodorowicz-Zaniewska D, Tyrak KE & Okon K The composition of T cell infiltrates varies in primary invasive breast cancer of different molecular subtypes as well as according to tumor size and nodal status. Virchows Arch 475, 13– 23 (2019). 31016433 PMC6611750
  33. [38]Li H Dysfunctional CD8 T Cells Form a Proliferative, Dynamically Regulated Compartment within Human Melanoma. Cell 176, 775– 789 e18 (2019). 30595452 PMC7253294
  34. [39]Yamada S, Shinozaki K & Agematsu K Involvement of CD27/CD70 interactions in antigen-specific cytotoxic T-lymphocyte (CTL) activity by perforin-mediated cytotoxicity. Clin Exp Immunol 130, 424– 30 (2002). 12452832 PMC1906551
  35. [41]Slyper M A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors. Nature Medicine 26, 792– 802 (2020). PMC7220853 32405060
  36. [42]Ruffell B Leukocyte composition of human breast cancer. Proc Natl Acad Sci U S A 109, 2796– 801 (2012). 21825174 PMC3287000
  37. [45]Chen J CCL18 from tumor-associated macrophages promotes breast cancer metastasis via PITPNM3. Cancer Cell 19, 541– 55 (2011). 21481794 PMC3107500
  38. [46]Qiu X Reversed graph embedding resolves complex single-cell trajectories. Nat Methods 14, 979– 982 (2017). 28825705 PMC5764547
  39. [47]Kumar A Specification and Diversification of Pericytes and Smooth Muscle Cells from Mesenchymoangioblasts. Cell Rep 19, 1902– 1916 (2017). 28564607 PMC6428685
  40. [48]Thiriot A Differential DARC/ACKR1 expression distinguishes venular from non-venular endothelial cells in murine tissues. BMC Biology 15, 45 (2017). 28526034 PMC5438556
  41. [49]Mailhos C Delta4, an endothelial specific notch ligand expressed at sites of physiological and tumor angiogenesis. Differentiation 69, 135– 44 (2001). 11798067
  42. [50]Ubezio B Synchronization of endothelial Dll4-Notch dynamics switch blood vessels from branching to expansion. Elife 5(2016). PMC4894757 27074663
  43. [51]Kryczek I CXCL12 and vascular endothelial growth factor synergistically induce neoangiogenesis in human ovarian cancers. Cancer Res 65, 465– 72 (2005). 15695388
  44. [52]Blanco R & Gerhardt H VEGF and Notch in tip and stalk cell selection. Cold Spring Harb Perspect Med 3, a006569 (2013). 23085847 PMC3530037
  45. [53]Jakobsson L Endothelial cells dynamically compete for the tip cell position during angiogenic sprouting. Nat Cell Biol 12, 943– 53 (2010). 20871601
  46. [54]Andersson A Single-cell and spatial transcriptomics enables probabilistic inference of cell type topography. Communications Biology 3, 565 (2020). 33037292 PMC7547664
  47. [55]Lakins MA, Ghorani E, Munir H, Martins CP & Shields JD Cancer-associated fibroblasts induce antigen-specific deletion of CD8 (+) T Cells to protect tumour cells. Nat Commun 9, 948 (2018). 29507342 PMC5838096
  48. [56]Newman AM Determining cell type abundance and expression from bulk tissues with digital cytometry. Nature Biotechnology 37, 773– 782 (2019). PMC6610714 31061481
  49. [57]Tsoucas D Accurate estimation of cell-type composition from gene expression data. Nat Commun 10, 2975 (2019). 31278265 PMC6611906
  50. [58]Tsai JH, Donaher JL, Murphy DA, Chau S & Yang J Spatiotemporal regulation of epithelial-mesenchymal transition is essential for squamous cell carcinoma metastasis. Cancer Cell 22, 725– 36 (2012). 23201165 PMC3522773
  51. [59]Katzenelenbogen Y Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer. Cell (2020). 32783915
  52. [61]Wu SZ Cryopreservation of human cancers conserves tumour heterogeneity for single-cell multi-omics analysis. Genome Medicine 13, 81 (2021). 33971952 PMC8111910
  53. [62]Lun ATL EmptyDrops: distinguishing cells from empty droplets in droplet-based single-cell RNA sequencing data. Genome Biol 20, 63 (2019). 30902100 PMC6431044
  54. [63]Zhao X, Rodland EA, Tibshirani R & Plevritis S Molecular subtyping for clinically defined breast cancer subgroups. Breast Cancer Res 17, 29 (2015). 25849221 PMC4365540
  55. [64]Patro R, Duggal G, Love MI, Irizarry RA & Kingsford C Salmon provides fast and bias-aware quantification of transcript expression. Nature Methods 14, 417– 419 (2017). 28263959 PMC5600148
  56. [65]Finak G MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biol 16, 278 (2015). 26653891 PMC4676162
  57. [67]Yu G, Wang LG, Han Y & He QY clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284– 7 (2012). 22455463 PMC3339379
  58. [68]Andersson A Spatial Deconvolution of HER2-positive Breast Tumors Reveals Novel Intercellular Relationships. bioRxiv, 2020.07.14.200600 (2020).
  59. [69]Efremova M, Vento-Tormo M, Teichmann SA & Vento-Tormo R CellPhoneDB: inferring cell–cell communication from combined expression of multi-subunit ligand–receptor complexes. Nature Protocols 15, 1484– 1506 (2020). 32103204

Cited by

Sorted by citations

Authors
Wu, S. Z., Al-Eryani, G., Roden, D. L., Junankar, S., Harvey, K., Andersson, A., Thennavan, A., Wang, C., Torpy, J. R., Bartonicek, N., Wang, T., Larsson, L., Kaczorowski, D., Weisenfeld, N. I., Uytingco, C. R., Chew, J. G., Bent, Z. W., Chan, C.-L., Gnanasambandapillai, V., Dutertre, C.-A., Gluch, L., Hui, M. N., Beith, J., Parker, A., Robbins, E., Segara, D., Cooper, C., Mak, C., Chan, B., Warrier, S., Ginhoux, F., Millar, E., Powell, J. E., Williams, S. R., Liu, X. S., O’Toole, S., Lim, E., Lundeberg, J., Perou, C. M., & Swarbrick, A.
Affiliations
  • Garvan Institute of Medical Research (AU)
  • Science for Life Laboratory (SE)
  • University of North Carolina at Chapel Hill (US)
  • Harvard University (US)
  • 10X Genomics (United States) (US)
  • Inserm (FR)
  • Chris O’Brien Lifehouse (AU)
  • UNSW Sydney (AU)
  • Royal Prince Alfred Hospital (AU)
  • St Vincent's Hospital Sydney (AU)
  • The University of Queensland (AU)
  • The University of Sydney (AU)
  • Agency for Science, Technology and Research (SG)
Corresponding author
Swarbrick, Alexander
Venue
Nature Genetics
Published
2021-09
Volume
53
Issue
9
Pages
1334-1347
ISSN
  • 1061-4036 (ISSN-L)
  • 1546-1718
OpenAlex
W3198511950
Cited by
1,775
References
11
Citations by year
2026: 459 · 2025: 430 · 2024: 424 · 2023: 316 · 2022: 140 · 2021: 6
FWCI
72.13
Type
journal article
Language
en
Open access
Locations
  • Nature Genetics (publishedVersion)
  • PubMed (publishedVersion)
  • PubMed Central (acceptedVersion)
  • Carolina Digital Repository (University of North Carolina at Chapel Hill) (submittedVersion)
  • UNC Libraries
Funders
  • National Breast Cancer Foundation (IIRS-19-106, IIRS-19-084, PRAC 16-006)
  • Petre Foundation
  • Department of Health | National Health and Medical Research Council
  • The White Butterfly Foundation The Sydney Breast Cancer Foundation John McMurtrie (AM) and Deborah McMurtrie
  • Breast Cancer Research Foundation
  • The Sydney Breast Cancer Foundation The Family and Friends of Michael O'Sullivan
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cancer Research › Breast Cancer Treatment Studies
  • Health Sciences › Medicine › Oncology › Cancer Immunotherapy and Biomarkers
Keywords
breast cancer, tumor heterogeneity, tumor microenvironment, macrophages, mesenchymal cells, CITE-seq, spatial transcriptomics, single-cell RNA sequencing, prognosis
MeSH terms
B-Lymphocytes, Breast Neoplasms, Female, Humans, Macrophages, Membrane Proteins, Biomarkers, Tumor, Gene Expression Regulation, Neoplastic, Sequence Analysis, RNA, CD8-Positive T-Lymphocytes, Myeloid Cells, Tumor Suppressor Proteins, Endothelial Cells, Single-Cell Analysis, Tumor Microenvironment, Transcriptome, B7-H1 Antigen

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review