Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis

Looking back on 30+ years of the Johnston County Osteoarthritis Project while looking forward with the Johnston County Health Study: A narrative review

Over the last 30 years, knowledge of the epidemiology of osteoarthritis (OA) has dramatically advanced, and Osteoarthritis and Cartilage has been on the forefront of disseminating research findings from large OA cohort studies, including the Johnston County OA Project (JoCoOA). The JoCoOA is a population-based, prospective longitudinal cohort that began roughly 30 years ago with a key focus on understanding prevalence, incidence, and progression of OA, as well as its risk factors, in a predominantly rural population of Black and White adults 45+ years old in a county in the southeastern United States. Selected OA results that will be discussed in this review include racial differences, lifetime risk, biomarkers, mortality, and OA risk factors. The new Johnston County Health Study will also be introduced. This new cohort study of OA and comorbid conditions builds upon current OA knowledge and JoCoOA infrastructure and is designed to reflect changes in demographics and urbanization in the county and the region.

Genomic Determinants of Knee Joint Biomechanics: An Exploration into the Molecular Basis of Locomotor Function, a Narrative Review

In recent years, the nexus between genetics and biomechanics has garnered significant attention, elucidating the role of genomic determinants in shaping the biomechanical attributes of human joints, specifically the knee. This review seeks to provide a comprehensive exploration of the molecular basis underlying knee joint locomotor function. Leveraging advancements in genomic sequencing, we identified specific genetic markers and polymorphisms tied to key biomechanical features of the knee, such as ligament elasticity, meniscal resilience, and cartilage health. Particular attention was devoted to collagen genes like COL1A1 and COL5A1 and their influence on ligamentous strength and injury susceptibility. We further investigated the genetic underpinnings of knee osteoarthritis onset and progression, as well as the potential for personalized rehabilitation strategies tailored to an individual's genetic profile. We reviewed the impact of genetic factors on knee biomechanics and highlighted the importance of personalized orthopedic interventions. The results hold significant implications for injury prevention, treatment optimization, and the future of regenerative medicine, targeting not only knee joint health but joint health in general.

Genetic Risk Factors for End-Stage Hip Osteoarthritis Treated With Total Hip Arthroplasty: A Genome-wide Association Study

BACKGROUND

Although a genetic component to hip osteoarthritis (OA) has been described, focused evaluation of the genetic components of end-stage disease is limited. We present a genomewide association study for patients undergoing total hip arthroplasty (THA) to characterize the genetic risk factors associated with end-stage hip osteoarthritis (ESHO), defined as utilization of the procedure.

METHODS

Patients who underwent primary THA for hip OA were identified in a national patient data repository using administrative codes. Fifteen thousand three hundred and fifty-five patients with ESHO and 374,193 control patients were identified. Whole genome regression of genotypic data for patients who underwent primary THA for hip OA corrected for age, sex, and body mass index (BMI) was performed. Multivariate logistic regression models were used to evaluate the composite genetic risk from the identified genetic variants.

RESULTS

There were 13 significant genes identified. Composite genetic factors resulted in an odds ratio 1.04 for ESHO (P < .001). The effect of genetics was lower than that of age (Odds Ratio (OR): 2.38; P < .001) and BMI (1.81; P < .001).

CONCLUSION

Multiple genetic variants, including 5 novel loci, were associated with end-stage hip OA treated with primary THA. Age and BMI were associated with greater odds of developing end-stage disease when compared to genetic factors.

The identification of distinct protective and susceptibility mechanisms for hip osteoarthritis: findings from a genome-wide association study meta-analysis of minimum joint space width and Mendelian randomisation cluster analyses

BACKGROUND

Hip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest.

METHODS

GWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA.

FINDINGS

50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82-1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life.

INTERPRETATIONS

One group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism.

FUNDING

Primarily funded by the Medical Research Council and Wellcome Trust.

Epigenetics as a mediator of genetic risk in osteoarthritis: role during development, homeostasis, aging, and disease progression

The identification of genomic loci that are associated with osteoarthritis (OA) has provided a starting point for understanding how genetic variation activates catabolic processes in the joint. However, genetic variants can only alter gene expression and cellular function when the epigenetic environment is permissive to these effects. In this review, we provide examples of how epigenetic shifts at distinct life stages can alter the risk for OA, which we posit is critical for the proper interpretation of genome-wide association studies (GWAS). During development, intensive work on the growth and differentiation factor 5 (GDF5) locus has revealed the importance of tissue-specific enhancer activity in controlling both joint development and the subsequent risk for OA. During homeostasis in adults, underlying genetic risk factors may help establish beneficial or catabolic "set points" that dictate tissue function, with a strong cumulative effect on OA risk. During aging, methylation changes and the reorganization of chromatin can "unmask" the effects of genetic variants. The destructive function of variants that alter aging would only mediate effects after reproductive competence and thus avoid any evolutionary selection pressure, as consistent with larger frameworks of biological aging and its relationship to disease. A similar "unmasking" may occur during OA progression, which is supported by the finding of distinct expression quantitative trait loci (eQTLs) in chondrocytes depending on the degree of tissue degradation. Finally, we propose that massively parallel reporter assays (MPRAs) will be a valuable tool to test the function of putative OA GWAS variants in chondrocytes from different life stages.

Analysis of RNA Polyadenylation in Healthy and Osteoarthritic Human Articular Cartilage

Polyadenylation (polyA) defines the 3' boundary of a transcript's genetic information. Its position can vary and alternative polyadenylation (APA) transcripts can exist for a gene. This causes variance in 3' regulatory domains and can affect coding sequence if intronic events occur. The distribution of polyA sites on articular chondrocyte transcripts has not been studied so we aimed to define their transcriptome-wide location in age-matched healthy and osteoarthritic knee articular cartilage. Total RNA was isolated from frozen tissue samples and analysed using the QuantSeq-Reverse 3' RNA sequencing approach, where each read runs 3' to 5' from within the polyA tail into the transcript and contains a distinct polyA site. Differential expression of transcripts was significant altered between healthy and osteoarthritic samples with enrichment for functionalities that were strongly associated with joint pathology. Subsequent examination of polyA site data allowed us to define the extent of site usage across all the samples. When comparing healthy and osteoarthritic samples, we found that differential use of polyadenylation sites was modest. However, in the genes affected, there was potential for the APA to have functional relevance. We have characterised the polyadenylation landscape of human knee articular chondrocytes and conclude that osteoarthritis does not elicit a widespread change in their polyadenylation site usage. This finding differentiates knee osteoarthritis from pathologies such as cancer where APA is more commonly observed.

Activating EGFR Signaling Attenuates Osteoarthritis Development Following Loading Injury in Mice

