SCODiA Report by MLBI Lab

Single-Cell Dissection of the Breast Cancer Microenvironment Reveals Distinct Genomic, Cellular, and Intercellular Communication Landscapes

This single-cell RNA-sequencing analysis comprehensively characterizes breast tissue, contrasting normal with primary tumor conditions. We identify widespread aneuploidy in tumor epithelial cells, profound remodeling of the tumor microenvironment with activated fibroblasts and altered immune cell populations, and distinct cell-cell interaction networks. These findings reveal critical genomic instability, metabolic reprogramming, and intercellular signaling hubs that drive breast cancer progression.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell RNA-seq Data by Condition, Sample, Cell Type, and Ploidy in Breast Tissue
  3. UMAP Visualization of Key Cell Type Markers and Minor Cell Type Annotations
  4. Overall Celltype_subset Marker Expression Profile
  5. Copy Number Variation Analysis in Tumor-Origin and Unassigned Cells Grouped by Sample
  6. CNV-based UMAP Visualization of Single-Cell RNA-seq Data
  7. Minor Cell Type Population Analysis in Breast Tissue: Normal vs. Primary Tumor
  8. T Cell and Innate Lymphoid Cell Subpopulation Analysis in Breast Cancer
  9. Macrophage Subset Population Analysis in Breast Tissue: Normal vs. Primary Tumor
  10. Analysis of T Cell Subset Proportions in Breast Tissue
  11. Ploidy Population Analysis of Epithelial and Unassigned Cells in Breast Tissue
  12. Cell-Cell Interaction Patterns in Primary Breast Tumors
  13. Primary Breast Tumor Cell-Cell Interaction Analysis
  14. Cell-Cell Interaction Analysis of Immune Checkpoint and Cell Cycle Related Genes in Breast Tissue
  15. Condition-Specific Cell-Cell Interaction Patterns in Breast Cancer
  16. Epithelial Cell Condition-Specific Surfaceome Markers in Breast Tissue
  17. Macrophage Condition-Specific Surfaceome Markers in Breast Tissue
  18. Fibroblast Condition-Specific Surfaceome Markers in Breast Tissue
  19. Differential Expression of Cell Cycle-Related Gene YWHAZ in Breast Epithelial Cells
  20. Gene Ontology (GSA) Analysis of Epithelial Cells in Breast Tissue
  21. Gene Set Enrichment Analysis (GSEA) of Breast Tissue Cell Types in Primary Tumor vs. Normal Conditions
  22. Discussion
  23. Query List

0. Dataset overview

Dataset Summary

Total Cells: 88,707 cells

Total Genes: 25,535 genes

Species: Human

Tissue: Breast

Conditions: primary_tumor, normal

1. UMAP Visualization of Single-Cell RNA-seq Data by Condition, Sample, Cell Type, and Ploidy in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis provides a comprehensive visualization of single-cell RNA-sequencing data from breast tissue on a Uniform Manifold Approximation and Projection (UMAP) embedding. The UMAP plots are colored by various cellular and sample-specific attributes: disease condition (primary_tumor vs. normal), individual sample origin, major cell type, minor cell type, inferred ploidy status, and detailed cell type subsets. These visualizations are crucial for assessing data quality, cell type annotation fidelity, potential batch effects, and identifying biologically meaningful patterns related to disease state and cellular characteristics.

Visual Summary

Condition and Sample Distribution

Cell Type Hierarchy

Ploidy Status

Biological Interpretation

The UMAP visualizations reveal several key biological insights into the breast tissue scRNA-seq dataset:

  1. Disease-Associated Cellular Changes: The clear distinction between normal and primary_tumor cells on the UMAP highlights significant shifts in cellular composition, transcriptional states, or both, within the tumor microenvironment compared to healthy tissue. The existence of tumor-specific clusters strongly suggests the presence of unique cell populations or altered states critical to tumor biology.
  2. Robust Cell Type Identification: The consistent clustering of cells by major, minor, and subset annotations demonstrates the high quality and specificity of the cell type assignments. The hierarchical nature of the clustering, where finer subsets reside within broader categories, supports a biologically meaningful organization of cellular identities.
  3. Aneuploidy as a Tumor Hallmark: The pronounced co-localization of Aneuploid cells with the Epithelial cell clusters is a highly significant finding. Given that "Epithelial cell" is specified as the Tumor origin celltype, this observation strongly suggests that the malignant epithelial cells within the primary tumor exhibit aneuploidy. Aneuploidy, the presence of an abnormal number of chromosomes, is a well-established hallmark of cancer cells and is associated with tumor development and progression. This finding validates the genomic instability commonly observed in breast cancer and confirms the successful identification of the neoplastic compartment within the dataset PubMed: aneuploidy cancer review.
  4. Tumor Microenvironment Complexity: The detailed celltype_subset UMAP showcases the intricate cellular heterogeneity within the breast tissue, particularly within the tumor microenvironment. The identification of various immune cell subtypes (e.g., T cell (Treg), Macrophage (M1/M2 subtypes)), stromal cells (e.g., Fibroblast), and different epithelial populations provides a rich foundation for investigating cell-cell interactions and their roles in tumor progression or immune response.

Annotation Notes

The UMAP plots serve as an excellent quality control and validation step for the dataset.

2. UMAP Visualization of Key Cell Type Markers and Minor Cell Type Annotations

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the expression levels of a panel of known cell type marker genes across the entire single-cell RNA-seq dataset, projected onto a UMAP embedding. Alongside these gene expression plots, a UMAP plot colored by celltype_minor annotations is provided as a reference. The purpose is to assess the consistency and quality of the predefined celltype_minor annotations against the expression patterns of canonical marker genes.

Visual Summary

The UMAP plots clearly delineate distinct clusters of cells. Each gene expression plot highlights specific regions of the UMAP with elevated expression (represented by warmer colors), corresponding to particular cell populations. The celltype_minor UMAP shows well-separated clusters for each annotated cell type, such as B cells, T cells (CD4+ and CD8+), Macrophages, Epithelial cells, Fibroblasts, Endothelial cells, etc.

Key observations from the visual comparison include:

Biological Interpretation

The expression patterns of the selected genes align remarkably well with their known biological roles as specific cell type markers, providing strong support for the accuracy of the celltype_minor annotations.

T Cell Markers:

B Cell & Plasma Cell Markers:

Myeloid Cell Markers:

Stromal Cell Markers:

Epithelial Cell Markers:

Endothelial Cell Marker:

Annotation Notes

The high concordance between the gene expression patterns and the celltype_minor annotations suggests that the cell type assignments are robust and accurate for the major cell populations represented in this breast tissue single-cell RNA-seq dataset. The clear separation of cell types on the UMAP and the specific localization of marker gene expression reinforce confidence in the underlying clustering and annotation process. This foundational check is critical for subsequent downstream analyses, such as differential gene expression or cell-cell interaction studies, ensuring that these analyses are performed on correctly identified cell populations.

3. Overall Celltype_subset Marker Expression Profile

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the expression of marker genes across different celltype_subset populations identified in the single-cell RNA-seq dataset from human breast tissue. The dot plot displays the mean expression level (color intensity) and the fraction of cells expressing a gene (dot size) for selected marker genes within each cell type. The primary goal is to assess the quality of the celltype_subset annotations by examining whether the expression patterns of these markers align with known cell type identities.

Visual Summary

The dot plot demonstrates a highly organized and distinct pattern of gene expression, with most marker genes showing strong specificity for particular celltype_subset groups.

