Abstract:Neoadjuvant chemotherapy (NAC) is one of the standard preoperative therapeutic approaches for breast cancer to downsize tumours, but its clinical efficacy is diverse among patients.1 To understand the causes of disparity in patient outcomes, it is critical to unveil how NAC remodels circulatory immune cells and tumour microenvironment (TME). In this study, we performed single-cell RNA sequencing of 132 823 cells isolated from primary breast tumours with different molecular subtypes, peripheral blood mononuclear cells (PBMCs) collected at three timepoints during patients undergoing NAC and post-NAC tumours respectively (Figure 1A, Table S1). First, we dissected the TME of 17 primary breast carcinomas by separating total cell populations into tumour cells, stromal cells and tumour-infiltrating immune cells. Independent analysis on the amount and gene features of each cell type was then conducted. Immune cells and stromal cells from different samples can be clustered together based on their cell types, while epithelial cells originating from different tumours were separately distributed on the t-SNE plot, suggesting a high degree of intertumour heterogeneity (Figures 1B, C and S1A). We isolated malignant cells from epithelial cells through assessing the copy number variation (CNV) inference (Figures 1D and S1B). Most patients exclusively expressed higher levels of one or more oncogenesis-associated genes (Figure 1E). We then compared the epithelial features of the malignant cells with published normal breast epithelia data. Our results showed that malignant cells from a single tumour could harbour the gene expression profiles derived from multiple subsets of normal epithelial cells (Figures 1F and S2). When analysed by the PAM50 classifier, tumour cells from an individual patient could be classified into multiple molecular subtypes, which demonstrated the presence of heterogeneity within a breast tumour (Figure 1G and S3). In the subsequent analysis of 14 functional states of malignant cells, we observed that epithelial–mesenchymal transition (EMT), invasion and metastasis signatures were enriched in samples from patients with triple-negative breast cancer (TNBC) (Figure 1H), which is consistent with previous findings indicating a more aggressive subtype of breast cancer.2 We then analysed the immune compartment and identified subclusters of T cells, B cells and macrophages. The results were consistent with those of previous studies, in which cytotoxic T cells, regulatory T cells, exhausted T cells, germinal centre (GC) B cells, plasma B cells, dendritic cells (DCs) and plasmacytoid-DCs were identified according to the expression of relevant marker genes (Figure 2A–D).3 Our cell–cell communication analysis revealed that T cells had more positive interaction with tumour cells (Figure 2E, F). Immunosuppressive interactions between tumour cells and T cells were identified including the ligand-receptor pairs NECTIN1–CD96 and NECTIN2–TIGIT (Figure 2G). This finding therefore suggested a plausible approach for breast carcinomas to escape from immune supervision and that TIGIT is a potential immune checkpoint of CD8+ T cells in breast cancer. We next analysed the TME and PBMCs of patients treated with NAC (Figures 3A, B and 4A, B). In paired tumours (pre- and post-treatment) from patient BC06, the malignant cells were isolated by inferCNV. We observed multiple enriched pathways associated with malignancy including angiogenesis, EMT and metastasis in post-treatment tumours (Figure 3E). Through PAM50 analysis, the compositions of tumour epithelial phenotypes were also found to be altered by NAC (Figure 3F). The fractions of T cells and B cells were decreased in the post-treatment tumour compared to the pre-treatment tumour of patient BC18 (Figure 3C) with diminished copy number aberrances (Figure 3D). The IFNα and IFNγ signalling pathways were consistently downregulated in immune cell populations of the post-treatment tumour (Figure 3H). In the analysis of paired tumours from patients BC06 and BC18, we demonstrated that NAC had not only reshaped the heterogeneity and functional features of tumour cells but also changed the proportions of tumour-infiltrating lymphocytes and inhibited the immune activation. Consistently, the fraction of B cells was significantly reduced during and after NAC in PBMCs, suggesting that the chemotherapy could affect the humoral immune response (Figure 4A–C). CXCR4, which is involved in maintaining the B cell population and function,4 was significantly downregulated by NAC as well (Figure 4D, E). The declined levels of STAT1, JUN and NFKBIZ in B cells, cytotoxic T cells, NK cells and DC2 cells respectively from pre-treatment stage to mid- and post-treatment stages indicated that the activation, proliferation and differentiation of these cells could be inhibited by NAC (Figures 4D, E and S4A, B, C).5, 6 We also noticed a slight increase in LAG3 and TIGIT mRNA levels in cytotoxic T cells after NAC (Figure S4A), which implies the occurrence of immune exhaustion induced by NAC.7 The GSVA results showed that the IFNα and IFNγ responses were downregulated in most immune cells after NAC in blood (Figure 4F). Moreover, the pathways which were responsible for transducing type I IFN signals, such as the JAK-STAT1 pathway, were also found to be decreased in T cells, NK cells and myeloid cells (Figure S4D).8 Collectively, these results suggested that the circulatory and local immune cells could be inhibited during or after NAC. The results were opposite to some previous findings that immune activation were commonly observed after NAC.9 However, a recent study revealed that the immune suppression in TME after NAC was linked to patients' resistance to the therapy and indicated poor prognosis, whereas patients with immune restores could benefit.10 In summary, our work provides a comprehensive single-cell transcriptome atlas of the microenvironment in breast tumour including malignant, stromal and immune cells. These results uncover detailed intratumour heterogeneity, immune diversity and the complex communication network of cancer cells and other components. In addition, the single-cell analysis on paired breast tumour samples as well as PBMCs collected during NAC proved that the therapy could reshape the systemic and local immune microenvironments, and the tumours in our study may enter the equilibrium and escape phases after NAC through decreased IFNα and IFNγ responses. In near future, prospective studies with more paired tumour biopsies from patients with pre- and post-treatment would be conducted to further confirm our findings and perform extensive in-depth investigation. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Neoadjuvant chemotherapy-induced remodeling of human hormonal receptor-positive breast cancer revealed by single-cell RNA sequencing

