Osteoarthritis (OA) is mainly caused by ageing, strain, trauma, and congenital joint abnormalities, resulting in articular cartilage degeneration. During the pathogenesis of OA, the changes in subchondral bone (SB) are not only secondary manifestations of OA, but also an active part of the disease, and are closely associated with the severity of OA. In different stages of OA, there were microstructural changes in SB. Osteocytes, osteoblasts, and osteoclasts in SB are important in the pathogenesis of OA. The signal transduction mechanism in SB is necessary to maintain the balance of a stable phenotype, extracellular matrix (ECM) synthesis, and bone remodelling between articular cartilage and SB. An imbalance in signal transduction can lead to reduced cartilage quality and SB thickening, which leads to the progression of OA. By understanding changes in SB in OA, researchers are exploring drugs that can regulate these changes, which will help to provide new ideas for the treatment of OA.

Cite this article: Bone Joint Res 2023;12(9):536–545.

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. Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. Methods. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Results. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p. discovery GWAS. = 1.21 × 10. -25. , p. replication GWAS. = 1.80 × 10. -12. ), CCDC170 (p. discovery GWAS. = 1.23 × 10. -11. , p. replication GWAS. = 3.22 × 10. -9. ), and SOX6 (p. discovery GWAS. = 4.41 × 10. -15. , p. replication GWAS. = 6.57 × 10. -14. ). Then, for the 36 common target genes, multiple gene ontology (GO) terms were detected for BMD such as positive regulation of cartilage development (p = 9.27 × 10. -3. ) and positive regulation of chondrocyte differentiation (p = 9.27 × 10. -3. ). Conclusion. We explored the potential roles of rSNP in the genetic mechanisms of BMD and identified multiple candidate genes. Our study results support the implication of regulatory genetic variants in the development of OP. Cite this article: Bone Joint Res 2023;12(2):147–154

The February 2023 Research Roundup³⁶⁰ looks at: Clinical and epidemiological features of scaphoid fracture nonunion; Routine sterile glove and instrument change at the time of abdominal wound closure to prevent surgical site infection (ChEETAh); Characterization of genetic risk of end-stage knee osteoarthritis treated with total knee arthroplasty; Platelet-rich plasma or autologous blood injection for plantar fasciitis; Volume and outcomes of joint arthroplasty; The hazards of absolute belief in the p-value laid bare.

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

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

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.

Adolescent idiopathic scoliosis (AIS), defined by an age at presentation of 11 to 18 years, has a prevalence of 0.47% and accounts for approximately 90% of all cases of idiopathic scoliosis. Despite decades of research, the exact aetiology of AIS remains unknown. It is becoming evident that it is the result of a complex interplay of genetic, internal, and environmental factors. It has been hypothesized that genetic variants act as the initial trigger that allow epigenetic factors to propagate AIS, which could also explain the wide phenotypic variation in the presentation of the disorder. A better understanding of the underlying aetiological mechanisms could help to establish the diagnosis earlier and allow a more accurate prediction of deformity progression. This, in turn, would prompt imaging and therapeutic intervention at the appropriate time, thereby achieving the best clinical outcome for this group of patients.

Cite this article: Bone Joint J 2022;104-B(8):915–921.

Aims. The aim of this study was to explore the genetic correlation and causal relationship between blood plasma proteins and rheumatoid arthritis (RA). Methods. Based on the genome-wide association studies (GWAS) summary statistics of RA from European descent and the GWAS summary datasets of 3,622 plasma proteins, we explored the relationship between RA and plasma proteins from three aspects. First, linkage disequilibrium score regression (LD score regression) was applied to detect the genetic correlation between RA and plasma proteins. Mendelian randomization (MR) analysis was then used to evaluate the causal association between RA and plasma proteins. Finally, GEO2R was used to screen the differentially expressed genes (DEGs) between patients with RA and healthy controls. Results. We found that seven kinds of plasma proteins had genetic correlations with RA, such as Soluble Receptor for Advanced Glycation End Products (sRAGE) (correlation coefficient = 0.2582, p = 0.049), vesicle transport protein USE1 (correlation coefficient = 0.1337, p = 0.018), and spermatogenesis-associated protein 20 (correlation coefficient = 0.3706, p = 0.018). There was a significant causal relationship between sRAGE and RA. By comparing the genes encoding seven plasma proteins, we found that only USE1 was a DEG associated with RA. Conclusion. Our study identified a set of candidate plasma proteins that showed signals correlated with RA. Since the results of this study need further experimental verification, they should be interpreted with caution. However, we hope that this paper will provide new insights for the discovery of pathogenic genes and RA pathogenesis in the future. Cite this article: Bone Joint Res 2022;11(2):134–142

