Integrated transcriptomic analysis systematically reveals the heterogeneity and molecular characterization of cancer-associated fibroblasts in osteosarcoma

Abstract

BackgroundOsteosarcoma (OS), with a peak incidence during the adolescent growth spurt, is correlated with poor prognosis for its high malignancy. The tumor microenvironment (TME) is highly complicated, with frequent interactions between tumor and stromal cells. The cancer-associated fibroblasts (CAFs) in the TME have been considered to actively involve in the progression, metastasis, and drug resistance of OS. This study aimed to characterize cellular heterogeneity and molecular characterization in CAFs subtypes and explore the potential targeting therapeutic strategies to improve the prognosis of OS patients. MethodsThe single-cell atlas of human OS tumor lesions were constructed from the GEO database. Then significant marker genes and potential biological functions for each CAFs subtype were identified and explored using the Seurat R package. Next, by performing the survival analyses and constructing the risk scores for CAFs subtypes, we aimed to identify and characterize the prognostic values of specific marker genes and different CAFs subtypes. Furthermore, we explored the therapeutic targets and innovative drugs targeting different CAFs subtypes based on the GDSC database. Finally, prognoses related CAFs subtypes were further validated through immunohistochemistry (IHC) on clinical OS specimens. ResultsOverall, nine main cell clusters and five subtypes of CAFs were identified. The differentially expressed marker genes for each CAFs clusters were then identified. Moreover, through Gene Ontology (GO) enrichment analysis, we defined the CAFs_2 (upregulated CXCL14 and C3), which was closely related to leukocyte migration and chemotaxis, as inflammatory CAFs (iCAFs). Likewise, we defined the CAFs_4 (upregulated CD74, HLA-DRA and HLA-DRB1), which was closely related to antigen process and presentation, as antigen-presenting CAFs (apCAFs). Furthermore, Kaplan-Meier analyses showed that CAFs_2 and CAFs_4 were correlated with poor clinical prognosis of OS patients. Meanwhile, therapeutic drugs targeting CAFs_2 and CAFs_4, such as 17-AAG/Docetaxel/Bleomycin and PHA-793887/NG-25/KIN001-102, were also explored, respectively. Finally, IHC assay confirmed the abundant CAFs_2 and CAFs_4 subtypes infiltration in the OS microenvironment compared with adjacent tissues. ConclusionOur study revealed the diversity, complexity, and heterogeneity of CAFs in OS, and complemented the single-cell atlas in OS TME.