SCODiA Report by MLBI Lab

Single-Cell Multi-Omics Analysis Reveals Tumor Microenvironment Remodeling, Genetic Instability, and Pathway Dysregulation in Primary Breast Cancer

This single-cell RNA sequencing report details profound biological shifts distinguishing normal breast tissue from primary tumors. UMAP visualizations show clear segregation of normal and malignant cell populations, with tumor epithelial cells exhibiting prevalent aneuploidy and widespread genomic instability. The tumor microenvironment undergoes extensive remodeling, marked by significant alterations in cell type proportions, enhanced cell-cell interactions, and global pathway dysregulation in tumor epithelial cells, cancer-associated fibroblasts, and tumor-associated macrophages. These findings collectively highlight the complex cellular and molecular landscape driving breast cancer progression.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell Transcriptomic Data Colored by Key Metadata
  3. Major Cell Type Score and Ploidy Visualization on UMAP
  4. Celltype_subset Marker Expression Analysis
  5. Copy Number Variation Analysis in Breast Cancer Tumor and Unassigned Cells Grouped by Sample
  6. CNV-based UMAP Visualization of Cell Type, Ploidy, Condition, and Sample
  7. Minor Cell Type Population Analysis in Breast Tissue
  8. Immune Lymphoid Cell Subset Composition in Normal and Primary Breast Tumor Tissues
  9. Primary Breast Tumor Condition Shows Altered T Cell Subset Proportions
  10. Macrophage Cell Population Validation Across Samples
  11. Ploidy Analysis of Epithelial and Unassigned Cells in Breast Tissue
  12. 유방암 미세환경 내 세포-세포 상호작용 패턴 분석
  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 Tissue
  16. Epithelial Cell Condition-Specific Surfaceome Markers in Breast Cancer
  17. Macrophage Condition-Specific Surfaceome Markers in Breast Tissue
  18. Fibroblast Condition-Specific Surfaceome Markers in Breast Tissue
  19. Cell-Type-Specific Surfaceome Marker Discovery in Breast Tissue
  20. Differential Expression of Cell Cycle-Related Genes in Breast Cancer Epithelial Cells: YWHAZ Upregulation
  21. Epithelial Cell Gene Ontology (GSA) Analysis: Condition and Ploidy-Specific Pathway Enrichment in Breast Tissue
  22. Gene Set Enrichment Analysis Reveals Widespread Pathway Dysregulation in Breast Cancer Microenvironment
  23. Discussion
  24. Query List

0. Dataset overview

Dataset Summary

1. UMAP Visualization of Single-Cell Transcriptomic Data Colored by Key Metadata

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents UMAP (Uniform Manifold Approximation and Projection) plots generated from single-cell RNA-seq data, providing a two-dimensional visualization of cellular relationships based on gene expression profiles. The UMAPs are colored by various metadata features, including condition (normal vs. primary_tumor), sample (individual patient samples), celltype_major, celltype_minor, celltype_subset (different levels of cell type annotation), and ploidy_dec (ploidy inference: Aneuploid/Diploid). These visualizations are crucial for assessing the overall structure of the dataset, evaluating the quality of cell type annotations, checking for batch effects, and understanding the distribution of different biological states across the cell populations.

Visual Summary

Condition UMAP

The UMAP colored by condition shows a clear separation between cells originating from normal tissue and primary_tumor tissue. While there are distinct clusters largely dominated by one condition, some regions exhibit a mixture of both normal and tumor cells, particularly in areas likely representing immune or stromal infiltrates within the tumor microenvironment. The bulk of the primary_tumor cells form several large, distinct clusters.

Sample UMAP

The sample UMAP reveals that cells from different patient samples are generally well-mixed within their respective cell type clusters. This indicates successful integration of data from multiple samples, suggesting that batch effects due to individual patient variability or experimental processing are largely mitigated, allowing for robust comparisons across samples.

Celltype_major, Celltype_minor, and Celltype_subset UMAPs

These three UMAPs progressively display more granular cell type annotations:

Ploidy_dec UMAP

The ploidy_dec UMAP shows Aneuploid cells predominantly localized within specific clusters, which largely overlap with the primary_tumor condition in the 'condition' UMAP and specifically with the Epithelial cell populations in the cell type UMAPs. Diploid cells are widespread across most clusters, particularly those associated with normal tissue and non-malignant cell types (e.g., stromal, immune, endothelial cells) within both normal and tumor contexts. A small number of Unclear ploidy cells are also present, mostly intermingled with diploid cells.

Biological Interpretation

The UMAP visualizations provide strong evidence for distinct cellular landscapes in normal breast tissue versus primary breast tumors. The clear separation of primary_tumor cells from normal cells on the condition UMAP highlights the profound transcriptomic changes associated with cancer development. The tumor microenvironment is complex, and the observed mixing of normal and primary_tumor cells in certain regions likely represents the presence of infiltrating immune cells, stromal cells, and endothelial cells that are part of the tumor-associated stroma, even if originating from normal tissue.

The robust clustering of cells by celltype_major, celltype_minor, and celltype_subset confirms the quality of the cell type annotations and the effectiveness of the dimensionality reduction method in capturing biological distinctions. The low number of unassigned cells across all annotation levels suggests comprehensive cell type identification within the dataset. The specific cell types identified, such as Luminal epithelial cell (a common subtype in breast cancer), various Macrophage polarizations (M1, M2A, M2B, M2C, M2D known for their roles in tumor immunity), and diverse T cell subsets (e.g., Treg, Cytotoxic, Th17), indicate the richness of the dataset for studying tumor immunology and stromal interactions.

A key finding from the ploidy_dec UMAP is the strong association of Aneuploid cells with the primary_tumor condition and, more specifically, with the Epithelial cell clusters. Given that Epithelial cell is designated as the tumor origin cell type in this breast tissue dataset, this observation is highly consistent with the known genomic instability and aneuploidy characteristic of epithelial-derived solid tumors like breast cancer. Aneuploidy, or the presence of an abnormal number of chromosomes, is a hallmark of cancer and contributes to tumor heterogeneity and progression PubMed search: aneuploidy cancer hallmarks. The detection of aneuploidy predominantly within the epithelial compartment reinforces these cells as the malignant population. The widespread presence of Diploid cells in other clusters further confirms the health and non-malignant nature of stromal, immune, and endothelial populations, which typically maintain a diploid state.

Annotation Notes

The consistency of clustering patterns across different annotation levels (major, minor, subset) and the clear biological distinctions observed (e.g., condition-specific cell populations, aneuploidy in epithelial tumor cells) suggest high-quality cell type annotation and data processing. The UMAPs effectively illustrate the cellular heterogeneity within the breast tissue samples and provide a reliable foundation for subsequent downstream analyses, such as differential gene expression or cell-cell interaction studies, within well-defined cell populations. The successful integration of samples indicates that comparisons across individuals can be made with confidence.

2. Major Cell Type Score and Ploidy Visualization on UMAP

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes single-cell RNA-seq data on a UMAP embedding, displaying major cell type scores, inferred ploidy status, and the assigned major cell type annotations. The goal is to understand the spatial distribution of different cell types and their ploidy state within the dataset, particularly in the context of breast tissue with primary tumor and normal conditions.

Visual Summary

The UMAP plots reveal several distinct clusters representing various cell populations.

Cell Type Score Distribution:

Biological Interpretation