Biological Interpretation

The observed marker gene expression patterns are largely consistent with established biological knowledge of cell types found in human breast tissue, affirming the quality of the celltype_subset annotations.

Annotation Notes

The comprehensive and specific marker expression patterns observed across the celltype_subset populations strongly validate the current cell type annotations within the AnnData object. The ability to identify both surface proteins and key intracellular transcription factors further enhances confidence in the fine-grained resolution of these cell identities. This robust annotation forms a solid foundation for subsequent in-depth analyses, such as differential gene expression or cell-cell interaction studies, ensuring that these analyses are performed on accurately defined cell populations.

References

  1. POU2F2 (OCT2): GeneCards (GeneCards)
  2. EBF1: GeneCards (GeneCards)
  3. ACKR1 (DARC): GeneCards (GeneCards)
  4. PROX1: GeneCards (GeneCards)
  5. FAP: GeneCards (GeneCards)
  6. ILC transcription factors: PubMed Search: "ILC differentiation transcription factors" (PubMed Search)
  7. CD68: GeneCards (GeneCards)
  8. SPP1 (Osteopontin): GeneCards (GeneCards)
  9. Keratins (KRT8, KRT18): GeneCards (GeneCards), (GeneCards)
  10. PIP: GeneCards (GeneCards)
  11. TPSAB1 (Tryptase): GeneCards (GeneCards)
  12. KIT: GeneCards (GeneCards)
  13. XBP1: GeneCards (GeneCards)
  14. ACTA2: GeneCards (GeneCards)
  15. CD8A: GeneCards (GeneCards)
  16. FOXP3: GeneCards (GeneCards)
  17. CTLA4: GeneCards (GeneCards)
  18. T cell subsets: PubMed Search: "T cell subset markers" (PubMed Search)

4. Copy Number Variation Analysis in Tumor-Origin and Unassigned Cells Grouped by Sample

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates copy number variations (CNVs) in Epithelial cells (identified as the tumor-origin cell type) and 'unassigned' cells, grouped by individual sample. The goal is to visualize CNV patterns across the genome for each sample and provide a summary of significantly amplified or deleted regions. This helps in understanding the genomic instability characteristic of cancer cells and validating the ploidy status of the samples.

Visual Summary

The visualization consists of two main parts: a CNV heatmap and a summary of significant CNV regions.

CNV Heatmap

Summary of Significantly Amplified Copy Number Regions

1q21.3:1q23.2 (Frequency: 0.64)

1q23.3:1q24.1 (Frequency: 0.45)

11q13.4:11q21 (Frequency: 0.27)

Biological Interpretation

  1. Confirmation of Tumor Status: The stark contrast in CNV burden between "Diploid Patient" and "Patient" samples strongly supports the accurate classification of these samples based on their ploidy status. The "Patient" samples, exhibiting widespread CNVs in Epithelial cells (the designated tumor-origin cell type) and 'unassigned' cells, are highly indicative of malignant tumor tissue. The 'unassigned' cells showing similar CNV patterns might represent tumor cells that could not be assigned to a specific epithelial subtype or are part of the tumor microenvironment with genomic alterations.
  2. Genomic Landscape of Breast Cancer: The identified recurrent CNVs, such as amplifications on 1q, 8q, 11q, and 20q, are well-established genomic alterations in breast cancer PMID: 29038234. These regions often harbor oncogenes whose increased copy number contributes to tumor initiation and progression.
  1. Potential Driver Genes: The specific genes highlighted (*NFASC*, *EIF3E*) within recurrently amplified regions warrant further investigation. *NFASC* (Neurofascin) is involved in cell adhesion and neural development, but its role in breast cancer, particularly upon amplification, could suggest altered cell-cell or cell-extracellular matrix interactions that contribute to tumor invasion and metastasis. *EIF3E* (Eukaryotic Translation Initiation Factor 3 Subunit E) plays a role in initiating protein synthesis, and its overexpression due to amplification could lead to increased translation of oncogenic proteins, promoting cell growth and survival.

Clinical or Translational Implications

5. CNV-based UMAP Visualization of Single-Cell RNA-seq Data

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a Uniform Manifold Approximation and Projection (UMAP) visualization where the embedding (X_cnv_umap1, X_cnv_umap2) is specifically computed using copy number variation (CNV) estimates from single-cell RNA-seq data. The resulting UMAP plots illustrate how cells group based on their CNV profiles, overlaid with annotations for major cell types, minor cell types, ploidy status (Aneuploid/Diploid), biological condition (normal/primary_tumor), and individual sample origin. This allows for an assessment of the relationship between genomic alterations (CNVs/ploidy) and cellular identity, disease state, and inter-sample variability.

Visual Summary

The five UMAP plots collectively reveal distinct patterns driven by cellular CNV profiles:

Biological Interpretation

The UMAP visualization, based on CNV estimates, provides critical insights into the genomic landscape of the breast tissue samples:

Annotation Notes

6. Minor Cell Type Population Analysis in Breast Tissue: Normal vs. Primary Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis provides a comprehensive view of the cellular composition, specifically focusing on minor cell type proportions, across different samples of normal breast tissue and primary breast tumors. The stacked bar plots normalize the cell counts within each sample to 100%, allowing for direct comparison of relative cell type abundances between conditions and across individual patients. This is crucial for understanding the overall cellular ecosystem in health and disease.

Visual Summary

The stacked bar plots display the relative proportions of various minor cell types for each sample, categorized by 'normal' and 'primary_tumor' conditions.

Biological Interpretation

The observed shifts in minor cell type populations between normal breast tissue and primary tumors highlight significant alterations in the cellular landscape during breast cancer progression.

  1. Tumor-associated Stromal and Immune Remodeling: The consistent dominance of Epithelial cells and Fibroblasts in tumors, along with an increase in Macrophages, reflects the well-established features of the tumor microenvironment (TME). Tumor-associated fibroblasts (CAFs) and tumor-associated macrophages (TAMs) are critical components of the TME, promoting tumor growth, invasion, angiogenesis, and immunosuppression [1]. The data shows an upregulation of macrophage presence in tumor samples, suggesting an active role for these immune cells in the disease context.
  2. Heterogeneity of Anti-tumor Immune Response: The variable presence of T cells (CD4+ and CD8+) in primary tumor samples points to inter-patient heterogeneity in immune infiltration. Some tumors appear "hot" (more infiltrated by T cells), while others are "cold" (less infiltrated). This heterogeneity is a known factor influencing patient response to immunotherapy in breast cancer [2].
  3. Significance of 'Unassigned' Cells: The high proportion of 'unassigned' cells in certain primary tumor samples is noteworthy. These cells could represent:

Clinical or Translational Implications

Understanding the precise cellular composition of breast tumors has several clinical and translational implications:

---

References:

[1] Macrophage Biology and Breast Cancer: NCBI (PubMed Central for a review on TAMs)

[2] Tumor-infiltrating lymphocytes in breast cancer: NCBI (PubMed Central for a review on TILs in breast cancer)

[3] Prognostic value of tumor-infiltrating lymphocytes in breast cancer: PubMed Search (PubMed search for "prognostic value tumor-infiltrating lymphocytes breast cancer")

[4] Tumor-associated macrophages in breast cancer prognosis: PubMed Search (PubMed search for "tumor-associated macrophages breast cancer prognosis")

