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. 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. 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. 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. 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 aim of this study was to investigate the global and local impact of fat on bone in obesity by using the diet-induced obese (DIO) mouse model. Methods. In this study, we generated a diet-induced mouse model of obesity to conduct lipidomic and 3D imaging assessments of bone marrow fat, and evaluated the correlated bone adaptation indices and bone mechanical properties. Results. Our results indicated that bone mass was reduced and bone mechanical properties were impaired in DIO mice. Lipidomic sequencing and bioinformatic analysis identified 373 differential lipids, 176 of which were upregulated and 197 downregulated. Functional enrichment analysis revealed a significant downregulation of the pathways: fat digestion and absorption (ko04975) and lipolysis regulation in adipocytes (ko04923) in DIO mice, leading to local fat accumulation. The use of 3D imaging confirmed the increase in fat accumulation within the bone marrow cavity of obese mice. Conclusion. Our study sheds light on the intricate interplay between fat and bone, and provides a non-toxic and non-invasive method for measuring marrow adipose tissue. Cite this article: Bone Joint Res 2023;12(9):580–589

Aims. The metabolic variations between the cartilage of osteoarthritis (OA) and Kashin-Beck disease (KBD) remain largely unknown. Our study aimed to address this by conducting a comparative analysis of the metabolic profiles present in the cartilage of KBD and OA. Methods. Cartilage samples from patients with KBD (n = 10) and patients with OA (n = 10) were collected during total knee arthroplasty surgery. An untargeted metabolomics approach using liquid chromatography coupled with mass spectrometry (LC-MS) was conducted to investigate the metabolomics profiles of KBD and OA. LC-MS raw data files were converted into mzXML format and then processed by the XCMS, CAMERA, and metaX toolbox implemented with R software. The online Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate the metabolites by matching the exact molecular mass data of samples with those from the database. Results. A total of 807 ion features were identified for KBD and OA, including 577 positive (240 for upregulated and 337 for downregulated) and 230 negative (107 for upregulated and 123 for downregulated) ions. After annotation, LC-MS identified significant expressions of ten upregulated and eight downregulated second-level metabolites, and 183 upregulated and 162 downregulated first-level metabolites between KBD and OA. We identified differentially expressed second-level metabolites that are highly associated with cartilage damage, including dimethyl sulfoxide, uric acid, and betaine. These metabolites exist in sulphur metabolism, purine metabolism, and glycine, serine, and threonine metabolism. Conclusion. This comprehensive comparative analysis of metabolism in OA and KBD cartilage provides new evidence of differences in the pathogenetic mechanisms underlying cartilage damage in these two conditions. Cite this article: Bone Joint Res 2024;13(7):362–371

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. 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. 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. 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. 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. 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 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 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. 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. 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. 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. 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. 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. 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. To assess the effect of physical exercise (PE) on the histological and transcriptional characteristics of proteoglycan-induced arthritis (PGIA) in BALB/c mice. Methods. Following PGIA, mice were subjected to treadmill PE for ten weeks. The tarsal joints were used for histological and genetic analysis through microarray technology. The genes differentially expressed by PE in the arthritic mice were obtained from the microarray experiments. Bioinformatic analysis in the DAVID, STRING, and Cytoscape bioinformatic resources allowed the association of these genes in biological processes and signalling pathways. Results. Arthritic mice improved their physical fitness by 42.5% after PE intervention; it induced the differential expression of 2,554 genes. The bioinformatic analysis showed that the downregulated genes (n = 1,371) were significantly associated with cellular processes that mediate the inflammation, including Janus kinase-signal transducer and activator of transcription proteins (JAK-STAT), Notch, and cytokine receptor interaction signalling pathways. Moreover, the protein interaction network showed that the downregulated inflammatory mediators interleukin (IL) 4, IL5, IL2 receptor alpha (IL2rα), IL2 receptor beta (IL2rβ), chemokine ligand (CXCL) 9, and CXCL12 were interacting in several pathways associated with the pathogenesis of arthritis. The upregulated genes (n = 1,183) were associated with processes involved in the remodelling of the extracellular matrix and bone mineralization, as well as with the processes of aerobic metabolism. At the histological level, PE attenuated joint inflammatory infiltrate and cartilage erosion. Conclusion. Physical exercise influences parameters intimately linked to inflammatory arthropathies. Research on the effect of PE on the pathogenesis process of arthritis is still necessary for animal and human models. Cite this article:Bone Joint Res. 2020;9(1):36–48

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