The visualizations provide strong evidence for the distinct cellular composition of the breast tissue samples.

  1. Robust Cell Type Identification: The major cell type scores align very well with the assigned celltype_major annotations. For example, the regions with high "HiCAT_major_score: T cell" precisely match the T cell cluster in the celltype_major plot, and similarly for Epithelial, Stromal, Myeloid, and Endothelial cells. This consistency suggests that the cell type annotations are robust and well-supported by the underlying gene expression profiles.
  2. Identification of Tumor Cells: A critical observation is the strong correlation between Epithelial cell annotations and Aneuploid status. Given that "Tumor origin celltype: Epithelial cell" and the tissue is Breast (often implying breast carcinoma), the large Aneuploid cluster, predominantly composed of Epithelial cells, is highly indicative of the malignant tumor cell population. Aneuploidy is a hallmark of cancer, and its enrichment within the epithelial compartment strongly supports its identification as the tumor cells [1]. The normal epithelial cells, likely present in the normal samples, would be expected to be diploid and might cluster separately or as a minor diploid component within the epithelial population.
  3. Immune and Stromal Compartments: The UMAP clearly delineates distinct clusters for various immune cell types (T cells, Myeloid cells, B cells, Mast cells) and stromal cells (Fibroblasts, Smooth muscle cells, Endothelial cells, etc., derived from celltype_minor context). These non-malignant cells are predominantly Diploid, as expected for healthy host cells. Their spatial separation from the aneuploid epithelial tumor cells highlights the complex cellular ecosystem of the tumor microenvironment and surrounding normal tissue.

Annotation Notes

The consistency between the calculated major cell type scores and the assigned celltype_major labels is excellent, reinforcing the quality and accuracy of the cell type annotations in this dataset. The distinct clustering of different cell types and the clear segregation of aneuploid vs. diploid populations further validate the biological integrity of the data representation. The small number of "unassigned" cells or "unclear" ploidy cells do not significantly detract from the overall clarity of the major populations.

---

References:

  1. Aneuploidy as a hallmark of cancer:

PubMed search: "aneuploidy cancer hallmark"

3. Celltype_subset Marker Expression Analysis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a dot plot illustrating the expression patterns of identified marker genes across various celltype_subset populations derived from single-cell RNA-seq data of human breast tissue. The purpose is to visually confirm the distinct molecular identities of these cell subsets based on their unique marker gene expression profiles. The markers were selected to include surfaceome-only genes, which is beneficial for potential downstream experimental validation methods like flow cytometry.

Visual Summary

The dot plot displays celltype_subset populations on the y-axis and marker genes on the x-axis. Each dot represents the expression of a specific gene in a given cell subset. The size of the dot correlates with the fraction of cells within that group expressing the gene (non-zero expression), while the color intensity indicates the mean expression level of the gene in that group. On the right side, a bar chart shows the total number of cells for each celltype_subset.

The plot reveals a clear, diagonal-like pattern of highly expressed and specific marker genes for most annotated celltype_subset populations. This pattern, emphasized by the red boxes, indicates strong specificity for the selected markers and good separation between cell types. The majority of celltype_subset groups show distinct clusters of genes with high expression (large, dark red dots) largely restricted to their assigned cell type, with minimal off-target expression in other groups.

Biological Interpretation

The marker gene expression patterns strongly support the assigned biological identities of the various celltype_subset populations in the human breast tissue dataset:

Endothelial Cell Lineage

Fibroblasts

Fibroblasts exhibit a robust and specific signature defined by numerous extracellular matrix (ECM) components and related genes, including DCN (decorin), LUM (lumican), multiple COL genes (COL1A1, COL1A2, COL3A1, COL5A1, COL6A2), FAP (Fibroblast Activation Protein), FBLN1, FBLN2, PDGFRA, and LOXL1. FAP is particularly notable as a marker for activated fibroblasts, including cancer-associated fibroblasts (CAFs), which play critical roles in tumor microenvironment remodeling [PubMed Search: "FAP cancer-associated fibroblasts"].

Innate Lymphoid Cells (ILCs)

The different ILC subsets show distinct but related expression profiles:

Epithelial Cells

Macrophages

Macrophage subsets display varying expression of general and specific markers:

Mast Cells

Mast cells are clearly identified by canonical markers KIT (CD117), a receptor tyrosine kinase essential for mast cell development and survival, and tryptases TPSAB1 and TPSB2, which are specific mast cell proteases [GeneCards: KIT], [GeneCards: TPSAB1]. SRGN (serglycin) is also expressed, involved in granule storage [GeneCards: SRGN].

Smooth Muscle Cells

Smooth muscle cells are well-defined by key contractile proteins and associated genes, including ACTA2 (alpha-smooth muscle actin), CALD1 (caldesmon 1), TPM2 (tropomyosin 2), TAGLN (transgelin), MYL9 (myosin light chain 9), and MYH11 (myosin heavy chain 11) [GeneCards: ACTA2].

T Cell Lineage

The T cell subsets exhibit highly specific marker expression consistent with their functional roles:

Annotation Notes

The comprehensive marker gene expression patterns observed in the dot plot provide strong evidence supporting the quality and accuracy of the celltype_subset annotations within this AnnData object. The clear, non-overlapping expression of specific markers for the vast majority of cell types indicates that these subsets are transcriptionally distinct and well-resolved. The inclusion of surfaceome-only markers further strengthens the utility of these annotations for downstream applications requiring cell isolation or targeted imaging. This level of granularity and distinct marker expression suggests a reliable foundation for further biological investigations into the roles of these specific cell subsets in breast tissue, across different conditions (primary_tumor vs. normal), or in response to various perturbations.

4. Copy Number Variation Analysis in Breast Cancer Tumor and Unassigned Cells Grouped by Sample

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates copy number variations (CNVs) in single-cell RNA-seq data focusing on 'Epithelial cell' (the defined tumor origin cell type) and 'unassigned' cells from breast tissue samples. Cells were grouped by sample, and CNVs were estimated and visualized using a heatmap of log2(CNR) (Copy Number Ratio). A summary heatmap and bar chart of significantly amplified cytogenetic bands across these samples were also generated to highlight recurrent amplifications. The primary goal is to identify common chromosomal alterations in tumor-associated cells and distinguish them from samples identified as diploid.

Visual Summary

CNV Heatmap (log2(CNR))

The first heatmap visualizes log2(CNR) values across genomic spots for individual samples.

Summary of Significantly Amplified Regions

The second figure provides a detailed summary of significant amplifications by cytogenetic band.

1p34.3:1p34.2 (Frequency ~0.67)

1q21.1:1q23.2 (Frequency ~0.50)

1q32.1:1q32.1 (NFASC) (Frequency ~0.83)

1q41:1q42 (Frequency ~0.42)

7p22.2:7p21.1 (Frequency ~0.70)

7p11.2:7q21.1 (Frequency ~0.50)

7q21.3:7q22.1 (Frequency ~0.42)

8p11.21:8q12.1 (Frequency ~0.33)

8q22.1:8q23.1 (EIF3E) (Frequency ~0.42)

8q23.2:8q24.13 (Frequency ~0.50)

8q24.3:9p24.2 (Frequency ~0.70)

10q26.3:11p15.5 (Frequency ~0.50)

11q13.4:11q21 (Frequency ~0.50)

15q26.3:16p13.3 (Frequency ~0.50)

16p13.3:16p13.3 (Frequency ~0.50)

17q25.3:18p11.31 (Frequency ~0.33)

Biological Interpretation

The analysis of CNVs in 'Epithelial cell' and 'unassigned' populations provides critical insights into the genomic landscape of breast cancer.

Clinical or Translational Implications

