Changes in a single-cell expression atlas of gastric cancer patients following NACT
To explore the effects on NACT in gastric cancer, we generated scRNA-seq profiles of pathologically confirmed gastric tumors by evaluating 4 samples received NACT, 4 samples without NACT and 3 pre-post treatment pairs. The clinicopathological characteristics of the patients are shown in Additional file 1: Table S1, included age, sex, tumor location, grade, Lauren classification and neoadjuvant response. After initial quality filtering, approximately 0.28 billion unique transcripts were obtained from 33,589 cells isolated. We classified these cells into clear transcriptional subtypes using clustering-distributed stochastic neighbor embedding (t-SNE) based on informative principal components (n = 21) (Additional file 3: Fig. S1A–C). Based on known classic markers, 11,535 cells (34%) were classified as epithelial cells, while the remaining cells (66%) included immune cells (i.e., myeloid, T, and B cells), endothelial cells and fibroblasts (Additional file 3: Fig. S1D).
We observed significantly fewer T cells in the post-treatment samples than in the pre-treatment samples (P = 0.010), while the other cell types showed no significant change (Additional file 3: Fig. S1E). To further verify this, we collected 40 patients with gastric cancer (20 post-treatment samples and 20 pre-treatment samples) and performed immunohistochemical staining (Additional file 3: Fig. S1F). The results showed similar downward trends in CD4+ and CD8+ T cell accumulation, indicating that NACT has important impair effects on immune cells (Additional file 3: Fig. S1G).
Next, we performed further canonical correlation analysis on non-epithelial cells in post-treatment samples and pre-treatment samples, which revealed a complex TME that contained 14 different cell types (plasma cells, B cells, CD4+ T cells, CD4+ regulatory T (Treg) cells, CD8+ T cells, cycling T cells, dendritic cell 3 (DC 3), dendritic cell 2 (DC 2), endothelial cells, fibroblasts, myofibroblasts, mast cells, monocytes, and macrophages) (Fig. 1A–D). We observed significantly more endothelial cells and fibroblasts in the post-treatment samples than in the pre-treatment samples (P = 0.034 and P = 0.005) (Fig. 1E). In addition, the proportions of CD4+ T cells and CD8+ T cells were significantly declined during the NACT (P = 0.018 and P = 0.010). However, there was no significant change in CD4+ Tregs cells, cycling T cells, B cells, plasma cells, macrophages, monocytes, dendritic cells, and mast cells (P > 0.05). The detailed changes in immune cell types between pre and post NACT were showed in Additional file 2: Table S2.
Altered expression of the nonnegative matrix factorization (NMF) 1 and NMF2 transcriptional programs improves the prognosis of patients treated with NACT
To investigate how NACT affects tumor cells, we compared the pathway activation of tumor cells in all post-treatment and pre-treatment samples by gene set variation analysis (GSVA). The results showed that pathways associated with cell proliferation (MYC targets, G2M checkpoint and oxidative phosphorylation) and protumor pathways (epithelial-mesenchymal transition (EMT) and angiogenesis) were upregulated in the post-treatment samples compared with the pre-treatment samples (Fig. 2A), indicating that surviving tumor cells not killed by NACT may have higher proliferative activities and a greater ability to metastasize. The up-regulation of the oxidative phosphorylation pathway after treatment suggests that the mitochondrial functional status of tumor cells may be altered. Therefore, we analyzed the changes of mitochondrial-related genes before and after treatment, and found that a series of mitochondrial genes were up-regulated after treatment (COX4I1, COX6C, COX7B, COX8A, NDUFA4, NDUFC2, PPA1 and UQCRH) (Additional file 4: Fig. S2A). These mitochondrial genes are mainly related to the ATP metabolic process and the electron transport chain [23], thus up-regulation of these genes indicates that energy metabolism of tumor cells was enhanced by NACT, which may promote the growth of tumor cells.
Moreover, immune-associated pathways (the interferon alpha (IFN-α) response, interferon gamma (INF-γ) response, PD1 signaling and MHC pathways) were downregulated after treatment (Fig. 2A), indicating that the TME developed a more immunosuppressive status. To analyze the toxicological responses of various types of cells to NACT, we simultaneously compared the changes of toxicology and drug-related pathways in tumor cells and other cell types before and after treatment. Interestingly, our results showed that the changes in these pathways were distinct for different cell types. The activities of multiple pathways (GENOTOXIC_DAMAGE, DRUG_ADME, DRUG_METABOLISM_CYTOCHROME_P4501) in tumor cells were significantly increased after treatment, while the levels of these pathways in various immune cells were down-regulated or did not change significantly after treatment (Additional file 4: Fig. S2B). These suggests that the toxic responses of different cells to NACT may be cell type-specific.
Then, to analyze changes in tumor cell gene expression between post-treatment and pre-treatment samples, we analyzed the differential gene expression of tumor cells. IFI27, LCN2, GAST, and SERPINA1 were downregulated after NACT in tumor cells, while OLFM4, MRPS35, MMP1 and MED21 were upregulated in tumor cells (Fig. 2B). IFI27 is involved in signaling pathways of apoptosis and type-I interferon [24], while the downregulation of IFI27 may enhance the ability of tumor cells to survive. LCN2 encodes a protein of the lipocalin family, which plays roles in innate immunity by binding bacterial siderophores [25]. LNC2 is also found to be overexpressed in multiple type of cancer and is involved in tumor growth and metastasis. Gastrin, a GAST-encoded protein, is required for stomach to secret digestive enzymes and hydrochloric acid [26]. Down-regulation of GAST suggest normal function of stomach may be affected, which may be associated with adverse effects of NACT. OLFM4 is reported to be upregulated in gastric cancer [27], and is reversely correlated with NF-κB/IL-8 pathway in gastric cancer [28]. Upregulation of OLFM4 by NACT may affect NF-κB mediated immune responses. MMP1 encodes a member of matrix metalloproteinases, and was associated with cancer invasion and metastasis by activating PAR-1 [28]. In addition, MMP1 may also regulate monocyte recruitment by inducing expression of CCL-2 [29]. The roles of MRPS35 and MED21 in tumor and immune regulation are still unknown, so further research is needed to explore the mechanism for the up-regulation of these genes and their implications for treatment.
As intratumoral heterogeneity is an important feature of tumors and may play an important role in chemotherapy resistance, we next analyzed the characteristics of tumor heterogeneity in gastric tumor cells in all samples to further study the intra-tumor heterogeneity dynamics during NACT. Unsupervised nonnegative matrix factorization (NMF) analysis revealed that there were three transcriptional programs (NMF1, NMF2 and NMF3) shared among these samples (Fig. 2C). NMF1 was defined by cell junction organization markers (CLDN4, ITGB1, JUP and LAMB3) and epithelial cell differentiation markers. NMF2 involved markers of cell cycling (for example, TOP2A, CDK1 and RRM2), and NMF3 was enriched for ribosomal genes, translation initiation and elongation factors (EIF3E and EEF1A1), and oxidative phosphorylation markers (Fig. 2D). We concluded that these programs reflected cell cycle activity (NMF2), undifferentiated progenitors (NMF3) and more differentiated epithelial cell programs (NMF1).
We further compared the difference in these three transcriptional programs between post-treatment and pre-treatment samples. The signature score of NMF1 was significantly downregulated after treatment, and gene expression analysis also showed that several cell junctions associated genes (CLDN4, CTNND1, ITGB1 and JUP) were significantly down-regulated after treatment (Additional file 4: Fig. S2C). In tissues, cell junctions interconnect cells and maintain homeostasis by regulating tissue barrier function, migration and proliferation of cells [30]. Intercellular communication through cell junctions may be involved in the development of diseases [31], defects of which would induce imbalance of tissue homeostasis and are widespread in cancers. Downregulation of cell junction genes post NACT suggest that function of cell–cell junctions might be altered by NACT, which may promote tumor invasion and migration. While the NMF2 signature was significantly upregulated after treatment, the NMF3 in post-treatment samples showed comparable expression (Fig. 2E).
We then determine if these changes following NACT would affect prognosis. The downregulated NMF1 and upregulated NMF2 signatures were both associated with improved overall survival outcomes based on The Cancer Genome Atlas (TCGA) database (Fig. 2F). These data indicated that gastric cancer patients treated with NACT following gastrectomy had a better prognosis than gastric cancer patients not treated with NACT following gastrectomy.
The changes of immune cells after NACT may be associated with immunosuppression
As the decreased proportion of T cells indicated that chemotherapy killed immune cells, subsequently forming a more immunosuppressive environment, we further investigated the changes of pathway and gene markers among T cells. Compared to the T cells in the pre-treatment samples, all T cells in the post-treatment samples showed lower cytotoxicity and proliferative marker expression patterns (Fig. 3A, top and middle panels) than those in pre-treatment tumors. This founding has also been reported in cervical cancer after radiotherapy [16]. The exhaustion levels of CD4+ T cells and cycling T cells were significantly higher in the pre-treatment samples than in the post-treatment samples, while CD8+ T cells and CD4+ Tregs showed no significant difference (Fig. 3A, bottom panel). We observed pervasive changes in the pathway activities of T cells between the post-treatment samples and the pretreatment samples. For CD4+ T cells, CD8+ T cells and CD4+ Treg cells, multiple immune-associated pathways were downregulated, including allograft rejection, the CTLA4 pathway, and the INF-α and INF-γ responses (Fig. 3B). Allograft rejection activities in this setting are presumably related to cells showing higher reactivity to malignancy cell-encoded neoepitopes [32]. Targets of MYC, which is involved in a proliferation-related pathway, showed decreased expression after NACT treatment, which was consistent with the observed decreased T cell proportions (Figs. 1E and 3B). After NACT, T cells showed decreased expression of metal-binding and proinflammatory genes, such as GNLY, ANXA1, CCL4, IL17A, CCL4L2 and XCL2 in CD4+ T cells and MT1G, CCL20, GNLY, and CCL3L3 in CD8+ T cells (Fig. 3C).
We further explored the characteristics of different myeloid cell subtypes comparing post-treatment versus pre-treatment samples. Monocytes and macrophages were identified based on known markers and the gene profiles of these myeloid cells were generated. Although the percentages of these cells kept stable after NACT (Fig. 1E), macrophages and monocytes exhibited varying levels of downregulation of MYC targets, the INF-α response, the INF-γ and IL6-JAK-STAT3 signaling (Additional file 5: Fig. S3A). These observations were consistent with CXCL10 downregulation in macrophages (Additional file 5: Fig. S3B). We furtherexamined the expression of known functional signatures in macrophages and monocytes, including M1-like and M2-like signatures. Unlike pre-treatment samples, the post-treatment samples showed a decreased expression level of M1-like signature (Additional file 5: Fig. S3C).
Stromal cells after NACT had increased angiogenesis activity
Then we investigated endothelial cells and fibroblasts at the pathway and gene level. First, post-treatment endothelial cells were shown to have MYC target upregulation (Fig. 3D, left panel), indicating that these cells were actively proliferating, which was consistent with the increased proportion of endothelial cells observed after NACT (Fig. 1E). Angiogenesis had increased activity in endothelial cells after chemotherapy (Fig. 3D, left panel), which is consistent with the function of residual tumor cells after NACT. It is well known abnormal blood vessel growth is a major factor in tumor occurrence and progression [7]. Recent evidence showed that tumor-associated endothelial cells are key players in cancer cell evasion of immune surveillance and enhanced chemoresistance [33], which was consistent with our results. Next, we examined fibroblasts, as fibroblasts are thought to be a highly plastic cell population in the TME and how fibroblast cells change with NACT treatment has not been evaluated. Compared to the pre-treatment samples, the post-treatment samples exhibited a higher percentage of fibroblast cells (Fig. 1E). Moreover, fibroblasts and myofibroblasts showed decreases in immune pathways, including the IFN-α response and IFN-γ responses (Fig. 3D, middle and right panels). CD69 gene was highly expressed in fibroblasts in the post-treatment samples (Fig. 3E). Some chemokine-encoding genes, such as CXCL10 and CXCL14, were downregulated in fibroblasts in post-treatment samples (Fig. 3E). Therefore, there is a close relationship between tumor-related fibroblasts and immune cells. Costa et al. indicated that fibroblasts contributed to immunosuppression in breast cancer [34]. Another study showed that proliferating T cells produced less IFN-γ and TNF-α when fibroblast cells were present in pancreatic cancer, thus contributing to diminished immune function [35].
Gastric cancer is characterized by rewired cell–cell interaction networks after NACT
After identifying different cell subtypes present in the gastric cancer TME, we further studied the associations between different transcriptional profiles in the TME and used CellPhone DB to identify the ligand-receptor pairs and molecular interactions between the major cell types (Fig. 4A). In general, the number of predicted interactions between tumor cells and immune cells was obviously reduced after treatment, while the interactions between endothelial cells and tumor cells were enhanced (Fig. 4B and Additional file 6: Fig. S4A). These findings suggest that a new immune microenvironment balance emerged after chemotherapy.
Next, analysis of secreted ligands for the detected cognate receptors demonstrated that extensive communication occurred among different cell types in the TME (Fig. 4C–F, Additional file 6: Fig. S4B–D). For tumor cells, interactions between tumor cells and T cell subsets via CXCL16-CXCR6, CCL20-CXCR3, PVR-CD226, and NECTIN2-CD226 decreased after treatment, indicating that recruitment and co-stimulation interactions between tumor cells and T cells were reduced by treatment, which is similar with the decreased activities in recruitment and stimulatory interactions were seen for macrophages. In contrast, the interaction scores of DLL4-NOTCH3, DLL1-NOTCH1, TGFBR3-TGFB1, and IGF1R-IGF1 between tumor cells and endothelial cells increased after treatment (Fig. 4D). Communication between CD8+ T cells and other immune cells via CCL20-CXCR3, CCL20-CCR6, CD52-SIGLEC10 and NECTIN2-CD226 was decreased after treatment (Fig. 4E, F), these results indicate that T cell associated activities of activation, proliferation and recruitment declined [36, 37]. These findings were consistent with the phenomenon of CD4+ T cells, CD4+ T cells saw similar drops in activation, proliferation and recruitment.