Osteoarthritis (OA), one of the most common motor system disorders, is a degenerative disease involving progressive joint destruction caused by a variety of factors. At present, OA has become the fourth most common cause of disability in the world. However, the pathogenesis of OA is complex and has not yet been clarified. Long non-coding RNA (lncRNA) refers to a group of RNAs more than 200 nucleotides in length with limited protein-coding potential, which have a wide range of biological functions including regulating transcriptional patterns and protein activity, as well as binding to form endogenous small interference RNAs (siRNAs) and natural microRNA (miRNA) molecular sponges. In recent years, a large number of lncRNAs have been found to be differentially expressed in a variety of pathological processes of OA, including extracellular matrix (ECM) degradation, synovial inflammation, chondrocyte apoptosis, and angiogenesis. Obviously, lncRNAs play important roles in regulating gene expression, maintaining the phenotype of cartilage and synovial cells, and the stability of the intra-articular environment. This article reviews the results of the latest research into the role of lncRNAs in a variety of pathological processes of OA, in order to provide a new direction for the study of OA pathogenesis and a new target for prevention and treatment. Cite this article: Bone Joint Res 2021;10(2):122–133

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. Osteoporosis is a metabolic disease resulting in progressive loss of bone mass as measured by bone mineral density (BMD). Physical exercise has a positive effect on increasing or maintaining BMD in postmenopausal women. The contribution of exercise to the regulation of osteogenesis in osteoblasts remains unclear. We therefore investigated the effect of exercise on osteoblasts in ovariectomized mice. Methods. We compared the activity of differentially expressed genes of osteoblasts in ovariectomized mice that undertook exercise (OVX+T) with those that did not (OVX), using microarray and bioinformatics. Results. Many inflammatory pathways were significantly downregulated in the osteoblasts after exercise. Meanwhile, IBSP and SLc13A5 gene expressions were upregulated in the OVX+T group. Furthermore, in in vitro assay, IBSP and SLc13A5 mRNAs were also upregulated during the osteogenic differentiation of MC3T3-E1 and 7F2 cells. Conclusion. These findings suggest that exercise may not only reduce the inflammatory environment in ovariectomized mice, indirectly suppressing the overactivated osteoclasts, but may also directly activate osteogenesis-related genes in osteoblasts. Exercise may thus prevent the bone loss caused by oestrogen deficiency through mediating the imbalance between the bone resorptive activity of osteoclasts and the bone formation activity of osteoblasts. Cite this article: W-B. Hsu, W-H. Hsu, J-S. Hung, W-J. Shen, R. W-W. Hsu. Transcriptome analysis of osteoblasts in an ovariectomized mouse model in response to physical exercise. Bone Joint Res 2018;7:601–608. DOI: 10.1302/2046-3758.711.BJR-2018-0075.R2

Objectives. Prosthetic joint infection (PJI) diagnosis is a major challenge in orthopaedics, and no reliable parameters have been established for accurate, preoperative predictions in the differential diagnosis of aseptic loosening or PJI. This study surveyed factors in synovial fluid (SF) for improving PJI diagnosis. Methods. We enrolled 48 patients (including 39 PJI and nine aseptic loosening cases) who required knee/hip revision surgery between January 2016 and December 2017. The PJI diagnosis was established according to the Musculoskeletal Infection Society (MSIS) criteria. SF was used to survey factors by protein array and enzyme-linked immunosorbent assay to compare protein expression patterns in SF among three groups (aseptic loosening and first- and second-stage surgery). We compared routine clinical test data, such as C-reactive protein level and leucocyte number, with potential biomarker data to assess the diagnostic ability for PJI within the same patient groups. Results. Cut-off values of 1473 pg/ml, 359 pg/ml, and 8.45 pg/ml were established for interleukin (IL)-16, IL-18, and cysteine-rich with EGF-like domains 2 (CRELD2), respectively. Receiver operating characteristic curve analysis showed that these factors exhibited an accuracy of 1 as predictors of PJI. These factors represent potential biomarkers for decisions associated with prosthesis reimplantation based on their ability to return to baseline values following the completion of debridement. Conclusion. IL-16, IL-18, and CRELD2 were found to be potential biomarkers for PJI diagnosis, with SF tests outperforming blood tests in accuracy. These factors could be useful for assessing successful debridement based on their ability to return to baseline values following the completion of debridement. Cite this article: M-F. Chen, C-H. Chang, L-Y. Yang, P-H. Hsieh, H-N. Shih, S. W. N. Ueng, Y. Chang. Synovial fluid interleukin-16, interleukin-18, and CRELD2 as novel biomarkers of prosthetic joint infections. Bone Joint Res 2019;8:179–188. DOI: 10.1302/2046-3758.84.BJR-2018-0291.R1

Aims

Rheumatoid arthritis (RA) is a common chronic immune disease. Berberine, as its main active ingredient, was also contained in a variety of medicinal plants such as Berberaceae, Buttercup, and Rutaceae, which are widely used in digestive system diseases in traditional Chinese medicine with anti-inflammatory and antibacterial effects. The aims of this article were to explore the therapeutic effect and mechanism of berberine on rheumatoid arthritis.

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

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.

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.

Cite this article: Bone Joint Res 2023;12(6):372–374.

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

The aim of this study was to determine the fracture haematoma (fxH) proteome after multiple trauma using label-free proteomics, comparing two different fracture treatment strategies.

Aims

