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. 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. 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. 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. To investigate the effects of senescent osteocytes on bone homeostasis in the progress of age-related osteoporosis and explore the underlying mechanism. Methods. In a series of in vitro experiments, we used tert-Butyl hydroperoxide (TBHP) to induce senescence of MLO-Y4 cells successfully, and collected conditioned medium (CM) and senescent MLO-Y4 cell-derived exosomes, which were then applied to MC3T3-E1 cells, separately, to evaluate their effects on osteogenic differentiation. Furthermore, we identified differentially expressed microRNAs (miRNAs) between exosomes from senescent and normal MLO-Y4 cells by high-throughput RNA sequencing. Based on the key miRNAs that were discovered, the underlying mechanism by which senescent osteocytes regulate osteogenic differentiation was explored. Lastly, in the in vivo experiments, the effects of senescent MLO-Y4 cell-derived exosomes on age-related bone loss were evaluated in male SAMP6 mice, which excluded the effects of oestrogen, and the underlying mechanism was confirmed. Results. The CM and exosomes collected from senescent MLO-Y4 cells inhibited osteogenic differentiation of MC3T3-E1 cells. RNA sequencing detected significantly lower expression of miR-494-3p in senescent MLO-Y4 cell-derived exosomes compared with normal exosomes. The upregulation of exosomal miR-494-3p by miRNA mimics attenuated the effects of senescent MLO-Y4 cell-derived exosomes on osteogenic differentiation. Luciferase reporter assay demonstrated that miR-494-3p targeted phosphatase and tensin homolog (PTEN), which is a negative regulator of the phosphoinositide 3-kinase (PI3K)/AKT pathway. Overexpression of PTEN or inhibition of the PI3K/AKT pathway blocked the functions of exosomal miR-494-3p. In SAMP6 mice, senescent MLO-Y4 cell-derived exosomes accelerated bone loss, which was rescued by upregulation of exosomal miR-494-3p. Conclusion. Reduced expression of miR-494-3p in senescent osteocyte-derived exosomes inhibits osteogenic differentiation and accelerates age-related bone loss via PTEN/PI3K/AKT pathway. Cite this article: Bone Joint Res 2024;13(2):52–65

Aims. To perform an incremental cost-utility analysis and assess the impact of differential costs and case volume on the cost-effectiveness of robotic arm-assisted unicompartmental knee arthroplasty (rUKA) compared to manual (mUKA). Methods. This was a five-year follow-up study of patients who were randomized to rUKA (n = 64) or mUKA (n = 65). Patients completed the EuroQol five-dimension questionnaire (EQ-5D) preoperatively, and at three months and one, two, and five years postoperatively, which was used to calculate quality-adjusted life years (QALYs) gained. Costs for the primary and additional surgery and healthcare costs were calculated. Results. rUKA was associated with a relative 0.012 QALY gain at five years, which was associated with an incremental cost per QALY of £13,078 for a unit undertaking 400 cases per year. A cost per QALY of less than £20,000 was achieved when ≥ 300 cases were performed per year. However, on removal of the cost for a revision for presumed infection (mUKA group, n = 1) the cost per QALY was greater than £38,000, which was in part due to the increased intraoperative consumable costs associated with rUKA (£626 per patient). When the absolute cost difference (operative and revision costs) was less than £240, a cost per QALY of less than £20,000 was achieved. On removing the cost of the revision for infection, rUKA was cost-neutral when more than 900 cases per year were undertaken and when the consumable costs were zero. Conclusion. rUKA was a cost-effective intervention with an incremental cost per QALY of £13,078 at five years, however when removing the revision for presumed infection, which was arguably a random event, this was no longer the case. The absolute cost difference had to be less than £240 to be cost-effective, which could be achieved by reducing the perioperative costs of rUKA or if there were increased revision costs associated with mUKA with longer follow-up. Cite this article: Bone Jt Open 2023;4(11):889–898

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. 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. Sagittal plane imbalance (SPI), or asymmetry between extension and flexion gaps, is an important issue in total knee arthroplasty (TKA). The purpose of this study was to compare SPI between kinematic alignment (KA), mechanical alignment (MA), and functional alignment (FA) strategies. Methods. In 137 robotic-assisted TKAs, extension and flexion stressed gap laxities and bone resections were measured. The primary outcome was the proportion and magnitude of medial and lateral SPI (gap differential > 2.0 mm) for KA, MA, and FA. Secondary outcomes were the proportion of knees with severe (> 4.0 mm) SPI, and resection thicknesses for each technique, with KA as reference. Results. FA showed significantly lower rates of medial and lateral SPI (2.9% and 2.2%) compared to KA (45.3%; p < 0.001, and 25.5%; p < 0.001) and compared to MA (52.6%; p < 0.001 and 29.9%; p < 0.001). There was no difference in medial and lateral SPI between KA and MA (p = 0.228 and p = 0.417, respectively). FA showed significantly lower rates of severe medial and lateral SPI (0 and 0%) compared to KA (8.0%; p < 0.001 and 7.3%; p = 0.001) and compared to MA (10.2%; p < 0.001 and 4.4%; p = 0.013). There was no difference in severe medial and lateral SPI between KA and MA (p = 0.527 and p = 0.307, respectively). MA resulted in thinner resections than KA in medial extension (mean difference (MD) 1.4 mm, SD 1.9; p < 0.001), medial flexion (MD 1.5 mm, SD 1.8; p < 0.001), and lateral extension (MD 1.1 mm, SD 1.9; p < 0.001). FA resulted in thinner resections than KA in medial extension (MD 1.6 mm, SD 1.4; p < 0.001) and lateral extension (MD 2.0 mm, SD 1.6; p < 0.001), but in thicker medial flexion resections (MD 0.8 mm, SD 1.4; p < 0.001). Conclusion. Mechanical and kinematic alignment (measured resection techniques) result in high rates of SPI. Pre-resection angular and translational adjustments with functional alignment, with typically smaller distal than posterior femoral resection, address this issue. Cite this article: Bone Jt Open 2024;5(8):681–687

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