Posttraumatic osteoarthritis (PTOA) results in joint pain, loss of joint function, and impaired quality of daily life in patients with limited treatment options. We previously demonstrated that epidermal growth factor receptor (EGFR) signaling is essential for maintaining chondroprogenitors during articular cartilage development and homeostasis. Here, we used a nonsurgical, loading-induced PTOA mouse model to investigate the protective action of EGFR signaling. A single bout of cyclic tibial loading at a peak force of 6 N injured cartilage at the posterior aspect of lateral femoral condyle. Similar loading at a peak force of 9 N ruptured the anterior cruciate ligament, causing additional cartilage damage at the medial compartment and ectopic cartilage formation in meniscus and synovium. Constitutively overexpression of an EGFR ligand, heparin binding EGF-like growth factor (HBEGF), in chondrocytes significantly reduced cartilage injury length, synovitis, and pain after 6 N loading and mitigated medial side cartilage damage and ectopic cartilage formation after 9 N loading. Mechanistically, overactivation of EGFR signaling protected chondrocytes from loading-induced apoptosis and loss of proliferative ability and lubricant synthesis. Overexpressing HBEGF in adult cartilage starting right before 6 N loading had similar beneficial effects. In contrast, inactivating EGFR in adult cartilage led to accelerated PTOA progression with elevated cartilage Mankin score and synovitis score and increased ectopic cartilage formation. As a therapeutic approach, we constructed a nanoparticle conjugated with the EGFR ligand TGFα. Intra-articular injections of this nanoconstruct once every 3 weeks for 12 weeks partially mitigated PTOA symptoms in cartilage and synovium after 6 N loading. Our findings demonstrate the anabolic actions of EGFR signaling in maintaining articular cartilage during PTOA development and shed light on developing a novel nanomedicine for PTOA. © 2022 American Society for Bone and Mineral Research (ASBMR).

Targeted Resequencing of Otosclerosis Patients from Different Populations Replicates Results from a Previous Genome-Wide Association Study

Otosclerosis is one of the most common causes of hearing loss in young adults. It has a prevalence of 0.3-0.4% in the European population. Clinical symptoms usually occur between the second and fifth decade of life. Different studies have been performed to unravel the genetic architecture of the disease. Recently, a genome-wide association study (GWAS) identified 15 novel risk loci and replicated the regions of three previously reported candidate genes. In this study, seven candidate genes from the GWAS were resequenced using single molecule molecular inversion probes (smMIPs). smMIPs were used to capture the exonic regions and the 3' and 5' untranslated regions (UTR). Discovered variants were tested for association with the disease using single variant and gene-based association analysis. The single variant results showed that 13 significant variants were associated with otosclerosis. Associated variants were found in five of the seven genes studied here, including AHSG, LINC01482, MARK3, SUPT3H and RELN. Conversely, burden testing did not show a major role of rare variants in the disease. In conclusion, this study was able to replicate five out of seven candidate genes reported in the previous GWAS. This association is likely mainly driven by common variants.

Single nucleotide polymorphism genes and mitochondrial DNA haplogroups as biomarkers for early prediction of knee osteoarthritis structural progressors: use of supervised machine learning classifiers

BACKGROUND

Knee osteoarthritis is the most prevalent chronic musculoskeletal debilitating disease. Current treatments are only symptomatic, and to improve this, we need a robust prediction model to stratify patients at an early stage according to the risk of joint structure disease progression. Some genetic factors, including single nucleotide polymorphism (SNP) genes and mitochondrial (mt)DNA haplogroups/clusters, have been linked to this disease. For the first time, we aim to determine, by using machine learning, whether some SNP genes and mtDNA haplogroups/clusters alone or combined could predict early knee osteoarthritis structural progressors.

METHODS

Participants (901) were first classified for the probability of being structural progressors. Genotyping included SNP genes TP63, FTO, GNL3, DUS4L, GDF5, SUPT3H, MCF2L, and TGFA; mtDNA haplogroups H, J, T, Uk, and others; and clusters HV, TJ, KU, and C-others. They were considered for prediction with major risk factors of osteoarthritis, namely, age and body mass index (BMI). Seven supervised machine learning methodologies were evaluated. The support vector machine was used to generate gender-based models. The best input combination was assessed using sensitivity and synergy analyses. Validation was performed using tenfold cross-validation and an external cohort (TASOAC).

RESULTS

From 277 models, two were defined. Both used age and BMI in addition for the first one of the SNP genes TP63, DUS4L, GDF5, and FTO with an accuracy of 85.0%; the second profits from the association of mtDNA haplogroups and SNP genes FTO and SUPT3H with 82.5% accuracy. The highest impact was associated with the haplogroup H, the presence of CT alleles for rs8044769 at FTO, and the absence of AA for rs10948172 at SUPT3H. Validation accuracy with the cross-validation (about 95%) and the external cohort (90.5%, 85.7%, respectively) was excellent for both models.

CONCLUSIONS

This study introduces a novel source of decision support in precision medicine in which, for the first time, two models were developed consisting of (i) age, BMI, TP63, DUS4L, GDF5, and FTO and (ii) the optimum one as it has one less variable: age, BMI, mtDNA haplogroup, FTO, and SUPT3H. Such a framework is translational and would benefit patients at risk of structural progressive knee osteoarthritis.

Machine Learning-Derived Acetabular Dysplasia and Cam Morphology Are Features of Severe Hip Osteoarthritis: Findings From UK Biobank

The contribution of shape changes to hip osteoarthritis (HOA) remains unclear, as is the extent to which these vary according to HOA severity. In the present study, we used statistical shape modeling (SSM) to evaluate relationships between hip shape and HOA of different severities using UK Biobank DXA images. We performed a cross-sectional study in individuals with left hip dual-energy X-ray absorptiometry (DXA) scans. Statistical shape modeling (SSM) was used to quantify hip shape. Radiographic HOA (rHOA) was classified using osteophyte size and number and joint space narrowing. HOA outcomes ranged in severity from moderate (grade 2) to severe (grade ≥3) rHOA, hospital-diagnosed HOA, and subsequent total hip replacement (THR). Confounder-adjusted logistic regression between the top 10 hip shape modes (HSMs) and OA outcomes was performed. Further models adjusted for alpha angle (AA) and lateral center-edge angle (LCEA), reflecting acetabular dysplasia and cam morphology, respectively. Composite HSM figures were produced combining HSMs associated with separate OA outcomes. A total of 40,311 individuals were included (mean 63.7 years, 47.8% male), of whom 5.7% had grade 2 rHOA, 1.7% grade ≥3 rHOA, 1.3% hospital-diagnosed HOA, and 0.6% underwent THR. Composite HSM figures for grade 2 rHOA revealed femoral neck widening, increased acetabular coverage, and enlarged lesser and greater trochanters. In contrast, grade ≥3 rHOA, hospital-diagnosed HOA, and THR were suggestive of cam morphology and reduced acetabular coverage. Associations between HSMs depicting cam morphology and reduced acetabular coverage and more severe HOA were attenuated by AA and LCEA adjustment, respectively. Relationships between hip shape and HOA differed according to severity. Notably, cam morphology and acetabular dysplasia were features of severe HOA, but unrelated to moderate disease, suggesting possible prognostic utility. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