5. CNV-based UMAP Visualization of Cell Type, Ploidy, Condition, and Sample

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes single-cell RNA-sequencing data on a Uniform Manifold Approximation and Projection (UMAP) embedding, specifically constructed using Copy Number Variation (CNV) estimates. The UMAP plots are colored by various metadata features: major cell type, minor cell type, ploidy status (ploidy_dec), condition (primary_tumor vs normal), and individual sample. This allows for an assessment of how these biological and technical factors contribute to the CNV landscape and cell population structure within the dataset.

Visual Summary

celltype_major and celltype_minor plots:

ploidy_dec plot:

condition plot:

sample plot:

Biological Interpretation

The UMAP embedding, specifically configured to emphasize CNV patterns, provides critical insights into the genomic landscape of the analyzed breast tissue samples.

  1. Tumor-specific Aneuploidy: The most striking observation is the clear segregation of cells based on their ploidy status, with "Aneuploid" cells forming distinct clusters that are almost exclusively derived from the "primary_tumor" condition. This aligns perfectly with the understanding that cancer cells, particularly epithelial cells (the designated Tumor origin celltype), often exhibit widespread chromosomal instability and aneuploidy as a hallmark of malignancy.
  2. Cell Type-Specific Genomic States: Epithelial cells, identified as the tumor origin cell type, predominantly co-localize with the aneuploid clusters and the primary tumor condition. In contrast, stromal cells, T cells, endothelial cells, and myeloid cells mainly occupy the diploid regions, consistent with their role as components of the tumor microenvironment or normal tissue, where they typically maintain a diploid genomic state.
  3. Data Quality and Annotation Validation: The robust separation of aneuploid/diploid cells and primary tumor/normal conditions on the CNV-based UMAP strongly validates the quality of the CNV estimates and the accuracy of the ploidy_dec and condition annotations. It also reinforces the biological relevance of the cell type assignments, especially distinguishing malignant epithelial cells from non-malignant cells.

Annotation Notes

The UMAP projections effectively capture significant biological variance related to copy number variations. The distinct clustering of cells based on ploidy, condition, and expected cell types (e.g., epithelial cells with aneuploidy in tumors) provides strong evidence for the quality of the CNV estimation and the biological integrity of the dataset. The minimal presence of "unassigned" cells and "unclear" ploidy status further supports the robustness of the current annotations.

6. Minor Cell Type Population Analysis in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a stacked bar plot showing the proportional distribution of minor cell types across individual samples for both normal breast tissue and primary breast tumors. The visualization allows for a direct comparison of cellular composition shifts associated with tumor development.

Visual Summary

The plot displays the relative proportions of 13 minor cell types (Dendritic cell, Endothelial cell, Epithelial cell, Fibroblast, ILC, Macrophage, Mast cell, NK cell, Smooth muscle cell, T cell CD4+, T cell CD8+, unassigned) for 4 normal samples and 10 primary tumor samples.

Key Observations:

Biological Interpretation

The observed shifts in minor cell type populations between normal breast tissue and primary tumors provide critical insights into the remodeling of the tumor microenvironment (TME) during breast cancer progression.

Clinical or Translational Implications

7. Immune Lymphoid Cell Subset Composition in Normal and Primary Breast Tumor Tissues

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a stacked bar plot visualizing the relative proportions of different lymphoid cell subsets, including T cells, Innate Lymphoid Cells (ILCs), and Natural Killer (NK) cells, across individual samples from normal breast tissue and primary breast tumors. The celltype_major category "T cell" was initially targeted, but the resulting celltype_subset breakdown broadly encompasses various lymphoid populations, reflecting a comprehensive view of the immune lymphoid compartment. Each bar represents 100% of the lymphoid cells identified within a given sample, allowing for a direct comparison of the internal composition of these immune cells between the two conditions.

Visual Summary

The visualization clearly separates samples into two conditions: 'normal' (4 samples) and 'primary_tumor' (10 samples).

Biological Interpretation

The observed shifts in lymphoid cell populations between normal breast tissue and primary tumors provide insights into the immune response dynamics within the tumor microenvironment (TME).

Clinical or Translational Implications

The compositional changes in the lymphoid compartment in primary breast cancer have significant clinical and translational implications:

8. Primary Breast Tumor Condition Shows Altered T Cell Subset Proportions

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates statistically significant differences in the proportions of T cell subsets between primary breast tumor tissue and normal breast tissue samples. The plot_box_for_celltype_population_with_signif_difference tool was used to compare cell type proportions (specifically at the 'subset' taxonomic level within the 'T cell' major population) between the 'primary_tumor' and 'normal' conditions, identifying subsets with a p-value less than 0.1 and a log2 fold change greater than 0.1.

Visual Summary

The boxplots illustrate the proportions of three T cell subsets that showed statistically significant differences between conditions:

Biological Interpretation

The observed shifts in T cell subset populations in the primary breast tumor microenvironment suggest a complex interplay of immune responses:

Collectively, these findings suggest that the breast tumor microenvironment is characterized by alterations in specific T cell and innate lymphoid populations, favoring an immune landscape that may suppress effective anti-tumor immunity (e.g., higher Th2) and potentially impair immune organization (e.g., lower LTI), while also showing an active but potentially complex B cell-supporting response (e.g., higher Tfh).

Clinical or Translational Implications

These differential proportions of T cell subsets between normal and primary tumor breast tissue have several potential clinical and translational implications:

9. Macrophage Cell Population Validation Across Samples

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the distribution of cells identified as 'Macrophage' (based on the celltype_minor annotation) across individual samples, separated by 'normal' and 'primary_tumor' conditions. The plot_celltype_population tool was used with a targets parameter specifically set to 'Macrophage', indicating a focus on this particular cell type.

Visual Summary

The bar plots display the selected 'Macrophage' population for each sample. Under the 'normal' condition, four samples are shown, and under the 'primary_tumor' condition, ten samples are displayed. In both conditions, for every single sample, the bar representing 'Macrophage' reaches 100% on the y-axis. This indicates that all cells considered for this specific visualization in each sample are annotated as 'Macrophage'.

Biological Interpretation

This visualization primarily serves as a confirmation of the presence and accurate selection of 'Macrophage' cells within the dataset. For each sample included in these plots, 100% of the cells being displayed are indeed identified as 'Macrophage' according to the celltype_minor annotation.

It is crucial to understand that this plot does not represent the overall proportion of macrophages relative to all other cell types within each tissue sample. Instead, it confirms that when the dataset is filtered or focused specifically on 'Macrophage' cells, those selected cells are consistently identified as such across all analyzed samples and conditions.

Biologically, the presence of macrophages in both normal breast tissue and primary tumor samples is highly expected. Macrophages are integral components of the immune system and play diverse roles in tissue homeostasis, inflammation, and disease, including cancer. In the context of breast cancer, tumor-associated macrophages (TAMs) are well-known to infiltrate the tumor microenvironment and can significantly influence tumor progression, metastasis, and response to therapy PubMed search: tumor associated macrophages breast cancer.

Annotation Notes / Limitations

While this plot validates the successful identification and selection of the 'Macrophage' population, it does not provide insights into key biological questions such as:

10. Ploidy 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) of epithelial and unassigned cells across various normal and primary breast tumor samples. Ploidy, the number of complete sets of chromosomes in a cell, is a critical feature often altered in cancer, with aneuploidy (abnormal chromosome number) being a hallmark of many solid tumors, including breast cancer. Given that epithelial cells are identified as the tumor origin cell type, examining their ploidy provides insight into genomic instability associated with tumorigenesis.