Objectives. The aim of this study was to provide a comprehensive understanding of alterations in messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs) in cartilage affected by osteoarthritis (OA). Methods. The expression profiles of mRNAs, lncRNAs, and circRNAs in OA cartilage were assessed using whole-transcriptome sequencing. Bioinformatics analyses included prediction and reannotation of novel lncRNAs and circRNAs, their classification, and their placement into subgroups. Gene ontology and pathway analysis were performed to identify differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed circRNAs (DECs). We focused on the overlap of DEGs and targets of DELs previously identified in seven high-throughput studies. The top ten DELs were verified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) in articular chondrocytes, both in vitro and in vivo. Results. In total, 739 mRNAs, 1152 lncRNAs, and 42 circRNAs were found to be differentially expressed in OA cartilage tissue. Among these, we identified 18 overlapping DEGs and targets of DELs, and the top ten DELs were screened by expression profile analysis as candidate OA-related genes. WISP2, ATF3, and CHI3L1 were significantly increased in both normal versus OA tissues and normal versus interleukin (IL)-1β-induced OA-like cell models, while ADAM12, PRELP, and ASPN were shown to be significantly decreased. Among the identified DELs, we observed higher expression of ENST00000453554 and MSTRG.99593.3, and lower expression of MSTRG.44186.2 and NONHSAT186094.1 in normal versus OA cells and tissues. Conclusion. This study revealed expression patterns of coding and noncoding RNAs in OA cartilage, which added sets of genes and noncoding RNAs to the list of candidate diagnostic biomarkers and therapeutic agents for OA patients. Cite this article: H. Li, H. H. Yang, Z. G. Sun, H. B. Tang, J. K. Min. Whole-transcriptome sequencing of knee joint cartilage from osteoarthritis patients. Bone Joint Res 2019;8:290–303. DOI: 10.1302/2046-3758.87.BJR-2018-0297.R1

