Aims. We aimed to develop a gene signature that predicts the occurrence of postmenopausal osteoporosis (PMOP) by studying its genetic mechanism. Methods. Five datasets were obtained from the Gene Expression Omnibus database. Unsupervised consensus cluster analysis was used to determine new PMOP subtypes. To determine the central genes and the core modules related to PMOP, the weighted gene co-expression network analysis (WCGNA) was applied. Gene Ontology enrichment analysis was used to explore the biological processes underlying key genes. Logistic regression univariate analysis was used to screen for statistically significant variables. Two algorithms were used to select important PMOP-related genes. A logistic regression model was used to construct the PMOP-related gene profile. The receiver operating characteristic area under the curve, Harrell’s concordance index, a calibration chart, and decision curve analysis were used to characterize PMOP-related genes. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the expression of the PMOP-related genes in the gene signature. Results. We identified three PMOP-related subtypes and four core modules. The muscle system process, muscle contraction, and actin filament-based movement were more active in the hub genes. We obtained five feature genes related to PMOP. Our analysis verified that the gene signature had good predictive power and applicability. The outcomes of the GSE56815 cohort were found to be consistent with the results of the earlier studies. qRT-PCR results showed that RAB2A and FYCO1 were amplified in clinical samples. Conclusion. The PMOP-related gene signature we developed and verified can accurately predict the risk of PMOP in patients. These results can elucidate the molecular mechanism of RAB2A and FYCO1 underlying PMOP, and yield new and improved treatment strategies, ultimately helping PMOP monitoring. Cite this article: Bone Joint Res 2022;11(8):548–560

Aims. Lumbar spinal stenosis (LSS) is a common skeletal system disease that has been partly attributed to genetic variation. However, the correlation between genetic variation and pathological changes in LSS is insufficient, and it is difficult to provide a reference for the early diagnosis and treatment of the disease. Methods. We conducted a transcriptome-wide association study (TWAS) of spinal canal stenosis by integrating genome-wide association study summary statistics (including 661 cases and 178,065 controls) derived from Biobank Japan, and pre-computed gene expression weights of skeletal muscle and whole blood implemented in FUSION software. To verify the TWAS results, the candidate genes were furthered compared with messenger RNA (mRNA) expression profiles of LSS to screen for common genes. Finally, Metascape software was used to perform enrichment analysis of the candidate genes and common genes. Results. TWAS identified 295 genes with permutation p-values < 0.05 for skeletal muscle and 79 genes associated for the whole blood, such as RCHY1 (PTWAS = 0.001). Those genes were enriched in 112 gene ontology (GO) terms and five Kyoto Encyclopedia of Genes and Genomes pathways, such as ‘chemical carcinogenesis - reactive oxygen species’ (LogP value = −2.139). Further comparing the TWAS significant genes with the differentially expressed genes identified by mRNA expression profiles of LSS found 18 overlapped genes, such as interleukin 15 receptor subunit alpha (IL15RA) (PTWAS = 0.040, PmRNA = 0.010). Moreover, 71 common GO terms were detected for the enrichment results of TWAS and mRNA expression profiles, such as negative regulation of cell differentiation (LogP value = −2.811). Conclusion. This study revealed the genetic mechanism behind the pathological changes in LSS, and may provide novel insights for the early diagnosis and intervention of LSS. Cite this article: Bone Joint Res 2023;12(6):387–396

Aims. Rotator cuff tear (RCT) is the leading cause of shoulder pain, primarily associated with age-related tendon degeneration. This study aimed to elucidate the potential differential gene expressions in tendons across different age groups, and to investigate their roles in tendon degeneration. Methods. Linear regression and differential expression (DE) analyses were performed on two transcriptome profiling datasets of torn supraspinatus tendons to identify age-related genes. Subsequent functional analyses were conducted on these candidate genes to explore their potential roles in tendon ageing. Additionally, a secondary DE analysis was performed on candidate genes by comparing their expressions between lesioned and normal tendons to explore their correlations with RCTs. Results. We identified 49 genes in torn supraspinatus tendons associated with advancing age. Among them, five age-related genes showed DE in lesioned tendons compared to normal tendons. Functional analyses and previous studies have highlighted their specific enrichments in biological functions, such as muscle development (e.g. myosin heavy chain 3 (MYH3)), transcription regulation (e.g. CCAAT enhancer binding brotein delta (CEBPD)), and metal ion homeostasis (e.g. metallothionein 1X (MT1X)). Conclusion. This study uncovered molecular aspects of tendon ageing and their potential links to RCT development, offering insights for targeted interventions. These findings enhance our understanding of the mechanisms of tendon degeneration, allowing potential strategies to be made for reducing the incidence of RCT. Cite this article: Bone Joint Res 2024;13(9):474–484

Aims. Degenerative cervical spondylosis (DCS) is a common musculoskeletal disease that encompasses a wide range of progressive degenerative changes and affects all components of the cervical spine. DCS imposes very large social and economic burdens. However, its genetic basis remains elusive. Methods. Predicted whole-blood and skeletal muscle gene expression and genome-wide association study (GWAS) data from a DCS database were integrated, and functional summary-based imputation (FUSION) software was used on the integrated data. A transcriptome-wide association study (TWAS) was conducted using FUSION software to assess the association between predicted gene expression and DCS risk. The TWAS-identified genes were verified via comparison with differentially expressed genes (DEGs) in DCS RNA expression profiles in the Gene Expression Omnibus (GEO) (Accession Number: GSE153761). The Functional Mapping and Annotation (FUMA) tool for genome-wide association studies and Meta tools were used for gene functional enrichment and annotation analysis. Results. The TWAS detected 420 DCS genes with p < 0.05 in skeletal muscle, such as ribosomal protein S15A (RPS15A) (PTWAS = 0.001), and 110 genes in whole blood, such as selectin L (SELL) (PTWAS = 0.001). Comparison with the DCS RNA expression profile identified 12 common genes, including Apelin Receptor (APLNR) (PTWAS = 0.001, PDEG = 0.025). In total, 148 DCS-enriched Gene Ontology (GO) terms were identified, such as mast cell degranulation (GO:0043303); 15 DCS-enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, such as the sphingolipid signalling pathway (ko04071). Nine terms, such as degradation of the extracellular matrix (R-HSA-1474228), were common to the TWAS enrichment results and the RNA expression profile. Conclusion. Our results identify putative susceptibility genes; these findings provide new ideas for exploration of the genetic mechanism of DCS development and new targets for preclinical intervention and clinical treatment. Cite this article: Bone Joint Res 2023;12(1):80–90

Aims. This study aimed to explore the biological and clinical importance of dysregulated key genes in osteoarthritis (OA) patients at the cartilage level to find potential biomarkers and targets for diagnosing and treating OA. Methods. Six sets of gene expression profiles were obtained from the Gene Expression Omnibus database. Differential expression analysis, weighted gene coexpression network analysis (WGCNA), and multiple machine-learning algorithms were used to screen crucial genes in osteoarthritic cartilage, and genome enrichment and functional annotation analyses were used to decipher the related categories of gene function. Single-sample gene set enrichment analysis was performed to analyze immune cell infiltration. Correlation analysis was used to explore the relationship among the hub genes and immune cells, as well as markers related to articular cartilage degradation and bone mineralization. Results. A total of 46 genes were obtained from the intersection of significantly upregulated genes in osteoarthritic cartilage and the key module genes screened by WGCNA. Functional annotation analysis revealed that these genes were closely related to pathological responses associated with OA, such as inflammation and immunity. Four key dysregulated genes (cartilage acidic protein 1 (CRTAC1), iodothyronine deiodinase 2 (DIO2), angiopoietin-related protein 2 (ANGPTL2), and MAGE family member D1 (MAGED1)) were identified after using machine-learning algorithms. These genes had high diagnostic value in both the training cohort and external validation cohort (receiver operating characteristic > 0.8). The upregulated expression of these hub genes in osteoarthritic cartilage signified higher levels of immune infiltration as well as the expression of metalloproteinases and mineralization markers, suggesting harmful biological alterations and indicating that these hub genes play an important role in the pathogenesis of OA. A competing endogenous RNA network was constructed to reveal the underlying post-transcriptional regulatory mechanisms. Conclusion. The current study explores and validates a dysregulated key gene set in osteoarthritic cartilage that is capable of accurately diagnosing OA and characterizing the biological alterations in osteoarthritic cartilage; this may become a promising indicator in clinical decision-making. This study indicates that dysregulated key genes play an important role in the development and progression of OA, and may be potential therapeutic targets. Cite this article: Bone Joint Res 2024;13(2):66–82

Aims. Knee osteoarthritis (OA) involves a variety of tissues in the joint. Gene expression profiles in different tissues are of great importance in order to understand OA. Methods. First, we obtained gene expression profiles of cartilage, synovium, subchondral bone, and meniscus from the Gene Expression Omnibus (GEO). Several datasets were standardized by merging and removing batch effects. Then, we used unsupervised clustering to divide OA into three subtypes. The gene ontology and pathway enrichment of three subtypes were analyzed. CIBERSORT was used to evaluate the infiltration of immune cells in different subtypes. Finally, OA-related genes were obtained from the Molecular Signatures Database for validation, and diagnostic markers were screened according to clinical characteristics. Quantitative reverse transcription polymerase chain reaction (qRT‐PCR) was used to verify the effectiveness of markers. Results. C1 subtype is mainly concentrated in the development of skeletal muscle organs, C2 lies in metabolic process and immune response, and C3 in pyroptosis and cell death process. Therefore, we divided OA into three subtypes: bone remodelling subtype (C1), immune metabolism subtype (C2), and cartilage degradation subtype (C3). The number of macrophage M0 and activated mast cells of C2 subtype was significantly higher than those of the other two subtypes. COL2A1 has significant differences in different subtypes. The expression of COL2A1 is related to age, and trafficking protein particle complex subunit 2 is related to the sex of OA patients. Conclusion. This study linked different tissues with gene expression profiles, revealing different molecular subtypes of patients with knee OA. The relationship between clinical characteristics and OA-related genes was also studied, which provides a new concept for the diagnosis and treatment of OA. Cite this article: Bone Joint Res 2023;12(12):702–711