Visual Summary

The bar plots display the proportion of aneuploid, diploid, and unclear cells for epithelial and unassigned cell populations within each sample, segregated by 'normal' and 'primary_tumor' conditions.

Biological Interpretation

The observed shift from predominantly diploid cells in normal tissue to a high and variable prevalence of aneuploid cells in primary tumors is a strong indicator of neoplastic transformation and genomic instability, a hallmark of cancer progression.

Clinical or Translational Implications

The detection of aneuploidy in tumor-origin epithelial cells has several clinical implications:

---

References:

  1. Aneuploidy as a hallmark of cancer:

PubMed Search: Aneuploidy cancer hallmark

  1. Tumor Heterogeneity:

PubMed Search: Breast cancer tumor heterogeneity aneuploidy

  1. Prognostic value of aneuploidy:

PubMed Search: Aneuploidy breast cancer prognosis

  1. Chromosomal Instability in Cancer:

PubMed Search: Chromosomal instability cancer prognosis

11. 유방암 미세환경 내 세포-세포 상호작용 패턴 분석

Report figure

[Analysis Visualization Results]...

Analysis Overview

본 분석은 단일 세포 RNA 시퀀싱 데이터를 기반으로 AnnData 객체에서 추출된 세포-세포 상호작용(Cell-Cell Interaction, CCI) 패턴을 normal 및 primary_tumor 조건별로 시각화한 결과입니다. 특히, Epithelial cell, Fibroblast, Macrophage, T cell CD4+, T cell CD8+ 등 주요 세포 유형 간의 리간드-수용체 상호작용에 초점을 맞추었습니다. CellPhoneDB를 사용하여 계산된 상호작용 중 각 조건에서 최대 80개의 유의미한 상호작용이 점도표(dot plot) 형태로 제시되었으며, 점의 크기는 상호작용의 통계적 유의성(-log10(p-value))을, 색상은 상호작용 리간드 및 수용체의 평균 발현 강도(log2(mean))를 나타냅니다. 특히 primary_tumor 조건에서는 종양 세포의 특성을 반영하는 이수성(Aneuploid) 상피세포(Epi)가 포함되어 분석되었습니다.

Visual Summary

Normal 조건에서의 세포-세포 상호작용:

Primary Tumor 조건에서의 세포-세포 상호작용:

ECM 리모델링 및 종양 관련 신호:

Biological Interpretation

이러한 세포-세포 상호작용의 변화는 유방암 미세환경의 재구성을 명확히 보여줍니다.

  1. ECM 리모델링 및 기질 활성화: normal 조직에서는 콜라겐-인테그린 상호작용이 조직 항상성 유지에 기여하지만, primary_tumor에서는 암세포(Aneuploid Epi)가 섬유아세포와 함께 ECM 리모델링에 적극적으로 참여하여 종양 세포의 침윤 및 전이를 촉진하는 환경을 조성합니다. 이는 암 관련 섬유아세포(CAFs)의 활성화와 관련이 깊습니다.
  2. 면역 미세환경의 변화: primary_tumor에서 대식세포와 T 세포 간의 상호작용이 증가하고, CXCL12-CXCR4 축과 같은 면역 조절 관련 상호작용이 나타나는 것은 종양 면역 환경이 활성화되거나 면역 억제적 방향으로 재편되고 있음을 시사합니다. 특히 Mac|T CD4+와 같은 상호작용은 종양 관련 대식세포(TAMs)가 T 세포 반응을 조절하는 메커니즘을 반영할 수 있습니다.
  3. 종양 성장 및 진행 촉진 신호: primary_tumor에서 SPP1-integrin, CXCL12-CXCR4, VEGF-VEGFR, HGF-MET, IL6-IL6R, WNT 등 다양한 성장 인자 및 신호 전달 경로가 활발하게 나타나는 것은 종양 세포의 증식, 생존, 혈관신생 및 전이를 직접적으로 지원하는 핵심적인 메커니즘을 나타냅니다. 예를 들어, SPP1 (Osteopontin)은 종양 진행, 면역 억제 및 항암 치료 저항성에 중요한 역할을 하는 것으로 알려져 있습니다 PubMed search: SPP1 cancer therapy. CXCL12-CXCR4 축은 암세포의 이동과 전이, 그리고 면역 세포 유입에 관여합니다 GeneCards: CXCR4. VEGF는 혈관신생의 주요 조절인자입니다 PubMed search: VEGF angiogenesis cancer.
  4. 이수성 상피세포의 역할: Aneuploid Epi는 유방암의 주요 종양 세포 집단으로, 이들이 섬유아세포, 대식세포 및 다른 이수성 상피세포와 광범위하게 상호작용하는 것은 종양 세포 자체가 미세환경 조성 및 종양 진행의 주도적인 역할을 수행함을 보여줍니다.

Clinical or Translational Implications

이러한 세포-세포 상호작용 패턴의 차이는 유방암의 진단 및 치료 전략 개발에 중요한 시사점을 제공합니다.

  1. 치료 표적 발굴:
  1. 바이오마커 개발:
  1. 병용 요법의 근거:
  1. 실험적 검증:

12. Primary Breast Tumor Cell-Cell Interaction Analysis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns in breast tissue under 'normal' and 'primary_tumor' conditions, utilizing single-cell RNA sequencing data. The plot_cci_dots tool was used to visualize the most significant and strong ligand-receptor interactions between different cell types for each condition. The dot size represents the negative log10 of the p-value (significance), and the color intensity indicates the mean expression level (strength) of the interaction, with warmer colors (yellow/green) signifying higher mean expression. We focused on the top 80 interactions by default. This comparison highlights how cell communication networks are altered in the tumor microenvironment (TME) compared to normal tissue, providing insights into potential mechanisms of tumor progression and identifying therapeutic targets.

Visual Summary

CCI for Normal Condition

CCI for Primary Tumor Condition

Biological Interpretation

The striking differences in CCI between normal and primary tumor breast tissue reflect profound changes in the cellular ecosystem during tumorigenesis and progression.

  1. Tumor-Stromal Remodeling: The extensive interactions of Aneuploid Epithelial cells with Fibroblasts, Endothelial cells, and Macrophages highlight the critical role of the tumor microenvironment (TME) in supporting tumor growth and invasion. Fibroblasts are likely reprogrammed into cancer-associated fibroblasts (CAFs), which are known to secrete ECM components and growth factors that promote tumor cell proliferation, survival, and metastasis [GeneCards: FAP].
  2. Increased ECM-mediated Interactions: The prominent and often upregulated integrin-mediated interactions (e.g., COL1A1, COL1A2, COL3A1, COL6A1, FN1, FBN1, LAMC1 integrin complexes) in the tumor condition signify extensive remodeling of the extracellular matrix. This ECM provides physical support, signaling cues, and scaffolds for cell migration and invasion, critical processes in cancer progression. Integrin signaling can also promote tumor cell survival and therapeutic resistance [PubMed: Integrin signaling cancer].
  3. Angiogenesis and Vascular Remodeling: The strong and numerous interactions involving Endothelial cells with Aneuploid Epithelial cells and Fibroblasts, particularly through growth factor pathways like VEGFA-NRP1, PGF-FLT1, and ANGPT2-TEK, indicate active angiogenesis. This process is essential for supplying nutrients and oxygen to the rapidly growing tumor, and for facilitating metastasis [PubMed: Angiogenesis cancer].
  4. Immune Cell Recruitment and Modulation: Enhanced chemokine signaling, specifically CXCL12-CXCR4 and CCL2-CCR2, suggests active recruitment of various immune cells, including T cells and macrophages, into the TME. While CXCL12-CXCR4 can attract both pro-tumorigenic cells (e.g., myeloid-derived suppressor cells, some regulatory T cells) and anti-tumorigenic cells, its overexpression is often associated with poor prognosis and metastasis in breast cancer [GeneCards: CXCL12]. Macrophages, particularly M2-like subsets observed in the TME, can switch their phenotype to become pro-tumorigenic, promoting angiogenesis, immune suppression, and metastasis [PubMed: Tumor associated macrophages].
  5. Ploidy-specific Interactions: The clear distinction between Diploid Epithelial cell and Aneuploid Epithelial cell interactions in the tumor plot underscores that the malignant, aneuploid cells are the primary drivers of these altered communication networks. This confirms the specific involvement of transformed cells in shaping the tumor microenvironment.

