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

Introduction: In attempting to unravel the complex cellular responses leading to prosthetic loosening investigators have been limited to studying gene expression of extracellular molecules about which most is known whereas new microarray technology allows simultaneous expression profiling of thousands of genes from a complex sample such as the membrane formed around loosened hip prostheses. Methods: Two groups of 8 patients were recruited who have undergone primary total hip arthroplasty for osteoarthritis and subsequently developed either septic or aseptic loosening +/− osteolysis. The control group consisted of one group of 5 patients with the same initial diagnosis who had undergone identical procedures, developed no clinical or radiological signs of aseptic or septic loosening, but had come to revision surgery for other complications as defined by the Swedish Hip register: fracture without previous osteolysis, dislocation, technical error, implant fracture, polyethylene wear or pain. Periprosthetic membrane was harvested at the time of revision surgery and subjected to RNA extraction. cDNA was then synthesized and hybridised to a Human Genome u95 Genechip ® array which contains a complete set of known human genes. Data normalisation, data filtering and pattern identification was performed using Genechip®3.1 software (Affymetrix, Santa Clara, CA). Results: This has revealed the involvement of a large number of genes coding for transcriptional regulators upstream from the extracellular and cell-cell signalling molecules already known to be involved in osteolysis and deep infection and which may ultimately control the responses to wear particles and bacterial challenge. Differential expression of genes involved in cell survival and death, cell growth regulation, cell metabolism, inflammation and immune response was found. Most interestingly pathways for control of local bone resorption and inflammatory response have been shown to be highly activated. Conclusions: The identification of these new pathogenetic mechanisms of total hip replacement failure make new indicators of disease susceptibility and prognosis plus new drug targets direct possibilities

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

The concept of translational research is always hampered by the problem that most of the disease phenotypes do not have a mono causal origin. Therefore most treatment schemes based on one to three drugs are not really productive for most of the patients even if the patients are carefully selected from the responder group. Here the array techniques has inspired many research groups to develop algorithms deriving interaction networks or regulatory networks from this type of data to better get rid of the complexity of the biochemical interactions. The challenge is to find networks and to select the group of master nodes which might be good targets for a balanced multi-drug treatment. This means not only to measure one data type with array techniques but to join array data from multiple platforms and different data levels. Our goal is to integrate these data types to form networks with a predictive character for osteosarcomas.

The existing web platform CAPweb/VAMP from the Institute Curie is based on a Java web-client and R. This platform is focused on array data analysis and visualisation, can be extended by additional R modules and is therefore an excellent choice to implement further algorithms for data integration and network prediction. We are now establishing algorithms beyond a pure association of effects like permutation procedures for optimal rank orders of effects in a given subset of 16 factors which can be assembled to bigger units and selection procedures of gene expression signals by gene dosage concepts.

The presented approach is sustainable because the platform can be constantly extended and improved. On the other hand this platform is end-user suitable. This is the best way to bring theoretical concepts to the bench scientist. As a consequence translational research will become more real and complex systems more feasible.

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

The pathophysiological basis of alterations in trabecular bone of patients with osteonecrosis of the femoral head (ONFH) remains unclear. ONFH has classically been considered a vascular disease with secondary changes in the subchondral bone. However, there is increasing evidence suggesting that ONFH could be a bone disease, since alterations in the functionality of bone tissue distant from the necrotic lesion have been observed. We comparatively studied the transcriptomic profile of trabecular bone obtained from the intertrochanteric region of patients with ONFH without an obvious aetiological factor, and patients with osteoarthritis (OA) undergoing total hip replacement in our Institution. To explore the biological processes that could be affected by ONFH, we compared the transcriptomic profile of trabecular bone from the intertrochanteric region and the femoral head of patients affected by this condition. Differential gene expression was studied using an Affymetrix microarray platform. Transcriptome analysis showed a differential signature in trabecular bone from the intertrochanteric region between patients with ONFH and those with OA. The gene ontology analyses of the genes overexpressed in bone tissue of patients with ONFH revealed a range of enriched biological processes related to cell adhesion and migration and angiogenesis. In contrast, most downregulated transcripts were involved in cell division. Trabecular bone in the intertrochanteric region and in the femoral head also exhibited a differential expression profile. Among the genes differentially expressed, we highlighted those related with cytokine production and immune response. This study identified a set of differently expressed genes in trabecular bone of patients with idiopathic ONFH, which might underlie the pathophysiology of this condition. Acknowledgements: This work was supported by grants PI18/00643 and PI22/00939 from ISCIII-FEDER, Ministerio de Ciencia, Innovación y Universidades (MICINN)-AES