Aims. This study explored the shared genetic traits and molecular interactions between postmenopausal osteoporosis (POMP) and sarcopenia, both of which substantially degrade elderly health and quality of life. We hypothesized that these motor system diseases overlap in pathophysiology and regulatory mechanisms. Methods. We analyzed microarray data from the Gene Expression Omnibus (GEO) database using weighted gene co-expression network analysis (WGCNA), machine learning, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to identify common genetic factors between POMP and sarcopenia. Further validation was done via differential gene expression in a new cohort. Single-cell analysis identified high expression cell subsets, with mononuclear macrophages in osteoporosis and muscle stem cells in sarcopenia, among others. A competitive endogenous RNA network suggested regulatory elements for these genes. Results. Signal transducer and activator of transcription 3 (STAT3) was notably expressed in both conditions. Single-cell analysis pinpointed specific cells with high STAT3 expression, and microRNA (miRNA)-125a-5p emerged as a potential regulator. Experiments confirmed the crucial role of STAT3 in osteoclast differentiation and muscle proliferation. Conclusion. STAT3 has emerged as a key gene in both POMP and sarcopenia. This insight positions STAT3 as a potential common therapeutic target, possibly improving management strategies for these age-related diseases. Cite this article: Bone Joint Res 2024;13(8):411–426

Aims. Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. Methods. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Results. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p. discovery GWAS. = 1.21 × 10. -25. , p. replication GWAS. = 1.80 × 10. -12. ), CCDC170 (p. discovery GWAS. = 1.23 × 10. -11. , p. replication GWAS. = 3.22 × 10. -9. ), and SOX6 (p. discovery GWAS. = 4.41 × 10. -15. , p. replication GWAS. = 6.57 × 10. -14. ). Then, for the 36 common target genes, multiple gene ontology (GO) terms were detected for BMD such as positive regulation of cartilage development (p = 9.27 × 10. -3. ) and positive regulation of chondrocyte differentiation (p = 9.27 × 10. -3. ). Conclusion. We explored the potential roles of rSNP in the genetic mechanisms of BMD and identified multiple candidate genes. Our study results support the implication of regulatory genetic variants in the development of OP. Cite this article: Bone Joint Res 2023;12(2):147–154

Aims. Impaired fracture repair in patients with type 2 diabetes mellitus (T2DM) is not fully understood. In this study, we aimed to characterize the local changes in gene expression (GE) associated with diabetic fracture. We used an unbiased approach to compare GE in the fracture callus of Zucker diabetic fatty (ZDF) rats relative to wild-type (WT) littermates at three weeks following femoral osteotomy. Methods. Zucker rats, WT and homozygous for leptin receptor mutation (ZDF), were fed a moderately high-fat diet to induce T2DM only in the ZDF animals. At ten weeks of age, open femoral fractures were simulated using a unilateral osteotomy stabilized with an external fixator. At three weeks post-surgery, the fractured femur from each animal was retrieved for analysis. Callus formation and the extent of healing were assessed by radiograph and histology. Bone tissue was processed for total RNA extraction and messenger RNA (mRNA) sequencing (mRNA-Seq). Results. Radiographs and histology demonstrated impaired fracture healing in ZDF rats with incomplete bony bridge formation and an influx of intramedullary inflammatory tissue. In comparison, near-complete bridging between cortices was observed in Sham WT animals. Of 13,160 genes, mRNA-Seq analysis identified 13 that were differentially expressed in ZDF rat callus, using a false discovery rate (FDR) threshold of 10%. Seven genes were upregulated with high confidence (FDR = 0.05) in ZDF fracture callus, most with known roles in inflammation. Conclusion. These findings suggest that elevated or prolonged inflammation contributes to delayed fracture healing in T2DM. The identified genes may be used as biomarkers to monitor and treat delayed fracture healing in diabetic patients. Cite this article: Bone Joint Res 2023;12(10):657–666

Aims. The aim of this study was to screen the entire genome for genetic markers associated with risk for anterior cruciate ligament (ACL) and posterior cruciate ligament (PCL) injury. Methods. Genome-wide association (GWA) analyses were performed using data from the Kaiser Permanente Research Board (KPRB) and the UK Biobank. ACL and PCL injury cases were identified based on electronic health records from KPRB and the UK Biobank. GWA analyses from both cohorts were tested for ACL and PCL injury using a logistic regression model adjusting for sex, height, weight, age at enrolment, and race/ethnicity using allele counts for single nucleotide polymorphisms (SNPs). The data from the two GWA studies were combined in a meta-analysis. Candidate genes previously reported to show an association with ACL injury in athletes were also tested for association from the meta-analysis data from the KPRB and the UK Biobank GWA studies. Results. There was a total of 2,214 cases of ACL and PCL injury and 519,869 controls within the two cohorts, with three loci demonstrating a genome-wide significant association in the meta-analysis: INHBA, AEBP2, and LOC101927869. Of the eight candidate genes previously studied in the literature, six were present in the current dataset, and only COL3A1 (rs1800255) showed a significant association (p = 0.006). Conclusion. Genetic markers in three novel loci in this study and one previously-studied candidate gene were identified as potential risk factors for ACL and PCL injury and deserve further validation and investigation of molecular mechanisms. Cite this article: Bone Jt Open 2021;2(6):414–421

Objectives. Osteoporosis is a chronic disease. The aim of this study was to identify key genes in osteoporosis. Methods. Microarray data sets GSE56815 and GSE56814, comprising 67 osteoporosis blood samples and 62 control blood samples, were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified in osteoporosis using Limma package (3.2.1) and Meta-MA packages. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to identify biological functions. Furthermore, the transcriptional regulatory network was established between the top 20 DEGs and transcriptional factors using the UCSC ENCODE Genome Browser. Receiver operating characteristic (ROC) analysis was applied to investigate the diagnostic value of several DEGs. Results. A total of 1320 DEGs were obtained, of which 855 were up-regulated and 465 were down-regulated. These differentially expressed genes were enriched in Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways, mainly associated with gene expression and osteoclast differentiation. In the transcriptional regulatory network, there were 6038 interactions pairs involving 88 transcriptional factors. In addition, the quantitative reverse transcriptase-polymerase chain reaction result validated the expression of several genes (VPS35, FCGR2A, TBCA, HIRA, TYROBP, and JUND). Finally, ROC analyses showed that VPS35, HIRA, PHF20 and NFKB2 had a significant diagnostic value for osteoporosis. Conclusion. Genes such as VPS35, FCGR2A, TBCA, HIRA, TYROBP, JUND, PHF20, NFKB2, RPL35A and BICD2 may be considered to be potential pathogenic genes of osteoporosis and may be useful for further study of the mechanisms underlying osteoporosis. Cite this article: B. Xia, Y. Li, J. Zhou, B. Tian, L. Feng. Identification of potential pathogenic genes associated with osteoporosis. Bone Joint Res 2017;6:640–648. DOI: 10.1302/2046-3758.612.BJR-2017-0102.R1

Aims. Osteoarthritis (OA) is the most prevalent joint disease. However, the specific and definitive genetic mechanisms of OA are still unclear. Methods. Tissue-related transcriptome-wide association studies (TWAS) of hip OA and knee OA were performed utilizing the genome-wide association study (GWAS) data of hip OA and knee OA (including 2,396 hospital-diagnosed hip OA patients versus 9,593 controls, and 4,462 hospital-diagnosed knee OA patients versus 17,885 controls) and gene expression reference to skeletal muscle and blood. The OA-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the messenger RNA (mRNA) expression profiles of hip OA and knee OA. Functional enrichment and annotation analysis of identified genes was performed by the DAVID and FUMAGWAS tools. Results. We detected 33 common genes, eight common gene ontology (GO) terms, and one common pathway for hip OA, such as calcium and integrin-binding protein 1 (CIB1) (PTWAS = 0.025, FCmRNA = -1.575 for skeletal muscle), adrenomedullin (ADM) (PTWAS = 0.022, FCmRNA = -4.644 for blood), Golgi apparatus (PTWAS <0.001, PmRNA = 0.012 for blood), and phosphatidylinositol 3' -kinase-protein kinase B (PI3K-Akt) signalling pathway (PTWAS = 0.033, PmRNA = 0.005 for blood). For knee OA, we detected 24 common genes, eight common GO terms, and two common pathways, such as histocompatibility complex, class II, DR beta 1 (HLA-DRB1) (PTWAS = 0.040, FCmRNA = 4.062 for skeletal muscle), Follistatin-like 1 (FSTL1) (PTWAS = 0.048, FCmRNA = 3.000 for blood), cytoplasm (PTWAS < 0.001, PmRNA = 0.005 for blood), and complement and coagulation cascades (PTWAS = 0.017, PmRNA = 0.001 for skeletal muscle). Conclusion. We identified a group of OA-associated genes and pathways, providing novel clues for understanding the genetic mechanism of OA. Cite this article:Bone Joint Res. 2020;9(3):130–138