Clinical or Translational Implications

The identified differential cell-cell interactions offer compelling insights for therapeutic development and patient stratification in breast cancer:

  1. Targeting Integrin-mediated Adhesion: The upregulation of various integrin complexes in the tumor, particularly those involving Aneuploid Epithelial cells and Fibroblasts, suggests that targeting specific integrin receptors could disrupt tumor cell adhesion, migration, and invasion, potentially hindering metastasis. Small molecule inhibitors or antibodies against specific integrins (e.g., αvβ3, α5β1) are under investigation in cancer therapies [PubMed: Integrin inhibitors cancer].
  2. Anti-angiogenic Strategies: The strong activation of VEGFA-NRP1, PGF-FLT1, and ANGPT2-TEK pathways confirms the importance of angiogenesis in breast cancer. Existing anti-angiogenic drugs, such as bevacizumab (targeting VEGFA), could be further optimized or combined with other therapies, or novel agents targeting NRP1, FLT1, or TEK could be explored [PubMed: Anti-angiogenic therapy breast cancer].
  3. Modulating Chemokine Axes: The CXCL12-CXCR4 and CCL2-CCR2 axes are critical for immune cell trafficking and can promote an immunosuppressive TME. Disrupting these axes, for example, with CXCR4 antagonists or CCR2 inhibitors, could impede the recruitment of pro-tumorigenic immune cells (e.g., macrophages, Tregs) and enhance anti-tumor immune responses [PubMed: Chemokine receptor inhibitors cancer].
  4. Targeting CAF-mediated Interactions: The prominent role of Fibroblasts in tumor interactions suggests that therapies aimed at reprogramming CAFs or disrupting their communication with tumor cells could be beneficial. This might involve targeting specific CAF-secreted factors or receptors involved in their activation [PubMed: CAF targeted therapies].
  5. Biomarker Identification: Specific ligand-receptor pairs with high activity in the primary tumor, especially those involving Aneuploid Epithelial cells, could serve as diagnostic or prognostic biomarkers, or predictors of response to targeted therapies. Further validation studies are warranted to explore their clinical utility.

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) mediated by a curated set of genes relevant to immune checkpoint and cell cycle pathways. Using CellPhoneDB results derived from single-cell RNA-seq data of breast tissue, we compare interaction patterns between normal tissue and primary tumor conditions. The plot_cci_dots tool visualizes significant ligand-receptor pairs, their interaction strength (mean expression), and statistical significance (p-value) across various interacting cell types. This allows for a deeper understanding of altered cellular communication in the tumor microenvironment.

Visual Summary

The provided dot plots illustrate significant cell-cell interactions for selected ligand-receptor pairs in both "normal" and "primary_tumor" conditions.

Biological Interpretation

The analysis highlights critical changes in cellular communication pathways involving growth factors and their receptors, which are instrumental in both cell cycle regulation and immune modulation, between normal breast tissue and primary tumors.

  1. Dominance of TGF-beta Signaling in Both Normal and Tumor Contexts, with Specific Tumor-Associated Alterations:
  1. Elevated EGFR Signaling in Tumor-Associated Stromal Cells:
  1. Specific Cell-type Interactions in the Tumor Microenvironment:
  1. Implications for Immune Checkpoint and Cell Cycle Control:

Clinical or Translational Implications

The observed patterns of cell-cell interactions, particularly the differential engagement of TGF-beta and EGFR signaling in the primary tumor, present several clinical and translational opportunities:

  1. Therapeutic Target Prioritization:
  1. Biomarker Discovery and Patient Stratification:
  1. Rational Combination Therapies:
  1. Experimental Validation:

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates statistically significant differences in cell-cell interactions (CCI) between 'normal' and 'primary_tumor' conditions in human breast tissue, focusing on interactions involving T cells, Myeloid cells, Stromal cells (Fibroblast, Smooth muscle cell), Endothelial cells, and Mast cells. The results are visualized as a dot plot, highlighting the top 25 most significant differential interactions for each condition.

Visual Summary

The dot plot effectively illustrates distinct patterns of cell-cell communication between normal and primary tumor breast tissue samples.

Biological Interpretation

The contrasting CCI profiles underscore fundamental biological differences between healthy and cancerous breast tissue microenvironments.

Clinical or Translational Implications

These findings point to specific cell-cell interaction axes that could serve as potential therapeutic targets or prognostic markers in breast cancer.

References

  1. Prostaglandin E2: PubMed search for "Prostaglandin E2 cancer immune regulation" https://pubmed.ncbi.nlm.nih.gov/?term=Prostaglandin+E2+cancer+immune+regulation
  2. SIRPA-CD47: GeneCards entry for CD47 https://www.genecards.org/cgi-bin/carddisp.pl?gene=CD47
  3. CCL2-ACKR1: PubMed search for "CCL2 ACKR1 cancer" https://pubmed.ncbi.nlm.nih.gov/?term=CCL2+ACKR1+cancer
  4. TNFSF10 (TRAIL): UniProt entry for TRAIL https://www.uniprot.org/uniprotkb/O14790/entry
  5. Collagen in Cancer: PubMed search for "collagen tumor microenvironment" https://pubmed.ncbi.nlm.nih.gov/?term=collagen+tumor+microenvironment
  6. Fibronectin in Cancer: GeneCards entry for FN1 https://www.genecards.org/cgi-bin/carddisp.pl?gene=FN1
  7. Integrin Inhibitors: PubMed search for "integrin inhibitors cancer therapy" https://pubmed.ncbi.nlm.nih.gov/?term=integrin+inhibitors+cancer+therapy

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

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates condition-specific surfaceome markers within Epithelial cells, identified as the tumor-origin cell type, from single-cell RNA-seq data of human breast tissue. The aim is to distinguish between normal and primary tumor conditions based on surface protein expression patterns. The visualization presents genes expressed on the cell surface, showing both the mean expression level and the fraction of cells expressing each gene across different patient samples, grouped by condition (normal vs. primary tumor) and ploidy status (diploid vs. likely aneuploid).

Visual Summary

The dot plot effectively illustrates clear differences in surfaceome marker expression between normal and primary tumor epithelial cells:

Biological Interpretation

The identified surfaceome markers provide crucial insights into the biological changes occurring in breast epithelial cells during tumorigenesis.

Clinical or Translational Implications

The identified condition-specific surfaceome markers in breast epithelial cells hold significant promise for clinical applications. Given their surface localization, these proteins are highly accessible for diagnostic, prognostic, and therapeutic interventions.

Biomarker Development

