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 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. 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. 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. 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