In osteoarthritis, chondrocytes acquire a hypertrophic phenotype that contributes to matrix degradation. Inflammation is proposed as trigger for the shift to a hypertrophic phenotype. Using in vitro culture of human chondrocytes and cartilage explants we could not find evidence for a role of inflammatory signalling activation. We found, however, that tissue repair macrophages may contribute to the onset of hypertrophy (doi: 10.1177/19476035211021907) Intra-articularly injected triamcinolone acetonide to inhibit inflammation in a murine model of collagenase-induced osteoarthritis, increased synovial macrophage numbers and osteophytosis, confirming the role of macrophages in chondrocyte hypertrophy occurring in osteophyte formation (doi: 10.1111/bph.15780). In search of targets to inhibit chondrocyte hypertrophy, we combined existing microarray data of different cartilage layers of murine growth plate and murine articular cartilage after induction of collagenase-induced osteoarthritis. We identified common differentially expressed genes and selected those known to be associated to inflammation. This revealed EPHA2, a tyrosine kinase receptor, as a new target. Using in silico, in vitro and in vivo models we demonstrated that inhibition of EPHA2 might be a promising treatment for osteoarthritis. Recently, single cell RNA-seq. has revealed detailed information about different populations of chondrocytes in articular cartilage during osteoarthritis. We re-analysed a published scRNA-seq data set of healthy and osteoarthritic cartilage to obtain the differentially expressed genes in the population of hypertrophic chondrocytes compared to the other chondrocytes, applied pathway analyses and then used drug databases to search for upstream inhibitors of these pathways. This drug repurposing approach led to the selection of 6 drugs that were screened and tested using several in vitro models with human chondrocytes and cartilage explants. In this lecture I will present this sequence of studies to highlight different approaches and models that can be used in the quest for a disease modifying drug for osteoarthritis

Abstract. Introduction. Risk factors for osteoarthritis include raised BMI and female gender. Whether these two factors influenced synovial gene expression was investigated using a triangulation and modelling strategy which generated 12 datasets of gene expression in synovial tissue from three knee pathologies with matching BMI groups, obese and overweight, and gender distributions. Methodology. Intra-operative synovial biopsies were immersed in RNAlater at 4oC before storage at -80oC. Total RNA was extracted using RNAeasy with gDNA removal. Following RT- PCR and quality assessment, cDNA was applied to Affymetrix Clariom D microarray gene chips. Bioinformatics analyses were performed. Linear models were prepared in limma with gender and BMI factors incorporated sequentially for each pathology comparison, generating 12 models of probes differentially expressed at FDR p<0.05 and Bayes number, B>0. Data analysis of differently expressed genes utilized Ingenuity Pathway Analysis and Cytoscape with Cluego and Cytohubba plug-ins. Results. Expression of 453 synovial genes was influenced by BMI and gender, 360 encode proteins such as HIF-1a, HSF1, HSPA4, HSPA5. Top canonical pathways include Unfolded protein response, Protein Ubiquiitation and Clathrin mediated endocytosis signalling linked by modulation of heat shock proteins, comparable to pathology dependent regulation. In addition BMI and gender modulate gene expression in the NRF2-mediated oxidative stress response pathway with down regulation of Glutathione-S-transferases potentially down regulating antioxidant defences. Conclusion. The enhanced risk of osteoarthritis induced by an elevated BMI and female gender maybe include differential expression of heat shock proteins and genes in the NRF2 pathway