Objectives. In order to screen the altered gene expression profile in peripheral blood mononuclear cells of patients with osteoporosis, we performed an integrated analysis of the online microarray studies of osteoporosis. Methods. We searched the Gene Expression Omnibus (GEO) database for microarray studies of peripheral blood mononuclear cells in patients with osteoporosis. Subsequently, we integrated gene expression data sets from multiple microarray studies to obtain differentially expressed genes (DEGs) between patients with osteoporosis and normal controls. Gene function analysis was performed to uncover the functions of identified DEGs. Results. A total of three microarray studies were selected for integrated analysis. In all, 1125 genes were found to be significantly differentially expressed between osteoporosis patients and normal controls, with 373 upregulated and 752 downregulated genes. Positive regulation of the cellular amino metabolic process (gene ontology (GO): 0033240, false discovery rate (FDR) = 1.00E + 00) was significantly enriched under the GO category for biological processes, while for molecular functions, flavin adenine dinucleotide binding (GO: 0050660, FDR = 3.66E-01) and androgen receptor binding (GO: 0050681, FDR = 6.35E-01) were significantly enriched. DEGs were enriched in many osteoporosis-related signalling pathways, including those of mitogen-activated protein kinase (MAPK) and calcium. Protein-protein interaction (PPI) network analysis showed that the significant hub proteins contained ubiquitin specific peptidase 9, X-linked (Degree = 99), ubiquitin specific peptidase 19 (Degree = 57) and ubiquitin conjugating enzyme E2 B (Degree = 57). Conclusion. Analysis of gene function of identified differentially expressed genes may expand our understanding of fundamental mechanisms leading to osteoporosis. Moreover, significantly enriched pathways, such as MAPK and calcium, may involve in osteoporosis through osteoblastic differentiation and bone formation. Cite this article: J. J. Li, B. Q. Wang, Q. Fei, Y. Yang, D. Li. Identification of candidate genes in osteoporosis by integrated microarray analysis. Bone Joint Res 2016;5:594–601. DOI: 10.1302/2046-3758.512.BJR-2016-0073.R1

Objectives. The molecular mechanism of rheumatoid arthritis (RA) remains elusive. We conducted a protein-protein interaction network-based integrative analysis of genome-wide association studies (GWAS) and gene expression profiles of RA. Methods. We first performed a dense search of RA-associated gene modules by integrating a large GWAS meta-analysis dataset (containing 5539 RA patients and 20 169 healthy controls), protein interaction network and gene expression profiles of RA synovium and peripheral blood mononuclear cells (PBMCs). Gene ontology (GO) enrichment analysis was conducted by DAVID. The protein association networks of gene modules were generated by STRING. Results. For RA synovium, the top-ranked gene module is HLA-A, containing TAP2, HLA-A, HLA-C, TAPBP and LILRB1 genes. For RA PBMCs, the top-ranked gene module is GRB7, consisting of HLA-DRB5, HLA-DRA, GRB7, CD63 and KIT genes. Functional enrichment analysis identified three significant GO terms for RA synovium, including antigen processing and presentation of peptide antigen via major histocompatibility complex class I (false discovery rate (FDR) = 4.86 × 10 – 4), antigen processing and presentation of peptide antigen (FDR = 2.33 × 10 – 3) and eukaryotic translation initiation factor 4F complex (FDR = 2.52 × 10 – 2). Conclusion. This study reported several RA-associated gene modules and their functional association networks. Cite this article: X. Xiao, J. Hao, Y. Wen, W. Wang, X. Guo, F. Zhang. Genome-wide association studies and gene expression profiles of rheumatoid arthritis: an analysis. Bone Joint Res 2016;5:314–319. DOI: 10.1302/2046-3758.57.2000502

Objectives. The aim of this study was to identify key pathological genes in osteoarthritis (OA). Methods. We searched and downloaded mRNA expression data from the Gene Expression Omnibus database to identify differentially expressed genes (DEGs) of joint synovial tissues from OA and normal individuals. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analyses were used to assess the function of identified DEGs. The protein-protein interaction (PPI) network and transcriptional factors (TFs) regulatory network were used to further explore the function of identified DEGs. The quantitative real-time polymerase chain reaction (qRT-PCR) was applied to validate the result of bioinformatics analysis. Electronic validation was performed to verify the expression of selected DEGs. The diagnosis value of identified DEGs was accessed by receiver operating characteristic (ROC) analysis. Results. A total of 1085 DEGs were identified. KEGG pathway analysis displayed that Wnt was a significantly enriched signalling pathway. Some hub genes with high interactions such as USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, and ZNF423 were identified in the PPI and TFs network. The results of qRT-PCR showed that GADD45B, ADAMTS1, and TFAM were down-regulated in joint synovial tissues of OA, which was consistent with the bioinformatics analysis. The expression levels of USP46, CPVL, FOSL2, and PTGS1 in electronic validation were compatible with the bio-informatics result. CPVL and TFAM had a potential diagnostic value for OA based on the ROC analysis. Conclusion. The deregulated genes including USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, ZNF423, ADAMTS1, and TFAM might be involved in the pathology of OA. Cite this article: X. Zhang, Y. Bu, B. Zhu, Q. Zhao, Z. Lv, B. Li, J. Liu. Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis. Bone Joint Res 2018;7:298–307. DOI: 10.1302/2046-3758.74.BJR-2017-0245.R1

Objectives. Adult mice lacking the transcription factor NFAT1 exhibit osteoarthritis (OA). The precise molecular mechanism for NFAT1 deficiency-induced osteoarthritic cartilage degradation remains to be clarified. This study aimed to investigate if NFAT1 protects articular cartilage (AC) against OA by directly regulating the transcription of specific catabolic and anabolic genes in articular chondrocytes. Methods. Through a combined approach of gene expression analysis and web-based searching of NFAT1 binding sequences, 25 candidate target genes that displayed aberrant expression in Nfat1. -/-. AC at the initiation stage of OA, and possessed at least four NFAT1 binding sites in the promoter of each gene, were selected and tested for NFAT1 transcriptional activities by chromatin immunoprecipitation (ChIP) and promoter luciferase reporter assays using chondrocytes isolated from the AC of three- to four-month-old wild-type mice or Nfat1. -/-. mice with early OA phenotype. Results. Chromatin immunoprecipitation assays revealed that NFAT1 bound directly to the promoter of 21 of the 25 tested genes encoding cartilage-matrix proteins, growth factors, inflammatory cytokines, matrix-degrading proteinases, and specific transcription factors. Promoter luciferase reporter assays of representative anabolic and catabolic genes demonstrated that NFAT1-DNA binding functionally regulated the luciferase activity of specific target genes in wild-type chondrocytes, but not in Nfat1. -/-. chondrocytes or in wild-type chondrocytes transfected with plasmids containing mutated NFAT1 binding sequences. Conclusion. NFAT1 protects AC against degradation by directly regulating the transcription of target genes in articular chondrocytes. NFAT1 deficiency causes defective transcription of specific anabolic and catabolic genes in articular chondrocytes, leading to increased matrix catabolism and osteoarthritic cartilage degradation. Cite this article: M. Zhang, Q. Lu, T. Budden, J. Wang. NFAT1 protects articular cartilage against osteoarthritic degradation by directly regulating transcription of specific anabolic and catabolic genes. Bone Joint Res 2019;8:90–100. DOI: 10.1302/2046-3758.82.BJR-2018-0114.R1

Objectives. We have previously investigated an association between the genome copy number variation (CNV) and acetabular dysplasia (AD). Hip osteoarthritis is associated with a genetic polymorphism in the aspartic acid repeat in the N-terminal region of the asporin (ASPN) gene; therefore, the present study aimed to investigate whether the CNV of ASPN is involved in the pathogenesis of AD. Methods. Acetabular coverage of all subjects was evaluated using radiological findings (Sharp angle, centre-edge (CE) angle, acetabular roof obliquity (ARO) angle, and minimum joint space width). Genomic DNA was extracted from peripheral blood leukocytes. Agilent’s region-targeted high-density oligonucleotide tiling microarray was used to analyse 64 female AD patients and 32 female control subjects. All statistical analyses were performed using EZR software (Fisher’s exact probability test, Pearson’s correlation test, and Student’s t-test). Results. CNV analysis of the ASPN gene revealed a copy number loss in significantly more AD patients (9/64) than control subjects (0/32; p = 0.0212). This loss occurred within a 60 kb region on 9q22.31, which harbours the gene for ASPN. The mean radiological parameters of these AD patients were significantly worse than those of the other subjects (Sharp angle, p = 0.0056; CE angle, p = 0.0076; ARO angle, p = 0.0065), and all nine patients required operative therapy such as total hip arthroplasty or pelvic osteotomy. Moreover, six of these nine patients had a history of operative or conservative therapy for developmental dysplasia of the hip. Conclusions. Copy number loss within the region harbouring the ASPN gene on 9q22.31 is associated with severe AD. A copy number loss in the ASPN gene region may play a role in the aetiology of severe AD. Cite this article: T. Sekimoto, M. Ishii, M. Emi, S. Kurogi, T. Funamoto, Y. Yonezawa, T. Tajima, T. Sakamoto, H. Hamada, E. Chosa. Copy number loss in the region of the ASPN gene in patients with acetabular dysplasia: ASPN CNV in acetabular dysplasia. Bone Joint Res 2017;6:439–445. DOI: 10.1302/2046-3758.67.BJR-2016-0094.R1

Objectives. X-linked hypophosphataemic rickets (XLHR) is a disease of impaired bone mineralization characterized by hypophosphataemia caused by renal phosphate wasting. The main clinical manifestations of the disorder are O-shaped legs, X-shaped legs, delayed growth, and bone pain. XLHR is the most common inheritable form of rickets, with an incidence of 1/20 000 in humans. It accounts for approximately 80% of familial cases of hypophosphataemia and serves as the prototype of defective tubular phosphate (PO4. 3+. ) transport, due to extra renal defects resulting in unregulated FGF23 activity. XLHR is caused by loss-of-function mutations in the PHEX gene. The aim of this research was to identify the genetic defect responsible for familial hypophosphataemic rickets in a four-generation Chinese Han pedigree and to analyze the function of this mutation. Methods. The genome DNA samples of all members in the pedigree were extracted from whole blood. We sequenced all exons of the PHEX and FGF23 genes, as well as the adjacent splice site sequence with Sanger sequencing. Next, we analyzed the de novo mutation c.1692 del A of the PHEX gene with an online digital service and investigated the mutant PHEX with SWISS-MODEL, immunofluorescence, and protein stability detection. Results. Through Sanger sequencing, we found a de novo mutation, c.1692 del A, in exon 16 of the PHEX gene in this pedigree. This mutation can make the PHEX protein become unstable and decay rapidly, which results in familial XLHR. Conclusion. We have found a de novo loss-of-function mutation, c.1692 del A, in exon 16 of the PHEX gene that can cause XLHR. Cite this article: J. Huang, X. Bao, W. Xia, L. Zhu, J. Zhang, J. Ma, N. Jiang, J. Yang, Q. Chen, T. Jing, J. Liu, D. Ma, G. Xu. Functional analysis of a de novo mutation c.1692 del A of the PHEX gene in a Chinese family with X-linked hypophosphataemic rickets. Bone Joint Res 2019;8:405–413. DOI: 10.1302/2046-3758.88.BJR-2018-0276.R1