EGFR Signaling Is Required for Maintaining Adult Cartilage Homeostasis and Attenuating Osteoarthritis Progression

The uppermost superficial zone of articular cartilage is the first line of defense against the initiation of osteoarthritis (OA). We previously used Col2-Cre to demonstrate that epidermal growth factor receptor (EGFR), a tyrosine kinase receptor, plays an essential role in maintaining superficial chondrocytes during articular cartilage development. Here, we showed that EGFR activity in the articular cartilage decreased as mice age. In mouse and human OA samples, EGFR activity was initially reduced at the superficial layer and then resurged in cell clusters within the middle and deep zone in late OA. To investigate the role of EGFR signaling in postnatal and adult cartilage, we constructed an inducible mouse model with cartilage-specific EGFR inactivation (Aggrecan-CreER Egfr^Wa5/flox , Egfr iCKO). EdU incorporation revealed that postnatal Egfr iCKO mice contained fewer slow-cycling cells than controls. EGFR deficiency induced at 3 months of age reduced cartilage thickness and diminished superficial chondrocytes, in parallel to alterations in lubricin production, cell proliferation, and survival. Furthermore, male Egfr iCKO mice developed much more severe OA phenotypes, including cartilage erosion, subchondral bone plate thickening, cartilage degeneration at the lateral site, and mechanical allodynia, after receiving destabilization of the medial meniscus (DMM) surgery. Similar OA phenotypes were also observed in female iCKO mice. Moreover, tamoxifen injections of iCKO mice at 1 month post-surgery accelerated OA development 2 months later. In summary, our data demonstrated that chondrogenic EGFR signaling maintains postnatal slow-cycling cells and plays a critical role in adult cartilage homeostasis and OA progression. © 2022 American Society for Bone and Mineral Research (ASBMR).

Genetics of osteoarthritis

Osteoarthritis genetics has been transformed in the past decade through the application of large-scale genome-wide association scans. So far, over 100 polymorphic DNA variants have been associated with this common and complex disease. These genetic risk variants account for over 20% of osteoarthritis heritability and the vast majority map to non-protein coding regions of the genome where they are presumed to act by regulating the expression of target genes. Statistical fine mapping, in silico analyses of genomics data, and laboratory-based functional studies have enabled the identification of some of these targets, which encode proteins with diverse roles, including extracellular signaling molecules, intracellular enzymes, transcription factors, and cytoskeletal proteins. A large number of the risk variants correlate with epigenetic factors, in particular cartilage DNA methylation changes in cis, implying that epigenetics may be a conduit through which genetic effects on gene expression are mediated. Some of the variants also appear to have been selected as humans adapted to bipedalism, suggesting that a proportion of osteoarthritis genetic susceptibility results from antagonistic pleiotropy, with risk variants having a positive role in joint formation but a negative role in the long-term health of the joint. Although data from an osteoarthritis genetic study has not yet directly led to a novel treatment, some of the osteoarthritis associated genes code for proteins that have available therapeutics. Genetic investigations are therefore revealing fascinating fundamental insights into osteoarthritis and can expose options for translational intervention.

Osteoarthritis Pathophysiology: Therapeutic Target Discovery may Require a Multifaceted Approach

Molecular understanding of osteoarthritis (OA) has greatly increased through careful analysis of tissue samples, preclinical models, and large-scale agnostic "-omic" studies. There is broad acceptance that systemic and biomechanical signals affect multiple tissues of the joint, each of which could potentially be targeted to improve patient outcomes. In this review six experts in different aspects of OA pathogenesis provide their independent view on what they believe to be good tractable approaches to OA target discovery. We conclude that molecular discovery has been high but future transformative studies require a multidisciplinary holistic approach to develop therapeutic strategies with high clinical efficacy.

Associations between adipokines gene polymorphisms and knee osteoarthritis: a meta-analysis

Background

Adipokines gene polymorphisms are speculated to be associated with the risk of knee osteoarthritis (OA), but evidence remains conflicting. This study therefore aimed to examine whether associations exist between adipokines gene polymorphisms and knee OA by considering the evidence collected from eligible studies through a meta-analysis.

Methods

A systematic search was performed on PubMed, Embase, Web of Science, China National Knowledge Infrastructure (CNKI), and Wanfang up to March 31, 2020. Meta-analysis was carried out by focusing on the associations between adipokines gene polymorphisms and knee OA with the allele model, dominant model, and recessive model.

Results

The present meta-analysis included 5 eligible studies for ADIPOQ rs1501299 with 1,021 cases and 1,097 controls, 3 eligible studies for ADIPOQ rs2241766 with 549 cases and 544 controls, 3 eligible studies for LEPR rs1137101 with 808 cases and 856 controls, 2 eligible studies for VISFATIN rs4730153 with 339 cases and 680 controls and 2 eligible studies for VISFATIN rs16872158 with 339 cases and 680 controls. Significant association was observed between LEPR rs1137101 and knee OA in the overall population (recessive: OR = 0.40, 95% CI 0.21–0.79). Limited data revealed that associations may exist between ADIPOQ rs2241766 and knee OA in Asians (dominant: OR = 1.35, 95% CI 1.03–1.78), between VISFATIN rs4730153 and knee OA in Asians (allele: OR = 0.58, 95% CI 0.41–0.83; dominant: OR = 0.57, 95% CI 0.39–0.83), and between VISFATIN rs16872158 and knee OA in Asians (allele: OR = 1.84, 95% CI 1.26–2.68; dominant: OR = 1.94, 95% CI 1.31–2.89).

Conclusions

Adipokines gene polymorphisms may be associated with knee OA. The association was observed in LEPR rs1137101 in the present study. In addition, limited data revealed that associations may also exist in ADIPOQ rs2241766, VISFATIN rs4730153 and VISFATIN rs16872158.

Prospero registration

CRD42020187664.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05111-4.

Hypoxia induces DOT1L in articular cartilage to protect against osteoarthritis

Osteoarthritis is the most prevalent joint disease worldwide, and it is a leading source of pain and disability. To date, this disease lacks curative treatment, as underlying molecular mechanisms remain largely unknown. The histone methyltransferase DOT1L protects against osteoarthritis, and DOT1L-mediated H3K79 methylation is reduced in human and mouse osteoarthritic joints. Thus, restoring DOT1L function seems to be critical to preserve joint health. However, DOT1L-regulating molecules and networks remain elusive, in the joint and beyond. Here, we identified transcription factors and networks that regulate DOT1L gene expression using a potentially novel bioinformatics pipeline. Thereby, we unraveled a possibly undiscovered link between the hypoxia pathway and DOT1L. We provide evidence that hypoxia enhanced DOT1L expression and H3K79 methylation via hypoxia-inducible factor-1 α (HIF1A). Importantly, we demonstrate that DOT1L contributed to the protective effects of hypoxia in articular cartilage and osteoarthritis. Intra-articular treatment with a selective hypoxia mimetic in mice after surgical induction of osteoarthritis restored DOT1L function and stalled disease progression. Collectively, our data unravel a molecular mechanism that protects against osteoarthritis with hypoxia inducing DOT1L transcription in cartilage. Local treatment with a selective hypoxia mimetic in the joint restores DOT1L function and could be an attractive therapeutic strategy for osteoarthritis.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Body Mass Index Differentially Moderates Heritability of Total Joint Replacement Due to Hip and Knee Osteoarthritis: A Cohort Study of 29,893 Swedish Twin Pairs