7. T Cell and Innate Lymphoid Cell Subpopulation Analysis in Breast Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis utilizes single-cell RNA sequencing data to characterize the cellular composition of the 'T cell' major population, further broken down into various T cell subsets and innate lymphoid cells (ILCs), across normal breast tissue and primary breast tumor samples. The plot_celltype_population tool was used to visualize the relative proportions of these immune cell subsets within each individual sample, categorized by 'condition' (normal vs. primary_tumor). This provides insights into the immune landscape shifts associated with breast tumorigenesis.

Visual Summary

The stacked bar plot effectively illustrates the relative proportions of T cell and ILC subsets for each sample, grouped by condition.

Normal vs. Primary Tumor Differences:

Biological Interpretation

The observed shifts in T cell and ILC subset populations between normal breast tissue and primary tumors provide key biological insights into the immune microenvironment of breast cancer.

Clinical or Translational Implications

8. Macrophage Subset Population Analysis in Breast Tissue: Normal vs. Primary Tumor

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis utilizes single-cell RNA sequencing data to visualize the proportional distribution of different macrophage subsets (M1, M2A, M2B, M2C, M2D) within individual samples, comparing normal breast tissue to primary breast tumor samples. The plot_celltype_population tool was used to specifically focus on the 'Macrophage' celltype_minor and further stratify them by their 'subset' classifications.

Visual Summary

The stacked bar plots display the relative proportions of five distinct macrophage subsets (M1, M2A, M2B, M2C, M2D) across individual samples grouped by condition ('normal' and 'primary_tumor').

Biological Interpretation

Macrophages are critical components of the tumor microenvironment (TME) and are highly plastic, polarizing into various phenotypes, broadly categorized as M1 (pro-inflammatory, anti-tumorigenic) and M2 (anti-inflammatory, pro-tumorigenic). The observed shifts in macrophage subsets provide insights into the immune landscape of breast cancer.

Clinical or Translational Implications

Understanding the balance and distribution of macrophage subsets has significant implications for breast cancer diagnosis, prognosis, and therapeutic strategies.

References:

  1. M2 Macrophage Roles: PubMed search for "M2 macrophages tissue repair immune regulation" https://pubmed.ncbi.nlm.nih.gov/?term=M2+macrophages+tissue+repair+immune+regulation
  2. M1 Macrophage Functions: PubMed search for "M1 macrophages anti-tumor immunity" https://pubmed.ncbi.nlm.nih.gov/?term=M1+macrophages+anti-tumor+immunity
  3. M2 and TAMs in Cancer: PubMed search for "M2 macrophages tumor associated macrophages cancer" https://pubmed.ncbi.nlm.nih.gov/?term=M2+macrophages+tumor+associated+macrophages+cancer
  4. M1/M2 Ratio Prognosis: PubMed search for "M1 M2 ratio prognosis cancer" https://pubmed.ncbi.nlm.nih.gov/?term=M1+M2+ratio+prognosis+cancer
  5. Macrophage Repolarization Therapy: PubMed search for "macrophage repolarization cancer therapy" https://pubmed.ncbi.nlm.nih.gov/?term=macrophage+repolarization+cancer+therapy

9. Analysis of T Cell Subset Proportions in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportions of specific T cell subset populations (T follicular helper cells (Tfh), Regulatory T cells (Treg), and Lymphoid Tissue Inducer cells (LTI)) within single-cell RNA-seq data from human breast tissue. The aim is to identify statistically significant differences in these cell type proportions between 'primary_tumor' and 'normal' conditions, providing insights into the immune microenvironment shifts associated with breast cancer.

Visual Summary

The boxplots illustrate the celltype proportion (as a percentage of total cells) for three T cell subsets across primary tumor and normal breast tissue samples. Black dots represent individual sample measurements, while the boxes indicate the interquartile range (IQR), with the line inside representing the median. Statistical significance between conditions is indicated by p-values.

Biological Interpretation

The observed shifts in T cell subset proportions between primary breast tumors and normal tissue highlight significant alterations in the immune landscape during oncogenesis.

Clinical or Translational Implications

The differential proportions of T cell subsets between primary tumors and normal breast tissue carry important clinical and translational implications:

10. Ploidy Population Analysis of Epithelial and Unassigned Cells in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the ploidy status (Aneuploid, Diploid, or Unclear) within the population of 'Epithelial cell' and 'unassigned' cells across individual samples from both normal and primary tumor breast tissue conditions. Given that 'Epithelial cell' is identified as the tumor origin cell type, this analysis provides insights into the genomic stability of these crucial cell populations in the context of breast cancer development.

Visual Summary

The stacked bar plot effectively illustrates the ploidy distribution within the combined "Epithelial cell and unassigned" population for each sample, separated by condition.

Biological Interpretation

The observed ploidy patterns strongly correlate with the pathological state of the tissue, providing critical biological insights:

Clinical or Translational Implications

The findings from this ploidy analysis have several important clinical and translational implications for breast cancer:

11. Cell-Cell Interaction Patterns in Primary Breast Tumors

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns within the primary breast tumor microenvironment. Using single-cell RNA-seq data, CellPhoneDB was applied to identify significant ligand-receptor interactions between various cell types, including the tumor-origin Epithelial cells (further delineated by ploidy status into Diploid and Aneuploid), Fibroblasts, Macrophages, and T cells (CD4+ and CD8+). The visualization highlights the top 80 significant interactions based on p-value and interaction strength (mean expression) specifically within the primary tumor condition, providing insights into the complex communication network driving tumor biology.

Visual Summary

The dot plot displays a complex landscape of cell-cell interactions within primary breast tumors.

Prominent Interacting Cell Types:

Biological Interpretation

The observed CCI patterns reveal critical communication axes within the breast tumor microenvironment that likely contribute to tumor progression, immune evasion, and stromal remodeling.

Immune Cell Recruitment and Modulation:

Other Key Pathways:

Clinical or Translational Implications

The identified cell-cell interaction patterns offer several avenues for clinical and translational applications in breast cancer.

Therapeutic Target Prioritization:

12. Primary Breast Tumor Cell-Cell Interaction Analysis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCI) within the primary breast tumor microenvironment using single-cell RNA sequencing data. The dot plot visualizes the most significant and highly expressed ligand-receptor pairs between different cell populations present in the primary_tumor condition. Interactions are filtered for a p-value cutoff of 0.05 and a mean expression cutoff of 0.01, displaying up to 80 prominent interactions. The size of the dots represents the statistical significance (as -log10(p-value)), and the color indicates the interaction strength (as log2(mean expression)). The analysis specifically includes interactions involving 'Aneuploid Epithelial' and 'Diploid Epithelial' cells, given that 'Epithelial cell' is identified as the tumor origin cell type and ploidy_dec distinguishes these populations.

Visual Summary

The visualization reveals a highly active and complex network of cell-cell communications within the primary breast tumor.

Key Ligand-Receptor Families:

Biological Interpretation