Aim. Osteoarthritis (OA) is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control the expression of genes and are likely to regulate the OA transcriptome. We performed integrative genomic analyses to define methylation-gene expression relationships in osteoarthritic cartilage. Patients and Methods. Genome-wide DNA methylation profiling of articular cartilage from five patients with OA of the knee and five healthy controls was conducted using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, California). Other independent genome-wide mRNA expression profiles of articular cartilage from three patients with OA and three healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Integrative pathway enrichment analysis of DNA methylation and mRNA expression profiles was performed using integrated analysis of cross-platform microarray and pathway software. Gene ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Results. We identified 1265 differentially methylated genes, of which 145 are associated with significant changes in gene expression, such as DLX5, NCOR2 and AXIN2 (all p-values of both DNA methylation and mRNA expression < 0.05). Pathway enrichment analysis identified 26 OA-associated pathways, such as mitogen-activated protein kinase (MAPK) signalling pathway (p = 6.25 × 10-4), phosphatidylinositol (PI) signalling system (p = 4.38 × 10-3), hypoxia-inducible factor 1 (HIF-1) signalling pathway (p = 8.63 × 10-3 pantothenate and coenzyme A (CoA) biosynthesis (p = 0.017), ErbB signalling pathway (p = 0.024), inositol phosphate (IP) metabolism (p = 0.025), and calcium signalling pathway (p = 0.032). Conclusion. We identified a group of genes and biological pathwayswhich were significantly different in both DNA methylation and mRNA expression profiles between patients with OA and controls. These results may provide new clues for clarifying the mechanisms involved in the development of OA. Cite this article: A. He, Y. Ning, Y. Wen, Y. Cai, K. Xu, Y. Cai, J. Han, L. Liu, Y. Du, X. Liang, P. Li, Q. Fan, J. Hao, X. Wang, X. Guo, T. Ma, F. Zhang. Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage. Bone Joint Res 2018;7:343–350. DOI: 10.1302/2046-3758.75.BJR-2017-0284.R1

Objectives. Several genome-wide association studies (GWAS) of bone mineral density (BMD) have successfully identified multiple susceptibility genes, yet isolated susceptibility genes are often difficult to interpret biologically. The aim of this study was to unravel the genetic background of BMD at pathway level, by integrating BMD GWAS data with genome-wide expression quantitative trait loci (eQTLs) and methylation quantitative trait loci (meQTLs) data. Method. We employed the GWAS datasets of BMD from the Genetic Factors for Osteoporosis Consortium (GEFOS), analysing patients’ BMD. The areas studied included 32 735 femoral necks, 28 498 lumbar spines, and 8143 forearms. Genome-wide eQTLs (containing 923 021 eQTLs) and meQTLs (containing 683 152 unique methylation sites with local meQTLs) data sets were collected from recently published studies. Gene scores were first calculated by summary data-based Mendelian randomisation (SMR) software and meQTL-aligned GWAS results. Gene set enrichment analysis (GSEA) was then applied to identify BMD-associated gene sets with a predefined significance level of 0.05. Results. We identified multiple gene sets associated with BMD in one or more regions, including relevant known biological gene sets such as the Reactome Circadian Clock (GSEA p-value = 1.0 × 10. -4. for LS and 2.7 × 10. -2. for femoral necks BMD in eQTLs-based GSEA) and insulin-like growth factor receptor binding (GSEA p-value = 5.0 × 10. -4. for femoral necks and 2.6 × 10. -2. for lumbar spines BMD in meQTLs-based GSEA). Conclusion. Our results provided novel clues for subsequent functional analysis of bone metabolism, and illustrated the benefit of integrating eQTLs and meQTLs data into pathway association analysis for genetic studies of complex human diseases. Cite this article: W. Wang, S. Huang, W. Hou, Y. Liu, Q. Fan, A. He, Y. Wen, J. Hao, X. Guo, F. Zhang. Integrative analysis of GWAS, eQTLs and meQTLs data suggests that multiple gene sets are associated with bone mineral density. Bone Joint Res 2017;6:572–576

Objectives. Previously, we reported the improved transfection efficiency of a plasmid DNA-chitosan (pDNA-CS) complex using a phosphorylatable nuclear localization signal-linked nucleic kinase substrate short peptide (pNNS) conjugated to chitosan (pNNS-CS). This study investigated the effects of pNNS-CS-mediated miR-140 and interleukin-1 receptor antagonist protein (IL-1Ra) gene transfection both in rabbit chondrocytes and a cartilage defect model. Methods. The pBudCE4.1-miR-140, pBudCE4.1-IL-1Ra, and negative control pBudCE4.1 plasmids were constructed and combined with pNNS-CS to form pDNA/pNNS-CS complexes. These complexes were transfected into chondrocytes or injected into the knee joint cavity. Results. High IL-1Ra and miR-140 expression levels were detected both in vitro and in vivo. In vitro, compared with the pBudCE4.1 group, the transgenic group presented with significantly increased chondrocyte proliferation and glycosaminoglycan (GAG) synthesis, as well as increased collagen type II alpha 1 chain (COL2A1), aggrecan (ACAN), and TIMP metallopeptidase inhibitor 1 (TIMP-1) levels. Nitric oxide (NO) synthesis was reduced, as were a disintegrin and metalloproteinase with thrombospondin type 1 motif 5 (ADAMTS-5) and matrix metalloproteinase (MMP)-13 levels. In vivo, the exogenous genes reduced the synovial fluid GAG and NO concentrations and the ADAMTS-5 and MMP-13 levels in cartilage. In contrast, COL2A1, ACAN, and TIMP-1 levels were increased, and the cartilage Mankin score was decreased in the transgenic group compared with the pBudCE4.1 group. Double gene combination produced greater efficacies than each single gene, both in vitro and in vivo. Conclusion. This study suggests that pNNS-CS is a good candidate for treating cartilage defects via gene therapy, and that IL-1Ra in combination with miR-140 produces promising biological effects on cartilage defects. Cite this article: R. Zhao, S. Wang, L. Jia, Q. Li, J. Qiao, X. Peng. Interleukin-1 receptor antagonist protein (IL-1Ra) and miR-140 overexpression via pNNS-conjugated chitosan-mediated gene transfer enhances the repair of full-thickness cartilage defects in a rabbit model. Bone Joint Res 2019;8:165–178. DOI: 10.1302/2046-3758.83.BJR-2018-0222.R1

Objectives. The objective of this study was to investigate the association of four single-nucleotide polymorphisms (SNPs) of the cannabinoid receptor 2 (CNR2) gene, gene-obesity interaction, and haplotype combination with osteoporosis (OP) susceptibility. Methods. Chinese patients with OP were recruited between March 2011 and December 2015 from our hospital. In this study, a total of 1267 post-menopausal female patients (631 OP patients and 636 control patients) were selected. The mean age of all subjects was 69.2 years (sd 15.8). A generalized multifactor dimensionality reduction (GMDR) model and logistic regression model were used to examine the interaction between SNP and obesity on OP. For OP patient-control haplotype analyses, the SHEsis online haplotype analysis software (. http://analysis.bio-x.cn/. ) was employed. Results. The logistic regression model revealed that the C allele of rs2501431 and the G allele of rs3003336 were associated with increased OP risk, compared with those with wild genotype. However, no significant correlations were found when analyzing the association of rs4237 and rs2229579 with OP risk. The GMDR analysis suggested that the interaction model composed of two factors, rs3003336 and abdominal obesity (AO), was the best model with statistical significance (p-value from sign test (P. sign. ) = 0.012), indicating a potential gene-environment interaction between rs3003336 and AO. Overall, the two-locus models had a cross-validation consistency of 10/10 and had a testing accuracy of 0.641. Abdominally obese subjects with the AG or GG genotype have the highest OP risk, compared with subjects with the AA genotype and normal waist circumference (WC) (odds ratio (OR) 2.23, 95% confidence interval (CI) 1.54 to 3.51). Haplotype analysis also indicated that the haplotype containing the rs3003336-G and rs2501431-C alleles was associated with a statistically increased OP risk. Conclusion. Our results suggested that the C allele of rs2501431 and the G allele of rs3003336 of the CNR2 gene, interaction between rs3003336 and AO, and the haplotype containing the rs3003336-G and rs2501431-C alleles were all associated with increased OP risk. Cite this article: Bone Joint Res 2019;8:544–549

Objectives. Previous genome-wide association studies (GWAS) have reported significant association of the single nucleotide polymorphism (SNP) rs8044769 in the fat mass and obesity-associated gene (FTO) with osteoarthritis (OA) risk in European populations. However, these findings have not been confirmed in Chinese populations. Methods. We systematically genotyped rs8044769 and evaluated the association between the genetic variants and OA risk in a case-controlled study including 196 OA cases and 442 controls in a northern Chinese population. Genotyping was performed using the Sequenom MassARRAY iPLEX platform. Results. We found that the variant T allele of rs8044769 showed no significant association of OA risk (p = 0.791), or association with body mass index (BMI) (pmeta = 0.786) in an additive genetic model. However, we detected a significant interaction between rs8044769 genotypes and BMI on OA risk (p = 0.037), as well as a borderline interaction between rs8044769 genotypes and age on OA risk (p = 0.062). Conclusions. Our findings indicate that rs8044769 in the FTO gene may not modify individual susceptibility to OA or increased BMI in the Chinese population. Further studies are warranted to validate and extend our findings. Cite this article: Prof L. Jiang. No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: A population-based study in China. Bone Joint Res 2016;5:169–174. DOI: 10.1302/2046-3758.55.2000589

Objectives. Effects of insulin-like growth factor 1 (IGF1), fibroblast growth factor 2 (FGF2) and bone morphogenetic protein 2 (BMP2) on the expression of genes involved in the proliferation and differentiation of osteoblasts in culture were analysed. The best sequence of growth factor addition that induces expansion of cells before their differentiation was sought. Methods. Primary human osteoblasts in in vitro culture were treated with IGF1, BMP2 or FGF2 (10 ng/ml) for 24 hours (IGF1) or 48 hours (BMP2 and FGF2). Experiments were performed during the exponential growth phase with approximately 1e7 cells per 75 cm. 2. flask. mRNA was reverse transcribed directly and analysed using RT-PCR Taqman assays. Expression levels of key genes involved in cell growth and differentiation (CDH11, TNFRSF11B, RUNX2, POSTN, ALP, WNT5A, LEF1, HSPA5, FOS, p21) were monitored using RT-PCR with gene-specific Taqman probes. . Results. Autocrine expression of BMP2 is stimulated by FGF2 and BMP2 itself. BMP2 and FGF2 act as proliferative factors as indicated by reduced expression of ALP and POSTN, whereas IGF1 exhibits a more subtle picture: the Wingless und Int-1 (Wnt) signalling pathway and the Smad pathway, but not p38 mitogen-activated protein (MAP) kinase signalling, were shown to be activated by IGF1, leading to proliferation and differentiation of the cells. . Conclusions. For future use of autologous bone cells in the management of bony defects, new treatment options take advantage of growth factors and differentiation factors. Thus, our results might help to guide the timely application of these factors for the expansion and subsequent differentiation of osteoblastic cells in culture. Cite this article: Bone Joint Res 2014;3:236–40

Objectives. This study aimed to investigate time-dependent gene expression of injured human anterior cruciate ligament (ACL), and to evaluate the histological changes of the ACL remnant in terms of cellular characterisation. Methods. Injured human ACL tissues were harvested from 105 patients undergoing primary ACL reconstruction and divided into four phases based on the period from injury to surgery. Phase I was < three weeks, phase II was three to eight weeks, phase III was eight to 20 weeks, and phase IV was ≥ 21 weeks. Gene expressions of these tissues were analysed in each phase by quantitative real-time polymerase chain reaction using selected markers (collagen types 1 and 3, biglycan, decorin, α-smooth muscle actin, IL-6, TGF-β1, MMP-1, MMP-2 and TIMP-1). Immunohistochemical staining was also performed using primary antibodies against CD68, CD55, Stat3 and phosphorylated-Stat3 (P-Stat3). . Results. Expression of IL-6 was mainly seen in phases I, II and III, collagen type 1 in phase II, MMP-1, 2 in phase III, and decorin, TGF-β1 and α-smooth muscle actin in phase IV. Histologically, degradation and scar formation were seen in the ACL remnant after phase III. The numbers of CD55 and P-Stat3 positive cells were elevated from phase II to phase III. . Conclusions. Elevated cell numbers including P-Stat3 positive cells were not related to collagens but to MMPs’ expressions

Objectives. Degenerative disc disease (DDD) and osteoarthritis (OA) are relatively frequent causes of disability amongst the elderly; they constitute serious socioeconomic costs and significantly impair quality of life. Previous studies to date have found that aggrecan variable number of tandem repeats (VNTR) contributes both to DDD and OA. However, current data are not consistent across studies. The purpose of this study was to evaluate systematically the relationship between aggrecan VNTR, and DDD and/or OA. Methods. This study used a highly sensitive search strategy to identify all published studies related to the relationship between aggrecan VNTR and both DDD and OA in multiple databases from January 1996 to December 2016. All identified studies were systematically evaluated using specific inclusion and exclusion criteria. Cochrane methodology was also applied to the results of this study. Results. The final selection of seven studies was comprehensively evaluated and includes results for 2928 alleles. The most frequent allele among all the studies was allele 27. After comparing the distributions of each allele with others, statistically significant differences have been found in the distribution of the alleles by the two groups, with an over-representation of allele (A)21 (disease: 3.22%, control: 0.44%). Thus, carrying A21 increased the risk of DDD. Such an association was not found to be statistically significant when considering the risk of OA. Conclusions. The findings suggest that VNTR A21 seems to be associated with higher risk to DDD, however, such an association may not be statistically significant regarding the risk of OA. Cite this article: L. Cong, G. Tu, D. Liang. A systematic review of the relationship between the distributions of aggrecan gene VNTR polymorphism and degenerative disc disease/osteoarthritis. Bone Joint Res 2018;7:308–317. DOI: 10.1302/2046-3758.74.BJR-2017-0207.R1

Aims. This study examined the relationship between obesity (OB) and osteoporosis (OP), aiming to identify shared genetic markers and molecular mechanisms to facilitate the development of therapies that target both conditions simultaneously. Methods. Using weighted gene co-expression network analysis (WGCNA), we analyzed datasets from the Gene Expression Omnibus (GEO) database to identify co-expressed gene modules in OB and OP. These modules underwent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and protein-protein interaction analysis to discover Hub genes. Machine learning refined the gene selection, with further validation using additional datasets. Single-cell analysis emphasized specific cell subpopulations, and enzyme-linked immunosorbent assay (ELISA), protein blotting, and cellular staining were used to investigate key genes. Results. WGCNA revealed critical gene modules for OB and OP, identifying the Toll-like receptor (TLR) signalling pathway as a common factor. TLR2 was the most significant gene, with a pronounced expression in macrophages. Elevated TLR2 expression correlated with increased adipose accumulation, inflammation, and osteoclast differentiation, linking it to OP development. Conclusion. Our study underscores the pivotal role of TLR2 in connecting OP and OB. It highlights the influence of TLR2 in macrophages, driving both diseases through a pro-inflammatory mechanism. These insights propose TLR2 as a potential dual therapeutic target for treating OP and OB. Cite this article: Bone Joint Res 2024;13(10):573–587

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. To explore the novel molecular mechanisms of histone deacetylase 4 (HDAC4) in chondrocytes via RNA sequencing (RNA-seq) analysis. Methods. Empty adenovirus (EP) and a HDAC4 overexpression adenovirus were transfected into cultured human chondrocytes. The cell survival rate was examined by real-time cell analysis (RTCA) and EdU and flow cytometry assays. Cell biofunction was detected by Western blotting. The expression profiles of messenger RNAs (mRNAs) in the EP and HDAC4 transfection groups were assessed using whole-transcriptome sequencing (RNA-seq). Volcano plot, Gene Ontology, and pathway analyses were performed to identify differentially expressed genes (DEGs). For verification of the results, the A289E/S246/467/632 A sites of HDAC4 were mutated to enhance the function of HDAC4 by increasing HDAC4 expression in the nucleus. RNA-seq was performed to identify the molecular mechanism of HDAC4 in chondrocytes. Finally, the top ten DEGs associated with ribosomes were verified by quantitative polymerase chain reaction (QPCR) in chondrocytes, and the top gene was verified both in vitro and in vivo. Results. HDAC4 markedly improved the survival rate and biofunction of chondrocytes. RNA-seq analysis of the EP and HDAC4 groups showed that HDAC4 induced 2,668 significant gene expression changes in chondrocytes (1,483 genes upregulated and 1,185 genes downregulated, p < 0.05), and ribosomes exhibited especially large increases. The results were confirmed by RNA-seq of the EP versus mutated HDAC4 groups and the validations in vitro and in vivo. Conclusion. The enhanced ribosome pathway plays a key role in the mechanism by which HDAC4 improves the survival rate and biofunction of chondrocytes. Cite this article: Bone Joint Res 2023;12(7):433–446

Aims. The antidiabetic agent metformin inhibits fibrosis in various organs. This study aims to elucidate the effects of hyperglycaemia and metformin on knee joint capsule fibrosis in mice. Methods. Eight-week-old wild-type (WT) and type 2 diabetic (db/db) mice were divided into four groups without or with metformin treatment (WT met(-/+), Db met(-/+)). Mice received daily intraperitoneal administration of metformin and were killed at 12 and 14 weeks of age. Fibrosis morphology and its related genes and proteins were evaluated. Fibroblasts were extracted from the capsules of 14-week-old mice, and the expression of fibrosis-related genes in response to glucose and metformin was evaluated in vitro. Results. The expression of all fibrosis-related genes was higher in Db met(-) than in WT met(-) and was suppressed by metformin. Increased levels of fibrosis-related genes, posterior capsule thickness, and collagen density were observed in the capsules of db/db mice compared with those in WT mice; these effects were suppressed by metformin. Glucose addition increased fibrosis-related gene expression in both groups of mice in vitro. When glucose was added, metformin inhibited the expression of fibrosis-related genes other than cellular communication network factor 2 (Ccn2) in WT mouse cells. Conclusion. Hyperglycaemia promotes fibrosis in the mouse knee joint capsule, which is inhibited by metformin. These findings can help inform the development of novel strategies for treating knee joint capsule fibrosis. Cite this article: Bone Joint Res 2024;13(7):321–331

Aims. Currently, the effect of drug treatment for osteoporosis is relatively poor, and the side effects are numerous and serious. Melatonin is a potential drug to improve bone mass in postmenopausal women. Unfortunately, the mechanism by which melatonin improves bone metabolism remains unclear. The aim of this study was to further investigate the potential mechanism of melatonin in the treatment of osteoporosis. Methods. The effects of melatonin on mitochondrial apoptosis protein, bmal1 gene, and related pathway proteins of RAW264.7 (mouse mononuclear macrophage leukaemia cells) were analyzed by western blot. Cell Counting Kit-8 was used to evaluate the effect of melatonin on cell viability. Flow cytometry was used to evaluate the effect of melatonin on the apoptosis of RAW264.7 cells and mitochondrial membrane potential. A reactive oxygen species (ROS) detection kit was used to evaluate the level of ROS in osteoclast precursors. We used bmal1-small interfering RNAs (siRNAs) to downregulate the Bmal1 gene. We established a postmenopausal mouse model and verified the effect of melatonin on the bone mass of postmenopausal osteoporosis in mice via micro-CT. Bmal1 lentiviral activation particles were used to establish an in vitro model of overexpression of the bmal1 gene. Results. Melatonin promoted apoptosis of RAW264.7 cells and increased the expression of BMAL1 to inhibit the activation of ROS and phosphorylation of mitogen-activated protein kinase (MAPK)-p38. Silencing the bmal1 gene weakened the above effects of melatonin. After that, we used dehydrocorydaline (DHC) to enhance the activation of MAPK-p38, and the effects of melatonin on reducing ROS levels and promoting apoptosis of RAW264.7 cells were also blocked. Then, we constructed a mouse model of postmenopausal osteoporosis and administered melatonin. The results showed that melatonin improves bone loss in ovariectomized mice. Finally, we established a model of overexpression of the bmal1 gene, and these results suggest that the bmal1 gene can regulate ROS activity and change the level of the MAPK-p38 signalling pathway. Conclusion. Our study confirmed that melatonin promotes the apoptosis of RAW264.7 cells through BMAL1/ROS/MAPK-p38, and revealed the therapeutic effect and mechanism of melatonin in postmenopausal osteoporosis. This finding enriches BMAL1 as a potential target for the treatment of osteoporosis and the pathogenesis of postmenopausal osteoporosis. Cite this article: Bone Joint Res 2023;12(11):677–690

Aims. Osteoarthritis (OA) is a common degenerative joint disease worldwide, which is characterized by articular cartilage lesions. With more understanding of the disease, OA is considered to be a disorder of the whole joint. However, molecular communication within and between tissues during the disease process is still unclear. In this study, we used transcriptome data to reveal crosstalk between different tissues in OA. Methods. We used four groups of transcription profiles acquired from the Gene Expression Omnibus database, including articular cartilage, meniscus, synovium, and subchondral bone, to screen differentially expressed genes during OA. Potential crosstalk between tissues was depicted by ligand-receptor pairs. Results. During OA, there were 626, 97, 1,060, and 2,330 differentially expressed genes in articular cartilage, meniscus, synovium, and subchondral bone, respectively. Gene Ontology enrichment revealed that these genes were enriched in extracellular matrix and structure organization, ossification, neutrophil degranulation, and activation at different degrees. Through ligand-receptor pairing and proteome of OA synovial fluid, we predicted ligand-receptor interactions and constructed a crosstalk atlas of the whole joint. Several interactions were reproduced by transwell experiment in chondrocytes and synovial cells, including TNC-NT5E, TNC-SDC4, FN1-ITGA5, and FN1-NT5E. After lipopolysaccharide (LPS) or interleukin (IL)-1β stimulation, the ligand expression of chondrocytes and synovial cells was upregulated, and corresponding receptors of co-culture cells were also upregulated. Conclusion. Each tissue displayed a different expression pattern in transcriptome, demonstrating their specific roles in OA. We highlighted tissue molecular crosstalk through ligand-receptor pairs in OA pathophysiology, and generated a crosstalk atlas. Strategies to interfere with these candidate ligands and receptors may help to discover molecular targets for future OA therapy. Cite this article: Bone Joint Res 2022;11(12):862–872

Aims. To assess the alterations in cell-specific DNA methylation associated with chondroitin sulphate response using peripheral blood collected from Kashin-Beck disease (KBD) patients before initiation of chondroitin sulphate treatment. Methods. Peripheral blood samples were collected from KBD patients at baseline of chondroitin sulphate treatment. Methylation profiles were generated using reduced representation bisulphite sequencing (RRBS) from peripheral blood. Differentially methylated regions (DMRs) were identified using MethylKit, while DMR-related genes were defined as those annotated to the gene body or 2.2-kilobase upstream regions of DMRs. Selected DMR-related genes were further validated by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) to assess expression levels. Tensor composition analysis was performed to identify cell-specific differential DNA methylation from bulk tissue. Results. This study revealed 21,060 hypermethylated and 44,472 hypomethylated DMRs, and 13,194 hypermethylated and 22,448 hypomethylated CpG islands for differential global methylation for chondroitin sulphate treatment response. A total of 12,666 DMR-related genes containing DMRs were identified in their promoter regions, such as CHL1 (false discovery rate (FDR) = 2.11 × 10. -11. ), RIC8A (FDR = 7.05 × 10. -4. ), and SOX12 (FDR = 1.43 × 10. -3. ). Additionally, RIC8A and CHL1 were hypermethylated in responders, while SOX12 was hypomethylated in responders, all showing decreased gene expression. The patterns of cell-specific differential global methylation associated with chondroitin sulphate response were observed. Specifically, we found that DMRs located in TESPA1 and ATP11A exhibited differential DNA methylation between responders and non-responders in granulocytes, monocytes, and B cells. Conclusion. Our study identified cell-specific changes in DNA methylation associated with chondroitin sulphate response in KBD patients. Cite this article: Bone Joint Res 2024;13(5):237–246

Aims. This study aimed, through bioinformatics analysis, to identify the potential diagnostic markers of osteoarthritis, and analyze the role of immune infiltration in synovial tissue. Methods. The gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified by R software. Functional enrichment analyses were performed and protein-protein interaction networks (PPI) were constructed. Then the hub genes were screened. Biomarkers with high value for the diagnosis of early osteoarthritis (OA) were validated by GEO datasets. Finally, the CIBERSORT algorithm was used to evaluate the immune infiltration between early-stage OA and end-stage OA, and the correlation between the diagnostic marker and infiltrating immune cells was analyzed. Results. A total of 88 DEGs were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DEGs were significantly enriched in leucocyte migration and interleukin (IL)-17 signalling pathways. Disease ontology (DO) indicated that DEGs were mostly enriched in rheumatoid arthritis. Six hub genes including FosB proto-oncogene, AP-1 transcription factor subunit (FOSB); C-X-C motif chemokine ligand 2 (CXCL2); CXCL8; IL-6; Jun proto-oncogene, AP-1 transcription factor subunit (JUN); and Activating transcription factor 3 (ATF3) were identified and verified by GEO datasets. ATF3 (area under the curve = 0.975) turned out to be a potential biomarker for the diagnosis of early OA. Several infiltrating immune cells varied significantly between early-stage OA and end-stage OA, such as resting NK cells (p = 0.016), resting dendritic cells (p = 0.043), and plasma cells (p = 0.043). Additionally, ATF3 was significantly correlated with resting NK cells (p = 0.034), resting dendritic cells (p = 0.026), and regulatory T cells (Tregs, p = 0.018). Conclusion. ATF3 may be a potential diagnostic marker for early diagnosis and treatment of OA, and immune cell infiltration provides new perspectives for understanding the mechanism during OA progression. Cite this article: Bone Joint Res 2022;11(9):679–689

Aims. This study aimed, through bioinformatics analysis and in vitro experiment validation, to identify the key extracellular proteins of intervertebral disc degeneration (IDD). Methods. The gene expression profile of GSE23130 was downloaded from the Gene Expression Omnibus (GEO) database. Extracellular protein-differentially expressed genes (EP-DEGs) were screened by protein annotation databases, and we used Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) to analyze the functions and pathways of EP-DEGs. STRING and Cytoscape were used to construct protein-protein interaction (PPI) networks and identify hub EP-DEGs. NetworkAnalyst was used to analyze transcription factors (TFs) and microRNAs (miRNAs) that regulate hub EP-DEGs. A search of the Drug Signatures Database (DSigDB) for hub EP-DEGs revealed multiple drug molecules and drug-target interactions. Results. A total of 56 EP-DEGs were identified in the differential expression analysis. EP-DEGs were enriched in the extracellular structure organization, ageing, collagen-activated signalling pathway, PI3K-Akt signalling pathway, and AGE-RAGE signalling pathway. PPI network analysis showed that the top ten hub EP-DEGs are closely related to IDD. Correlation analysis also demonstrated a significant correlation between the ten hub EP-DEGs (p＜0.05), which were selected to construct TF–gene interaction and TF–miRNA coregulatory networks. In addition, ten candidate drugs were screened for the treatment of IDD. Conclusion. The findings clarify the roles of extracellular proteins in IDD and highlight their potential as promising novel therapeutic targets. Cite this article: Bone Joint Res 2023;12(9):522–535