BACKGROUND

Osteoarthritis and obesity are diseases with high prevalence, and they share common etiologies. We investigated the sex-specific genetic susceptibility to hip and knee osteoarthritis necessitating total joint replacement (TJR), and how body mass index (BMI) moderated the heritability of these osteoarthritis phenotypes.

METHODS

We linked 29,893 twin pairs with information on BMI in the Swedish Twin Registry with the Swedish National Patient Register to identify twins who underwent primary TJR of the hip or knee combined with a concomitant diagnosis of primary osteoarthritis of these joints. Structural equation modeling was used to calculate the heritability of hip and knee osteoarthritis treated with TJR, with estimates adjusted for the first available BMI, birth year, and sex. We also investigated how heritability varied with BMI treated as a continuous variable.

RESULTS

Similar heritability estimates for hip replacement (0.65 [95% confidence interval (CI), 0.59 to 0.70]) and knee replacement (0.57 [95% CI, 0.50 to 0.64]) were found. Heritability decreased with higher BMI in both sexes for hip replacement and in men for knee replacement. In contrast, heritability for knee replacement increased with higher BMI in women; the estimate was 0.37 (90% likelihood interval [LI], 0.25 to 0.49) for a BMI of 20 kg/m2 and 0.87 (90% LI, 0.68 to 0.94) for a BMI of 35 kg/m2.

CONCLUSIONS

In our population, heritability explained, on average, about half of the susceptibility to undergo primary TJR of the hip or knee with the indication of primary osteoarthritis, but it varied with BMI and sex. We demonstrated substantial heritability for knee replacement in obese women.

LEVEL OF EVIDENCE

Prognostic Level III. See Instructions for Authors for a complete description of levels of evidence.

The role of circRNA derived from RUNX2 in the serum of osteoarthritis and its clinical value

Background

Circular RNA (circRNA) has been shown to affect the pathological process of osteoarthritis (OA) and is expected to become a potential marker for disease diagnosis. This study aimed to investigate the association between circRNA derived from the gene of runt‐related transcription factor 2 (RUNX2) and OA risk.

Methods

The expression profile of RUNX2‐derived circRNAs in serum of OA patients was detected. Then, the cytological localization of screened differential circRNAs was studied. Luciferase (LUC) reporter assay was used to identify the microRNA (miRNA) sponge capacity of the circRNAs. Bioinformatics analysis was used to construct the functional pathway of this circRNA‐miRNAs network. And then, the diagnostic value of RUNX2‐derived circRNAs in OA was evaluated.

Results

RUNX2‐derived hsa_circ_0005526 (circ_RUNX2) is significantly highly expressed in OA serum and mainly located in the cytoplasm within the cartilage cell by sponging multiple miRNAs (miR‐498, miR‐924, miR‐361‐3p, and miR‐665). Bioinformatics analysis showed ECM‐receptor interaction pathway ranked the most significant pathway of circ_RUNX2‐miRNAs regulatory network in KEGG database. The ROC curve showed that there may be good diagnostic value of serum circ_RUNX2 in OA.

Conclusion

RUNX2‐derived circ_RUNX2 may be involved in OA development via ECM‐receptor interaction pathways and may be used as potential clinical indicator of OA.

RUNX2 promotes malignant progression in gastric cancer by regulating COL1A1

BACKGROUND

Runt-related transcription factor 2 (RUNX2) is an important gene that has been implicated in the progression of human cancer. Aberrant expression of RUNX2 predicts gastric cancer (GC) metastasis. However, the molecular mechanism of RUNX2 remains unknown.

OBJECTIVE

We hypothesize that RUNX2 promotes GC metastasis by regulating the extracellular matrix component collagen type I alpha 1 (COL1A1).

METHODS

The GEPIA database and immunohistochemical staining of 60 GC tissues were used to analyse the correlations between RUNX2 or COL1A1 expression and clinicopathological features, and the Kaplan-Meier method was used to evaluate survival. RT-PCR, western blotting and immunofluorescence were used to detect RUNX2 and COL1A1 expression in GC cells. Migration and invasion assays were performed to assess the influence of RUNX2 and COL1A1 on metastasis.

RESULTS

RUNX2 and COL1A1 were highly expressed at both the gene and protein levels in GC, and patients who were positive for RUNX2 and COL1A1 had shorter survival. RUNX2 and COL1A1 expression linearly correlated with each other (r= 0.15, p< 0.01) and with clinical stage and lymph node metastasis (p< 0.05). Overexpressing RUNX2in vitro enhanced COL1A1 expression and promoted GC cell invasion and migration, whereas COL1A1 knockdown inhibited the increase in cell metastatic capacity promoted by RUNX2. In vivo, GC cells overexpressing RUNX2 promoted lung metastasis, and the downregulation of COL1A1 reduced the metastasis promoted by RUNX2.

CONCLUSIONS

RUNX2 may promote GC metastasis by regulating COL1A1. RUNX2/COL1A1 can be employed as a novel target for therapy in GC.

WNT Signalling in Osteoarthritis and Its Pharmacological Targeting

Osteoarthritis (OA) is a highly disabling musculoskeletal condition affecting millions of people worldwide. OA is characterised by progressive destruction and irreversible morphological changes of joint tissues and architecture. At molecular level, de-regulation of several pathways contributes to the disruption of tissue homeostasis in the joint. Overactivation of the WNT/β-catenin signalling pathway has been associated with degenerative processes in OA. However, the multiple layers of complexity in the modulation of the signalling and the still insufficient knowledge of the specific molecular drivers of pathogenetic mechanisms have made difficult the pharmacological targeting of this pathway for therapeutic purposes. This review aims to provide an overview of the WNT/β-catenin signalling in OA with a particular focus on its role in the articular cartilage. Pathway components whose targeting showed therapeutic potential will be highlighted and described. A specific section will be dedicated to Lorecivivint, the first inhibitor of the β-catenin-dependent pathway currently in phase III clinical trial as OA-modifying agent.

Causal Discovery in Radiographic Markers of Knee Osteoarthritis and Prediction for Knee Osteoarthritis Severity With Attention-Long Short-Term Memory