The observed cell-cell interactions provide critical insights into the biological processes driving primary breast tumor progression and shaping the tumor microenvironment (TME).

  1. ECM Remodeling and Tumor Cell-Stroma Adhesion (Integrins): The widespread integrin-mediated interactions underscore the profound importance of extracellular matrix (ECM) remodeling and cell adhesion in breast cancer. Fibroblasts, acting as cancer-associated fibroblasts (CAFs), secrete collagen and fibronectin, which then interact with integrins on both tumor epithelial cells and other stromal cells (endothelial cells). This interaction network facilitates tumor cell migration, invasion, survival, and resistance to therapy, and is fundamental to the structural and functional integrity of the TME. GeneCards: Integrin alpha family, GeneCards: Collagen family
  2. Angiogenesis (PGF-FLT1, VEGF-VEGFR2): The strong PGF-FLT1 (and PGF-NRP1) and VEGF-VEGFR2 interactions, particularly between Fibroblasts/Endothelial cells and Epithelial cells, are hallmarks of active angiogenesis. Placental growth factor (PGF) and Vascular Endothelial Growth Factor (VEGF) are potent pro-angiogenic factors that promote the formation of new blood vessels, essential for providing nutrients and oxygen to the rapidly growing tumor, and for facilitating metastasis. This robust signaling indicates a highly pro-angiogenic TME in primary breast tumors. PubMed Search: PGF angiogenesis cancer, PubMed Search: VEGF angiogenesis cancer
  3. Immune Cell Recruitment and Stromal Interaction (Chemokines):
  1. Pro-inflammatory and Pro-tumorigenic Signaling (IL6-IL6R): The presence of IL6-IL6 receptor interactions indicates an active inflammatory component within the primary tumor microenvironment. IL-6 is a pleiotropic cytokine that can drive tumor cell proliferation, survival, angiogenesis, and immune evasion in breast cancer. GeneCards: IL6
  2. Aneuploid Epithelial Cells as Drivers: The active participation of Aneuploid Epithelial cells in a broad range of strong interactions, particularly with Fibroblasts and Endothelial cells, supports their role as the primary malignant population actively co-opting the microenvironment for growth and survival. The interactions with stromal cells (e.g., integrins, PGF-FLT1) are critical for their malignant phenotype.

Clinical or Translational Implications

The identified cell-cell interactions offer compelling insights for therapeutic targeting and experimental validation in breast cancer.

  1. Targeting Integrin-Mediated Adhesion: Given the pervasive role of integrins in mediating interactions between tumor cells (Aneuploid Epithelial) and the TME (Fibroblasts, Endothelial cells), integrin inhibitors represent a promising class of therapeutic agents to disrupt tumor cell adhesion, migration, and invasion. Experimental validation could involve *in vitro* assays of tumor cell migration/invasion in 3D co-culture models with CAFs, and *in vivo* studies using integrin-blocking antibodies to assess their impact on primary tumor growth and metastasis.
  2. Anti-Angiogenic Strategies: The strong PGF-FLT1 and VEGF-VEGFR2 signaling pathways highlight the importance of angiogenesis in primary breast tumors. Existing anti-VEGF therapies (e.g., bevacizumab) or novel inhibitors targeting PGF or FLT1 could be investigated, potentially in combination with other agents, to effectively curb tumor angiogenesis and growth. Clinical trials could evaluate their efficacy in breast cancer patients with similar CCI profiles. PubMed Search: Anti-angiogenic therapy breast cancer
  3. Modulating Immune and Stromal Communications: The CXCL12-CXCR4 axis represents a potential target to alter the immunosuppressive and pro-metastatic TME. CXCR4 antagonists could be explored to reduce recruitment of immunosuppressive cells and inhibit metastasis. Similarly, targeting CCL2-CCR1 could modulate macrophage infiltration and polarization. Experimental validation would involve assessing the impact of such interventions on immune cell composition and function within the TME, as well as tumor growth and metastasis.
  4. IL-6 Pathway Inhibition: Given the pro-tumorigenic role of IL-6, inhibitors of IL-6 or its receptor could be explored, particularly in patients exhibiting high IL-6 signaling.
  5. Biomarker Development: The specific ligand-receptor pairs identified as highly active could serve as potential biomarkers for patient stratification, predicting response to targeted therapies, or identifying patients at higher risk of aggressive disease. Further research is needed to correlate these CCI patterns with clinical outcomes.

13. Cell-Cell Interaction Analysis of Immune Checkpoint and Cell Cycle Related Genes in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interactions (CCI) using CellPhoneDB results derived from single-cell RNA-seq data of human breast tissue, comparing normal and primary tumor conditions. The analysis specifically filters for a predefined list of genes associated with immune checkpoint pathways and cell cycle regulation. The plot_cci_dots tool visualizes significant ligand-receptor interactions between different cell types, with dot size representing the statistical significance (-log10(p-value)) and dot color indicating the interaction strength (log2(mean expression)). A key aspect of this analysis is the differentiation between Diploid Epi (diploid epithelial cells) and Aneuploid Epi (aneuploid epithelial cells, representing the likely malignant population in the tumor context, as Epithelial cell is the tumor origin celltype).

Visual Summary

Comparison of Normal vs. Primary Tumor:

Specific Observations by Condition:

Normal Tissue CCI:

Primary Tumor Tissue CCI:

Biological Interpretation

The analysis reveals a profound shift in cellular communication within the breast tissue transitioning from a normal to a primary tumor state, particularly when focusing on genes related to immune checkpoints and cell cycle regulation.

  1. TGF-beta Signaling as a Central Driver of Tumor Microenvironment Remodeling: The consistent and strong activation of TGF-beta signaling across numerous cell-cell interfaces in the primary tumor is a critical finding. In the normal tissue, TGF-beta primarily mediates homeostatic epithelial-endothelial and stromal interactions. However, in the tumor, Aneuploid Epi cells become key producers and responders to TGF-beta, engaging in autocrine loops and paracrine interactions with Endothelial cells, Fibroblasts, and Macrophages. This is highly relevant as TGF-beta is a pleiotropic cytokine with dual roles: initially tumor-suppressive, but later promoting tumor progression by fostering epithelial-mesenchymal transition, enhancing angiogenesis, suppressing anti-tumor immunity, and driving desmoplasia [PMID: 32669614]. The strong Mac|Aneuploid Epi interactions via TGF-beta suggest potential M2-like macrophage polarization, contributing to an immunosuppressive TME.
  2. Altered EGFR Signaling in Tumor Progression: While EGFR signaling is present in normal breast tissue, its ligand repertoire and interacting partners change in the primary tumor. The emergence of EREG_EGFR interactions and its involvement with Aneuploid Epi suggests a specific activation of the EGFR pathway that supports tumor cell proliferation and survival. EGFR is a well-established oncogenic driver, and its altered signaling profile in Aneuploid Epi highlights its role in malignant transformation and progression [GeneCards: EGFR].
  3. Complex Immune Cell Crosstalk in the TME: The primary tumor environment shows a significant increase in interactions involving various immune cell types (T cells, Macrophages, Mast cells) with tumor cells and stromal components. While direct PD-1/PD-L1 interactions were not highlighted as top signals in this specific analysis, the presence of CD86_CD28 interactions, even if primarily stromal, indicates costimulatory molecule engagement. LCK_CD8_receptor and CD83_IFNGR1 interactions point to T cell receptor signaling and interferon pathway engagement, respectively, at the tumor-stromal-immune interfaces. The overall increase in immune cell interactions, particularly those involving Macrophages and T cells with Aneuploid Epi cells often through TGFB signaling, suggests a dynamic, yet potentially immunosuppressive, immune landscape characteristic of many solid tumors.
  4. Role of Cell Cycle Regulators (Indirectly via Growth Factors): Many of the filtered "cell cycle" genes are intracellular effectors (CDKs, cyclins, MCMs). Their influence on cell-cell communication is primarily indirect, by modulating the expression or activity of growth factor receptors or their ligands. The observed upregulation of EGFR and TGF-beta signaling directly impacts cell proliferation and cell cycle progression in the tumor cells, illustrating how external signals from the TME can drive aberrant cell cycle activity in cancer.

Clinical or Translational Implications

