Abstract
INTRODUCTION: Regeneration of fibrochondrocytes is essential for the healing of the tendon-bone interface (TBI), which is similar to the formation of neurogenic heterotopic ossification (HO). Through single-cell integrative analysis, this study explored the homogeneity of HO cells and fibrochondrocytes. METHODS: This study integrated six datasets, namely, GSE94683, GSE144306, GSE168153, GSE138515, GSE102929, and GSE110993. The differentiation trajectory and key transcription factors (TFs) for HO occurrence were systematically analyzed by integrating single-cell RNA (scRNA) sequencing, bulk RNA sequencing, and assay of transposase accessible chromatin seq. The differential expression and enrichment pathways of TFs in heterotopically ossified tissues were identified. RESULTS: HO that mimicked pathological cells was classified into HO1 and HO2 cell subsets. Results of the pseudo-temporal sequence analysis suggested that HO2 is a differentiated precursor cell of HO1. The analysis of integrated scRNA data revealed that ectopically ossified cells have similar transcriptional characteristics to cells in the fibrocartilaginous zone of tendons. The modified SCENIC method was used to identify specific transcriptional regulators associated with ectopic ossification. Xbp1 was defined as a common key transcriptional regulator of ectopically ossified tissues and the fibrocartilaginous zone of tendons. Subsequently, the CellPhoneDB database was completed for the cellular ligand-receptor analysis. With further pathway screening, this study is the first to propose that Xbp1 may upregulate the Notch signaling pathway through Jag1 transcription. Twenty-four microRNAs were screened and were found to be potentially associated with upregulation of XBP1 expression after acute ischemic stroke. CONCLUSION: A systematic analysis of the differentiation landscape and cellular homogeneity facilitated a molecular understanding of the phenotypic similarities between cells in the fibrocartilaginous region of tendon and HO cells. Furthermore, by identifying Xbp1 as a hub regulator and by conducting a ligand-receptor analysis, we propose a potential Xbp1/Jag1/Notch signaling pathway.