The goal of this study is to build a prognostic model to predict the severity of radiographic knee osteoarthritis (KOA) and to identify long-term disease progression risk factors for early intervention and treatment. We designed a long short-term memory (LSTM) model with an attention mechanism to predict Kellgren/Lawrence (KL) grade for knee osteoarthritis patients. The attention scores reveal a time-associated impact of different variables on KL grades. We also employed a fast causal inference (FCI) algorithm to estimate the causal relation of key variables, which will aid in clinical interpretability. Based on the clinical information of current visits, we accurately predicted the KL grade of the patient's next visits with 90% accuracy. We found that joint space narrowing was a major contributor to KOA progression. Furthermore, our causal structure model indicated that knee alignments may lead to joint space narrowing, while symptoms (swelling, grinding, catching, and limited mobility) have little impact on KOA progression. This study evaluated a broad spectrum of potential risk factors from clinical data, questionnaires, and radiographic markers that are rarely considered in previous studies. Using our statistical model, providers are able to predict the risk of the future progression of KOA, which will provide a basis for selecting proper interventions, such as proceeding to joint arthroplasty for patients. Our causal model suggests that knee alignment should be considered in the primary treatment and KOA progression was independent of clinical symptoms.

Association of rs2862851 in TGFA Gene with Peripheral TGFA Levels and the Severity of Knee Osteoarthritis in the Han Chinese Population

Background: Osteoarthritis (OA) is a complex joint disorder characterized by sclerosis of subchondral bone. The knee is one of the most commonly affected joints. Given that the genetic mechanisms underlying knee OA remain elusive, our study aims were to first confirm the association of the TGFA gene alleles with the risk of knee OA and, second, to evaluate the relationship between peripheral TGFA concentrations and knee OA in an independent Han Chinese population. Materials and Methods: We performed a case-control study consisting of 392 knee OA patients and 808 unrelated healthy controls. Single-marker-based association analyses and haplotype-based analyses using 3 single nucleotide polymorphisms (SNPs) were performed to confirm the association of TGFA gene alleles with the risk of knee OA. Furthermore, we used enzyme-linked immunosorbent assay (ELISA) kits to detect the peripheral blood TGFA concentrations in patients and healthy controls and then evaluated the relationships between the TGFA alleles and genotypes with serum TGFA levels. Results: We replicated the genetic association of the rs2862851 T allele with the risk of knee OA (p = 1.68 × 10^-4, OR = 1.41). Moreover, we observed that the peripheral TGFA concentrations were higher in knee OA patients than in healthy controls (p = 8.15 × 10^-13). The peripheral TGFA concentrations were significantly different among the various rs2862851 genotypes for both cases (p = 4.16 × 10^-16) and controls (p = 7.24 × 10^-19). The individuals with the TT genotype in both cases and controls, had the highest peripheral TGFA concentrations. Moreover, with the increase in knee OA grade, peripheral TGFA concentration also increased (p = 1.36 × 10^-72). Conclusion: Our study confirmed the association of the TGFA gene with the risk of knee OA and identified a positive correlation between peripheral TGFA levels and the severity of knee OA in the Han Chinese population, providing clues for understanding the etiology of knee OA.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene

Background

Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology.

Methods

We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees.

Results

We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci.

Conclusions

We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

Of mice and men: converging on a common molecular understanding of osteoarthritis

Despite an increasing burden of osteoarthritis in developed societies, target discovery has been slow and there are currently no approved disease-modifying osteoarthritis drugs. This lack of progress is due in part to a series of misconceptions over the years: that osteoarthritis is an inevitable consequence of ageing, that damaged articular cartilage cannot heal itself, and that osteoarthritis is driven by synovial inflammation similar to that seen in rheumatoid arthritis. Molecular interrogation of disease through ex-vivo tissue analysis, in-vitro studies, and preclinical models have radically reshaped the knowledge landscape. Inflammation in osteoarthritis appears to be distinct from that seen in rheumatoid arthritis. Recent randomised controlled trials, using treatments repurposed from rheumatoid arthritis, have largely been unsuccessful. Genome-wide studies point to defects in repair pathways, which accords well with recent promise using growth factor therapies or Wnt pathway antagonism. Nerve growth factor has emerged as a robust target in osteoarthritis pain in phase 2-3 trials. These studies, both positive and negative, align well with those in preclinical surgical models of osteoarthritis, indicating that pathogenic mechanisms identified in mice can lead researchers to valid human targets. Several novel candidate pathways are emerging from preclinical studies that offer hope of future translational impact. Enhancing trust between industry, basic, and clinical scientists will optimise our collective chance of success.

MiR-155 promotes interleukin-1β-induced chondrocyte apoptosis and catabolic activity by targeting PIK3R1-mediated PI3K/Akt pathway

Osteoarthritis (OA) is a common joint disease characterized by progressive cartilage degradation, in which elevated chondrocyte apoptosis and catabolic activity play an important role. MicroRNA-155 (miR-155) has recently been shown to regulate apoptosis and catabolic activity in some pathological circumstances, yet, whether and how miR-155 is associated with OA pathology remain unexplored. We report here that miR-155 level is significantly up-regulated in human OA cartilage biopsies and also in primary chondrocytes stimulated by interleukin-1β (IL-1β), a pivotal pro-catabolic factor promoting cartilage degradation. Moreover, miR-155 inhibition attenuates and its overexpression promotes IL-1β-induced apoptosis and catabolic activity in chondrocytes in vitro. We also demonstrate that the PIK3R1 (p85α regulatory subunit of phosphoinositide 3-kinase (PI3K)) is a target of miR-155 in chondrocytes, and more importantly, PIK3R1 restoration abrogates miR-155 effects on chondrocyte apoptosis and catabolic activity. Mechanistically, PIK3R1 positively regulates the transduction of PI3K/Akt pathway, and a specific Akt inhibitor reverses miR-155 effects on promoting chondrocyte apoptosis and catabolic activity, phenocopying the results obtained via PIK3R1 knockdown, hence establishing that miR-155 promotes chondrocyte apoptosis and catabolic activity through targeting PIK3R1-mediated PI3K/Akt pathway activation. Altogether, our study discovers novel roles and mechanisms of miR-155 in regulating chondrocyte apoptosis and catabolic activity, providing an implication for therapeutically intervening cartilage degradation and OA progression.

Can metabolic profiling provide a new description of osteoarthritis and enable a personalised medicine approach?

Osteoarthritis (OA) is a multifactorial disease contributing to significant disability and economic burden in Western populations. The aetiology of OA remains poorly understood, but is thought to involve genetic, mechanical and environmental factors. Currently, the diagnosis of OA relies predominantly on clinical assessment and plain radiographic changes long after the disease has been initiated. Recent advances suggest that there are changes in joint fluid metabolites that are associated with OA development. If this is the case, biochemical and metabolic biomarkers of OA could help determine prognosis, monitor disease progression and identify potential therapeutic targets. Moreover, for focussed management and personalised medicine, novel biomarkers could sub-stratify patients into OA phenotypes, differentiating metabolic OA from post-traumatic, age-related and genetic OA. To date, OA biomarkers have concentrated on cytokine action and protein signalling with some progress. However, these remain to be adopted into routine clinical practice. In this review, we outline the emerging metabolic links to OA pathogenesis and how an elucidation of the metabolic changes in this condition may provide future, more descriptive biomarkers to differentiate OA subtypes.