The observed patterns of cell-cell interactions have significant clinical and translational implications, particularly for therapeutic targeting and understanding tumor biology.

  1. TGF-beta Pathway as a Promising Therapeutic Target: The widespread and robust TGF-beta signaling in the primary tumor, particularly involving Aneuploid Epi cells, stromal, and immune cells, positions the TGF-beta pathway as a high-priority therapeutic target. Inhibiting TGF-beta signaling could disrupt multiple pro-tumorigenic processes, including tumor cell proliferation, immune evasion, and fibrosis. Various agents targeting TGF-beta, such as receptor kinase inhibitors or ligand traps, are currently under investigation for cancer treatment [PubMed search: TGF-beta inhibitor cancer therapy].
  2. EGFR Pathway Targeting in Aneuploid Epithelial Cells: The strong EGFR signaling, especially involving EREG_EGFR with Aneuploid Epi cells, suggests that these malignant cells may be susceptible to EGFR inhibitors. Further validation could determine if specific EGFR inhibitors (e.g., gefitinib, erlotinib, or afatinib) could be effective in a subset of breast cancer patients characterized by this signaling profile [GeneCards: EGFR].
  3. Modulating the Tumor-Immune Microenvironment: The extensive immune-stromal-tumor interactions, particularly the role of TGF-beta in macrophage-tumor crosstalk, highlight the potential for combination therapies that target both tumor cell intrinsic pathways and the immunosuppressive TME. For instance, combining TGF-beta inhibitors with immune checkpoint blockade could potentially overcome resistance mechanisms and enhance anti-tumor immunity by reprogramming the TME and reversing immune suppression [PubMed search: TGF-beta immune checkpoint blockade combination].
  4. Biomarker Discovery: The specific ligand-receptor pairs identified as highly active in the primary tumor, such as TGFB1_TGFbeta_receptor1 or EREG_EGFR on Aneuploid Epi cells, could serve as potential biomarkers for patient stratification or response prediction to targeted therapies. Further studies could investigate their expression levels and functional relevance in patient samples.

14. Condition-Specific Cell-Cell Interaction Patterns in Breast Cancer

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates statistically significant differences in cell-cell interactions (CCI) between normal breast tissue and primary breast tumors, focusing on major immune and stromal cells. The dot plot visualizes the strength (standardized mean, color intensity) and statistical significance (-log10(p), dot size) of various ligand-receptor pairs (CCI indices) across individual patient samples within each condition. This approach helps identify specific intercellular communication pathways that are dysregulated in the tumor microenvironment.

Visual Summary

The dot plot clearly differentiates CCI patterns between normal and primary tumor samples:

Biological Interpretation

The observed differences highlight a profound remodeling of the cellular communication landscape in primary breast tumors compared to normal tissue.

  1. Tumor Microenvironment Remodeling by ECM-Integrin Axis: The most prominent finding is the massive upregulation of interactions centered around collagens, fibronectin, and integrins in primary tumors.
  1. Normal Tissue Homeostasis and Surveillance: The interactions observed in normal tissue likely reflect physiological processes.

Clinical or Translational Implications

15. Epithelial Cell Condition-Specific Surfaceome Markers in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers for Epithelial cells from human breast tissue, comparing cells derived from "Diploid" samples (representing normal tissue, given the reference condition for DEG is 'normal' and the sample grouping) versus "primary_tumor" samples. The dot plot visualizes the expression of up to 50 top surfaceome markers per condition, showing both the fraction of cells expressing each marker (dot size) and the mean expression level within expressing cells (dot color intensity). The goal is to identify unique surface markers that distinguish normal epithelial cells from tumor epithelial cells, which can have implications for diagnosis and targeted therapies.

Visual Summary

The dot plot is organized to show samples (rows) grouped first by ploidy inference ("Diploid" for the upper cluster, which represents normal tissue epithelial cells in this context) and then by condition ("primary_tumor" for the lower cluster, representing tumor epithelial cells). The columns represent individual surfaceome genes.

Biological Interpretation

The analysis highlights significant differences in the surface proteome of epithelial cells between normal-like (diploid) and primary tumor conditions in breast tissue. These differences reflect the altered biological state and functions of cancer cells.

Clinical or Translational Implications

The identified condition-specific surfaceome markers in epithelial cells hold significant clinical and translational potential for breast cancer.

16. Macrophage Condition-Specific Surfaceome Markers in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies and visualizes condition-specific surfaceome markers in Macrophage cells isolated from human breast tissue, comparing normal tissue samples with primary_tumor samples. The dot plot illustrates the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) across different patient samples within each condition. Focusing on surfaceome markers is particularly relevant for therapeutic targeting and cell-type characterization.

Visual Summary

The dot plot clearly segregates macrophage surface markers into two distinct groups based on their expression patterns in normal versus primary_tumor conditions.

Biological Interpretation

The distinct sets of surface markers highlight the profound functional reprogramming and adaptation of macrophages in response to the tumor microenvironment compared to their roles in normal tissue homeostasis.

Normal Tissue Macrophages:

These markers collectively suggest a macrophage phenotype oriented towards immune surveillance, tissue maintenance, and resolution of inflammation in a healthy state.

Primary Tumor-Associated Macrophages (TAMs):

The broad upregulation of various surface markers in TAMs reflects their multifaceted roles in tumor progression:

This analysis underscores the highly adaptable nature of macrophages, which acquire distinct phenotypes and functions in the tumor microenvironment to promote cancer progression.

Clinical or Translational Implications

The identified condition-specific surfaceome markers in macrophages offer significant clinical and translational potential:

Modulating the activity or depletion of TAMs could reprogram the tumor microenvironment, making it less permissive for tumor growth and potentially more susceptible to other immunotherapies.

17. Fibroblast Condition-Specific Surfaceome Markers in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify surfaceome markers that are specifically expressed in Fibroblasts under different conditions (normal vs. primary tumor) in human breast tissue. Using single-cell RNA-seq data, the plot_markers_and_expression_dot tool generated a dot plot visualizing the expression patterns of these markers across individual patient samples. The analysis focused on genes encoding surface proteins to identify potential therapeutic targets or biomarkers accessible on the cell surface.

Visual Summary

The dot plot effectively illustrates the differential expression of surfaceome markers in Fibroblasts across various patient samples, categorized by their condition (normal or primary tumor).

Biological Interpretation

The identified surfaceome markers provide strong biological insights into the distinct roles and activation states of Fibroblasts in the breast tumor microenvironment.

Key CAF-associated Surface Markers:

Clinical or Translational Implications

The identification of these condition-specific surfaceome markers in Fibroblasts has significant clinical and translational potential.

18. Differential Expression of Cell Cycle-Related Gene YWHAZ in Breast Epithelial Cells

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the differential gene expression of a pre-defined set of cell cycle-related genes between primary tumor and normal conditions, specifically within Epithelial cells (identified as the major disease-related cell type and tumor origin cell type in this dataset). The aim was to identify genes with statistically significant expression differences (p-value < 0.1 and absolute log2 Fold Change > 0.1) and visualize their expression patterns. The provided boxplot highlights the expression of YWHAZ, a key cell cycle regulator, showing its distribution across the two conditions.

Visual Summary

The boxplot displays the gene expression levels of YWHAZ (Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Zeta) in Epithelial cells, comparing primary_tumor and normal conditions.

Biological Interpretation

The significant upregulation of YWHAZ in primary tumor Epithelial cells suggests its potential involvement in breast cancer progression, particularly through its role in cell cycle regulation.

Clinical or Translational Implications

