Cancer stem cells orchestrate immune evasion through extracellular vesicle-mediated non-canonical signaling pathways
Data files
Mar 24, 2026 version files 202.76 MB
-
Figure1B_raw_data-CD3_T_tsne.acs
2.21 MB
-
Figure1B_raw_data-concat_T8_high_T8_low_CD3_T_download_16000_cells.downloading
2.31 MB
-
Figure1B_raw_data-T8_high_HSY_1156553_TUMOR_raw_data_010.fcs
54.18 MB
-
Figure1B_raw_data-T8_Low_GCH_1071022_TUMOR_raw_data_013.fcs
51.04 MB
-
Figure2B_raw_data-20250804_Figure_2B_tsne.acs
2.79 MB
-
Figure2B_raw_data-Ts-KD1.fcs
2.10 MB
-
Figure2B_raw_data-Ts-KD2.fcs
1.99 MB
-
Figure2B_raw_data-Ts-KD3.fcs
17.21 MB
-
Figure2B_raw_data-Ts-NC1.fcs
2.11 MB
-
Figure2B_raw_data-Ts-NC2.fcs
33.05 MB
-
Figure2B_raw_data-Ts-NC3.fcs
32.95 MB
-
Figure2G_raw_data-17-Aug-2025.acs
273.79 KB
-
Figure2G_raw_data-concat2_1.fcs
267.82 KB
-
Processed_data.pdf
264.44 KB
-
README.md
5.91 KB
Abstract
Cancer stem cells (CSCs) within tumors exhibit a remarkable capacity to evade immune surveillance, thereby escaping elimination by the host immune system, yet the molecular mechanisms underlying this immune-evasive plasticity remain incompletely defined. Compared to other breast cancer subtypes, triple-negative breast cancer (TNBC) is particularly enriched with CSCs exhibiting basal-like molecular features characterized by elevated stemness indices, positioning it as an exemplary model for investigating bidirectional CSCs-TME crosstalk. Here, we mapped the spatial immunometabolic landscape of TNBC using high-dimensional flow cytometry to simultaneously interrogate 50 TME-associated proteins at single-cell resolution, including CSC markers tetraspanin-8 (TSPAN8) and CD44. A total of 161 treatment-naïve TNBC specimens were analyzed, stratified into CSC-high (CSCHi) and CSC-low (CSCLo) cohorts based on TSPAN8 and CD44 protein expression thresholds. TNBCs with elevated stemness indices drive effector-to-regulatory T cell (Treg) phenotypic conversion, establishing immunosuppressive niches. Blocking EV-associated TSPAN8 using an anti-TSPAN8 monoclonal antibody synergized with PD-1 checkpoint inhibition, providing preclinical validation for dual targeting of CSC plasticity and immune checkpoint pathways. The raw and processed data from high-dimensional single-cell proteomics analysis performed on the BD FACSymphony Flow Cytometry system, as shown in Figures 1B, 2B, and 2G, have been deposited.
Dataset DOI: 10.5061/dryad.05qfttfhx
Description of the data and file structure
To characterize the local breast cancer tumor microenvironment (BC-TME), we performed high-dimensional single-cell proteomic profiling on 161 fresh tumor specimens obtained from 65 treatment-naïve patients with triple-negative breast cancer (TNBC). A 50-parameter antibody panel was designed to assess TME-associated proteins, including major lineage markers for distinguishing immune cell subsets (e.g., CD4, CD8, TCRγδ, CD69, CD19, CD14, CD15, CD16, CD203, CD206, CD56, CD11b, CD66b, CD123, CD11c, HLA-DR), immune regulatory molecules (e.g., PD-1, CTLA-4, CD25, CD127, TIM-3, CD204, CD47), epithelial markers (e.g., HER2, Mucin-1), stromal markers (e.g., FAP), and cancer stem cell (CSC) markers (TSPAN8, CD44, and CD24). Specimens were analyzed by spectral flow cytometry and stratified into CSC-high (CSCHigh) and CSC-low (CSCLow) groups based on median TSPAN8/CD44 co-expression levels within the CD45-FAP-CD24- population. High-dimensional data were processed using nonlinear dimensionality reduction (t-SNE) and unsupervised PhenoGraph clustering to resolve the heterogeneity of tumor-infiltrating leukocytes in both individual and aggregated analyses (Figure 1B in Processed_data).
The raw data for Figure 1B consist of four files, which are described as follows:
1. Figure1B_raw_data-CD3_T_tsne.acs: Flow cytometry data for Figure 1B underwent nonlinear dimensionality reduction (t-SNE) analysis and unsupervised Phenograph clustering to analyze T cell subsets in CSC-High and CSC-Low tumors.
2. Figure1B_raw_data-concat_T8_high_T8_low_CD3_T_download_16000_cells.downloading: This flow cytometry dataset consists of T cells from a CSC-high tumor sample and a CSC-low tumor sample, which were downsampled to equal numbers and concatenated for subsequent t-SNE analysis.
3. Figure1B_raw_data-T8_high_HSY_1156553_TUMOR_raw_data_010.fcs: This is a representative raw flow cytometry data from a CSC-high tumor sample.
4. Figure1B_raw_data-T8_Low_GCH_1071022_TUMOR_raw_data_013.fcs: This is a representative raw flow cytometry data from a CSC-low tumor sample.
We hypothesized that the increased Treg infiltration in situ is attributable to a paracrine or systemic effect exerted by CSCs. To test this, an in vitro co-culture system using peripheral blood mononuclear cells (PBMCs) from healthy donors was established. Conditioned medium (CM) collected from fluorescence-activated cell sorting (FACS)-sorted TSPAN8+ human TNBC MDA-MB-231 or BC MCF-7 cells (CM-TS+) and the same cells with TSPAN8 depletion (CM-TS+(kd)) were supplemented into the culture system. As anticipated, co-culturing with CM-TS+ led to a significant increase in the frequency of CD4+FOXP3+ Tregs compared to co-culturing with CM-TS+(kd), as determined by FCM, without altering the frequency of naïve CD4+ T cells (Figure 2B in Processed_data).
The raw data for Figure 2B consist of seven files, which are described as follows:
1. Figure2B_raw_data-20250804_Figure_2B_tsne.acs: This is the data for Figure 2B obtained from flow cytometry, which underwent t-SNE analysis and unsupervised Phenograph clustering. The data were used to generate the Figure 2B image, which includes analysis of T cell subsets in both CM-TS+ and CM-TS+(kd) groups of PBMCs.
2-4. Figure2B_raw_data-Ts-KD1.fcs, Figure2B_raw_data-Ts-KD2.fcs, and Figure2B_raw_data-Ts-KD3.fcs: The raw flow cytometry data from three CM-TS+(kd) groups of PBMCs.
5-7. Figure2B_raw_data-Ts-NC1.fcs, Figure2B_raw_data-Ts-NC2.fcs, and Figure2B_raw_data-Ts-NC3.fcs: The raw flow cytometry data from three CM-TS+ groups of PBMCs.
Using an indirect co-culture system, we treated PBMCs with indicated EVs. Incubation with EVs-TS+ resulted in a marked increase in CD4+ T cells, particularly CD4+FOXP3+ Tregs, and a decrease in CD8+ T cells (especially the CD8+FOXP3− subset) in PBMCs compared to all control conditions (Figure 2G in Processed_data). These data, consistent with prior observations in patient tumor specimens, indicate that EVs derived from TSPAN8+ CSCs mediate the enhanced abundance of Tregs.
The raw data for Figure 2G consist of four files, which are described as follows:
1. Figure2G_raw_data-17-Aug-2025.acs: This is the data for Figure 2G obtained from flow cytometry, which underwent t-SNE analysis and unsupervised Phenograph clustering. The data were used to generate the Figure 2G image, which includes analysis of T cell subsets in EV-depletion, EVs-TS+ and EVs-TS+(kd) groups of PBMCs.
2. Figure2G_raw_data-concat2_1.fcs: This flow cytometry data contains T cells from EV-depletion, EVs-TS+ and EVs-TS+(kd) groups of PBMCs that were downsampled to equal numbers and concatenated for subsequent t-SNE analysis.
The processed data of Figure 1B, 2B, and 2G are presented in Processed_data.pdf.
Code/software
The raw data for Figure 1B, 2B, and 2G can be viewed and analyzed by using FlowJo (version 10.5.3).
Access information
Other publicly accessible locations of the data:
- n/a
Data was derived from the following sources:
- n/a
Human subjects data
Approval for experiments involving human tissue samples was obtained from the Human Ethics Committee at Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China (Approval Letter No. 2019SQ220). Tumor and matched adjacent normal tissue specimens were identified and confirmed by pathologists at the same hospital.
Tumor tissues and matched adjacent normal specimens were collected from 65 patients with primary BC. Following surgical resection, fresh tissue samples were minced with surgical scalpels in ice-cold PBS and then digested with type IV collagenase (Sigma-Aldrich) for 1 h at 37°C on a shaker. The resulting cell suspension was filtered through a 70 μm strainer and rinsed with PBS for cell counting. Cells were subsequently stained with fluorescently conjugated primary antibodies (Table 1). Flow cytometric analysis was performed on a BD FACSymphony instrument, and the data were analyzed using FlowJo software (version 10.5.3). The gating strategy: cells were initially gated based on forward scatter (FSC-A) and side scatter (SSC-A). Single cells were then selected using FSC-A versus FSC-H. Viable cells were identified by excluding those stained with a viability dye. Among the live cells, CD45+^^ cells were defined as immune cells and further subdivided into functional immune phenotypes using established lineage markers. CD45-^^ cells were classified as non-immune. Based on TSPAN8 and CD44 co-expression within the CD45-FAP-CD24- population, samples were categorized into either CSClow or CSChigh groups.