Aims. The effects of remnant preservation on the anterior cruciate ligament (ACL) and its relationship with the tendon graft remain unclear. We hypothesized that the co-culture of remnant cells and bone marrow stromal cells (BMSCs) decreases apoptosis and enhances the activity of the hamstring tendons and tenocytes, thus aiding ACL reconstruction. Methods. The ACL remnant, bone marrow, and hamstring tendons were surgically harvested from rabbits. The apoptosis rate, cell proliferation, and expression of types I and III collagen, transforming growth factor-β (TGF-β), vascular endothelial growth factor (VEGF), and tenogenic genes (scleraxis (SCX), tenascin C (TNC), and tenomodulin (TNMD)) of the hamstring tendons were compared between the co-culture medium (ACL remnant cells (ACLRCs) and BMSCs co-culture) and control medium (BMSCs-only culture). We also evaluated the apoptosis, cell proliferation, migration, and gene expression of hamstring tenocytes with exposure to co-culture and control media. Results. Compared to BMSCs-only culture medium, the co-culture medium showed substantially decreased early and late apoptosis rates, attenuation of intrinsic and extrinsic apoptotic pathways, and enhanced proliferation of the hamstring tendons and tenocytes. In addition, the expression of collagen synthesis, TGF-β, VEGF, and tenogenic genes in the hamstring tendons and tenocytes significantly increased in the co-culture medium compared to that in the control medium. Conclusion. In the presence of ACLRCs and BMSCs, the hamstring tendons and tenocytes significantly attenuated apoptosis and enhanced the expression of collagen synthesis, TGF-β, VEGF, and tenogenic genes. This in vitro study suggests that the ACLRCs mixed with BMSCs could aid regeneration of the hamstring tendon graft during ACL reconstruction. Cite this article: Bone Joint Res 2023;12(1):9–21