Abstract. Introduction. It is increasingly evident that synovium may play a larger role in the aetiology of osteoarthritis. We compared gene expression in whole tissue synovial biopsies from end-stage knee osteoarthritis and knee trauma patients with that of their paired explant cultures to determine how accurately cultured cells represent holistic synovial function. Methodology. Synovial tissue biopsies were taken from 16 arthroplasty patients and 8 tibial plateau fracture patients with no osteoarthritis. Pairs of whole tissue fragments were either immediately immersed in RNAlater Stabilisation Solution at 4o C before transfer to -80o C storage until RNA extraction; or weighed, minced and cultured at 500mg tissues/5ml media in a humidified incubator at 37oC, 5% CO2. After sub-culturing total RNA was extracted using RNAeasy Plus Mini Kit with gDNA removal. Following RT-PCR and quality assessment, cDNA was applied to Affymetrix Clariom D microarray gene chips. Bioinformatics analyses were performed. Results. PCA analysis illustrates the clear separation of expression array data from cultured cells compared with their parental whole tissues and no segregation between cells derived from osteoarthritic or trauma tissues. A differentially expressed gene heat map demonstrated the hierarchical independence of cultured cells from their paired sample parental tissues. The biological pathways enriched by these gene expression differences emphasise the activities of macrophages and lymphocytes lost from culture. Conclusion. Adherent synovial cells grown from different knee pathologies lose the expression patterns characteristic of their originating pathology. Interpretation of data needs caution as the cells are not representative of whole synovium

Abstract. Introduction. Synovitis impacts osteoarthritis symptomatology and progression. The transcription factors controlling synovial gene expression have not been described. This study analyses gene expression in synovium samples from 16 patients with osteoarthritis with 9 undergoing arthroscopic and 8 knee trauma surgery for non-arthritic pathologies. Methodology. Intra-operative synovial biopsies were immersed in RNAlater at 4oC before storage at -80oC. Total RNA was extracted using RNAeasy. After purification, RT-PCR and quality assessment, cDNA was applied to Affymetrix Clariom D microarray gene chips. Bioinformatics analyses were performed. Linear models were prepared in limma with gender and BMI factors incorporated sequentially for each pathology comparison, generating 12 models of probes differentially expressed at FDR p<0.05 and Bayes number, B>0. Data analysis of differently expressed genes utilized Ingenuity Pathway Analysis and Cytoscape with Cluego and Cytohubba plug-ins. Results. Amongst the 2084 genes with significantly differential expression (DEG), 135 had transcription regulator capabilities and 121 a nuclear location. IPA analysis of OATKR and arthroscopic tissue comparison DEG identified 12 nuclear transcription factors linked to 31 DEG whose encoded proteins located within cytoplasmic and cell membrane compartments. All 12 were significantly up-regulated and acting in pathways up-regulating transcription of DNA and RNA, cell survival and angiogenesis while down-regulating senescence and apoptosis. NFE2L2, integral to the TGF-beta signalling pathway, was identified as a bottleneck gene. Conclusion. This analysis indicates the complexity of synovial gene expression regulation and offers target genes and pathways for evaluation during osteoarthritis pathogenesis