Phosphoinositide-3-kinase regulatory subunit 1 gene polymorphisms are associated with schizophrenia and bipolar disorder in the Han Chinese population

Schizophrenia (SCZ) and bipolar disorder (BD) are severe psychiatric disorders that share many genetic risk factors. This study aimed to investigate the association of phosphoinositide-3-kinase regulatory subunit1 (PIK3R1) gene rs3756668 and rs3730089 polymorphisms with SCZ and BD risks and determine the expression levels of PIK3R1. A total of 548 SCZ cases, 512 BD cases, and 598 healthy controls were included in this study. Single nucleotide polymorphisms (SNPs) were genotyped using the Sequenom MassARRAY platform, and quantitative reverse transcription polymerase chain reaction was conducted to examine the mRNA expression of PIK3R1. The genotypic distribution of rs3756668 in the BD group was significantly different from that in the healthy controls (P = 0.038). After adjustment for gender and age was made, rs3730089 was significantly associated with the risk of SCZ [AA/(AG + GG): OR = 2.25, P_adj = 0.040; AA/GG: OR = 2.27, P_adj = 0.038]. The SNP rs3756668 was associated with the susceptibility of BD (AA+GG/AG: OR = 0.73, P = 0.011) and the association remained after adjusting for gender and age. The mRNA level of PIK3R1 was significantly upregulated in patients with BD compared with that in the control group (P < 0.001). In terms of the diagnostic value of PIK3R1 for BD, the receiver operating characteristic curve analysis showed an area under the curve of 0.809 with 74.0% sensitivity and 73.9% specificity. PIK3R1 may be the shared susceptibility gene of SCZ and BD and may be a potential diagnostic biomarker for BD.

Overexpression of MIG-6 in the cartilage induces an osteoarthritis-like phenotype in mice

Background

Osteoarthritis (OA) is the most common form of arthritis and characterized by degeneration of the articular cartilage. Mitogen-inducible gene 6 (Mig-6) has been identified as a negative regulator of the epidermal growth factor receptor (EGFR). Cartilage-specific Mig-6 knockout (KO) mice display increased EGFR signaling, an anabolic buildup of the articular cartilage, and formation of chondro-osseous nodules. Since our understanding of the EGFR/Mig-6 network in the cartilage remains incomplete, we characterized mice with cartilage-specific overexpression of Mig-6 in this study.

Methods

Utilizing knee joints from cartilage-specific Mig-6-overexpressing (Mig-6^over/over) mice (at multiple time points), we evaluated the articular cartilage using histology, immunohistochemical staining, and semi-quantitative histopathological scoring (OARSI) at multiple ages. MicroCT analysis was employed to examine skeletal morphometry, body composition, and bone mineral density.

Results

Our data show that cartilage-specific Mig-6 overexpression did not cause any major developmental abnormalities in the articular cartilage, although Mig-6^over/over mice have slightly shorter long bones compared to the control group. Moreover, there was no significant difference in bone mineral density and body composition in any of the groups. However, our results indicate that Mig-6^over/over male mice show accelerated cartilage degeneration at 12 and 18 months of age. Immunohistochemistry for SOX9 demonstrated that the number of positively stained cells in Mig-6^over/over mice was decreased relative to controls. Immunostaining for MMP13 appeared increased in areas of cartilage degeneration in Mig-6^over/over mice. Moreover, staining for phospho-EGFR (Tyr-1173) and lubricin (PRG4) was decreased in the articular cartilage of Mig-6^over/over mice.

Conclusion

Overexpression of Mig-6 in the articular cartilage causes no major developmental phenotype; however, these mice develop earlier OA during aging. These data demonstrate that Mig-6/EGFR pathways are critical for joint homeostasis and might present a promising therapeutic target for OA.

Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5-UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Connections between 3' end processing and DNA damage response: Ten years later

Ten years ago we reviewed how the cellular DNA damage response (DDR) is controlled by changes in the functional and structural properties of nuclear proteins, resulting in a timely coordinated control of gene expression that allows DNA repair. Expression of genes that play a role in DDR is regulated not only at transcriptional level during mRNA biosynthesis but also by changing steady-state levels due to turnover of the transcripts. The 3' end processing machinery, which is important in the regulation of mRNA stability, is involved in these gene-specific responses to DNA damage. Here, we review the latest mechanistic connections described between 3' end processing and DDR, with a special emphasis on alternative polyadenylation, microRNA and RNA binding proteins-mediated deadenylation, and discuss the implications of deregulation of these steps in DDR and human disease. This article is categorized under: RNA Processing > 3' End Processing RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions RNA in Disease and Development > RNA in Disease.

Unpicking observational relationships between hip shape and osteoarthritis: hype or hope?

PURPOSE OF REVIEW

To review recent findings concerning the observational relationship between hip shape and hip osteoarthritis (HOA) and their shared genetic influences, and the potential for clinical application.

RECENT FINDINGS

Recent observational studies have strengthened the evidence that specific shape deformities, such as cam and acetabular dysplasia, are related to HOA. Statistical shape modelling has emerged as a method to measure hip shape holistically, with the added advantage that this can be applied to dual X-ray absorptiometry scan images. This has led to several additional aspects of hip shape variation being identified, such as a wider femoral neck and larger lesser trochanter, in association with HOA. Furthermore, this method has formed the basis of genetic studies identifying novel genetic influences on hip shape, several of which are shared with known genetic risk factors for HOA.

SUMMARY

Shared genetic influences of hip shape and HOA raise the possibility that hip shape plays a casual role in the development of HOA, justifying preventive approaches aiming to combat these adverse consequences.

Genome Engineering for Osteoarthritis: From Designer Cells to Disease-Modifying Drugs

BACKGROUND

Osteoarthritis (OA) is a highly prevalent degenerative joint disease involving joint cartilage and its surrounding tissues. OA is the leading cause of pain and disability worldwide. At present, there are no disease-modifying OA drugs, and the primary therapies include exercise and nonsteroidal anti-inflammatory drugs until total joint replacement at the end-stage of the disease.

METHODS

In this review, we summarized the current state of knowledge in genetic and epigenetic associations and risk factors for OA and their potential diagnostic and therapeutic applications.

RESULTS

Genome-wide association studies and analysis of epigenetic modifications (such as miRNA expression, DNA methylation and histone modifications) conducted across various populations support the notion that there is a genetic basis for certain subsets of OA pathogenesis.

CONCLUSION

With recent advances in the development of genome editing technologies such as the CRISPR-Cas9 system, these genetic and epigenetic alternations in OA can be used as platforms from which potential biomarkers for the diagnosis, prognosis, drug response, and development of potential personalized therapeutic targets for OA can be approached. Furthermore, genome editing has allowed the development of "designer" cells, whereby the receptors, gene regulatory networks, or transgenes can be modified as a basis for new cell-based therapies.