Aims. The aim of this study was to compare the early postoperative mortality and morbidity in older patients with a fracture of the femoral neck, between those who underwent total hip arthroplasty (THA) and those who underwent hemiarthroplasty. Methods. This nationwide, retrospective cohort study used data from the Japanese Diagnosis Procedure Combination database. We included older patients (aged ≥ 60 years) who underwent THA or hemiarthroplasty after a femoral neck fracture, between July 2010 and March 2022. A total of 165,123 patients were included. The THA group was younger (mean age 72.6 (SD 8.0) vs 80.7 years (SD 8.1)) and had fewer comorbidities than the hemiarthroplasty group. Patients with dementia or malignancy were excluded because they seldom undergo THA. The primary outcome measures were mortality and complications while in hospital, and secondary outcomes were readmission and reoperation within one and two years after discharge, and the costs of hospitalization. We conducted an instrumental variable analysis (IVA) using differential distance as a variable. Results. The IVA analysis showed that the THA group had a significantly higher rate of complications while in hospital (risk difference 6.3% (95% CI 2.0 to 10.6); p = 0.004) than the hemiarthroplasty group, but there was no significant difference in the rate of mortality while in hospital (risk difference 0.3% (95% CI -1.7 to 2.2); p = 0.774). There was no significant difference in the rate of readmission (within one year: risk difference 1.3% (95% CI -1.9 to 4.5); p = 0.443; within two years: risk difference 0.1% (95% CI -3.2 to 3.4); p = 0.950) and reoperation (within one year: risk difference 0.3% (95% CI -0.6 to 1.1); p = 0.557; within two years: risk difference 0.1% (95% CI -0.4 to 0.7); p = 0.632) after discharge. The costs of hospitalization were significantly higher in the THA group than in the hemiarthroplasty group (difference $2,634 (95% CI $2,496 to $2,772); p < 0.001). Conclusion. Among older patients undergoing surgery for a femoral neck fracture, the risk of early complications was higher after THA than after hemiarthroplasty. Our findings should aid in clinical decision-making in these patients. Cite this article: Bone Joint J 2024;106-B(9):986–993