Aims. The decrease in the number of satellite cells (SCs), contributing to myofibre formation and reconstitution, and their proliferative capacity, leads to muscle loss, a condition known as sarcopenia. Resistance training can prevent muscle loss; however, the underlying mechanisms of resistance training effects on SCs are not well understood. We therefore conducted a comprehensive transcriptome analysis of SCs in a mouse model. Methods. We compared the differentially expressed genes of SCs in young mice (eight weeks old), middle-aged (48-week-old) mice with resistance training intervention (MID+ T), and mice without exercise (MID) using next-generation sequencing and bioinformatics. Results. After the bioinformatic analysis, the PI3K-Akt signalling pathway and the regulation of actin cytoskeleton in particular were highlighted among the top ten pathways with the most differentially expressed genes involved in the young/MID and MID+ T/MID groups. The expression of Gng5, Atf2, and Rtor in the PI3K-Akt signalling pathway was higher in the young and MID+ T groups compared with the MID group. Similarly, Limk1, Arhgef12, and Araf in the regulation of the actin cytoskeleton pathway had a similar bias. Moreover, the protein expression profiles of Atf2, Rptor, and Ccnd3 in each group were paralleled with the results of NGS. Conclusion. Our results revealed that age-induced muscle loss might result from age-influenced genes that contribute to muscle development in SCs. After resistance training, age-impaired genes were reactivated, and age-induced genes were depressed. The change fold in these genes in the young/MID mice resembled those in the MID + T/MID group, suggesting that resistance training can rejuvenate the self-renewing ability of SCs by recovering age-influenced genes to prevent sarcopenia. Cite this article: Bone Joint Res 2022;11(2):121–133