This analysis provides a robust foundation for further functional studies and validation of these promising surfaceome markers as potential diagnostic tools and therapeutic targets in 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 for Macrophage cells in breast tissue, comparing normal tissue with primary tumors. The dot plot displays the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) across individual patient samples, grouped by condition. The focus is on surfaceome proteins, which are particularly relevant for cell-surface-mediated interactions and potential therapeutic targeting.

Visual Summary

The dot plot effectively distinguishes macrophage surfaceome profiles between normal and primary tumor conditions, as well as highlighting heterogeneity within primary tumors.

Biological Interpretation

The observed differential expression of surfaceome markers points to a significant phenotypic shift in macrophages within the primary breast tumor microenvironment compared to normal breast tissue.

The marked contrast between the normal and tumor-associated macrophage surfaceome profiles highlights the profound reprogramming macrophages undergo in the presence of cancer, adapting functions that often support tumor progression.

Clinical or Translational Implications

The identification of condition-specific surfaceome markers on macrophages in breast cancer has several important clinical and translational implications:

These targets offer avenues for developing novel macrophage-targeted immunotherapies, such as antibodies or small molecules, to re-educate pro-tumorigenic macrophages or deplete them from the tumor microenvironment.

17. Fibroblast Condition-Specific Surfaceome Markers in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis identifies condition-specific surfaceome markers for Fibroblasts, comparing normal breast tissue with primary breast tumor samples. The dot plot visualizes the expression of these markers across individual patient samples, indicating both the mean expression level (color intensity) and the fraction of cells expressing the marker (dot size) within the Fibroblast population. This approach aims to uncover surface-accessible proteins that distinguish tumor-associated fibroblasts from those in normal tissue.

Visual Summary

The dot plot clearly segregates fibroblast populations based on their gene expression profiles, aligning with the "normal" and "primary_tumor" conditions.

Biological Interpretation

The observed upregulation of specific surfaceome markers in primary tumor fibroblasts points to their transformation into Cancer-Associated Fibroblasts (CAFs), a critical component of the tumor microenvironment (TME). CAFs play diverse roles in promoting tumor growth, invasion, metastasis, angiogenesis, and immune suppression.

Notable markers and their biological roles:

Clinical or Translational Implications

The identification of these highly expressed surfaceome markers on breast cancer-associated fibroblasts holds significant clinical and translational potential:

18. Cell-Type-Specific Surfaceome Marker Discovery in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify condition-specific surfaceome markers for T cell CD4+ populations and other cell types within breast tissue from single-cell RNA sequencing data. The plot_markers_and_expression_dot tool was used, configured to find up to 50 surfaceome markers. While the user query specifically asked for "condition-specific" markers for T cell CD4+, the generated dot plot displays markers that are highly *cell-type-specific* across various immune and stromal cell subsets (e.g., Endothelial cells, Fibroblasts, Macrophages, different T cell subsets, etc.). This suggests the marker discovery focused on distinguishing cell types from each other rather than comparing conditions (primary tumor vs. normal) within a single cell type. The plot visualizes the mean expression level (color intensity) and the fraction of cells expressing each marker (dot size) for each cell type.

Visual Summary

The dot plot effectively illustrates the distinct expression patterns of various surfaceome-enriched genes across 25 different cell subsets found in the breast tissue.

Biological Interpretation

The analysis revealed distinct surfaceome marker profiles for various cell types in the breast microenvironment, including several CD4+ T cell subsets.

Beyond T cells, the plot also effectively identified markers for other key cell types:

The inclusion of several intracellular markers (e.g., FOXP3, STAT1, GATA3, RORC, GZMB) in a "surfaceome only" analysis suggests that either the definition of surfaceome was broad or these highly canonical intracellular markers were specifically allowed in the marker finding process to ensure robust cell type identification.

Clinical or Translational Implications

The identification of cell-type-specific surfaceome markers has significant implications for both understanding breast biology and developing clinical strategies.

19. Differential Expression of Cell Cycle-Related Genes in Breast Cancer Epithelial Cells: YWHAZ Upregulation

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis aimed to identify and visualize statistically significant expression differences of selected cell cycle-related genes within specific cell types (Epithelial cell, Fibroblast, Macrophage) between 'primary_tumor' and 'normal' breast tissue conditions. The provided boxplot specifically illustrates the differential expression of the YWHAZ gene in Epithelial cells, which are the tumor-origin cells in breast tissue.

Visual Summary

The boxplot for the gene YWHAZ in Epithelial cells reveals a marked difference in its expression levels between the two conditions:

Biological Interpretation

The significant upregulation of YWHAZ (Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta) in Epithelial cells within the primary tumor environment is a key finding, especially considering that Epithelial cells are the cells of tumor origin in breast cancer.

Clinical or Translational Implications

20. Epithelial Cell Gene Ontology (GSA) Analysis: Condition and Ploidy-Specific Pathway Enrichment in Breast Tissue

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis utilizes Gene Ontology (GSA) to identify enriched biological pathways and processes in epithelial cells within the provided single-cell RNA-seq dataset. The goal is to understand how pathway activity differs across various cellular states, specifically comparing diploid versus other epithelial cells, normal versus other epithelial cells, and primary tumor versus other epithelial cells. The results are visualized as bar plots, showing the statistical significance (-log(p-val) and -log(q-val)) of enriched terms.

Visual Summary

The visualizations present three bar plots, each detailing the most significantly enriched Gene Ontology (GO) terms and other pathway database terms (e.g., from KEGG) for epithelial cells under different comparison contexts:

  1. Diploid_vs_others: Pathways enriched in diploid epithelial cells compared to all other (presumably aneuploid) epithelial cells.
  2. normal_vs_others: Pathways enriched in normal epithelial cells compared to all other (presumably primary tumor) epithelial cells.
  3. primary_tumor_vs_others: Pathways enriched in primary tumor epithelial cells compared to all other (presumably normal) epithelial cells.

Each plot displays the top enriched "Term" on the y-axis, ranked by statistical significance, with corresponding -log(p-val) and -log(q-val) on the x-axis. A higher value indicates greater statistical significance. The primary_tumor_vs_others comparison shows the highest significance levels for enriched terms, followed by normal_vs_others, and then Diploid_vs_others.

Biological Interpretation

Diploid Epithelial Cells (vs. others)

Epithelial cells maintaining diploidy, often representing a more stable genomic state, exhibit enrichment in pathways that suggest a baseline of cellular function and interaction with the microenvironment:

Normal Epithelial Cells (vs. others)

Normal epithelial cells, when compared to other epithelial cells (likely predominantly tumor cells), show distinct pathway enrichments:

Primary Tumor Epithelial Cells (vs. others)

Epithelial cells from primary tumors display a pronounced shift in their active pathways, reflecting cancer-associated biology:

Clinical or Translational Implications

The GSA results provide valuable insights into the biological underpinnings of breast cancer and offer potential avenues for therapeutic exploration:

21. Gene Set Enrichment Analysis Reveals Widespread Pathway Dysregulation in Breast Cancer Microenvironment

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results for six key cell types (Endothelial cell, Epithelial cell, Fibroblast, ILC, Macrophage, and Smooth muscle cell) within the breast tissue context. The comparisons are primarily between primary tumor and normal conditions for each cell type, with an additional comparison for Epithelial cells based on ploidy status (Diploid vs. others, likely Aneuploid). The dot plot visualizes the Normalized Enrichment Score (NES) and statistical significance (-log(P)) for a curated set of 80 gene sets, providing insight into activated (red dots) or suppressed (blue dots) biological pathways.

Visual Summary

