COL10A1 expression distinguishes a subset of cancer-associated fibroblasts present in the stroma of high-risk basal cell carcinoma

Basal cell carcinoma (BCC) is the most frequently diagnosed skin cancer and the most common malignancy in humans. Different morphological subtypes of BCC are associated with low- or high-risk of recurrence and aggressiveness, but the underlying biology of how the individual subtypes arise remains largely unknown. Because the majority of BCCs appear to arise from mutations in the same pathway, we hypothesized that BCC development, growth and invasive potential is also influenced by the tumor microenvironment and in particular by cancer-associated fibroblasts (CAFs) and their secreted factors. In this work, we aimed to characterize the stroma of the different BCC subtypes with a focus on CAF populations. To investigate the stromal features of the different BCC subtypes, we applied laser-capture microdissection (LCM) followed by RNA sequencing. A cohort of 15 BCC samples from 5 different “pure” subtypes (superficial, nodular, micronodular, sclerosing and basosquamous; n=3 each) were selected and included in the analysis. Healthy skin was used as a control (n=6). We show that the stroma of the different BCC subtypes have distinct gene expression signatures. Nodular and micronodular seem to have the most similar signatures, while superficial and sclerosing the most different. By comparing low- and high-risk BCC subtypes, we observed that Collagen 10A1 (COL10A1) is overexpressed in the stroma of sclerosing/infiltrative and basosquamous but not micronodular high-risk subtypes. Those findings were confirmed by immunohistochemistry in a cohort of 89 different BCC and 13 healthy skin samples. Moreover, scRNAseq analysis of BCCs of two independent datasets showed that the COL10A1-expressing population of cells is associated with the stroma adjacent to invasive BCC and shows extracellular matrix remodeling features.