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

Objectives

Several genome-wide association studies (GWAS) of bone mineral density (BMD) have successfully identified multiple susceptibility genes, yet isolated susceptibility genes are often difficult to interpret biologically. The aim of this study was to unravel the genetic background of BMD at pathway level, by integrating BMD GWAS data with genome-wide expression quantitative trait loci (eQTLs) and methylation quantitative trait loci (meQTLs) data

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

Objectives

Sustained intra-articular delivery of pharmacological agents is an attractive modality but requires use of a safe carrier that would not induce cartilage damage or fibrosis. Collagen scaffolds are widely available and could be used intra-articularly, but no investigation has looked at the safety of collagen scaffolds within synovial joints. The aim of this study was to determine the safety of collagen scaffold implantation in a validated in vivo animal model of knee arthrofibrosis.

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

Objectives

Studies which consider the molecular mechanisms of degeneration and regeneration of cartilaginous tissues are seriously hampered by problematic ribonucleic acid (RNA) isolations due to low cell density and the dense, proteoglycan-rich extracellular matrix of cartilage. Proteoglycans tend to co-purify with RNA, they can absorb the full spectrum of UV light and they are potent inhibitors of polymerase chain reaction (PCR). Therefore, the objective of the present study is to compare and optimise different homogenisation methods and RNA isolation kits for an array of cartilaginous tissues.

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