Aims

Deciphering the genetic relationships between major depressive disorder (MDD) and osteoarthritis (OA) may facilitate an understanding of their biological mechanisms, as well as inform more effective treatment regimens. We aim to investigate the mechanisms underlying relationships between MDD and OA in the context of common genetic variations.

Aims

Despite the interest in the association of gut microbiota with bone health, limited population-based studies of gut microbiota and bone mineral density (BMD) have been made. Our aim is to explore the possible association between gut microbiota and BMD.

Aims

The aim of this study was to screen the entire genome for genetic markers associated with risk for anterior cruciate ligament (ACL) and posterior cruciate ligament (PCL) injury.

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

Cite this article: Bone Joint Res 2021;10(2):122–133.

Aims

Rheumatoid arthritis (RA) is a systematic autoimmune disorder, characterized by synovial inflammation, bone and cartilage destruction, and disease involvement in multiple organs. Although numerous drugs are employed in RA treatment, some respond little and suffer from severe side effects. This study aimed to screen the candidate therapeutic targets and promising drugs in a novel method.

While total joint replacement (TJR) is considered as an effective intervention to relieve pain and restore joint function for end-stage osteoarthritis (OA) patients, a significant proportion of the patients are dissatisfied with their surgery outcomes. The aim of this study was to identify genetic factors that can predict patients who do or do not benefit from these surgical procedures by a genome-wide association study (GWAS). Study participants were derived from the Newfoundland Osteoarthritis Study (NFOAS) which consisted of 1086 TJR patients. Non-responders to TJR was defined as patients who did not reach the minimum clinically important difference (MCID) based on the self administered Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) in terms of pain reduction or function improvment. DNA was extracted from the blood samples of the study participants and genotyped by Illumina GWAS genotyping platform. Over two million single nucleotide polymorphisms (SNPs) across the genome were genotyped and tested for assocition with non-responders. 39 non-responders and 44 age, sex, and BMI matched responders were included in this study. Four chromosome regions on chromosomes 5, 7, 8, and 12 were suggested to be associated with non-responders with p < 1 0–5. The most promising one was on chromosome 5 with the lead SNP rs17118094 (p=1.7×10–6) which can classify 72% of non-responders accurately. The discriminatory power of this SNP alone is very promising as indicated by an area under the curve (AUC) of 0.72 with 95% confidence interval of 0.63 to 0.81, which is much better than any previously studied predictors mentioned above. All the patients who carry two copies of the G allele (minor allele) of rs17118094 were non-responders and 75% of those who carry one copy of the G allele were non-responders. The discriminatory ability of the lead SNPs on chromosomes 7 and 12 were comparable to the one on chromosome 5 with an AUC of 0.74, and 88% of patients who carry two copies of the A allele of rs10244798 on chromosome 7 were non-responders. Similarly, 88% of patients who carry two copies of the C allele of rs10773476 on chromosome 12 were non-responders. While the discriminatory ability of rs9643244 on chromosome 8 was poor with an AUC of 0.26, its strong association with non-responders warrants a further investigation in the region. The study identified four genomic regions harboring genetic factors for non-responders to TJR. The lead SNPs in those regions have great discriminatory ability to predict non-responders and could be used to create a genetic prediction model for clinical unitilty and application