The dot plot effectively illustrates distinct patterns of pathway enrichment and depletion across different cell types and conditions.

Biological Interpretation

Widespread Activation of Proliferation, Metabolic, and Signaling Pathways in the Tumor Microenvironment

Across Epithelial cells, Fibroblasts, Macrophages, Endothelial cells, and Smooth muscle cells from primary tumors, a consistent and significant enrichment (large red dots) is observed for pathways central to cancer progression:

Differential Pathway Activity in Aneuploid vs. Diploid Epithelial Cells

The comparison of Epithelial cell: Diploid_vs_others (likely Diploid vs. Aneuploid epithelial cells) provides critical insights into the molecular characteristics linked to ploidy status:

Cell-Type Specific Contributions to the Tumor Microenvironment

Clinical or Translational Implications

The pervasive activation of proliferation, metabolic, and key signaling pathways across multiple cell types in the breast tumor microenvironment has significant clinical implications:

This comprehensive GSEA provides a detailed molecular fingerprint of the breast tumor microenvironment, emphasizing the interconnectedness of different cell types in driving disease progression and identifying potential common vulnerabilities for therapeutic exploitation.

22. Discussion

The comprehensive single-cell analysis reveals a striking and multifaceted transformation from normal breast tissue to primary breast cancer, encompassing cellular composition, genomic stability, intercellular communication, and functional pathway activity. UMAP visualizations consistently delineate distinct cellular landscapes for normal and primary tumor conditions, with aneuploid epithelial cells being a hallmark feature predominantly confined to the tumor compartment. This genomic instability in epithelial cells, designated as the tumor origin, is further substantiated by extensive recurrent copy number variations (CNVs) in primary tumor samples, notably amplifications on chromosomes 1q, 7p/q, and 8q, implicating genes like NFASC and EIF3E in tumor progression.

The tumor microenvironment (TME) undergoes profound remodeling, evidenced by significant shifts in cell type proportions. While normal tissue is dominated by epithelial cells, primary tumors show a reduced epithelial fraction alongside a notable increase in stromal cells, particularly cancer-associated fibroblasts (CAFs), and diverse immune populations, including macrophages (tumor-associated macrophages, TAMs) and various T cell subsets. Functional analysis of these infiltrating immune cells indicates an altered immune landscape in tumors, with increased proportions of Th2 and T follicular helper (Tfh) cells and a decrease in Lymphoid Tissue inducer (LTI) cells, collectively suggesting a shift towards immunosuppression and potential disruption of organized lymphoid structures, while also reflecting active B cell responses. Macrophages within the TME also exhibit a distinct surfaceome profile, characterized by upregulation of pro-tumorigenic markers such as EREG, CXCR4, and MET, distinguishing them from normal tissue macrophages.

Cell-cell interaction (CCI) analyses unveil a dramatically denser and more complex communication network in primary tumors compared to normal tissue. Aneuploid epithelial cells actively engage in extensive crosstalk with CAFs, endothelial cells, and TAMs, underscoring their role in shaping the TME. A central finding is the widespread upregulation of extracellular matrix (ECM)-integrin interactions, involving various collagens and fibronectin, which signifies extensive ECM remodeling crucial for tumor cell adhesion, migration, and invasion. Angiogenic pathways (e.g., VEGF-VEGFR, PGF-FLT1) are highly activated in tumor endothelial cells, supporting neovascularization. Chemokine signaling, particularly the CXCL12-CXCR4 and CCL2-CCR2 axes, is enhanced, facilitating the recruitment of immune and stromal cells. Furthermore, TGF-beta signaling is pervasively activated in the primary tumor, notably through integrin αvβ6-mediated activation of latent TGF-beta, contributing to fibrosis, immune suppression, and epithelial-stromal crosstalk. In contrast, normal tissue CCIs are more sparse and reflect homeostatic processes and baseline immune surveillance, such as ProstaglandinE2 and SIRPA-CD47 signaling.

Functional pathway analyses (GSEA and GSA) corroborate these findings, revealing widespread pathway dysregulation across tumor-associated cell types. Primary tumor epithelial cells, CAFs, and TAMs consistently exhibit significant enrichment in pathways related to cellular proliferation (DNA replication, ribosome biogenesis), metabolic reprogramming (oxidative phosphorylation, purine/pyrimidine metabolism, glutathione metabolism), and key signaling cascades (MAPK, Hippo, mTOR, TGF-beta). This pervasive metabolic and signaling rewiring supports aggressive tumor growth and adaptation to the challenging TME. Notably, these pro-oncogenic programs are predominantly driven by the aneuploid epithelial cell population, as evidenced by their depletion in diploid epithelial cells. The significant upregulation of YWHAZ, a cell cycle regulator, in tumor epithelial cells further underscores the heightened proliferative drive. An intriguing and recurrent finding from Gene Ontology analysis is the enrichment of pathways associated with neurodegenerative diseases in both normal and tumor epithelial cells. While unexpected, this may point to shared fundamental cellular stress responses, mitochondrial dysfunction, or protein aggregation mechanisms that are potentially exacerbated in cancer. Collectively, these results highlight the intricate interplay between malignant cells and their microenvironment, driven by genomic instability, extensive cellular communication, and metabolic adaptation, which underpins breast cancer progression.

Hypotheses:

  1. Aneuploidy in breast epithelial cells is a primary driver of the observed tumor-specific transcriptional programs, cell-cell interaction patterns, and metabolic dysregulation within the tumor microenvironment.
  2. The extensive extracellular matrix remodeling and increased integrin-mediated interactions in primary breast tumors create a pro-tumorigenic niche that promotes cancer cell invasion, metastasis, and therapeutic resistance.
  3. The altered T cell subset composition (increased Th2 and Tfh, decreased LTI) and phenotypic shift of tumor-associated macrophages contribute to an immunosuppressive microenvironment, facilitating immune evasion and tumor progression.
  4. Hyperactivated TGF-beta signaling, particularly through integrin αvβ6 activation, plays a central role in promoting fibrosis, immune suppression, and epithelial-stromal crosstalk in breast cancer.
  5. The upregulation of YWHAZ in primary tumor epithelial cells contributes to uncontrolled cell cycle progression and survival, serving as a critical hub in the malignant transformation.