Osteoarthritis (OA) is the most common form of arthritis and one of the ten most disabling diseases in developed countries. Total joint replacement (TJR) is considered by far as the most effective treatment for end-stage OA patients. The majority of patients achieve symptomatic improvement following TJR. However, about 22% of the TJR patients either do not improve or deteriorate after surgery. Several potential non-genetic predictors for the TJR outcome have been investigated. However, the results were either inconclusive or had very limited predictive power. The aim of this study was to identify genetic variants for the poor outcome of TJR in primary OA patients by a genome-wide association study (GWAS). Study participants were total knee or hip replacement patients due to primary OA who were recruited to the Newfoundland Osteoarthritis Study (NFOAS) before 2017. The Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) was used to assess pain and functional impairment pre- and 3.99±1.38 years post-surgery. Two non-responder classification criteria were used in our study. One was defined by an absolute WOMAC change score. Participants with a change score less than 7/20 points for pain were considered as pain non-responders; and those with less than 22/68 points for function were classified as function non-responders. The second one was the Outcome Measures in Arthritis Clinical Trials and the Osteoarthritis Research Society International (OMERACT-OARSI) criteria. Blood DNA samples were genotyped using the Illumina GWAS microarrays genotyping platform. The quality control (QC) filtering was performed on GWAS data before the association of the genetic variants with non-responders to TJR was tested using the GenABEL package in R with adjustment for the relatedness of the study population and using the commonly accepted GWAS significance threshold p < 5*10. −8. to control multiple testing. In total, 316 knee and 122 hip OA patients (mean age 65.45±7.62 years, and 58% females) passed the QC check. These study participants included 368 responders and 56 non-responders to pain, and 364 responders and 68 non-responders to function based on the absolute WOMAC point score change classification. While 377 responders and 56 non-responders to pain, and 366 responders and 71 non-responders to function were identified by the OMERACT-OARSI classification criteria. Interestingly, the same results were obtained by both classification methods, and we found that the G allele of rs4797006 was significantly associated with pain non-responders with odds ratio (OR) of 5.12 (p<7.27×10. -10. ). This SNP is in intron one of the melanocortin receptor 5 (MC5R) gene on chr18. This gene plays central roles in immune response, pain sensitivity, and negative regulation of inflammatory response to antigenic stimulus. The A allele of rs200752023 was associated with function non-responders with OR of 4.41 (p<3.29×10. -8. ). The SNP is located in intron three of the RNA Binding Fox-1 Homolog 3 (RBFOX3) gene on chr17 which has been associated with numerous neurological disorders. Our data suggested that two chromosomal regions are associated with TJR poor outcomes and could be the novel targets for developing strategies to improve the outcome of the TJR

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

Periprosthetic joint infection (PJI) is a potentially devastating complication of joint replacement surgery. Osteocytes comprise 90–95% of all cells in hard bone tissue, are long-lived and are becoming increasingly recognised as a critical cell type in the regulation of bone and systemic physiology. The purpose of this study was to examine role of these cells in PJI pathophysiology and aetiology, with the rationale that their involvement could contribute to the difficulty in detecting and clearing PJI. This study examined the ability of human osteocytes to become infected by Staphylococcus aureus and the responses of both the host cell and pathogen in this scenario. Several S. aureus (MRSA) strains were tested for their ability to infect human primary osteocyte-like cells in vitro and human bone samples ex vivo. Bone biopsies were retrieved from patients undergoing revision total hip arthroplasty for either aseptic loosening associated with osteolysis, or for PJI. Retrieved bacterial colony number from cell lysates and colony morphology were determined. Gene expression was measured by microarray/bioinformatics analysis and/or real-time RT-PCR. Exposure to planktonic S. aureus (approx. 100 CFU/cell) resulted in intracellular infection of human osteocyte-like cells. We found no evidence of increased rates of osteocyte cell death in bacteria exposed cultures. Microarray analysis of osteocyte gene expression 24h following exposure revealed more than 1,500 differentially expressed genes (fold-change more than 2, false discovery rate p < 0.01). The gene expression patterns were consistent with a strong innate immune response and altered functionality of the osteocytes. Consistent patterns of host gene expression were observed between experimentally infected osteocyte-like cultures and human bone, and in PJI patient bone samples. Internalised bacteria switched to the quasi-dormant small colony variant (SCV) form over a period of 5d, and the ensuing infection appeared to reach a stable state. S. aureus infection of viable osteocytes was also identified in bone taken from PJI patients. We have demonstrated [1] that human osteocytes can become infected by S. aureus and respond robustly by producing immune mediators. The bony location of the infected osteocyte may render them refractory to clearance by immune cells, and osteocytes may therefore be an immune-privileged cell type. The phenotypic switch of S. aureus to SCV, a form less sensitive to most antibiotics and one associated with intracellular survival, suggests that infection of osteocytes may contribute to a chronic disease state. The osteocyte may therefore serve as a reservoir of bacteria for reinfection, perhaps explaining the high prevalence of infections that only become apparent after long periods of time or recur following surgical/medical treatment. Our findings also provide a biological rationale for the recognised need for aggressive bone debridement in the surgical management of PJI