Abstract P2-23-01: Patient-derived Xenografts Allow Deconvolution of Single Agent and Combination Chemotherapy Responses in Triple-Negative Breast Cancer

Effects of Core Needle Biopsy and Subsequent Neoadjuvant Chemotherapy on Molecular Alterations and Outcome in Breast Cancer

Genomic and transcriptomic analysis of breast cancer identifies novel signatures associated with response to neoadjuvant chemotherapy

The Impact of Neoadjuvant Chemotherapy on the Tumor Microenvironment in Advanced High-Grade Serous Carcinoma

Single-Cell RNA Sequencing Reveals Tumor Microenvironment Remodeling after Neoadjuvant Immunotherapy in Non-small Cell Lung Cancer

Unveiling Neoadjuvant Chemotherapy-Induced Immune Landscape Remodeling and Metabolic Reprogramming in Lung Adenocarcinoma by Scrna-Sequencing

Integrated single-cell and bulk RNA sequencing analysis identifies a neoadjuvant chemotherapy-related gene signature for predicting survival and therapy in breast cancer

Unveiling chemotherapy-induced immune landscape remodeling and metabolic reprogramming in lung adenocarcinoma by scRNA-sequencing

Tumor microenvironment remodeling after neoadjuvant immunotherapy in non-small cell lung cancer revealed by single-cell RNA sequencing

The single cell immunogenomic landscape after neoadjuvant immunotherapy combined chemotherapy in esophageal squamous cell carcinoma

Multi-omic subtyping and molecular determinants of neoadjuvant therapeutic regimen for locally advanced HER2-positive and triple-negative breast cancer: Implications of novel therapy.

Neoadjuvant Endocrine Therapy: A Potential Way to Make Cold Hormone Receptor-Rich Breast Cancer Hot.

Chemotherapy induces dynamic immune responses in breast cancers that impact treatment outcome

Abstract 7607: Single-cell genomics study before and after neoadjuvant therapy for triple-negative breast cancer

Neoadjuvant Chemotherapy Alters Peripheral and Tumour‐infiltrating Immune Cells in Breast Cancer Revealed by Single‐cell RNA Sequencing

Single-cell profiling of immune cells associated with response to neoadjuvant chemoimmunotherapy in IIIA non-small cell lung cancer (NSCLC).

Data from Chemotherapy Induces Breast Cancer Stemness in Association with Dysregulated Monocytosis

Cancer-associated Fibroblasts Undergoing Neoadjuvant Chemotherapy Suppress Rectal Cancer Revealed by Single-Cell and Spatial Transcriptomics.

Dynamic Changes in the NK-, Neutrophil-, and B-cell Immunophenotypes Relevant in High Metastatic Risk Post Neoadjuvant Chemotherapy-Resistant Early Breast Cancers

Chemotherapy induces breast cancer stemness in association with dysregulated monocytosis.