Potential therapeutic targets:

  1. TGF-beta Pathway (Ligands/Receptors/Activation via Integrin αvβ6): Consistently identified as highly active and pervasive in the primary tumor microenvironment, especially involving aneuploid epithelial cells and fibroblasts. Integrin αvβ6-mediated activation of latent TGF-beta is also prominent, indicating robust local signaling that promotes tumor growth, immune evasion, and fibrosis. Evidence: CCI analyses (Sections 11, 12, 13) show strong TGFB1-TGFbeta_receptor1/2 and TGFB1_integrin_avB6_complex interactions in primary tumor. GSA and GSEA (Sections 20, 21) also show enrichment of TGF-beta signaling pathway. Validation: Test TGF-beta inhibitors or αvβ6-integrin blocking antibodies in *in vitro* co-culture models of tumor cells and CAFs, or *in vivo* PDX models. Assess impact on tumor growth, invasion, fibrosis, and immune infiltration.
  2. Integrin-mediated ECM interactions (e.g., Collagen-Integrin, Fibronectin-Integrin complexes): Extensive upregulation of various collagen-integrin and fibronectin-integrin interactions in primary tumors, often involving aneuploid epithelial cells and fibroblasts. These interactions are critical for ECM remodeling, cell adhesion, migration, and invasion, central to metastasis and tumor progression. Evidence: CCI analyses (Sections 11, 12, 14) show robust upregulation of COL6A1/A2/A3/18A1-integrin_a1b1/a2b1/a3b1/avb5/a5b1_complex and FN1-integrin_a5b1_complex in primary tumors. Validation: Develop or repurpose specific integrin inhibitors (e.g., targeting α1β1, α2β1, α5β1, αVβ5) and test their efficacy in blocking tumor cell invasion, migration, and metastasis in *in vitro* assays and *in vivo* models. Assess effects on tumor stiffness and ECM architecture.
  3. NECTIN4 (on Epithelial Cells): NECTIN4 is a surface protein highly and specifically upregulated on primary tumor epithelial cells, serving as a distinct tumor-associated antigen. It is implicated in increased invasiveness and its surface accessibility makes it an ideal therapeutic target. Evidence: Section 15 identifies NECTIN4 as a highly expressed and prevalent condition-specific surfaceome marker in primary tumor epithelial cells (dark red dots, large sizes). Validation: Evaluate NECTIN4-targeting antibody-drug conjugates (ADCs) or monoclonal antibodies in preclinical breast cancer models. Assess their ability to inhibit tumor growth and metastasis. Confirm NECTIN4 expression by immunohistochemistry in patient cohorts.
  4. CXCR4 (on Macrophages and Fibroblasts): CXCR4 is a chemokine receptor upregulated in tumor-associated macrophages and involved in significant cell-cell interactions (e.g., CXCL12-CXCR4) with fibroblasts and T cells in primary tumors. It plays roles in tumor growth, metastasis, and immune cell trafficking, contributing to an immunosuppressive TME. Evidence: CCI analyses (Sections 11, 12) show robust CXCL12-CXCR4 interactions in primary tumor. Section 16 identifies CXCR4 as a condition-specific surfaceome marker for macrophages in tumors. Validation: Use CXCR4 antagonists to evaluate their effect on TAM recruitment to the tumor, macrophage polarization, tumor growth, and metastasis in *in vivo* models. Combine with immune checkpoint inhibitors to assess synergistic effects.
  5. YWHAZ (on Epithelial Cells): YWHAZ, a member of the 14-3-3 protein family, is a key cell cycle regulator significantly upregulated in primary tumor epithelial cells (tumor origin cell type), indicative of its role in driving uncontrolled proliferation and survival, hallmarks of cancer. Evidence: Section 19 presents boxplots showing significantly elevated YWHAZ expression in primary tumor epithelial cells compared to normal (p ≤ 0.001). Validation: Perform *in vitro* functional studies using siRNA/CRISPR to deplete YWHAZ in breast cancer cell lines and assess effects on cell proliferation, apoptosis, and cell cycle progression. Validate *in vivo* in PDX models.

Follow-up validation ideas:

  1. Validate recurrent CNV regions (e.g., 1q, 8q amplifications, specific genes like NFASC, EIF3E) in larger cohorts of breast cancer patient samples using FISH or array CGH, correlating with patient outcomes and specific tumor subtypes.
  2. Apply spatial transcriptomics or multiplex immunofluorescence/proteomics to breast tumor tissues to validate the spatial localization of specific cell types (e.g., FAP+ CAFs, CD163+ TAMs, NECTIN4+ epithelial cells) and key ligand-receptor interactions (e.g., SPP1-integrin, CXCL12-CXCR4, TGFB1-integrin_avB6_complex) identified by CCI analysis.
  3. Quantify and functionally characterize T cell subsets (Th2, LTI, Tfh) and macrophage polarizations (M1/M2 markers like CD163, EREG, CXCR4, MET) in fresh tumor biopsies and matched normal tissue using flow cytometry or high-plex immunostaining, correlating with clinical outcomes.
  4. Conduct perturbation assays (in vitro/in vivo) such as co-culture experiments with patient-derived tumor epithelial cells, CAFs, and TAMs to functionally validate key CCIs (e.g., blocking TGF-beta signaling, integrin antagonists, CXCR4 inhibitors) and their impact on tumor cell proliferation, migration, and immune cell function. Use in vivo patient-derived xenograft (PDX) or syngeneic mouse models with gene editing (CRISPR/shRNA) to target upregulated genes (e.g., YWHAZ, FAP, NECTIN4, EGFR, ERBB3) in tumor or stromal cells, followed by assessing tumor growth, metastasis, and immune infiltration.
  5. Perform metabolic flux analysis on isolated tumor epithelial cells, CAFs, and TAMs to validate the predicted metabolic reprogramming (e.g., oxidative phosphorylation, glutathione metabolism) and identify specific metabolic vulnerabilities that can be therapeutically targeted.

Limitations:

This report is based on single-cell RNA sequencing data, providing transcriptomic insights which infer cellular function and interactions; direct functional consequences require experimental validation. CNV estimates derived from RNA-seq are inferential and would benefit from orthogonal genomic validation methods. Cell-cell interaction predictions from CellPhoneDB suggest potential communication but require experimental confirmation of physical contact and functional impact. The inherent heterogeneity among patient samples and the specific dataset size may limit the generalizability of some findings. While minor, 'unassigned' cell populations could harbor uncharacterized cell states. Lastly, observational data cannot establish causality, and further mechanistic studies are necessary to elucidate the precise roles of identified targets and pathways in breast cancer progression.

23. 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 major cell type scores on UMAP 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. Select tumor origin cells and unassigned cells, group them by sample, show a CNV heatmap, and also show a summary of significantly amplified copy number regions, and save.
  5. Show CNV patterns on UMAP, including major cell type, minor cell type, ploidy results, condition, and sample in 2 columns, and 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. If there are statistically significant differences between conditions in T cell subset populations, show boxplots and save. Set ncols appropriately based on the total number of panels.
  9. Show a subset population bar plot for Macrophages and save.
  10. Select tumor origin cells and unassigned cells, show their ploidy population as a bar plot, and save.
  11. Show cell-cell interaction patterns including Epithelial cell, Fibroblast, Macrophage, and T cell by condition, and save. For cell-cell interactions, select up to 80 interactions per condition.
  12. Show cell-cell interactions by condition and save the result. Select at most 80 cell-cell interactions for each condition.
  13. Select only genes related to immune checkpoint and cell cycle pathways, show cell-cell interactions for these genes, and save.
  14. Find statistically significant differences in cell-cell interactions between conditions for T cell, Myeloid cell, Stromal cell, Endothelial cell, and Mast cell, show them as a dot plot, and save. Set max_n_items_per_group = 25.
  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 Macrophage, show them as a dot plot, and save. Only surfaceome markers, up to 50 per condition.
  17. Extract condition-specific markers for Fibroblast, show them as a dot plot, and save. Only surfaceome markers, up to 50 per condition.
  18. Extract condition-specific markers for T cell CD4+, show them as a dot plot, and save. Only surfaceome markers, up to 50 per condition.
  19. For genes related to the Cell cycle pathway, for Epithelial cell, Fibroblast, and Macrophage, show boxplots of statistically significant expression differences between conditions, and save. Set max_n_items_to_plot = 24, and ncols appropriately so that the width:height ratio is about 2:3 based on the total number of panels.
  20. Show Gene Ontology (GSA) analysis results for epithelial cells as bar plots and save the result.
  21. Show a dot plot of Gene set enrichment analysis results for Endothelial cell, Epithelial cell, Fibroblast, ILC, Macrophage, and Smooth muscle cell, and save. Use RdBu_r as the color map, and set n_pws_to_show = 80.
Top