Osteoarthritis (OA) is the fastest growing global health problem, with a total joint replacement being the only effective treatment for patients with end stage OA. Many groups are examining the use of bone marrow or adipose derived mesenchymal stem cells (MSCs) to repair cartilage, or modulate inflammation to promote healing, however, little efficacy in promoting cartilage repair, or reducing patient symptoms over temporary treatments such as micro-fracture has been observed. There is a growing body of literature demonstrating that MSCs derived from the synovial lining of the joint are superior in terms of chondrogenic differentiation and while improvements in clinical outcome measures have been observed with synovial MSCs, results from clinical studies are still highly variable. Based on our results, we believe this variability in clinical studies with MSCs results in part from the isolation, expansion and re-injection of distinct MSCs subtypes in normal vs. OA tissues, each with differing regenerating potential. However, it remains unknown if this heterogeneity is natural (e.g. multiple MSC subtypes present) or if MSCs are influenced by factors in vivo (disease state/stage). Therefore, in this study, we undertook an ‘omics’ screening approach on MSCs from normal and OA knee synovial tissue. Specifically, we characterized their global proteome and genomic expression patterns to determine if multiple MSC from normal and OA joints are distinct at the protein/gene expression level and/if so, what proteins/genes are differentially expressed between MSCs derived from normal and OA synovial tissue. Synovium tissue was collected from OA patients undergoing joint replacement and normal cadaveric knees. The in vitro adipogenic, chondrogenic and osteogenic differentiation potential of the MSCs was analyzed via qPCR and histology. Fully characterized MSC populations where then analyzed through an unbiased shotgun proteomics, and microarray analysis. Synovial MSCs isolated from both OA and normal knees demonstrated similar multipotent differentiation capacity. Likewise, both OA and normal MSCs display the typical MSCs cell surface marker profile in vitro (CD90+, CD44+, CD73+, CD105+). Using shotgun proteomics, 7720 unique peptides corresponding to 2183 proteins were identified and quantified between normal and OA MSCs. Of these 2183 proteins, 994 were equally expressed in normal and OA, MSCs, 324 were upregulated in OA MSCs (with 50 proteins exclusively expressed in OA MSCs), 630 proteins were upregulated in normal MSCs (with 16 proteins exclusively expressed in normal MSCs). Microarray analysis of normal and OA MSCs demonstrated a similar result in where, 967 genes were differentially expressed between normal and OA MSCs, with 423 genes upregulated in OA, and 544 genes upregulated in normal MSCs. In this project, we have demonstrated that although normal and OA synovial derived MSCs demonstrate similar multipotent differentiation potential and cell surface markers expression, these cells demonstrated significant differences at the molecular level (protein and gene expression). Further research is required to determine if these differences influence functional differences in vitro and/or in vivo and what drives this dramatic change in the regulatory pathways within normal vs. OA synovial MSCs

