Abstract
BACKGROUND: With cancer-associated fibroblasts (CAFs) as the main cell type, the rich myxoid stromal components in chordoma tissues may likely contribute to its development and progression. METHODS: Single-cell RNA sequencing (scRNA-seq), spatial transcriptomics, bulk RNA-seq, and multiplexed quantitative immunofluorescence (QIF) were used to dissect the heterogeneity, spatial distribution, and clinical implication of CAFs in chordoma. RESULTS: We sequenced here 72 097 single cells from 3 primary and 3 recurrent tumor samples, as well as 3 nucleus pulposus samples as controls using scRNA-seq. We identified a unique cluster of CAF in recurrent tumors that highly expressed hypoxic genes and was functionally enriched in endoplasmic reticulum stress (ERS). Pseudotime trajectory and cell communication analyses showed that this ERS-CAF subpopulation originated from normal fibroblasts and widely interacted with tumoral and immune cells. Analyzing the bulk RNA-seq data from 126 patients, we found that the ERS-CAF signature score was associated with the invasion and poor prognosis of chordoma. By integrating the results of scRNA-seq with spatial transcriptomics, we demonstrated the existence of ERS-CAF in chordoma tissues and revealed that this CAF subtype displayed the most proximity to its surrounding tumor cells. In subsequent QIF validation involving 105 additional patients, we confirmed that ERS-CAF was abundant in the chordoma microenvironment and located close to tumor cells. Furthermore, both ERS-CAF density and its distance to tumor cells were correlated with tumor malignant phenotype and adverse patient outcomes. CONCLUSIONS: These findings depict the CAF landscape for chordoma and may provide insights into the development of novel treatment approaches.