Epigenetics as a New Frontier in Orthopedic Regenerative Medicine and Oncology

Skeletal regenerative medicine aims to repair or regenerate skeletal tissues using pharmacotherapies, cell-based treatments, and/or surgical interventions. The field is guided by biological principles active during development, wound healing, aging, and carcinogenesis. Skeletal development and tissue maintenance in adults represent highly intricate biological processes that require continuous adjustments in the expression of cell type-specific genes that generate, remodel, and repair the skeletal extracellular matrix. Errors in these processes can facilitate musculoskeletal disease including cancers or injury. The fundamental molecular mechanisms by which cell type-specific patterns in gene expression are established and retained during successive mitotic divisions require epigenetic control, which we review here. We focus on epigenetic regulatory proteins that control the mammalian epigenome at the level of chromatin with emphasis on proteins that are amenable to drug intervention to mitigate skeletal tissue degeneration (e.g., osteoarthritis and osteoporosis). We highlight recent findings on a number of druggable epigenetic regulators, including DNA methyltransferases (e.g., DNMT1, DNMT3A, and DNMT3B) and hydroxylases (e.g., TET1, TET2, and TET3), histone methyltransferases (e.g., EZH1, EZH2, and DOT1L) as well as histone deacetylases (e.g., HDAC3, HDAC4, and HDAC7) and histone acetyl readers (e.g., BRD4) in relation to the development of bone or cartilage regenerative drug therapies. We also review how histone mutations lead to epigenomic catastrophe and cause musculoskeletal tumors. The combined body of molecular and genetic studies focusing on epigenetic regulators indicates that these proteins are critical for normal skeletogenesis and viable candidate drug targets for short-term local pharmacological strategies to mitigate musculoskeletal tissue degeneration. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1465-1474, 2019.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci

OBJECTIVE

To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization.

METHODS

We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses.

RESULTS

We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10^-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage.

CONCLUSION

PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Functional testing of thousands of osteoarthritis-associated variants for regulatory activity

To date, genome-wide association studies have implicated at least 35 loci in osteoarthritis but, due to linkage disequilibrium, the specific variants underlying these associations and the mechanisms by which they contribute to disease risk have yet to be pinpointed. Here, we functionally test 1,605 single nucleotide variants associated with osteoarthritis for regulatory activity using a massively parallel reporter assay. We identify six single nucleotide polymorphisms (SNPs) with differential regulatory activity between the major and minor alleles. We show that the most significant SNP, rs4730222, exhibits differential nuclear protein binding in electrophoretic mobility shift assays and drives increased expression of an alternative isoform of HBP1 in a heterozygote chondrosarcoma cell line, in a CRISPR-edited osteosarcoma cell line, and in chondrocytes derived from osteoarthritis patients. This study provides a framework for prioritization of GWAS variants and highlights a role of HBP1 and Wnt signaling in osteoarthritis pathogenesis.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Wnt and RUNX2 mediate cartilage breakdown by osteoarthritis synovial fibroblast-derived ADAMTS-7 and -12

Failure of therapeutic approaches for the treatment of osteoarthritis (OA) based on the inhibition of metalloproteinases, might be because of their constitutive expression in homeostasis, together with their network complexity. The knowledge of this network would contribute to selective target pathological conditions. In this sense, blockade of mediators produced by neighbouring joint cells, such as synovial fibroblasts (SF), would prevent cartilage damage. Thus, we studied the contribution of ADAMTS‐7 and ‐12 from SF to cartilage oligomeric matrix protein (COMP) degradation, and the signalling pathways involved in their expression. We report for the first time in SF, the involvement of ERK‐Runx2 axis and Wnt/β‐catenin signalling in ADAMTS‐12 and ADAMTS‐7 expressions, respectively, with the subsequent consequences in COMP degradation from cartilage extracellular matrix. After stimulation with IL‐1β or fibronectin fragments, we showed that ERK inhibition decreased Runx2 activation and ADAMTS‐12 expression in OA‐SF, also reducing Fn‐fs‐induced COMP degradation. Blockage of Wnt signalling by DKK1 reduced ADAMTS‐7 and COMP degradation in OA‐SF as well. In addition, Wnt7B expression was induced by IL‐1β and by itself, also increasing ADAMTS‐7. Our results could contribute to the development of disease‐modifying OA drugs targeting ADAMTS‐7 and ‐12 for the prevention of extracellular matrix components degradation like COMP.

Increased susceptibility to develop spontaneous and post-traumatic osteoarthritis in Dot1l-deficient mice

OBJECTIVE

We earlier identified that the histone methyltransferase Disruptor of telomeric silencing 1-like (DOT1L) is as a master protector of cartilage health via limiting excessive activation of the Wnt pathway. However, cartilage-specific homozygous Dot1l knockout mice exhibited a severe growth phenotype and perinatal death, which hampered their use in induced or ageing models of osteoarthritis (OA). The aim of this study was to generate and examine haploinsufficient and inducible conditional Dot1l-deficient mouse models to evaluate the importance of DOT1L during post-traumatic or ageing-associated OA onset and progression.

METHOD

We used cartilage-specific heterozygous and postnatal tamoxifen-inducible Dot1l knockout mice and performed destabilization of the medial meniscus (DMM) and ageing as OA models. Mice were examined histologically using X-rays and micro-computed tomography (μCT), and cartilage damage and osteophyte formation were assessed based on OARSI guidelines. Immunohistochemistry of DOT1L, H3K79me2, TCF1 and COLX was performed.

RESULTS

Both Dot1l-deficient strains exhibit a phenotype characterized by joint remodeling with extensive osteophyte formation and ectopic ossification upon ageing, indicating accelerated development of spontaneous osteoarthritis. In the DMM-induced OA mouse model, absence of Dot1l resulted in increased cartilage damage. Wnt signalling hyper-activation and ectopic chondrocyte hypertrophy were observed in the articular cartilage of both Dot1l-deficient mice.

CONCLUSIONS

This study demonstrated the functional relevance of DOT1L in vivo during the development of OA using genetically modified mice. Thus, maintaining or enhancing DOT1L activity during ageing or after trauma might prevent OA onset and progression.

Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage

OBJECTIVE

Multiple single-nucleotide polymorphisms (SNPs) conferring susceptibility to osteoarthritis (OA) mark imbalanced expression of positional genes in articular cartilage, reflected by unequally expressed alleles among heterozygotes (allelic imbalance [AI]). We undertook this study to explore the articular cartilage transcriptome from OA patients for AI events to identify putative disease-driving genetic variation.

METHODS

AI was assessed in 42 preserved and 5 lesioned OA cartilage samples (from the Research Arthritis and Articular Cartilage study) for which RNA sequencing data were available. The count fraction of the alternative alleles among the alternative and reference alleles together (φ) was determined for heterozygous individuals. A meta-analysis was performed to generate a meta-φ and P value for each SNP with a false discovery rate (FDR) correction for multiple comparisons. To further validate AI events, we explored them as a function of multiple additional OA features.