Commonly used alterations of prosthetic surfaces include grit-blasting (GB), plasma-sprayed titanium (Ti) or hydroxyapatite (HA) coating. Systemic concentrations of cobalt (Co) and chromium (Cr) are elevated in patients with metal-on-metal hip replacement, but can occur for all modular hip replacements. Here, we use whole genome microarrays to assess differential gene expression in primary human osteoblasts grown in vitro and on these prosthesis surfaces following exposure to clinically relevant concentrations of Co and Cr. Mesenchymal cells obtained from bone-fragments of 3 patients undergoing joint replacement surgery were differentiated into osteoblasts. Subsequently, cells were cultured in vitro on tissue-culture plates (TCP), or on GB, Ti and HA surfaces (JRI Orthopaedics Ltd, Sheffield, UK). Following 24hr exposure to a combination of clinically equivalent concentrations of Co2+:Cr3+, RNA was extracted and hybridized to SurePrint-G3 Gene Expression Microarray. Probe signals were normalised using ‘Limma’ package on R-Bioconductor and differential gene expression assessed with empirical Bayes approach (Log2FC>1.00, P<0.001 for differentially expressed genes). For cells grown on TCP, 11 genes were upregulated with 500μg/L Co2+:Cr3+. Of these, 4 were associated to HIF-1 signalling based on KEGG pathway analysis (P=5.4e-5). Exposure to 1000μg/L Co2+:Cr3+ altered expression at 164 loci for HA surfaces, and a separate 50 loci for Ti surfaces compared to GB surfaces. Genes for osteoblast differentiation (BMP2 and RGS2) were downregulated on HA surfaces compared to GB, whilst genes for cell-adhesion (ESAM), vesicular trafficking (RAB37) and protection against oxidative damage (NRF2) were upregulated. Ti surfaces caused an upregulation in ERBB3 and CNTF, which are associated with inhibition of osteoblast differentiation and mineralisation, when compared to GB surfaces. This study confirms the role of HIF-1 signalling in response to prosthesis generated metal ions, and is the first to provide a comprehensive genome-wide insight into transcriptional response of osteoblasts at prosthesis surface to clinically equivalent metal exposure

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

It is supposed that disturbed vascularization is a major cause for the development of an atrophic non-union. However, an actual study revealed normal vessel formation in human non-union tissues [1]. An animal study using an atrophic non-union model should clarify the influence of the inhibition of angiogenesis by the inhibitor Fumagillin on bone healing and the underlying processes including inflammation, chondrogenesis, angiogenesis and osteogenesis. For each group and time point (3, 7, 14, 21 and 42 days) 5–6 adult female Sprague Dawley rats were analyzed. The tibia was osteotomized and stabilized intramedullary with a k-wire coated with the drug carrier PDLLA (control group) or PDLLA +10% Fumagillin (atrophy group). Microarrays: Total-RNA were pooled per group, labeled with the Agilent single-color Quick-Amp Labeling Kit Cy3 and hybridized on Agilent SurePrint G3 Rat Gene Expression microarrays. After feature extraction and quantile normalization, relevant biological processes were identified using GeneOntology. Genes with an expression value below the 25. percentile were excluded. Heatmaps were used for visualization. The analysis of inflammatory genes revealed an upregulation of monocyte/macrophage- relevant factors such as the chemokines Ccl2 and Ccl12 and the surface marker CD14. Other factors involved in the early inflammation process such as Il1a, Tnf and Il6 were not affected. Chondrogenic markers including Collagen Type II, -IX, -X, Mmp9, Mmp13, Hapln1, Ucma, Runx2, Sox5 and -9 were downregulated in this group. Furthermore, osteogenic factors were less regulated within the middle stage of healing (day 14–21). This gene panel included Bmps, Bmp antagonists, Bmp- and Tgfb receptors, integrines and matrix proteins. qPCR analysis of angiogenic genes showed an upregulation of Angpt2, Fgf1 and -2, but not for Vegfa over the later healing time points. We demonstrated in a previous study that inhibiting angiogenesis in an osteotomy model led to a reduction in vessel formation and to the development of an atrophic non-union phenotype [2]. The microarray analysis indicated no prolonged inflammatory reaction in the atrophy group. But the upregulation of chemokines together with a delay in hematoma degradation signs to a mismatch between recruitment and demand of macrophages from the vessel system. Furthermore, chondrogenesis was completely blocked, which was shown by a downregulation of chondrogenic but also osteogenic markers being involved in chondrogenic processes. A reduced recruitment of MSCs might be a possible explanation. Although, microarray data revealed only minor expression changes regarding angiogenic genes, validation by q-PCR showed an upregulation of Angpt2, Fgf1 and -2 over the later healing time points. Due to the heterogeneity of the callus tissue it might be that variations of gene expression of a single tissue type will be masked by the expression levels of other tissue types. This issue is even more pronounced when analyzing different time points and by pooling the samples