To examine how eukaryotic translation initiation factor 5A (eIF5A) regulates osteoarthritis (OA) during mechanical overload and the specific mechanism.

Objectives. Diabetes mellitus (DM) is known to impair fracture healing. Increasing evidence suggests that some microRNA (miRNA) is involved in the pathophysiology of diabetes and its complications. We hypothesized that the functions of miRNA and changes to their patterns of expression may be implicated in the pathogenesis of impaired fracture healing in DM. Methods. Closed transverse fractures were created in the femurs of 116 rats, with half assigned to the DM group and half assigned to the control group. Rats with DM were induced by a single intraperitoneal injection of streptozotocin. At post-fracture days five, seven, 11, 14, 21, and 28, miRNA was extracted from the newly generated tissue at the fracture site. Microarray analysis was performed with miRNA samples from each group on post-fracture days five and 11. For further analysis, real-time polymerase chain reaction (PCR) analysis was performed at each timepoint. Results. Microarray analysis showed that there were 14 miRNAs at day five and 17 miRNAs at day 11, with a greater than twofold change in the DM group compared with the control group. Among these types of miRNA, five were selected based on a comparative and extended literature review. Real-time PCR analysis revealed that five types of miRNA (miR-140-3p, miR-140-5p, miR-181a-1-3p, miR-210-3p, and miR-222-3p) were differentially expressed with changing patterns of expression during fracture healing in diabetic rats compared with controls. Conclusions. Our findings provide information to further understand the pathology of impaired fracture healing in a diabetic rat model. These results may allow the potential development of molecular therapy using miRNA for the treatment of impaired fracture healing in patients with DM. Cite this article: S. Takahara, S. Y. Lee, T. Iwakura, K. Oe, T. Fukui, E. Okumachi, T. Waki, M. Arakura, Y. Sakai, K. Nishida, R. Kuroda, T. Niikura. Altered expression of microRNA during fracture healing in diabetic rats. Bone Joint Res 2018;7:139–147. DOI: 10.1302/2046-3758.72.BJR-2017-0082.R1

Aims

To investigate the optimal thresholds and diagnostic efficacy of commonly used serological and synovial fluid detection indexes for diagnosing periprosthetic joint infection (PJI) in patients who have rheumatoid arthritis (RA).

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

Aims

Circular RNA (circRNA) is involved in the regulation of articular cartilage degeneration induced by inflammatory factors or oxidative stress. In a previous study, we found that the expression of circStrn3 was significantly reduced in chondrocytes of osteoarthritis (OA) patients and OA mice. Therefore, the aim of this paper was to explore the role and mechanism of circStrn3 in osteoarthritis.

Osteoarthritis (OA) is a chronic degenerative joint disease characterized by progressive cartilage degradation, synovial membrane inflammation, osteophyte formation, and subchondral bone sclerosis. Pathological changes in cartilage and subchondral bone are the main processes in OA. In recent decades, many studies have demonstrated that activin-like kinase 3 (ALK3), a bone morphogenetic protein receptor, is essential for cartilage formation, osteogenesis, and postnatal skeletal development. Although the role of bone morphogenetic protein (BMP) signalling in articular cartilage and bone has been extensively studied, many new discoveries have been made in recent years around ALK3 targets in articular cartilage, subchondral bone, and the interaction between the two, broadening the original knowledge of the relationship between ALK3 and OA. In this review, we focus on the roles of ALK3 in OA, including cartilage and subchondral bone and related cells. It may be helpful to seek more efficient drugs or treatments for OA based on ALK3 signalling in future.

Aims

Fracture-related infections (FRIs) are a major concern for patients and healthcare systems, yet their impact on mental health has been largely overlooked. This study aimed to assess the longitudinal impact of FRI on patients’ quality of life.

Aims

Metal particles detached from metal-on-metal hip prostheses (MoM-THA) have been shown to cause inflammation and destruction of tissues. To further explore this, we investigated the histopathology (aseptic lymphocyte-dominated vasculitis-associated lesions (ALVAL) score) and metal concentrations of the periprosthetic tissues obtained from patients who underwent revision knee arthroplasty. We also aimed to investigate whether accumulated metal debris was associated with ALVAL-type reactions in the synovium.

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.

Aims

The management of periprosthetic joint infection (PJI) remains a major challenge in orthopaedic surgery. In this study, we aimed to characterize the local bone microstructure and metabolism in a clinical cohort of patients with chronic PJI.

Aims

The involvement of long non-coding RNA (lncRNA) in bone marrow mesenchymal stem cell (MSC) osteogenic differentiation during osteoporosis (OP) development has attracted much attention. In this study, we aimed to disclose how LINC01089 functions in human mesenchymal stem cell (hMSC) osteogenic differentiation, and to study the mechanism by which LINC01089 regulates MSC osteogenesis.

Aims

This study aimed to investigate the role and mechanism of meniscal cell lysate (MCL) in fibroblast-like synoviocytes (FLSs) and osteoarthritis (OA).

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.