Aims. Osteoarthritis (OA) is the most prevalent joint disease. However, the specific and definitive genetic mechanisms of OA are still unclear. Methods. Tissue-related transcriptome-wide association studies (TWAS) of hip OA and knee OA were performed utilizing the genome-wide association study (GWAS) data of hip OA and knee OA (including 2,396 hospital-diagnosed hip OA patients versus 9,593 controls, and 4,462 hospital-diagnosed knee OA patients versus 17,885 controls) and gene expression reference to skeletal muscle and blood. The OA-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the messenger RNA (mRNA) expression profiles of hip OA and knee OA. Functional enrichment and annotation analysis of identified genes was performed by the DAVID and FUMAGWAS tools. Results. We detected 33 common genes, eight common gene ontology (GO) terms, and one common pathway for hip OA, such as calcium and integrin-binding protein 1 (CIB1) (PTWAS = 0.025, FCmRNA = -1.575 for skeletal muscle), adrenomedullin (ADM) (PTWAS = 0.022, FCmRNA = -4.644 for blood), Golgi apparatus (PTWAS <0.001, PmRNA = 0.012 for blood), and phosphatidylinositol 3' -kinase-protein kinase B (PI3K-Akt) signalling pathway (PTWAS = 0.033, PmRNA = 0.005 for blood). For knee OA, we detected 24 common genes, eight common GO terms, and two common pathways, such as histocompatibility complex, class II, DR beta 1 (HLA-DRB1) (PTWAS = 0.040, FCmRNA = 4.062 for skeletal muscle), Follistatin-like 1 (FSTL1) (PTWAS = 0.048, FCmRNA = 3.000 for blood), cytoplasm (PTWAS < 0.001, PmRNA = 0.005 for blood), and complement and coagulation cascades (PTWAS = 0.017, PmRNA = 0.001 for skeletal muscle). Conclusion. We identified a group of OA-associated genes and pathways, providing novel clues for understanding the genetic mechanism of OA. Cite this article:Bone Joint Res. 2020;9(3):130–138

Objectives

The objective of this study was to investigate the association of four single-nucleotide polymorphisms (SNPs) of the cannabinoid receptor 2 (CNR2) gene, gene-obesity interaction, and haplotype combination with osteoporosis (OP) susceptibility.

Objective. Low back pain (LBP) is a common debilitating condition with great socioeconomic impact. Identifying individuals at risk of LBP is challenging. We have shown IgG N-glycans are associated with LBP. Herewith, we used polygenic risk scores (PRS) from IgG-glycome to test predictability for LBP. Methods. Clusters of IgG-glycans were identified using weighted correlation network approach in TwinsUK (n = 4246). Genome-wide association studies were carried out for the clusters and top associated SNPs (p<5e-8) were extracted. Weighted PRS was calculated as the sum of the number of copies of effect allele from GWAS multiplied by their effect size using the UK Biobank data (n = 350000). The predictive capacity of the PRS for back pain in UK Biobank was estimated using logistic regression. Results. Multiple SNPs were found to be associated with the glycan clusters near genes known to be involved in glycosylation and the inflammatory response (e.g. ST6GAL1, B4GALT1, FUT8). A total of 175 SNPs was used to calculate weighted PRS. In UK Biobank the PRS was a statically significant, but poor, predictor of the risk of back pain (β = 0.126±0.050, p = 0.015, R. 2. = 2.6e-5). The SNPs on chromosome 14 in regulatory regions of FUT8 gene, one of the key governors of core fucosylation, were found to be significantly associated with back pain in UK Biobank (FDR-adjusted p-value < 0.05). Conclusions. These pilot data suggest that genetic component of glycosylation may be associated with the risk of LBP; however, its predictive ability is poor. Conflict of Interest: YSA is a co-owner of Maatschap PolyOmica. GL is a founder and CEO of Genos Glycoscience Research Laboratory. MBF and FMKW declare no conflict of interests. Sources of Funding: The research has been supported by the EC FP7 project PainOmics (grant agreement #602736) and conducted using the UK Biobank Resource (project # 18219)