Skip to main content

Spatial transcriptomics atlas reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment components in colorectal cancer

Abstract

Background

The tumor-promoting role of tumor microenvironment (TME) in colorectal cancer has been widely investigated in cancer biology. Cancer-associated fibroblasts (CAFs), as the main stromal component in TME, play an important role in promoting tumor progression and metastasis. Hence, we explored the crosstalk between CAFs and microenvironment in the pathogenesis of colorectal cancer in order to provide basis for precision therapy.

Methods

We integrated spatial transcriptomics (ST) and bulk-RNA sequencing datasets to explore the functions of CAFs in the microenvironment of CRC. In detail, single sample gene set enrichment analysis (ssGSEA), gene set variation analysis (GSVA), pseudotime analysis and cell proportion analysis were utilized to identify the cell types and functions of each cell cluster. Immunofluorescence and immunohistochemistry were applied to confirm the results based on bioinformatics analysis.

Results

We profiled the tumor heterogeneity landscape and identified two distinct types of CAFs, which myo-cancer-associated fibroblasts (mCAFs) is associated with myofibroblast-like cells and inflammatory-cancer-associated fibroblasts (iCAFs) is related to immune inflammation. When we carried out functional analysis of two types of CAFs, we uncovered an extensive crosstalk between iCAFs and stromal components in TME to promote tumor progression and metastasis. Noticeable, some anti-tumor immune cells such as NK cells, monocytes were significantly reduced in iCAFs-enriched cluster. Then, ssGSEA analysis results showed that iCAFs were related to EMT, lipid metabolism and bile acid metabolism etc. Besides, when we explored the relationship of chemotherapy and microenvironment, we detected that iCAFs influenced immunosuppressive cells and lipid metabolism reprogramming in patient who underwent chemotherapy. Additionally, we identified the clinical role of iCAFs through a public database and confirmed it were related to poor prognosis.

Conclusions

In summary, we identified two types of CAFs using integrated data and explored their functional significance in TME. This in-depth understanding of CAFs in microenvironment may help us to elucidate its cancer-promoting functions and offer hints for therapeutic studies.

Introduction

As one of the most common tumors in the world, it is estimated that more than 1.9 million new colorectal cancer (CRC) cases and more than 935,000 deaths in 2020 by the International Agency for Research on Cancer (IARC), ranking third in incidence and second in mortality among known cancers [1]. Though substantial progress has been made in cancer pathogenesis and drug therapy, the clinical prognosis of colorectal cancer is still poor due to tumor metastasis [2].

In recent years, the tumor microenvironment (TME) has become a research hotspot of mechanisms of tumor biology and drug development. The tumor microenvironment is a complex system that contains a variety of cellular and noncellular components, such as: immune cells, endothelial cells, cancer-associated fibroblasts (CAFs) and cytokines et al. [3, 4]. More and more evidences indicate that the crosstalk between various stromal components in TME and tumor cells are critical factors affecting tumor growth and metastasis [5]. CAFs, as one of the most abundant stromal cell types in TME, can remodel the extracellular matrix (ECM), and promote tumor progression through the interactions with tumor cells and immune cells by secreting various growth factors, chemokines and cytokines [3]. Some studies have identified different types of CAFs at the single-cell level according to unique gene signatures or functions. However, the types and functions of CAFs varies in different studies or tumor types, suggesting that the functional roles of CAFs in TME is complicated and yet has not been clearly explained, so it is worth for further exploration [6,7,8].

In this research, we analyzed previously published spatial transcriptomics (ST) data and profiled a spatial atlas of TME inside colorectal cancer tissues [9]. Our work highlights the role of CAFs in colorectal cancer. In detail, we combined spatial transcriptomics with ssGSEA to precisely identify different cell types. We found that there were subsets enriched with inflammatory-cancer-associated fibroblasts (iCAFs) and myo-cancer-associated fibroblasts (mCAFs) in CRC with bioinformatics analysis and experimental verification. Through functional enrichment analysis, the functional role of CAFs in TME were further investigated. In addition, we explored the correlation between CAFs and prognosis by analyzing bulk RNA-sequencing in public dataset. These results promote an in-depth understanding of the functions of CAFs in TME and provide a basis for CRC precision therapy.

Methods

Data sources

