Aims. Osteoarthritis (OA) is a common degenerative joint disease. The osteocyte transcriptome is highly relevant to osteocyte biology. This study aimed to explore the osteocyte transcriptome in subchondral bone affected by OA. Methods. Gene expression profiles of OA subchondral bone were used to identify disease-relevant genes and signalling pathways. RNA-sequencing data of a bone loading model were used to identify the loading-responsive gene set. Weighted gene co-expression network analysis (WGCNA) was employed to develop the osteocyte mechanics-responsive gene signature. Results. A group of 77 persistent genes that are highly relevant to extracellular matrix (ECM) biology and bone remodelling signalling were identified in OA subchondral lesions. A loading responsive gene set, including 446 principal genes, was highly enriched in OA medial tibial plateaus compared to lateral tibial plateaus. Of this gene set, a total of 223 genes were identified as the main contributors that were strongly associated with osteocyte functions and signalling pathways, such as ECM modelling, axon guidance, Hippo, Wnt, and transforming growth factor beta (TGF-β) signalling pathways. We limited the loading-responsive genes obtained via the osteocyte transcriptome signature to identify a subgroup of genes that are highly relevant to osteocytes, as the mechanics-responsive osteocyte signature in OA. Based on WGCNA, we found that this signature was highly co-expressed and identified three clusters, including early, late, and persistently responsive genes. Conclusion. In this study, we identified the mechanics-responsive osteocyte signature in OA-lesioned subchondral bone. Cite this article: Bone Joint Res 2022;11(6):362–370

Aims. Rheumatoid arthritis (RA) is a systematic autoimmune disorder, characterized by synovial inflammation, bone and cartilage destruction, and disease involvement in multiple organs. Although numerous drugs are employed in RA treatment, some respond little and suffer from severe side effects. This study aimed to screen the candidate therapeutic targets and promising drugs in a novel method. Methods. We developed a module-based and cumulatively scoring approach that is a deeper-layer application of weighted gene co-expression network (WGCNA) and connectivity map (CMap) based on the high-throughput datasets. Results. Four noteworthy RA-related modules were identified, revealing the immune- and infection-related biological processes and pathways involved in RA. HLA-DMA, HLA-DMB, HLA-DPA1, HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DRB1, BLNK, BTK, CD3D, CD4, IL2RG, INPP5D, LCK, PTPRC, RAC2, SYK, and VAV1 were recognized as the key hub genes with high connectivity in gene regulation networks and gene pathway networks. Moreover, the long noncoding RNAs (lncRNAs) in the RA-related modules, such as FAM30A and NEAT1, were identified as the indispensable interactors with the hub genes. Finally, candidate drugs were screened by developing a cumulatively scoring approach based on the selected modules. Niclosamide and the other compounds of T-type calcium channel blocker, IKK inhibitor, and PKC activator, HIF activator, and proteasome inhibitor, which harbour the similar gene signature with niclosamide, were promising drugs with high specificity and broad coverage for the RA-related modules. Conclusion. This study provides not only the promising targets and drugs for RA but also a novel methodological insight into the target and drug screening. Cite this article: Bone Joint Res 2020;9(8):501–514