Comprehensive analysis on subchondral bone marrow lesions of human osteoarthritis by integrating bulk and single-cell transcriptomes

Abstract

ObjectiveThis study aims to demonstrate the cellular composition and underlying mechanisms in subchondral bone marrow lesions (BMLs) of knee osteoarthritis (OA).MethodsBMLs were assessed by MRI Osteoarthritis Knee Score (MOAKS)≥2. Bulk RNA-sequencing (bulk-seq) and BML-specific differentially expressed genes (DEGs) analysis were performed among subchondral bone samples (including OA-BML=3, paired OA-NBML=3; non-OA=3). The hub genes of BMLs were identified by verifying in independent datasets and multiple bioinformatic analyses. To further estimate cell-type composition of subchondral bone, we utilized two newly developed deconvolution algorithms (MuSiC, MCP-counter) in transcriptomic datasets, based on signatures from open-accessed single-cell RNA sequencing (scRNA-seq). Finally, competing endogenous RNA (ceRNA) and transcription factor (TF) networks were constructed through multiple predictive databases, and validated by public non-coding RNA profiles.ResultsA total of 86 BML-specific DEGs (up 79, down 7) were identified. IL11 and VCAN were identified as core hub genes. The “has-miR-424-5p/lncRNA PVT1” was determined as crucial network, targeting IL11 and VCAN, respectively. More importantly, two deconvolution algorithms produced approximate estimations of cell-type composition, and the cluster of heterotopic-chondrocyte was discovered abundant in BMLs, and positively correlated with the expression of hub genes.ConclusionIL11 and VCAN were identified as the core hub genes of BMLs, and their molecular networks were determined as well. We profiled the characteristics of subchondral bone at single-cell level and determined that the heterotopic-chondrocyte was abundant in BMLs and was closely linked to IL11 and VCAN. Our study may provide new insights into the microenvironment and pathological molecular mechanism of BMLs, and could lead to novel therapeutic strategies.