The spatial transcriptomics dataset (10X Genomics) was downloaded from a spatial transcriptomics research website (http://www.cancerdiversity.asia/scCRLM/). The spatial transcriptomics data of two patients were used, one of whom did not receive neoadjuvant chemotherapy treatment (NACT) was named colon1 (ST-P1), and the other received neoadjuvant chemotherapy with partial response (PR) was named colon2 (ST-P3) [9]. Bulk RNA-sequencing dataset was derived from the COADREAD cohort of The Cancer Genome Atlas (TCGA) and downloaded from UCSC XENA (https://xena.ucsc.edu/).

Spatial transcriptomics data processing

We used the R package Seurat (v4.1.0) to process space transcriptomics data and used log-normalization to standardize data [10]. We used functions SelectIntegrationFeatures, FindIntegrationAnchors, and IntegrateData to remove batch effects and integrate Seurat object into a single ST dataset. To reduce the dimensionality of ST data, function RunPCA was performed, then functions FindNeighbors and FindClusters were used to cluster similar ST spots. Different clusters were preliminarily annotated based on hematoxylin–eosin staining (H&E) sections and unsupervised clustering analysis. When we used cell markers to annotate clusters, we found that some clusters highly expressed multiple cell markers, so ssGSEA algorithm [11] were performed on scoring common cell types based on the average expression matrix of different clusters, and studies have confirmed its robustness in ST [12, 13].

Enrichment analysis

Differentially expressed genes (DEGs) of cell clusters were identified by the function FindAllMarkers, and top 10 DEGs of each cluster that were ranked according to the log2FC were applied to data visualization. Function clusterProfiler (v3.18.1) [14] was used for KEGG pathway analysis. ssGSEA (Hallmark Gene sets from Molecular Signatures Database, MSigDB) was conducted with GSVA (v1.38.2) [15]. We set adj.p.val < 0.05, |avg. logFC|> 0.5 as cut-off criteria.

Bulk RNA-seq data analysis

We downloaded TCGA colorectal cancer cohort data (COADREAD) for transcriptome analysis. We selected the top 10 DEGs in the cluster as cluster-specific geneset, and calculated the scores of genesets in the transcriptomic data through GSVA, median was set as cut-off value. R packages survival (v3.2-10), survminer (0.4.9) were used for survival analysis. In addition, ImmuneScore and StromalScore were calculated using ESTIMATE algorithm, and their correlation with specific cell type was analyzed by correlation analysis, p value < 0.05 as statistically significant.

Pseudotime analysis

In order to analyze the cell type differentiations in the tumor microenvironment, monocle2 (v2.18.0) [16] was used for trajectory analysis to find the transitional relationships among different clusters. The plot_genes_in_pseudotime function was applied to discover transitional changes in gene expression levels among different clusters.

Immunofluorescence assay

Specific proteins expression of distinct cancer-associated fibroblasts within colorectal cancer were analyzed using Immunofluorescence staining. Paraffin-embedded colorectal cancer tissues (n = 2) were made slices, heat mediated antigen retrieval with TRIS–EDTA (pH = 8) for 20 min and then blocked with 5% BSA for 30 min at room temperature. After incubated with anti-PDGFRA (rabbit, 1:100, HuaBio, ET1702-49) and anti-RGS5 (rabbit, 1:300, Affinity, Cat. #: DF4417) separately at 4℃ overnight, the sections were washed with PBS for three times, then 488 conjugated goat anti-rabbit IgG (1:800 dilution, HA1121) and 594 conjugated goat anti-rabbit IgG (1:800 dilution, HA1122) were separately added and incubated for 1 h. Finally, DAPI was applied to stain cell nucleus. Sections without incubation with primary antibody or secondary antibody were used as control (without any light).

Immunohistochemical analysis

Immunohistochemical analysis was carried out to explore the specific protein expression in distinct tissues. Paraffin-embedded colorectal cancer or para-carcinoma tissues (n = 5) were applied in this project. Sections were firstly treated with boiling TRIS–EDTA (pH = 8) to repair antigen for 20 min. Then endogenous peroxidases were wiped off. Thereafter, tissues were incubated with 5% BSA for 30 min. Anti-PDGFRA (rabbit, 1:100, HuaBio, ET1702-49) was added in sections at 4 °C overnight. After washed with PBS for three times and probed with HRP conjugated compact polymer system for 30 min, DAB was used as the chromogenic agent. Hematoxylin was used to dye the cell nucleus.

Statistical analysis and visualization

All statistical analysis were based on R (v4.0.5), and data visualization was performed on R packages Seurat (v4.1.0), ggplot2 (v3.3.5), ggsignif (v0.6.1), pheatmap (v1.0.12) and ggstatsplot (v0.9.1) [17].

Results

Spatial transcriptomics to identify intra-tumor heterogeneity in CRC

In this study, we focused on two patients’ spatial transcriptomics data from a previously published CRC research (ST-P1, untreated, ST-P3, NACT with PR) [9]. Firstly, we used unsupervised clustering to cluster similar ST spots and divide into nine distinct clusters (Fig. 1A, B). Secondly, we annotated the clusters according to H&E sections and cell markers, then five morphological regions were identified (Fig. 1A, C, D). The results showed that unsupervised clustering analysis could effectively cluster ST spots with similar features into the same cluster, such as fibroblast and tumor regions (Fig. 1E). In addition, unsupervised clustering analysis could also subdivide tissues sections, which facilitates the discovery of tissue heterogeneity that is invisible to naked eye.

Fig. 1
figure 1

Spatial transcriptomics (ST) to identify intra-tumor heterogeneity in colorectal cancer. A Spatial images of unsupervised clustering results (left) and hematoxylin–eosin staining (H&E) sections annotations (right). B, C UMAP plot of cell types and PCA clustering result profiled in the presenting work. D The expression level of cell markers across different morphological regions through the Dot plot. E Spatial plots show the spatial expression pattern of fibroblast and tumor markers (COL1A1, EPCAM) in this study. F Heatmap of cell-type scores estimated by ssGSEA

Because each spot of spatial transcriptomics contains more than one cell, its accuracy is lower than that of single-cell sequencing. Therefore, we scored well-defined genesets of 12 cell types [7, 18,19,20,21,22,23,24] by ssGSEA algorithm to recognize cell types that contained in each cluster. The annotation results had demonstrated the effectiveness of using ssGSEA in identifying cell populations of each cluster. As shown in Fig. 1F, fibroblast, enterocyte and smooth muscle were precisely annotated, and the lamina propria is composed of loose connective tissue that contains a variety of cells such as: dendritic cells, NK cells, etc. [25], while the tumor region was significantly enriched with monocytes, NK cells, and epithelial cells, indicating the existence of immune-inflammatory microenvironment in tumor region, which was consistent with previous studies [4, 5]. In all, spatial transcriptome analysis combined with ssGSEA can accurately determine the cell types contained in cell populations and offset the insufficiency of resolution in ST.

CAFs-enriched subgroups are identified in CRC

To analyze the heterogeneity of tumor region, principal component analysis (PCA) was used to cluster tumor region into 5 subclusters (Tumor-subcluster 0 ~ 4) (Fig. 2A). To find the cell types and functions of each cluster, we analyzed cluster-specific differential expression genes (DEGs) and found that subcluster 0 highly expressed fibroblast markers, such as: COL1A2, COL1A1, SPARC and COL3A1, etc. (Fig. 2B). Therefore, we considered subcluster 0 as CAFs-enriched cluster. Li et al. [19] previously identified two types of CAFs at single-cell sequencing of CRC, in which CAF-A expressed genes related to ECM remodeling, while CAF-B expressed cell markers of myo-fibroblasts, such as ACTA2 and TAGLN. Elyada et al. [6] identified two distinct types of CAFs from human pancreatic cancer, named mCAFs and iCAFs, an emerging study [7] also confirmed the presence of two different fibroblasts subtypes in bladder cancer. As seen in Fig. 2C, to confirm the existence of CAFs with specific signatures, we applied immunofluorescence analysis and we detected the expression of PDGFRA (iCAFs marker) and RGS5 (mCAFs marker) in colorectal cancer [7]. Thereafter, we utilized ssGSEA to further investigate the cell subtype of subcluster 0. And the result showed that subcluster 0 enriched with mCAFs (mCAFs-enriched cluster) (Fig. 2D, F).

Fig. 2
figure 2

CAFs-enriched subgroups are identified in CRC. A UMAP plot of unsupervised clustering result of tumor subclusters. B Heatmap of differentially expressed genes (DEGs) of different Tumor-subclusters. C Immunofluorescence analysis of specific proteins expression of iCAFs and mCAFs. D Heatmap of tumor microenvironment (TME) cell-type scores calculated by ssGSEA in tumor region. E KEGG enrichment functions of up-regulated genes in Tumor-subcluster 0. F Violin plots show the expression of tumor, inflammatory-cancer-associated fibroblasts (iCAFs) and myo-cancer-associated fibroblasts (mCAFs) markers in tumor region. G Heatmap of tumor microenvironment (TME) cell-type scores calculated by ssGSEA in fibroblast region

To determine the functions of mCAFs-enriched cluster, we carried out KEGG pathway enrichment analysis. As shown in Fig. 2E, mCAFs-enriched cluster (Tumor-subcluster 0) was related to ECM-receptor interaction, focal adhesion, and proteoglycans in cancer et al., suggesting that mCAFs has the function of ECM remodeling in TME, which was consistent with previous studies [6, 7, 26]. It was proved again that ssGSEA was robust for spatial transcriptomics data.

In addition, we also conducted in-depth analysis on fibroblast region (Additional file 1: Fig. S1). The result showed that Fibro-subcluster 1 enriched iCAFs, named iCAFs-enriched cluster, and Fibro-subcluster 2 enriched mCAFs (Fig. 2G), because it highly expressed EPCAM, thus we identified it as mCAFs-enriched tumor cluster (Fig. 3D, E). Interestingly, when subclusters annotated by cell markers, we found that marker of CAFs, such as RGS5 and ACTA2 in mCAFs, PDGFRA and CXCL12 in iCAFs, were not specifically expressed in a certain subcluster (Fig. 2F, Additional file 1: Fig. S1D). Therefore, it further indicates that the accuracy of 10× spatial transcriptome data is limited, and spatial transcriptomics combined with multidimensional analysis (such as ssGSEA) can provide more detailed information. In detail, we found that anti-tumor immune cells such as T cells and dendritic cells were significantly enriched in mCAFs-enriched tumor cluster, iCAFs and macrophages were co-enriched in iCAFs-enriched cluster, however the anti-tumor immune cells, especially NK cells, reduced significantly in iCAFs-enriched cluster (Fig. 2G). It has been reported that CAFs can inhibit the immune system by releasing cytokines, chemokines and other compounds, thus leading to tumor metastasis [7, 27], indicating the important role of iCAFs in the immune microenvironment and its value as an anti-cancer drug target.

Fig. 3
figure 3

iCAFs promote oncogenesis by altering tumor microenvironment. A Heatmap shows ssGSEA result to explore the functional roles of different fibro-subclusters. B, C Trajectory of differentiation of different Fibro-subclusters predicted by monocle 2. D, E The differential expression of EPCAM in distinct Fibro-subclusters. F Gene expression level in single spot ordered along the pseudotime for EMT markers

iCAFs promote oncogenesis by altering tumor microenvironment

We applied ssGSEA enrichment analysis on fibroblast region to reveal the functions of iCAFs in TME, and our results indicated that iCAFs were associated with epithelial mesenchymal transition (EMT), cholesterol homeostasis, bile acid metabolism, and fatty acid metabolism (Fig. 3A). To further explore the transitional relationships of EMT phenotype in distinct clusters (Fig. 3B, C), we examined the expression level of EMT markers in different clusters through pseudotime analysis (Fig. 3F). Strikingly, the results demonstrated an increased expression of EMT signatures in iCAFs-enrich cluster, while lack of expression in tumor site (Fibro-subcluster 2). These results suggest that CAFs, rather than tumor epithelial cells, might be the main culprit in promoting EMT and leading to tumor metastasis [19]. Metabolic reprogramming also plays an important role in tumor proliferation and metastasis [28]. In this study, iCAFs were supposed to be correlated with lipid metabolism through the bioinformatics analysis (Fig. 3A), suggesting that the changes of lipid metabolic activity of iCAFs in TME might be a potential mechanism in promoting tumor tumorigenesis.

Chemotherapy alters the tumor microenvironment

When we integrated two patients’ datasets, we found a patient-specific ST expression pattern (Fig. 4B), so we used CCA [10] to remove the batch effects on ST data. It is noteworthy that there was still obvious heterogeneity between colon1 and colon2 (Fig. 4A), therefore we supposed that chemotherapeutic drugs may be one reason for explaining the differential expression profiling in TME. In detail, we analyzed the cell compositions of two patients, results showed that the putative proportion of iCAFs relatively increased in all clusters (Fig. 4C left, Fig. 4D and Additional file 1: Figs. S2, S3) of colon2 (NACT with PR), while the proportions of some anti-tumor immune cells decreased, such as: NK cells, monocytes, etc. (Fig. 4C left, Fig. 4D and Additional file 1: Figs. S2, S3). In addition, we also examined the metabolism pattern changes caused by chemotherapeutic drugs, and we found that the metabolic activity significantly decreased in colon2 (Fig. 4E). However, when we focused on iCAFs-enriched cluster, we detected that fatty acid metabolic activity did not decrease in colon2 (Fig. 4F), then we supposed this phenomenon may be associated with iCAFs. Perhaps, it is the decrease of these tumor-killing cells and metabolism patterns changes caused by iCAFs that explain the internal mechanisms of drug resistance.

Fig. 4
figure 4

Chemotherapy alters the tumor microenvironment. A, B UMAP plots show patient-specific spatial transcriptomics (ST) expression patterns between colon1 and colon2, A UMAP plot of ST data after CCA, B UMAP plot of ST data without CCA. C Cell compositions predicted in different patients by ssGSEA, predicted cell compositions in all clusters (left), predicted cell composition in fibroblast cluster (middle), predicted cell composition in tumor cluster (right). D Boxplots show the cell compositions of iCAFs, NK cells and monocytes in all clusters. E, F Box plots show the patient-specific metabolism patterns scored by ssGSEA. p value < 0.05 was considered as cut-off criteria

iCAFs is associated with clinical prognosis and immune infiltration

To correlate spatial transcriptomics data with public dataset, we evaluated the clinical significance of iCAFs in TCGA COADREAD cohort. Clinicopathological analysis showed that iCAFs were significantly associated with lymph node invasion (p < 0.05), the higher of iCAFs, the more possibility of lymph node invasion (Fig. 5A–D). In addition, univariate (Fig. 5E, F) and multivariate Cox proportional hazards regression (Fig. 5G) were used to analyze the association of the accumulation of iCAFs with prognosis. Multivariate analysis after adjusting confounding factors showed that the accumulation of iCAFs was an independent prognostic factor for OS (p = 0.035) (Fig. 5G). To explore the differential expression of iCAFs in distinct tissues, we specially detected the expression of iCAFs marker (PDGFRA) in tumor tissues and para-carcinoma tissues (n = 5) using immunohistochemical method, and our results demonstrated that the higher proportion of iCAFs in tumor tissues rather than para-carcinoma tissues (Fig. 5H). These results further indicated that iCAFs were associated with poor prognosis. When we focused on the relationship between iCAFs and immune inflammation, ESTIMATE algorithm was applied to calculate the immune score and stromal score for TCGA cohort, and the results showed that iCAFs was significantly correlated with immune score (r = 0.39, p < 0.05) (Fig. 5I, J), suggesting that iCAFs could interact with immunosuppressive cells [27, 29], thus inhibit the anti-tumor inflammatory response in TME.

Fig. 5
figure 5

iCAFs is associated with clinical prognosis and immune infiltration. AD Box plots show the correlations between iCAFs and clinicopathologic features. E, F K-M survival plots show that high iCAFs predicted poor prognosis in TCGA COADREAD cohort, OS: overall survival, DSS: disease specific survival. G Forest map shows that iCAFs is an independent prognostic factor for OS. H Immunohistochemical analysis of PDGFRA expression in colorectal cancer or para-carcinoma tissues. I, J The diagrams show the relationships between iCAFs and immune score as well as stromal score. p value < 0.05 was considered as cut-off criteria

Fig. 6
figure 6

Workflow of bioinformatic analysis. Two types of datasets were utilized in this project. ST data was used to explore the role of CAFs in TME within patients with colorectal cancer, while bulk RNA-seq dataset was applied to validate the clinical translational value of CAFs identified in ST data. Then immunofluorescence staining and immunohistochemical was applied to validate the results based on bioinformatics analysis and we confirmed the functions of two CAFs types in CRC

Discussion

Currently, targeted drugs combined with surgery has largely reduced the mortality of colorectal cancer, however tumor metastasis and drug resistance are still important factors leading to poor prognosis of CRC patients [8], so new therapeutic strategies are waiting to be discovered. Previous studies pay more attention to the biology of tumor cells, nevertheless, emerging research gradually noticed the important roles of tumor microenvironment in tumorigenesis, metastasis and drug resistance [5, 30], while CAFs as the main stromal components in the microenvironment, can affect the tumor growth and metastasis through various mechanisms [3]. In this research, we utilized spatial transcriptome data and public TCGA cohort dataset to elucidate CAFs functions and its interactions with the microenvironment, hoping to promote the development of drug treatment strategies, a workflow was drew to clearly show the overview of our project (Fig. 6).

Li et al. firstly divided CAFs into two distinct subpopulations (CAF-A and CAF-B) based on the single-cell sequencing data in CRC. CAF-A was related to ECM remodeling, while CAF-B was similar to myo-fibroblasts. However, the functional roles of these two types of CAFs remain unclear [19]. In recent years, studies have classified CAFs into iCAFs and mCAFs according to their functions in pancreatic cancer, prostate cancer and triple negative breast cancer [6, 7, 31]. We also identified the existence of these two types of CAFs in CRC with fundamental experiment. An increasing number of research reveal that CAFs through various approaches such as receptor activation, cytokine and cytotoxic production, to inhibit NK cells [3, 27, 32], and we also found that NK cells were decreased in iCAFs-enriched cluster, which was consistent with these studies. In addition, we detected that iCAFs and macrophages were co-enriched in the same subcluster by using bioinformatics analysis. Zhang et al. reported that CAFs can inhibit NK cells proliferation by regulating tumor-associated macrophages (TAM) [29], suggesting the extensively interactions between iCAFs and immune cells in TME. Besides, the relationships of chemotherapy drugs and tumor microenvironment were also examined, and we found that the proportion of iCAFs increased while the number of NK cells and monocytes decreased in patient who underwent neoadjuvant chemotherapy (colon2). Since iCAFs can inhibit immune cells in TME, the decrease of immune cells may be related to iCAFs [27]. Therefore, we speculated that the increase of iCAFs in patient who underwent chemotherapy might be an important mechanism for the occurrence of drug resistance. In the future, more fundamental experiments should be conducted to elucidate the underlying mechanism.

Although epithelial-mesenchymal transformation (EMT) has been widely studied, we still have limited understanding of the mechanism of this phenomenon [33]. Our study indicated that iCAFs up-regulated EMT signature, while tumor cells lacked of an EMT marker, suggesting that EMT may be driven by fibroblasts rather than tumor cells [19, 33, 34]. In addition, metabolic reprogramming is also an important mechanism for tumor progression. More and more studies have shown that CAFs can promote tumor metastasis through the metabolic crosstalk with tumor cells or other stromal cells in TME [28, 35]. Despite the known glycolysis and amino acid metabolism play important roles in TME, there are emerging evidences demonstrate that lipid metabolism is also essential for tumor development [36, 37]. This is consistent with our ssGSEA pathway analysis. Furthermore, when we investigated the relationship between chemotherapy and metabolism, we discovered opposite trends between fatty acid metabolism pattern and other metabolism patterns in iCAFs-enriched cluster, and we speculated this phenomenon might lead to drug resistance. Therefore, targeted therapy toward iCAFs as well as its metabolites in tumor microenvironment is likely to be an effective way to inhibit tumor metastasis.

By correlating our spatial transcriptomics results to a public database, we confirmed that iCAFs was associated with immune infiltration and clinical outcomes in CRC. As main stromal cells in tumor microenvironment, CAFs can be detected in almost all solid tumors, thus targeting iCAFs can include the vast majority of patients, which is an ideal choice for colorectal cancer treatment.

There are still some limitations in this project. In this study, spatial transcriptome data were used to profile the heterogeneity of TME. Although this emerging technology can provide spatial location information of cells and facilitate the identification of cell types, due to the fact that 10X Genomics platform contains 1–10 cells in each spatial spot [38], the accuracy of spatial transcriptome data is lower than that of single-cell sequencing, so more efforts needed to improve the resolution of ST. Therefore, aiming to annotate cell types of each cluster more precisely, ssGSEA algorithm were applied in this project. Compared with morphological regions annotation, each cluster had been accurately annotated, which once again confirmed that ssGSEA algorithm is feasible for ST. Hence, the combination of ST and ssGSEA can provide a more comprehensive understanding of the cell types contained in the subpopulation, so as to discover the co-enriched cell types and the potential interactions between themselves. In addition, the small sample in our study may also be a possible limitation. And our project concerned on the functions of iCAFs in TME, the value of mCAFs was not fully recognized, so ongoing efforts are required to explore the functions of mCAFs.

Conclusions

In summary, our studies identified two different types of CAFs in CRC and provided an in-depth perspective on the interactions between CAFs and the tumor microenvironment. This study provides a practicable idea for drug therapy strategies targeting CRC in the future.

Availability of data and materials

Colorectal cancer spatial transcriptomics dataset analyzed in this research can be found at: http://www.cancerdiversity.asia/scCRLM/. Bulk transcriptomics of colorectal cancer comes from TCGA cohort COADREAD at: https://xena.ucsc.edu/.

Abbreviations

ST:

Spatial transcriptomics

CAFs:

Cancer-associated fibroblasts

TME:

Tumor microenvironment

mCAFs:

Myo-cancer-associated fibroblasts

iCAFs:

Inflammatory-cancer-associated fibroblasts

CRC:

Colorectal cancer

ssGSEA:

Single sample gene set enrichment analysis

ECM:

Extracellular matrix

NACT:

Neoadjuvant chemotherapy treatment

TCGA:

The Cancer Genome Atlas

DEGs:

Differentially expressed genes

GSVA:

Gene Set Variation Analysis

PCA:

Principal component analysis

EMT:

Epithelial mesenchymal transition

PR:

Partial response

TAM:

Tumor-associated macrophages

H&E:

Hematoxylin–eosin staining

IF:

Immunofluorescence

IHC:

Immunohistochemistry

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  Google Scholar 

  2. Peng S, Chen D, Cai J, Yuan Z, Huang B, Li Y, et al. Enhancing cancer-associated fibroblast fatty acid catabolism within a metabolically challenging tumor microenvironment drives colon cancer peritoneal metastasis. Mol Oncol. 2021;15(5):1391–411.

    Article  Google Scholar 

  3. Chen X, Song E. Turning foes to friends: targeting cancer-associated fibroblasts. Nat Rev Drug Discov. 2019;18(2):99–115.

    CAS  Article  Google Scholar 

  4. Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res. 2019;79(18):4557–66.

    CAS  Article  Google Scholar 

  5. Jin MZ, Jin WL. The updated landscape of tumor microenvironment and drug repurposing. Signal Transduct Target Ther. 2020;5(1):166.

    Article  Google Scholar 

  6. Elyada E, Bolisetty M, Laise P, Flynn WF, Courtois ET, Burkhart RA, et al. Cross-species single-cell analysis of pancreatic ductal adenocarcinoma reveals antigen-presenting cancer-associated fibroblasts. Cancer Discov. 2019;9(8):1102–23.

    CAS  Article  Google Scholar 

  7. Chen Z, Zhou L, Liu L, Hou Y, Xiong M, Yang Y, et al. Single-cell RNA sequencing highlights the role of inflammatory cancer-associated fibroblasts in bladder urothelial carcinoma. Nat Commun. 2020;11(1):5077.

    CAS  Article  Google Scholar 

  8. Kamali Zonouzi S, Pezeshki PS, Razi S, Rezaei N. Cancer-associated fibroblasts in colorectal cancer. Clin Transl Oncol. 2021;24(5):757–69.

    Article  Google Scholar 

  9. Wu Y, Yang S, Ma J, Chen Z, Song G, Rao D, et al. Spatiotemporal immune landscape of colorectal cancer liver metastasis at single-cell level. Cancer Discov. 2022;12(1):134–53.

    CAS  Article  Google Scholar 

  10. Hao Y, Hao S, Andersen-Nissen E, Mauck WM 3rd, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184(13):3573–87.

    CAS  Article  Google Scholar 

  11. Foroutan M, Bhuva DD, Lyu R, Horan K, Cursons J, Davis MJ. Single sample scoring of molecular phenotypes. BMC Bioinform. 2018;19(1):404.

    CAS  Article  Google Scholar 

  12. Salmén F, Vickovic S, Larsson L, Stenbeck L, Vallon-Christersson J, Ehinger A, et al. Multidimensional transcriptomics provides detailed information about immune cell distribution and identity in HER2+ breast tumors. bioRxiv. 2018.

  13. Stur E, Corvigno S, Xu M, Chen K, Tan Y, Lee S, et al. Spatially resolved transcriptomics of high-grade serous ovarian carcinoma. iScience. 2022;25(3):103923.

    CAS  Article  Google Scholar 

  14. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    CAS  Article  Google Scholar 

  15. Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. 2013;14:7.

    Article  Google Scholar 

  16. Qiu X, Mao Q, Tang Y, Wang L, Chawla R, Pliner HA, et al. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods. 2017;14(10):979–82.

    CAS  Article  Google Scholar 

  17. Patil I. Visualizations with statistical details: the “ggstatsplot” approach. J Open Sour Softw. 2021;6(61):3167.

    Article  Google Scholar 

  18. Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.

    CAS  Article  Google Scholar 

  19. Li H, Courtois ET, Sengupta D, Tan Y, Chen KH, Goh JJL, et al. Reference component analysis of single-cell transcriptomes elucidates cellular heterogeneity in human colorectal tumors. Nat Genet. 2017;49(5):708–18.

    CAS  Article  Google Scholar 

  20. Wang Y, Song W, Wang J, Wang T, Xiong X, Qi Z, et al. Single-cell transcriptome analysis reveals differential nutrient absorption functions in human intestine. J Exp Med. 2020;217(2):e20191130.

    Article  Google Scholar 

  21. Zhang L, Li Z, Skrzypczynska KM, Fang Q, Zhang W, O’Brien SA, et al. Single-cell analyses inform mechanisms of myeloid-targeted therapies in colon cancer. Cell. 2020;181(2):442–59.

    CAS  Article  Google Scholar 

  22. Zhang Y, Song J, Zhao Z, Yang M, Chen M, Liu C, et al. Single-cell transcriptome analysis reveals tumor immune microenvironment heterogenicity and granulocytes enrichment in colorectal cancer liver metastases. Cancer Lett. 2020;470:84–94.

    CAS  Article  Google Scholar 

  23. Mei Y, Xiao W, Hu H, Lu G, Chen L, Sun Z, et al. Single-cell analyses reveal suppressive tumor microenvironment of human colorectal cancer. Clin Transl Med. 2021;11(6): e422.

    CAS  Article  Google Scholar 

  24. Wang H, Gong P, Chen T, Gao S, Wu Z, Wang X, et al. Colorectal cancer stem cell states uncovered by simultaneous single-cell analysis of transcriptome and telomeres. Adv Sci. 2021;8(8):2004320.

    CAS  Article  Google Scholar 

  25. Pelaseyed T, Bergstrom JH, Gustafsson JK, Ermund A, Birchenough GM, Schutte A, et al. The mucus and mucins of the goblet cells and enterocytes provide the first defense line of the gastrointestinal tract and interact with the immune system. Immunol Rev. 2014;260(1):8–20.

    CAS  Article  Google Scholar 

  26. Han C, Liu T, Yin R. Biomarkers for cancer-associated fibroblasts. Biomark Res. 2020;8(1):64.

    Article  Google Scholar 

  27. Mao X, Xu J, Wang W, Liang C, Hua J, Liu J, et al. Crosstalk between cancer-associated fibroblasts and immune cells in the tumor microenvironment: new findings and future perspectives. Mol Cancer. 2021;20(1):131.

    CAS  Article  Google Scholar 

  28. Martinez-Reyes I, Chandel NS. Cancer metabolism: looking forward. Nat Rev Cancer. 2021;21(10):669–80.

    CAS  Article  Google Scholar 

  29. Zhang R, Qi F, Zhao F, Li G, Shao S, Zhang X, et al. Cancer-associated fibroblasts enhance tumor-associated macrophages enrichment and suppress NK cells function in colorectal cancer. Cell Death Dis. 2019;10(4):273.

    Article  Google Scholar 

  30. Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.

    CAS  Article  Google Scholar 

  31. Sebastian A, Hum NR, Martin KA, Gilmore SF, Peran I, Byers SW, et al. Single-cell transcriptomic analysis of tumor-derived fibroblasts and normal tissue-resident fibroblasts reveals fibroblast heterogeneity in breast cancer. Cancers. 2020;12(5):1307.

    CAS  Article  Google Scholar 

  32. Turley SJ, Cremasco V, Astarita JL. Immunological hallmarks of stromal cells in the tumour microenvironment. Nat Rev Immunol. 2015;15(11):669–82.

    CAS  Article  Google Scholar 

  33. Mosa MH, Michels BE, Menche C, Nicolas AM, Darvishi T, Greten FR, et al. A Wnt-induced phenotypic switch in cancer-associated fibroblasts inhibits EMT in colorectal cancer. Cancer Res. 2020;80(24):5569–82.

    CAS  Article  Google Scholar 

  34. Hu JL, Wang W, Lan XL, Zeng ZC, Liang YS, Yan YR, et al. CAFs secreted exosomes promote metastasis and chemotherapy resistance by enhancing cell stemness and epithelial-mesenchymal transition in colorectal cancer. Mol Cancer. 2019;18(1):91.

    CAS  Article  Google Scholar 

  35. Demircioglu F, Wang J, Candido J, Costa ASH, Casado P, de Luxan DB, et al. Cancer associated fibroblast FAK regulates malignant cell metabolism. Nat Commun. 2020;11(1):1290.

    CAS  Article  Google Scholar 

  36. Wu D, Zhuo L, Wang X. Metabolic reprogramming of carcinoma-associated fibroblasts and its impact on metabolic heterogeneity of tumors. Semin Cell Dev Biol. 2017;64:125–31.

    CAS  Article  Google Scholar 

  37. Gong J, Lin Y, Zhang H, Liu C, Cheng Z, Yang X, et al. Reprogramming of lipid metabolism in cancer-associated fibroblasts potentiates migration of colorectal cancer cells. Cell Death Dis. 2020;11(4):267.

    CAS  Article  Google Scholar 

  38. Ji AL, Rubin AJ, Thrane K, Jiang S, Reynolds DL, Meyers RM, et al. Multimodal analysis of composition and spatial architecture in human squamous cell carcinoma. Cell. 2020;182(2):497–514.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We really appreciate the Wu et al. [9] offered spatial transcriptomics data and the public database TCGA for allowing access to their data.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

ZP designed the study and carried out the data analyses. ZP, MY, HD, and ZF interpreted the results. ZP and KH draft the manuscript, MY, ZF, HD revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kongwang Hu.

Ethics declarations

Ethics approval and consent to participate

The data analyzed during the study come from the relevant studies where written informed consent was received prior to participation. All the tissue samples involved in this study came from previous paraffin-embedded surgical specimens derived from patients with colorectal cancer, informed consent has been signed and received. To protect personal privacy, all the personal details were not obtained. This study was approved by Clinical Medical Research Ethics Committee of First Affiliated Hospital of Anhui Medical University.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Figures S1–S3, Tables S1–S9.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Peng, Z., Ye, M., Ding, H. et al. Spatial transcriptomics atlas reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment components in colorectal cancer. J Transl Med 20, 302 (2022). https://doi.org/10.1186/s12967-022-03510-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03510-8

Keywords

  • Spatial transcriptomics
  • Colorectal cancer
  • Cancer-associated fibroblasts
  • Tumor microenvironment