Aims. The present study investigated receptor activator of nuclear factor kappa-Β ligand (RANKL), osteoprotegerin (OPG), and Runt-related transcription factor 2 (RUNX2) gene expressions in giant cell tumour of bone (GCTB) patients in relationship with tumour recurrence. We also aimed to investigate the influence of CpG methylation on the transcriptional levels of RANKL and OPG. Methods. A total of 32 GCTB tissue samples were analyzed, and the expression of RANKL, OPG, and RUNX2 was evaluated by quantitative polymerase chain reaction (qPCR). The methylation status of RANKL and OPG was also evaluated by quantitative methylation-specific polymerase chain reaction (qMSP). Results. We found that RANKL and RUNX2 gene expression was upregulated more in recurrent than in non-recurrent GCTB tissues, while OPG gene expression was downregulated more in recurrent than in non-recurrent GCTB tissues. Additionally, we proved that changes in DNA methylation contribute to upregulating the expression of RANKL and downregulating the expression of OPG, which are critical for bone homeostasis and GCTB development. Conclusion. Our results suggest that the overexpression of RANKL/RUNX2 and the lower expression of OPG are associated with recurrence in GCTB patients. Cite this article: Bone Joint Res 2024;13(2):84–91

Aims. Dupuytren’s contracture is characterized by increased fibrosis of the palmar aponeurosis, with eventual replacement of the surrounding fatty tissue with palmar fascial fibromatosis. We hypothesized that adipocytokines produced by adipose tissue in contact with the palmar aponeurosis might promote fibrosis of the palmar aponeurosis. Methods. We compared the expression of the adipocytokines adiponectin and leptin in the adipose tissue surrounding the palmar aponeurosis of male patients with Dupuytren’s contracture, and of male patients with carpal tunnel syndrome (CTS) as the control group. We also examined the effects of adiponectin on fibrosis-related genes and proteins expressed by fibroblasts in the palmar aponeurosis of patients with Dupuytren’s contracture. Results. Adiponectin expression in the adipose tissue surrounding the palmar aponeurosis was significantly lower in patients with Dupuytren’s contracture than in those with CTS. The expression of fibrosis-related genes and proteins, such as types 1 and 3 collagen and α-smooth muscle actin, was suppressed in a concentration-dependent manner by adding AdipoRon, an adiponectin receptor agonist. The expression of fibrosis-related genes and proteins was also suppressed by AdipoRon in the in vitro model of Dupuytren’s contracture created by adding TGF-β to normal fibroblasts collected from patients with CTS. Conclusion. Fibrosis of the palmar aponeurosis in Dupuytren’s contracture in males may be associated with adiponectin expression in the adipose tissue surrounding the palmar aponeurosis. Although fibroblasts within the palmar aponeurosis are often the focus of attention when elucidating the pathogenesis of Dupuytren’s contracture, adiponectin expression in adipose tissues warrants closer attention in future research. Cite this article: Bone Joint Res 2023;12(8):486–493

Aims. This study aimed to demonstrate the promoting effect of elastic fixation on fracture, and further explore its mechanism at the gene and protein expression levels. Methods. A closed tibial fracture model was established using 12 male Japanese white rabbits, and divided into elastic and stiff fixation groups based on different fixation methods. Two weeks after the operation, a radiograph and pathological examination of callus tissue were used to evaluate fracture healing. Then, the differentially expressed proteins (DEPs) were examined in the callus using proteomics. Finally, in vitro cell experiments were conducted to investigate hub proteins involved in this process. Results. Mean callus volume was larger in the elastic fixation group (1,755 mm. 3. (standard error of the mean (SEM) 297)) than in the stiff fixation group (258 mm. 3. (SEM 65)). Pathological observation found that the expression levels of osterix (OSX), collagen, type I, alpha 1 (COL1α1), and alkaline phosphatase (ALP) in the callus of the elastic fixation group were higher than those of the stiff fixation group. The protein sequence of the callus revealed 199 DEPs, 124 of which were highly expressed in the elastic fixation group. In the in vitro study, it was observed that a stress of 200 g led to upregulation of thrombospondin 1 (THBS1) and osteoglycin (OGN) expression in bone marrow mesenchymal stem cells (BMSCs). Additionally, these genes were found to be upregulated during the osteogenic differentiation process of the BMSCs. Conclusion. Elastic fixation can promote fracture healing and osteoblast differentiation in callus, and the ability of elastic fixation to promote osteogenic differentiation of BMSCs may be achieved by upregulating genes such as THBS1 and OGN. Cite this article: Bone Joint Res 2024;13(10):559–572

Aims. This study aimed to explore the role of small colony variants (SCVs) of Staphylococcus aureus in intraosseous invasion and colonization in patients with periprosthetic joint infection (PJI). Methods. A PJI diagnosis was made according to the MusculoSkeletal Infection Society (MSIS) for PJI. Bone and tissue samples were collected intraoperatively and the intracellular invasion and intraosseous colonization were detected. Transcriptomics of PJI samples were analyzed and verified by polymerase chain reaction (PCR). Results. SCVs can be isolated from samples collected from chronic PJIs intraoperatively. Transmission electron microscopy (TEM) and immunofluorescence (IF) showed that there was more S. aureus in bone samples collected from chronic PJIs, but much less in bone samples from acute PJIs, providing a potential mechanism of PJI. Immunofluorescence results showed that SCVs of S. aureus were more likely to invade osteoblasts in vitro. Furthermore, TEM and IF also demonstrated that SCVs of S. aureus were more likely to invade and colonize in vivo. Cluster analysis and principal component analysis (PCA) showed that there were substantial differences in gene expression profiles between chronic and acute PJI. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these differentially expressed genes were enriched to chemokine-related signal pathways. PCR also verified these results. Conclusion. Our study has shown that the S. aureus SCVs have a greater ability to invade and colonize in bone, resulting in S. aureus remaining in bone tissues long-term, thus explaining the pathogenesis of chronic PJI. Cite this article: Bone Joint Res 2022;11(12):843–853

Aims. To test the hypothesis that reseeded anterior cruciate ligament (ACL)-derived cells have a better ability to survive and integrate into tendon extracellular matrix (ECM) and accelerate the ligamentization process, compared to adipose-derived mesenchymal stem cells (ADMSCs). Methods. Acellularized tibialis allograft tendons were used. Tendons were randomly reseeded with ACL-derived cells or ADMSCs. ACL-derived cells were harvested and isolated from remnants of ruptured ACLs during reconstruction surgery and cultured at passage three. Cell suspensions (200 µl) containing 2 × 10. 6. ACL-derived cells or ADMSCs were prepared for the purpose of reseeding. At days 1, 3, and 7 post-reseeding, graft composites were assessed for repopulation with histological and immunohistochemical analysis. Matrix protein contents and gene expression levels were analyzed. Results. In the graft reseeded with ACL-derived cells, a large number of elongated cells that integrated into the matrix were evident at day 3 and day 7. However, in the graft reseeded with ADMSCs, only a small number of elongated cells were found integrated into the matrix. Immunofluorescence for Ki-67 and type I collagen confirmed the pronounced production of type I collagen by Ki-67-positive ACL-derived cells integrated into the ECM. A messenger RNA (mRNA) expression assay demonstrated significantly higher gene expression levels of types I (p = 0.013) and III (p = 0.050) collagen in the composites reseeded with ACL-derived cells than ADMSCs. Conclusion. ACL-derived cells, when reseeded to acellularized tendon graft, demonstrated earlier better survival and integration in the tendon ECM and resulted in higher gene expression levels of collagen, which may be essential to the normal ligamentization process compared to ADMSCs. Cite this article: Bone Joint Res 2022;11(11):777–786