The finding that YWHAZ is significantly overexpressed in primary breast tumor Epithelial cells opens several potential clinical and translational avenues:

---

References:

  1. YWHAZ (14-3-3 zeta) gene information: GeneCards. https://www.genecards.org/cgi-bin/carddisp.pl?gene=YWHAZ
  2. 14-3-3 Proteins in Cell Cycle Regulation: PubMed search for "14-3-3 cell cycle regulation". https://pubmed.ncbi.nlm.nih.gov/?term=14-3-3+cell+cycle+regulation
  3. Role of YWHAZ in Cancer: PubMed search for "YWHAZ breast cancer proliferation". https://pubmed.ncbi.nlm.nih.gov/?term=YWHAZ+breast+cancer+proliferation

19. Gene Ontology (GSA) Analysis of Epithelial Cells in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Ontology (GO) enrichment results (using Gene Set Analysis, GSA) for epithelial cells, a key cell type in breast tissue and the identified tumor origin cell type. The analysis compares gene expression profiles in two distinct contexts:

  1. Diploid_vs_others: Epithelial cells classified as "Diploid" (likely representing less transformed or normal-like cells) versus all other epithelial cells.
  2. primary_tumor_vs_others: Epithelial cells from "primary_tumor" condition versus all other conditions (presumably normal tissue epithelial cells).

The bar plots visualize the most significantly enriched GO terms, ranked by their statistical significance (-log(p-val) and -log(q-val)). These enriched terms highlight the biological processes and pathways that are differentially active in the specified comparison groups within the epithelial cell population.

Visual Summary

The two bar plots display the top significantly enriched Gene Ontology terms. Each plot consists of two panels: one showing -log(p-val) and the other showing -log(q-val) for each term. Higher values indicate greater statistical significance.

GSA for Epithelial cell: Diploid_vs_others

This plot shows terms enriched in diploid epithelial cells compared to non-diploid (likely aneuploid) epithelial cells.

Highly significant terms (high -log(p-val) and -log(q-val)) include

"Coronavirus disease"

"Focal adhesion"

"Ribosome biogenesis in eukaryotes"

"Epstein-Barr virus infection"

"Kaposi sarcoma-associated herpesvirus infection"

"Ribosome"

"TNF signaling pathway"

"Proteoglycans in cancer"

"Cellular senescence"

"Pathways in cancer"

GSA for Epithelial cell: primary_tumor_vs_others

This plot shows terms enriched in epithelial cells from primary tumors compared to epithelial cells from normal tissue.

Highly significant terms (high -log(p-val) and -log(q-val)) include

"Thermogenesis"

"Oxidative phosphorylation"

"Huntington disease"

"Parkinson disease"

"Prion disease"

"Non-alcoholic fatty liver disease"

"Alzheimer disease"

"Diabetic cardiomyopathy"

"Pathways of neurodegeneration"

"Amyotrophic lateral sclerosis"

"Endocytosis"

"Protein processing in endoplasmic reticulum"

"Lysosome"

"Ubiquitin mediated proteolysis"

Biological Interpretation

Epithelial Cells: Diploid vs. Other Ploidy States

The enrichment of terms like "Focal adhesion" and "Ribosome biogenesis" in diploid epithelial cells suggests active baseline cellular functions related to cell structure, adhesion, and protein synthesis. The presence of "Cellular senescence" and "Pathways in cancer" in this context could indicate that diploid epithelial cells, especially if they are adjacent to tumor regions or under oncogenic stress, might be engaging in tumor-suppressive mechanisms or early cellular responses to transformation. "TNF signaling pathway" and pathways related to various viral infections (e.g., Epstein-Barr virus, Human papillomavirus) might reflect either general immune responses present in the tissue or a specific host-response pattern in these less transformed cells, potentially against oncogenic viruses known to be implicated in some cancers. The "p53 signaling pathway" and "NF-kappa B signaling pathway" are critical regulators of cell cycle, apoptosis, and inflammation, which are active in maintaining cellular homeostasis or initiating stress responses.

Epithelial Cells: Primary Tumor vs. Normal Tissue

Epithelial cells from primary tumors exhibit a distinct biological signature compared to normal epithelial cells. The most prominent enrichments point towards significant metabolic reprogramming and cellular stress responses, characteristic hallmarks of cancer cells.

In summary, diploid epithelial cells appear to be engaging in basal maintenance, stress responses, and potentially early tumor-suppressive mechanisms. In contrast, primary tumor epithelial cells exhibit profound metabolic shifts, heightened protein processing and degradation, and widespread cellular stress responses, consistent with their high proliferative rate and adaptation to the tumor microenvironment.

Clinical or Translational Implications

The distinct pathway enrichments in primary tumor epithelial cells offer potential avenues for therapeutic intervention and biomarker discovery in breast cancer:

Understanding these context-specific pathway activations provides valuable insight into the underlying biology of breast cancer progression and potential therapeutic vulnerabilities.

20. Gene Set Enrichment Analysis (GSEA) of Breast Tissue Cell Types in Primary Tumor vs. Normal Conditions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis utilizes Gene Set Enrichment Analysis (GSEA) to identify enriched or depleted biological pathways across key major cell types (Endothelial cell, Epithelial cell, Fibroblast, ILC, Macrophage, Smooth muscle cell) in breast tissue. The comparisons include cells from normal vs. others conditions, and primary_tumor vs. others conditions. For Epithelial cells, an additional comparison was performed for Diploid cells vs. others, providing insight into the role of ploidy. The dot plot visualizes the Normalized Enrichment Score (NES) and statistical significance (-log(p-val)) for 80 selected gene sets, highlighting significant shifts in cellular biology within the tumor microenvironment.

Visual Summary

The dot plot effectively displays the GSEA results, with cell type and comparison conditions on the X-axis and biological pathways (gene sets) on the Y-axis.

Biological Interpretation

Widespread Metabolic Reprogramming in the Tumor Microenvironment

A striking and consistent finding across almost all examined cell types (Endothelial, Epithelial, Fibroblast, ILC, Macrophage, Smooth muscle cell) in primary_tumor_vs_others comparisons is the strong enrichment of metabolic pathways such as Glycolysis / Gluconeogenesis and Pyruvate metabolism. This indicates a systemic shift towards aerobic glycolysis, often referred to as the Warburg effect, which provides energy and building blocks for rapid proliferation not just in cancer cells but also in associated stromal and immune cells. This metabolic reprogramming is a hallmark of cancer and profoundly influences the tumor microenvironment.

Enhanced ECM Remodeling and Cell Adhesion

Pathways related to extracellular matrix (ECM) interaction and cell adhesion, specifically ECM-receptor interaction and Focal adhesion, are significantly enriched in Endothelial cells, Epithelial cells, Fibroblasts, Macrophages, and Smooth muscle cells within primary_tumor_vs_others. This highlights the extensive remodeling of the tumor microenvironment, which is critical for tumor growth, invasion, and metastasis in breast cancer.

Reference: PubMed search: ECM remodeling cancer metastasis

Activated Pro-Tumorigenic Signaling Pathways

Several crucial signaling pathways are consistently enriched in primary tumor samples across multiple cell types:

Functional Reprogramming of Stromal and Immune Cells

Distinct Biology of Diploid Epithelial Cells

The comparison of Epithelial cell: Diploid_vs_others versus Epithelial cell: primary_tumor_vs_others reveals critical differences:

Clinical or Translational Implications

  1. Biomarker Discovery: The consistently enriched pathways (e.g., VEGF signaling, Wnt signaling, metabolic pathways, ECM interaction) across multiple cell types in primary tumors represent a rich source for potential diagnostic, prognostic, or predictive biomarkers in breast cancer.
  2. Therapeutic Targets: The observed widespread metabolic reprogramming (glycolysis) and activation of specific signaling pathways (VEGF, Wnt, AGE-RAGE) across cancer cells and the tumor microenvironment suggest that therapeutic strategies targeting these common vulnerabilities could have broad anti-tumor effects, influencing not only cancer cells but also their supportive stromal and immune components.
  1. Tumor Microenvironment-Focused Therapies: The active and integrated roles of stromal (fibroblasts, endothelial cells, smooth muscle cells) and immune cells (macrophages, ILCs) in supporting tumor growth emphasize the importance of developing therapies that consider the entire tumor ecosystem, rather than focusing solely on cancer cells.
  2. Ploidy as a Prognostic Factor: The distinct pathway activities associated with Diploid versus primary_tumor epithelial cells underscore the potential of ploidy status as a significant indicator of cellular state and likely tumor aggressiveness. Further investigation into the fate and response to therapy of diploid tumor cells or co-existing normal cells within the tumor could inform personalized treatment strategies.

21. Discussion

The comprehensive single-cell analysis of breast tissue provides profound insights into the distinct cellular and molecular features of the primary tumor microenvironment compared to normal tissue. A cornerstone finding is the widespread genomic instability, evidenced by extensive aneuploidy, predominantly observed in epithelial cells within primary tumors. UMAP visualizations based on CNV estimates clearly segregate aneuploid epithelial cells from diploid counterparts and other cell types, strongly supporting their malignant nature and highlighting genomic alteration as a primary driver of cellular heterogeneity in cancer. Recurrent amplifications on chromosomes 1q, 8q, 11q, and 20q, including genes like *NFASC* and *EIF3E*, are consistent with known genomic aberrations in breast cancer, underscoring their potential role in tumorigenesis and validating the ploidy inference.

The cellular composition shifts dramatically in the primary tumor. Epithelial cells remain dominant, but a significant proportion of 'unassigned' cells in tumors may represent highly aberrant or novel malignant cell states. Fibroblasts are consistently abundant, indicative of an activated stromal compartment, while macrophages show a noticeable increase, pointing towards heightened immune cell infiltration. Within immune populations, T cell subsets are highly variable, with some tumors exhibiting increased cytotoxic T cells. Notably, both T follicular helper (Tfh) and regulatory T (Treg) cells are significantly elevated in primary tumors, while Lymphoid Tissue Inducer (LTI) cells are reduced. Macrophages in tumors show an increased proportion of M1 phenotypes, suggesting an anti-tumor inflammatory response, yet M2-like subsets remain prominent, implying a complex and often immunosuppressive immune landscape. This co-existence of M1 and M2 phenotypes in tumors, rather than a clear M2 dominance often reported, suggests a more nuanced, dynamic balance of macrophage polarization that could vary greatly between patients and within different tumor regions.

Cell-cell interaction (CCI) analyses reveal extensive reprogramming of intercellular communication within the tumor. A pervasive feature is the highly active extracellular matrix (ECM)-integrin signaling, primarily between cancer-associated fibroblasts (CAFs), endothelial cells, and epithelial cells, involving collagens (e.g., COL6A1/2/3, COL18A1) and fibronectin (FN1) with various integrin complexes. This indicates significant ECM remodeling and mechanotransduction driving tumor invasion and angiogenesis. Other critical pathways include robust pro-angiogenic signaling (PGF-FLT1, VEGF-VEGFR2) and active TGF-beta signaling, particularly originating from aneuploid epithelial cells and extensively involving fibroblasts and macrophages, which is known to promote immune suppression and desmoplasia. Epithelial cells from primary tumors also show upregulation of surface markers like ERBB3, CA12, ADAM15, and ALCAM, reflecting their transformed state, while tumor-associated macrophages (TAMs) express unique markers such as EREG, ITGB8, CXCR4, and MET. CAFs themselves acquire distinct surfaceome markers including FAP, CD276, and MMP14. The finding that YWHAZ, a cell cycle regulator, is significantly upregulated in tumor epithelial cells further highlights altered proliferative control.

Gene Ontology (GSA) and Gene Set Enrichment Analysis (GSEA) reinforce these findings, showing widespread metabolic reprogramming (glycolysis, pyruvate metabolism) across nearly all cell types in the primary tumor microenvironment, not just the malignant cells. This generalized metabolic shift, coupled with enriched ECM-receptor interaction, focal adhesion, and VEGF/Wnt signaling pathways, paints a picture of a highly collaborative and metabolically active ecosystem supporting tumor growth. The distinct biological profiles of diploid versus aneuploid epithelial cells, with diploid cells enriched for basic cellular maintenance and the primary tumor cells exhibiting profound metabolic shifts and stress responses, further delineates the molecular consequences of genomic instability.

Hypotheses:

  1. Aneuploid epithelial cells, through aberrant WNT and TGF-beta signaling, directly orchestrate the activation and ECM-remodeling functions of cancer-associated fibroblasts (CAFs), driving tumor invasion and metastasis.
  2. The observed increase in T follicular helper (Tfh) and regulatory T (Treg) cells, alongside specific pro-tumorigenic macrophage (TAM) populations characterized by markers like CXCR4 and EREG, collectively contribute to an immunosuppressive tumor microenvironment, enabling immune evasion.
  3. Widespread metabolic reprogramming, including enhanced glycolysis and pyruvate metabolism across tumor epithelial cells, fibroblasts, and macrophages, creates a synergistic metabolic landscape that fuels tumor growth and suppresses anti-tumor immunity.
  4. The reduction of Lymphoid Tissue Inducer (LTI) cells in primary tumors leads to impaired formation or functionality of tertiary lymphoid structures (TLS), diminishing local anti-tumor immune responses.