Aims. This study aimed to describe practice variation in the use of total hip arthroplasty (THA) for older patients with femoral neck fracture and to determine the association between patient, surgeon, and institution factors and treatment with THA. Methods. We performed a cross-sectional analysis of 49,597 patients aged 60 years and older from Ontario, Canada, who underwent hemiarthroplasty or THA for femoral neck fracture between 2002 and 2017. This population-based study used routinely collected healthcare databases linked through ICES (formerly known as the Institute for Clinical Evaluative Sciences). Multilevel logistic regression modelling was used to quantify the association between patient, surgeon, and institution-level variables and whether patients were treated with THA. Variance partition coefficient and median odds ratios were used to estimate the variation attributable to higher-level variables and the magnitude of effect of higher-level variables, respectively. Results. Over the study period, 9.4% of patients (n = 4,638) were treated with THA. Patient factors associated with higher likelihood of treatment by THA included: younger age, male sex, and diagnosis with rheumatoid arthritis. Long-term care residence, use of home care services prior to hip fracture, diagnosis of dementia, higher comorbidity burden, and the most marginalized group were negatively associated with treatment by THA. Treating surgeon and institution accounted for 54.2% and 17.8% of the total variation in treatment with THA, respectively. Surgeon volume of THA procedures in the 365 days prior to surgery was the strongest higher-level predictor of treatment with THA. Specific treating surgeons and institutions still accounted for significant proportions of the variability in treatment with THA (40.3% and 19.5% of total observed variation, respectively) after controlling for available patient, surgeon, and institution-level variables. Conclusion. The strongest predictors for treatment of patients with femoral neck fracture with THA were patient age, treating surgeon, and treating institution. This practice variation highlights differential access to care for patients. Cite this article: Bone Joint J 2023;105-B(2):180–189

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. Adenosine, lidocaine, and Mg. 2+. (ALM) therapy exerts differential immuno-inflammatory responses in males and females early after anterior cruciate ligament (ACL) reconstruction (ACLR). Our aim was to investigate sex-specific effects of ALM therapy on joint tissue repair and recovery 28 days after surgery. Methods. Male (n = 21) and female (n = 21) adult Sprague-Dawley rats were randomly divided into ALM or Saline control treatment groups. Three days after ACL rupture, animals underwent ACLR. An ALM or saline intravenous infusion was commenced prior to skin incision, and continued for one hour. An intra-articular bolus of ALM or saline was also administered prior to skin closure. Animals were monitored to 28 days, and joint function, pain, inflammatory markers, histopathology, and tissue repair markers were assessed. Results. Despite comparable knee function, ALM-treated males had reduced systemic inflammation, synovial fluid angiogenic and pro-inflammatory mediators, synovitis, and fat pad fibrotic changes, compared to controls. Within the ACL graft, ALM-treated males had increased expression of tissue repair markers, decreased inflammation, increased collagen organization, and improved graft-bone healing. In contrast to males, females had no evidence of persistent systemic inflammation. Compared to controls, ALM-treated females had improved knee extension, gait biomechanics, and elevated synovial macrophage inflammatory protein-1 alpha (MIP-1α). Within the ACL graft, ALM-treated females had decreased inflammation, increased collagen organization, and improved graft-bone healing. In articular cartilage of ALM-treated animals, matrix metalloproteinase (MMP)-13 expression was blunted in males, while in females repair markers were increased. Conclusion. At 28 days, ALM therapy reduces inflammation, augments tissue repair patterns, and improves joint function in a sex-specific manner. The study supports transition to human safety trials. Cite this article: Bone Joint Res 2024;13(6):279–293

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. Exosomes (exo) are involved in the progression of osteoarthritis (OA). This study aimed to investigate the function of dysfunctional chondrocyte-derived exo (DC-exo) on OA in rats and rat macrophages. Methods. Rat-derived chondrocytes were isolated, and DCs induced with interleukin (IL)-1β were used for exo isolation. Rats with OA (n = 36) or macrophages were treated with DC-exo or phosphate-buffered saline (PBS). Macrophage polarization and autophagy, and degradation and chondrocyte activity of cartilage tissues, were examined. RNA sequencing was used to detect genes differentially expressed in DC-exo, followed by RNA pull-down and ribonucleoprotein immunoprecipitation (RIP). Long non-coding RNA osteoarthritis non-coding transcript (OANCT) and phosphoinositide-3-kinase regulatory subunit 5 (PIK3R5) were depleted in DC-exo-treated macrophages and OA rats, in order to observe macrophage polarization and cartilage degradation. The PI3K/AKT/mammalian target of rapamycin (mTOR) pathway activity in cells and tissues was measured using western blot. Results. DC-exo inhibited macrophage autophagy (p = 0.002) and promoted M1 macrophage polarization (p = 0.002). DC-exo at 20 μg/ml induced collagen degradation (p < 0.001) and inflammatory cell infiltration (p = 0.023) in rats. OANCT was elevated in DC (p < 0.001) and in cartilage tissues of OA patients (p < 0.001), and positively correlated with patients’ Kellgren-Lawrence grade (p < 0.001). PIK3R5 was increased in DC-exo-treated cartilage tissues (p < 0.001), and OANCT bound to fat mass and obesity-associated protein (FTO) (p < 0.001). FTO bound to PIK3R5 (p < 0.001) to inhibit the stability of PIK3R5 messenger RNA (mRNA) (p < 0.001) and disrupt the PI3K/AKT/mTOR pathway (p < 0.001). Conclusion. Exosomal OANCT from DC could bind to FTO protein, thereby maintaining the mRNA stability of PIK3R5, further activating the PI3K/AKT/mTOR pathway to exacerbate OA. Cite this article: Bone Joint Res 2022;11(9):652–668

Aims. Bone turnover markers (BTMs) follow distinct trends after fractures and limited evidence suggests differential levels in BTMs in patients with delayed healing. The effect of vitamin D, and other factors that influence BTMs and fracture healing, is important to elucidate the use of BTMs as surrogates of fracture healing. We sought to determine whether BTMs can be used as early markers of delayed fracture healing, and the effect of vitamin D on BTM response after fracture. Methods. A total of 102 participants aged 18 to 50 years (median 28 years (interquartile range 23 to 35)), receiving an intramedullary nail for a tibial or femoral shaft fracture, were enrolled in a randomized controlled trial comparing vitamin D. 3. supplementation to placebo. Serum C-terminal telopeptide of type I collagen (CTX; bone resorption marker) and N-terminal propeptide of type I procollagen (P1NP; bone formation marker) were measured at baseline, six weeks, and 12 weeks post-injury. Clinical and radiological fracture healing was assessed at three months. Results. CTX and P1NP concentrations peaked at six weeks in all groups. Elevated six-week CTX and P1NP were associated with radiological healing at 12 weeks post-injury (odds ratio (OR) 10.5; 95% confidence interval 2.71 to 53.5, p = 0.002). We found no association between CTX or P1NP and functional healing. Baseline serum 25(OH)D showed a weak inverse relationship with P1NP (p = 0.036) and CTX (p = 0.221) at 12 weeks, but we observed no association between vitamin D supplementation and either BTM. Conclusion. Given the association between six-week BTM concentrations and three-month radiological fracture healing, CTX and P1NP appear to be potential surrogate markers of fracture healing. Cite this article: Bone Joint Res 2022;11(4):239–250