RESULTS

We observed a total of 2,070 SNPs that consistently marked AI of 1,031 unique genes in articular cartilage. Of these genes, 49 were found to be significantly differentially expressed (fold change <0.5 or >2, FDR <0.05) between preserved and paired lesioned cartilage, and 18 had previously been reported to confer susceptibility to OA and/or related phenotypes. Moreover, we identified notable highly significant AI SNPs in the CRLF1, WWP2, and RPS3 genes that were related to multiple OA features.

CONCLUSION

We present a framework and resulting data set for researchers in the OA research field to probe for disease-relevant genetic variation that affects gene expression in pivotal disease-affected tissue. This likely includes putative novel compelling OA risk genes such as CRLF1, WWP2, and RPS3.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

EGFR Signaling: Friend or Foe for Cartilage?

Recent studies using genetically modified mice, pharmacological approaches, and human samples have highlighted an important role for the epidermal growth factor receptor (EGFR), selected ligands, and downstream components in endochondral bone formation and joint homeostasis. Although most data demonstrate an important function of this pathway in endochondral ossification and articular cartilage growth, conflicting results on its role in osteoarthritis have been reported. In some contexts, inactivation of EGFR signaling has been shown to protect joints from surgically induced osteoarthritis, whereas in others, similar manipulations worsened joint pathology. The current review summarizes recent studies of cartilage EGFR signaling in long bone development and diseases, provides potential explanations for the reported discrepancies, and suggests directions for future work to clarify the potential of this pathway as target for osteoarthritis treatment. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

Multitrait genome association analysis identifies new susceptibility genes for human anthropometric variation in the GCAT cohort

Background

Heritability estimates have revealed an important contribution of SNP variants for most common traits; however, SNP analysis by single-trait genome-wide association studies (GWAS) has failed to uncover their impact. In this study, we applied a multitrait GWAS approach to discover additional factor of the missing heritability of human anthropometric variation.

Methods

We analysed 205 traits, including diseases identified at baseline in the GCAT cohort (Genomes For Life- Cohort study of the Genomes of Catalonia) (n=4988), a Mediterranean adult population-based cohort study from the south of Europe. We estimated SNP heritability contribution and single-trait GWAS for all traits from 15 million SNP variants. Then, we applied a multitrait-related approach to study genome-wide association to anthropometric measures in a two-stage meta-analysis with the UK Biobank cohort (n=336 107).

Results

Heritability estimates (eg, skin colour, alcohol consumption, smoking habit, body mass index, educational level or height) revealed an important contribution of SNP variants, ranging from 18% to 77%. Single-trait analysis identified 1785 SNPs with genome-wide significance threshold. From these, several previously reported single-trait hits were confirmed in our sample with LINC01432 (p=1.9×10⁻⁹) variants associated with male baldness, LDLR variants with hyperlipidaemia (ICD-9:272) (p=9.4×10⁻¹⁰) and variants in IRF4 (p=2.8×10⁻⁵⁷), SLC45A2 (p=2.2×10⁻¹³⁰), HERC2 (p=2.8×10⁻¹⁷⁶), OCA2 (p=2.4×10⁻¹²¹) and MC1R (p=7.7×10⁻²²) associated with hair, eye and skin colour, freckling, tanning capacity and sun burning sensitivity and the Fitzpatrick phototype score, all highly correlated cross-phenotypes. Multitrait meta-analysis of anthropometric variation validated 27 loci in a two-stage meta-analysis with a large British ancestry cohort, six of which are newly reported here (p value threshold <5×10⁻⁹) at ZRANB2-AS2, PIK3R1, EPHA7, MAD1L1, CACUL1 and MAP3K9.

Conclusion

Considering multiple-related genetic phenotypes improve associated genome signal detection. These results indicate the potential value of data-driven multivariate phenotyping for genetic studies in large population-based cohorts to contribute to knowledge of complex traits.

What did we learn from 'omics' studies in osteoarthritis

PURPOSE OF REVIEW

'Omics' technologies developed for the massive analysis of the major biologically relevant molecules (genes, proteins, metabolites) have been applied to the study of osteoarthritis (OA) for more than a decade.

RECENT FINDINGS

'Omics' studies have undoubtedly contributed to increase the knowledge on pathogenic processes related with OA and have provided hundreds to thousands of molecules that might have a putative biomarker utility for this disease.

SUMMARY

This review describes the most recent 'omics' studies in OA research, their conclusions, and discuss those remaining challenges. Still many validation studies must be performed in large and well-characterized cohorts for the translation of the findings from 'omics' strategies to clinical applications. The development of tools for the intelligent integration of 'omics' data with clinical and imaging information is also mandatory to take full profit of the work that has been already performed.

Development of a prediction model for future risk of radiographic hip osteoarthritis

OBJECTIVE

To develop and validate a prognostic model for incident radiologic hip osteoarthritis (HOA) and determine the value of previously identified predictive factors.

DESIGN

We first validated previously reported predictive factors for HOA by performing univariate and multivariate analyses for all predictors in three large prospective cohorts (total sample size of 4548 with 653 incident cases). The prognostic model was developed in 2327 individuals followed for 10 years from the Rotterdam Study-I (RS-I) cohort. External validation of the model was tested on discrimination in two other cohorts: RS-II (n = 1435) and the Cohort Hip and Cohort Knee (CHECK) study (n = 786).

RESULTS

From the total number of 28 previously reported predictive factors, we were able to replicate 13 factors, while 15 factors were not significantly predictive in a meta-analysis of the three cohorts. The basic model including the demographic, questionnaire, and clinical examination variables (area under the receiver-operating characteristic curve (AUC) = 0.67) or genetic markers (AUC = 0.55) or urinary C-terminal cross-linked telopeptide of type II collagen (uCTX-II) levels (AUC = 0.67) alone were poor predictors of HOA in all cohorts. Imaging factors showed the highest predictive value for the development of HOA (AUC = 0.74). Addition of imaging variables to the basic model led to substantial improvement in the discriminative ability of the model (AUC = 0.78) compared with uCTX-II (AUC = 0.74) or genetic markers (AUC = 0.68). Applying external validation, similar results were observed in the RS-II and the CHECK cohort.

CONCLUSIONS

The developed prediction model included demographic, a limited number of questionnaire, and imaging risk factors seems promising for prediction of HOA.

Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis

Osteoarthritis is a common complex disease imposing a large public-health burden. Here, we performed a genome-wide association study for osteoarthritis, using data across 16.5 million variants from the UK Biobank resource. After performing replication and meta-analysis in up to 30,727 cases and 297,191 controls, we identified nine new osteoarthritis loci, in all of which the most likely causal variant was noncoding. For three loci, we detected association with biologically relevant radiographic endophenotypes, and in five signals we identified genes that were differentially expressed in degraded compared with intact articular cartilage from patients with osteoarthritis. We established causal effects on osteoarthritis for higher body mass index but not for triglyceride levels or genetic predisposition to type 2 diabetes.