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

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

The molecular mechanism of rheumatoid arthritis (RA) remains elusive. We conducted a protein-protein interaction network-based integrative analysis of genome-wide association studies (GWAS) and gene expression profiles of RA.

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

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

Objectives. MicroRNAs (miRNAs) have been reported as key regulators of bone formation, signalling, and repair. Fracture healing is a proliferative physiological process where the body facilitates the repair of a bone fracture. The aim of our study was to explore the effects of microRNA-186 (miR-186) on fracture healing through the bone morphogenetic protein (BMP) signalling pathway by binding to Smad family member 6 (SMAD6) in a mouse model of femoral fracture. Methods. Microarray analysis was adopted to identify the regulatory miR of SMAD6. 3D micro-CT was performed to assess the bone volume (BV), bone volume fraction (BVF, BV/TV), and bone mineral density (BMD), followed by a biomechanical test for maximum load, maximum radial degrees, elastic radial degrees, and rigidity of the femur. The positive expression of SMAD6 in fracture tissues was measured. Moreover, the miR-186 level, messenger RNA (mRNA) level, and protein levels of SMAD6, BMP-2, and BMP-7 were examined. Results. MicroRNA-186 was predicted to regulate SMAD6. Furthermore, SMAD6 was verified as a target gene of miR-186. Overexpressed miR-186 and SMAD6 silencing resulted in increased callus formation, BMD and BV/TV, as well as maximum load, maximum radial degrees, elastic radial degrees, and rigidity of the femur. In addition, the mRNA and protein levels of SMAD6 were decreased, while BMP-2 and BMP-7 levels were elevated in response to upregulated miR-186 and SMAD6 silencing. Conclusion. In conclusion, the study indicated that miR-186 could activate the BMP signalling pathway to promote fracture healing by inhibiting SMAD6 in a mouse model of femoral fracture. Cite this article: Bone Joint Res 2019;8:550–562

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

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

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

MicroRNAs (miRNAs ) are small non-coding RNAs that regulate gene expression. We hypothesised that the functions of certain miRNAs and changes to their patterns of expression may be crucial in the pathogenesis of nonunion. Healing fractures and atrophic nonunions produced by periosteal cauterisation were created in the femora of 94 rats, with 1:1 group allocation. At post-fracture days three, seven, ten, 14, 21 and 28, miRNAs were extracted from the newly generated tissue at the fracture site. Microarray and real-time polymerase chain reaction (PCR) analyses of day 14 samples revealed that five miRNAs, miR-31a-3p, miR-31a-5p, miR-146a-5p, miR-146b-5p and miR-223-3p, were highly upregulated in nonunion. Real-time PCR analysis further revealed that, in nonunion, the expression levels of all five of these miRNAs peaked on day 14 and declined thereafter. . Our results suggest that miR-31a-3p, miR-31a-5p, miR-146a-5p, miR-146b-5p and miR-223-3p may play an important role in the development of nonunion. These findings add to the understanding of the molecular mechanism for nonunion formation and may lead to the development of novel therapeutic strategies for its treatment. Cite this article: Bone Joint J 2015; 97-B:1144–51

Aims. Animal models have been developed that allow simulation of post-traumatic joint contracture. One such model involves contracture-forming surgery followed by surgical capsular release. This model allows testing of antifibrotic agents, such as rosiglitazone. Methods. A total of 20 rabbits underwent contracture-forming surgery. Eight weeks later, the animals underwent a surgical capsular release. Ten animals received rosiglitazone (intramuscular initially, then orally). The animals were sacrificed following 16 weeks of free cage mobilisation. The joints were tested biomechanically, and the posterior capsule was assessed histologically and via genetic microarray analysis. Results. There was no significant difference in post-traumatic contracture between the rosiglitazone and control groups (33° (standard deviation (. sd. ) 11) vs 37° (. sd. 14), respectively; p = 0.4). There was no difference in number or percentage of myofibroblasts. Importantly, there were ten genes and 17 pathways that were significantly modulated by rosiglitazone in the posterior capsule. Discussion. Rosiglitazone significantly altered the genetic expression of the posterior capsular tissue in a rabbit model, with ten genes and 17 pathways demonstrating significant modulation. However, there was no significant effect on biomechanical or histological properties. Cite this article: M. P. Abdel. Effectiveness of rosiglitazone in reducing flexion contracture in a rabbit model of arthrofibrosis with surgical capsular release: A biomechanical, histological, and genetic analysis. Bone Joint Res 2016;5:11–17. doi: 10.1302/2046-3758.51.2000593

Objectives. The goal of this study was to determine whether intra-articular administration of the potentially anti-fibrotic agent decorin influences the expression of genes involved in the fibrotic cascade, and ultimately leads to less contracture, in an animal model. Methods. A total of 18 rabbits underwent an operation on their right knees to form contractures. Six limbs in group 1 received four intra-articular injections of decorin; six limbs in group 2 received four intra-articular injections of bovine serum albumin (BSA) over eight days; six limbs in group 3 received no injections. The contracted limbs of rabbits in group 1 were biomechanically and genetically compared with the contracted limbs of rabbits in groups 2 and 3, with the use of a calibrated joint measuring device and custom microarray, respectively. Results. There was no statistical difference in the flexion contracture angles between those limbs that received intra-articular decorin versus those that received intra-articular BSA (66° vs 69°; p = 0.41). Likewise, there was no statistical difference between those limbs that received intra-articular decorin versus those who had no injection (66° vs 72°; p = 0.27). When compared with BSA, decorin led to a statistically significant increase in the mRNA expression of 12 genes (p < 0.01). In addition, there was a statistical change in the mRNA expression of three genes, when compared with those without injection. . Conclusions. In this model, when administered intra-articularly at eight weeks, 2 mg of decorin had no significant effect on joint contractures. However, our genetic analysis revealed a significant alteration in several fibrotic genes. Cite this article: Bone Joint Res 2014;3:82–8

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.

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

Long non-coding RNAs (lncRNAs) act as crucial regulators in osteoporosis (OP). Nonetheless, the effects and potential molecular mechanism of lncRNA PCBP1 Antisense RNA 1 (PCBP1-AS1) on OP remain largely unclear. The aim of this study was to explore the role of lncRNA PCBP1-AS1 in the pathogenesis of OP.

Aims

This study aimed, through bioinformatics analysis and in vitro experiment validation, to identify the key extracellular proteins of intervertebral disc degeneration (IDD).

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.

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.

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.