Aims. This study aimed to uncover the hub long non-coding RNAs (lncRNAs) differentially expressed in osteoarthritis (OA) cartilage using an integrated analysis of the competing endogenous RNA (ceRNA) network and co-expression network. Methods. Expression profiles data of ten OA and ten normal tissues of human knee cartilage were obtained from the Gene Expression Omnibus (GEO) database (GSE114007). The differentially expressed messenger RNAs (DEmRNAs) and lncRNAs (DElncRNAs) were identified using the edgeR package. We integrated human microRNA (miRNA)-lncRNA/mRNA interactions with DElncRNA/DEmRNA expression profiles to construct a ceRNA network. Likewise, lncRNA and mRNA expression profiles were used to build a co-expression network with the WGCNA package. Potential hub lncRNAs were identified based on an integrated analysis of the ceRNA network and co-expression network. StarBase and Multi Experiment Matrix databases were used to verify the lncRNAs. Results. We detected 1,212 DEmRNAs and 49 DElncRNAs in OA and normal knee cartilage. A total of 75 dysregulated lncRNA-miRNA interactions and 711 dysregulated miRNA-mRNA interactions were obtained in the ceRNA network, including ten DElncRNAs, 69 miRNAs, and 72 DEmRNAs. Similarly, 1,330 dysregulated lncRNA-mRNA interactions were used to construct the co-expression network, which included ten lncRNAs and 407 mRNAs. We finally identified seven hub lncRNAs, named MIR210HG, HCP5, LINC00313, LINC00654, LINC00839, TBC1D3P1-DHX40P1, and ISM1-AS1. Subsequent enrichment analysis elucidated that these lncRNAs regulated extracellular matrix organization and enriched in osteoclast differentiation, the FoxO signalling pathway, and the tumour necrosis factor (TNF) signalling pathway in the development of OA. Conclusion. The integrated analysis of the ceRNA network and co-expression network identified seven hub lncRNAs associated with OA. These lncRNAs may regulate extracellular matrix changes and chondrocyte homeostasis in OA progress. Cite this article:Bone Joint Res. 2020;9(3):90–98

Aims. This study investigates head-neck taper corrosion with varying head size in a novel hip simulator instrumented to measure corrosion related electrical activity under torsional loads. Methods. In all, six 28 mm and six 36 mm titanium stem-cobalt chrome head pairs with polyethylene sockets were tested in a novel instrumented hip simulator. Samples were tested using simulated gait data with incremental increasing loads to determine corrosion onset load and electrochemical activity. Half of each head size group were then cycled with simulated gait and the other half with gait compression only. Damage was measured by area and maximum linear wear depth. Results. Overall, 36 mm heads had lower corrosion onset load (p = 0.009) and change in open circuit potential (OCP) during simulated gait with (p = 0.006) and without joint movement (p = 0.004). Discontinuing gait’s joint movement decreased corrosion currents (p = 0.042); however, wear testing showed no significant effect of joint movement on taper damage. In addition, 36 mm heads had greater corrosion area (p = 0.050), but no significant difference was found for maximum linear wear depth (p = 0.155). Conclusion. Larger heads are more susceptible to taper corrosion; however, not due to frictional torque as hypothesized. An alternative hypothesis of taper flexural rigidity differential is proposed. Further studies are necessary to investigate the clinical significance and underlying mechanism of this finding. Cite this article: Bone Jt Open 2021;2(11):1004–1016