Potential therapeutic targets:

  1. TGF-beta pathway (e.g., TGFB1, TGF-beta receptors): TGF-beta signaling is a central driver of pro-tumorigenic processes in the breast cancer microenvironment, extensively mediating interactions between aneuploid epithelial cells, fibroblasts, and macrophages to promote immune evasion, fibrosis, and tumor progression. Evidence: CCI analyses (Sections 11, 13) show strong TGFB1/2/3 and TGF-beta receptor interactions across multiple cell types in primary tumors. GSEA (Section 20) indicates widespread enrichment of TGF-beta related pathways. Macrophage-epithelial interactions via TGFB ligands suggest potential M2-like macrophage polarization. Validation: Evaluate the efficacy of TGF-beta inhibitors (receptor kinases or ligand traps) in breast cancer preclinical models (e.g., PDX, organoids) to reduce tumor growth, metastasis, CAF activation (FAP, collagen deposition), and immune suppression (Treg/TAM modulation). Assess changes in target engagement via IHC/western blot.
  2. FAP (Fibroblast Activation Protein): FAP is a canonical surface marker highly and specifically expressed on activated Cancer-Associated Fibroblasts (CAFs), which are critical for ECM remodeling, immune suppression, and supporting tumor growth and invasion. Evidence: Fibroblast condition-specific markers (Section 17) show FAP as a top upregulated surfaceome marker in primary tumor fibroblasts. Validation: Develop or utilize FAP-targeting antibody-drug conjugates (ADCs) or CAR-T cell therapies in *in vivo* breast cancer models. Validate target specificity and anti-tumor efficacy, including reduction in tumor volume, metastasis, and changes in tumor stromal composition (e.g., collagen content, other CAF markers).
  3. ERBB3 (HER3): ERBB3 is a surface receptor specifically overexpressed on primary tumor epithelial cells, a key oncogenic driver in breast cancer, often associated with tumor growth, survival, and therapy resistance. Evidence: Epithelial cell condition-specific surfaceome markers (Section 15) clearly show ERBB3 as highly upregulated and specific to primary tumor epithelial cells. Validation: Test ERBB3-targeting agents (e.g., monoclonal antibodies, ADCs) in *in vitro* 3D cultures of patient-derived tumor epithelial cells and *in vivo* xenograft models. Assess impact on cell proliferation, survival, and sensitivity to standard-of-care therapies. Confirm ERBB3 protein expression in clinical samples via IHC.
  4. CXCR4: CXCR4 is a chemokine receptor highly expressed on tumor-associated macrophages and involved in fibroblast-fibroblast interactions, crucial for immune cell migration, tumor cell invasion, and creating an immunosuppressive/pro-metastatic microenvironment. Evidence: Macrophage condition-specific surfaceome markers (Section 16) show CXCR4 as a top upregulated marker in TAMs. CCI analysis (Section 12) identifies CXCL12-CXCR4 interactions particularly between fibroblasts and macrophages. Validation: Utilize CXCR4 antagonists in preclinical breast cancer models to evaluate their ability to reduce TAM and Treg infiltration, inhibit tumor cell invasion/metastasis, and enhance the efficacy of immunotherapies. Monitor changes in immune cell populations via flow cytometry from tumor biopsies.
  5. Integrin receptors (e.g., αVβ1, αVβ5, α5β1): Integrins are broadly involved in extensive ECM-receptor interactions and focal adhesion across epithelial, fibroblast, and endothelial cells in primary tumors, mediating tumor growth, invasion, angiogenesis, and resistance to therapy. Evidence: CCI analyses (Sections 11, 12, 14) consistently highlight diverse integrin complexes (e.g., FN1-integrin_a3b1, COL6A1-integrin_a1b1) as highly active, especially in primary tumor conditions and across tumor-stromal-endothelial cell interfaces. GSEA (Section 20) shows enrichment of 'ECM-receptor interaction' and 'Focal adhesion' pathways. Validation: Employ integrin-blocking antibodies or small molecule inhibitors (specific to αVβ1, αVβ5, α5β1) in *in vitro* assays for tumor cell adhesion, migration, and invasion on various ECM substrates. *In vivo*, assess their impact on primary tumor growth, angiogenesis (CD31 staining), and metastasis in relevant breast cancer models.

Follow-up validation ideas:

  1. Perform spatial transcriptomics or multiplexed immunostaining (e.g., CODEX, IMC) to validate the co-localization and direct interaction of aneuploid epithelial cells, CAFs (FAP+, ACTA2+), and immune cells (Tregs, TAMs) within tumor tissues, particularly focusing on WNT and TGF-beta ligand-receptor co-expression.
  2. Utilize *in vitro* co-culture models of tumor epithelial cells, fibroblasts, and macrophages, coupled with targeted perturbation assays (e.g., siRNA knockdown, small molecule inhibitors of WNT, TGF-beta, or CXCR4), to assess their impact on ECM deposition, cell migration, invasion, and immune cell polarization (e.g., flow cytometry for M1/M2 markers).
  3. Conduct metabolic tracing experiments (e.g., 13C-glucose/glutamine tracing) on sorted tumor epithelial cells, fibroblasts, and macrophages from both normal and tumor conditions to quantitatively validate the altered metabolic fluxes, particularly glycolysis and oxidative phosphorylation, identified by GSEA.
  4. Validate the functional impact of reduced LTI cells by comparing the presence and organization of tertiary lymphoid structures (TLS) in breast cancer patients with high versus low LTI cell proportions via immunohistochemistry. Perturbation of LTI cells in *in vivo* tumor models could assess its effect on anti-tumor immunity and tumor progression.

Limitations:

This report is based on single-cell RNA-sequencing data, which provides transcriptomic snapshots. While CNV estimates are included, they are inferred from RNA-seq and may not fully capture the complexity of genomic alterations compared to dedicated genomic assays. Cell-cell interaction predictions are computational and require experimental validation to confirm functional relevance. The 'unassigned' cell populations in tumor samples warrant further investigation. The dataset represents human breast tissue; however, inter-patient heterogeneity, as observed, suggests that generalizability to all breast cancer subtypes or patients may require broader validation cohorts. Causality of observed associations cannot be definitively established from this correlative analysis.

22. Query List

  1. Show UMAPs including condition, sample, major cell type, minor cell type, ploidy_dec, and celltype_subset in 2 columns and save.
  2. Show expression levels of CD3D, CD4, CD8A, CD79A, MS4A1, MZB1, CD14, LYZ, FBLN1, NOTCH3, EPCAM, MUC1, CD34 genes on UMAP along with minor cell type annotations. Set ncols=4 and save.
  3. Show the overall Celltype_subset marker expression dot plot and save it. Keep target_cell as None and keep the other arguments at their default values.
  4. Filter for tumor origin cells (Epithelial cell) and unassigned cells, group by sample, show CNV heatmap, and include a summary of significantly amplified copy number regions, then save.
  5. Show CNV patterns on UMAP, including major cell type, minor cell type, ploidy results, condition, and sample in 2 columns, then save.
  6. Show a population bar plot of minor cell types and save.
  7. Show a subset population bar plot for T cells and save.
  8. Show a subset population bar plot for Macrophages and save.
  9. For T cell subset populations, show a boxplot for statistically significant differences between conditions if any, and save. Determine ncols appropriately based on the total number of panels.
  10. Filter for tumor origin cells (Epithelial cell) and unassigned cells, show their ploidy population as a bar plot, and save.
  11. Show cell-cell interaction patterns per condition, including tumor origin cells (Epithelial cell), Fibroblast, Macrophage, and T cells. Select up to 80 cell-cell interactions for each condition and save.
  12. Show cell-cell interactions by condition and save the result. Select at most 80 cell-cell interactions for each condition.
  13. Filter for immune checkpoint and cell cycle pathway related genes, show cell-cell interactions for these genes, and save.
  14. For major immune and stromal cells, find statistically significant differences in cell-cell interactions between conditions and show them as a dot plot. Set max_n_items_per_group = 25 and save.
  15. Show the condition-specific markers for tumor-origin cells (Epithelial cell) as a dot plot and save the result. Use only surfaceome markers, up to 50 markers per condition.
  16. Extract condition-specific markers for Macrophages, show them as a dot plot, and save. Only show surfaceome markers, up to 50 per condition.
  17. Extract condition-specific markers for Fibroblasts, show them as a dot plot, and save. Only show surfaceome markers, up to 50 per condition.
  18. For cell cycle pathway related genes, select genes with statistically significant expression differences between conditions in major disease-related cells (Epithelial cell), show them as a boxplot, and save. Set max_n_items_to_plot = 24, and determine ncols such that the aspect ratio of width to height is approximately 2x3 based on the total number of panels.
  19. Show Gene Ontology (GSA) analysis results for epithelial cells as bar plots and save the result.
  20. Show a dot plot of Gene Set Enrichment Analysis results for major cell types (Endothelial cell, Epithelial cell, Fibroblast, ILC, Macrophage, Smooth muscle cell). Set color map to RdBu_r and n_pws_to_show = 80, then save.
↑ Top