Aims. The association of auraptene (AUR), a 7-geranyloxycoumarin, on osteoporosis and its potential pathway was predicted by network pharmacology and confirmed in experimental osteoporotic mice. Methods. The network of AUR was constructed and a potential pathway predicted by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) terms enrichment. Female ovariectomized (OVX) Institute of Cancer Research mice were intraperitoneally injected with 0.01, 0.1, and 1 mM AUR for four weeks. The bone mineral density (BMD) level was measured by dual-energy X-ray absorptiometry. The bone microstructure was determined by histomorphological changes in the femora. In addition, biochemical analysis of the serum and assessment of the messenger RNA (mRNA) levels of osteoclastic markers were performed. Results. In total, 65.93% of the genes of the AUR network matched with osteoporosis-related genes. Osteoclast differentiation was predicted to be a potential pathway of AUR in osteoporosis. Based on the network pharmacology, the BMD and bone mineral content levels were significantly (p < 0.05) increased in the whole body, femur, tibia, and lumbar spine by AUR. AUR normalized the bone microstructure and the serum alkaline phosphatase (ALP), bone-specific alkaline phosphatase (bALP), osteocalcin, and calcium in comparison with the OVX group. In addition, AUR treatment reduced TRAP-positive osteoclasts and receptor activator of nuclear factor kappa-B ligand (RANKL). +. nuclear factor of activated T cells 1 (NFATc1). +. expression in the femoral body. Moreover, the expressions of initiators for osteoclastic resorption and bone matrix degradation were significantly (p < 0.05) regulated by AUR in the lumbar spine of the osteoporotic mice. Conclusion. AUR ameliorated bone loss by downregulating the RANKL/NFATc1 pathway, resulting in improvement of osteoporosis. In conclusion, AUR might be an ameliorative cure that alleviates bone loss in osteoporosis via inhibition of osteoclastic activity. Cite this article: Bone Joint Res 2022;11(5):304–316

Aims. Developmental dysplasia of the hip (DDH) is a complex musculoskeletal disease that occurs mostly in children. This study aimed to investigate the molecular changes in the hip joint capsule of patients with DDH. Methods. High-throughput sequencing was used to identify genes that were differentially expressed in hip joint capsules between healthy controls and DDH patients. Biological assays including cell cycle, viability, apoptosis, immunofluorescence, reverse transcription polymerase chain reaction (RT-PCR), and western blotting were performed to determine the roles of the differentially expressed genes in DDH pathology. Results. More than 1,000 genes were differentially expressed in hip joint capsules between healthy controls and DDH. Both gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that extracellular matrix (ECM) modifications, muscle system processes, and cell proliferation were markedly influenced by the differentially expressed genes. Expression of Collagen Type I Alpha 1 Chain (COL1A1), COL3A1, matrix metalloproteinase-1 (MMP1), MMP3, MMP9, and MMP13 was downregulated in DDH, with the loss of collagen fibres in the joint capsule. Expression of transforming growth factor beta 1 (TGF-β1) was downregulated, while that of TGF-β2, Mothers against decapentaplegic homolog 3 (SMAD3), and WNT11 were upregulated in DDH, and alpha smooth muscle actin (αSMA), a key myofibroblast marker, showed marginal increase. In vitro studies showed that fibroblast proliferation was suppressed in DDH, which was associated with cell cycle arrest in G0/G1 and G2/M phases. Cell cycle regulators including Cyclin B1 (CCNB1), Cyclin E2 (CCNE2), Cyclin A2 (CCNA2), Cyclin-dependent kinase 1 (CDK1), E2F1, cell division cycle 6 (CDC6), and CDC7 were downregulated in DDH. Conclusion. DDH is associated with the loss of collagen fibres and fibroblasts, which may cause loose joint capsule formation. However, the degree of differentiation of fibroblasts to myofibroblasts needs further study. Cite this article: Bone Joint Res 2021;10(9):558–570

Aims. Transcription factor nuclear factor kappa B (NF-κB) plays a major role in the pathogenesis of chronic inflammatory diseases in all organ systems. Despite its importance, NF-κB targeted drug therapy to mitigate chronic inflammation has had limited success in preclinical studies. We hypothesized that sex differences affect the response to NF-κB treatment during chronic inflammation in bone. This study investigated the therapeutic effects of NF-κB decoy oligodeoxynucleotides (ODN) during chronic inflammation in male and female mice. Methods. We used a murine model of chronic inflammation induced by continuous intramedullary delivery of lipopolysaccharide-contaminated polyethylene particles (cPE) using an osmotic pump. Specimens were evaluated using micro-CT and histomorphometric analyses. Sex-specific osteogenic and osteoclastic differentiation potentials were also investigated in vitro, including alkaline phosphatase, Alizarin Red, tartrate-resistant acid phosphatase staining, and gene expression using reverse transcription polymerase chain reaction (RT-PCR). Results. Local delivery of NF-κB decoy ODN in vivo increased osteogenesis in males, but not females, in the presence of chronic inflammation induced by cPE. Bone resorption activity was decreased in both sexes. In vitro osteogenic and osteoclastic differentiation assays during inflammatory conditions did not reveal differences among the groups. Receptor activator of nuclear factor kappa Β ligand (Rankl) gene expression by osteoblasts was significantly decreased only in males when treated with ODN. Conclusion. We demonstrated that NF-κB decoy ODN increased osteogenesis in male mice and decreased bone resorption activity in both sexes in preclinical models of chronic inflammation. NF-κB signalling could be a therapeutic target for chronic inflammatory diseases involving bone, especially in males. Cite this article: Bone Joint Res 2024;13(1):28–39

Aims. The aim of this study was to explore the genetic correlation and causal relationship between blood plasma proteins and rheumatoid arthritis (RA). Methods. Based on the genome-wide association studies (GWAS) summary statistics of RA from European descent and the GWAS summary datasets of 3,622 plasma proteins, we explored the relationship between RA and plasma proteins from three aspects. First, linkage disequilibrium score regression (LD score regression) was applied to detect the genetic correlation between RA and plasma proteins. Mendelian randomization (MR) analysis was then used to evaluate the causal association between RA and plasma proteins. Finally, GEO2R was used to screen the differentially expressed genes (DEGs) between patients with RA and healthy controls. Results. We found that seven kinds of plasma proteins had genetic correlations with RA, such as Soluble Receptor for Advanced Glycation End Products (sRAGE) (correlation coefficient = 0.2582, p = 0.049), vesicle transport protein USE1 (correlation coefficient = 0.1337, p = 0.018), and spermatogenesis-associated protein 20 (correlation coefficient = 0.3706, p = 0.018). There was a significant causal relationship between sRAGE and RA. By comparing the genes encoding seven plasma proteins, we found that only USE1 was a DEG associated with RA. Conclusion. Our study identified a set of candidate plasma proteins that showed signals correlated with RA. Since the results of this study need further experimental verification, they should be interpreted with caution. However, we hope that this paper will provide new insights for the discovery of pathogenic genes and RA pathogenesis in the future. Cite this article: Bone Joint Res 2022;11(2):134–142

Aims. Glucose-insulin-potassium (GIK) is protective following cardiac myocyte ischaemia-reperfusion (IR) injury, however the role of GIK in protecting skeletal muscle from IR injury has not been evaluated. Given the similar mechanisms by which cardiac and skeletal muscle sustain an IR injury, we hypothesized that GIK would similarly protect skeletal muscle viability. Methods. A total of 20 C57BL/6 male mice (10 control, 10 GIK) sustained a hindlimb IR injury using a 2.5-hour rubber band tourniquet. Immediately prior to tourniquet placement, a subcutaneous osmotic pump was placed which infused control mice with saline (0.9% sodium chloride) and treated mice with GIK (40% glucose, 50 U/l insulin, 80 mEq/L KCl, pH 4.5) at a rate of 16 µl/hr for 26.5 hours. At 24 hours following tourniquet removal, bilateral (tourniqueted and non-tourniqueted) gastrocnemius muscles were triphenyltetrazolium chloride (TTC)-stained to quantify percentage muscle viability. Bilateral peroneal muscles were used for gene expression analysis, serum creatinine and creatine kinase activity were measured, and a validated murine ethogram was used to quantify pain before euthanasia. Results. GIK treatment resulted in a significant protection of skeletal muscle with increased viability (GIK 22.07% (SD 15.48%)) compared to saline control (control 3.14% (SD 3.29%)) (p = 0.005). Additionally, GIK led to a statistically significant reduction in gene expression markers of cell death (CASP3, p < 0.001) and inflammation (NOS2, p < 0.001; IGF1, p = 0.007; IL-1β, p = 0.002; TNFα, p = 0.012), and a significant reduction in serum creatine kinase (p = 0.004) and creatinine (p < 0.001). GIK led to a significant reduction in IR-related pain (p = 0.030). Conclusion. Systemic GIK infusion during and after limb ischaemia protects murine skeletal muscle from cell death, kidneys from reperfusion metabolites, and reduces pain by reducing post-ischaemic inflammation. Cite this article: Bone Joint Res 2023;12(3):212–218

Aims. To verify whether secretory leucocyte protease inhibitor (SLPI) can promote early tendon-to-bone healing after anterior cruciate ligament (ACL) reconstruction. Methods. In vitro: the mobility of the rat bone mesenchymal stem cells (BMSCs) treated with SLPI was evaluated by scratch assay. Then the expression levels of osteogenic differentiation-related genes were analyzed by real-time quantitative PCR (qPCR) to determine the osteogenic effect of SLPI on BMSCs. In vivo: a rat model of ACL reconstruction was used to verify the effect of SLPI on tendon-to-bone healing. All the animals of the SLPI group and the negative control (NC) group were euthanized for histological evaluation, micro-CT scanning, and biomechanical testing. Results. SLPI improved the migration ability of BMSCs and upregulated the expression of genes related to osteogenic differentiation of BMSCs in vitro. In vivo, the SLPI group had higher histological scores at the tendon-bone interface by histological evaluation. Micro-CT showed more new bone formation and bone ingrowth around the grafted tendon in the SLPI group. Evaluation of the healing strength of the tendon-bone connection showed that the SLPI group had a higher maximum failure force and stiffness. Conclusion. SLPI can effectively promote early tendon-to-bone healing after ACL reconstruction via enhancing the migration and osteogenic differentiation of BMSCs. Cite this article: Bone Joint Res 2022;11(7):503–512

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