Severe osteoarthritis of the hand associates with common variants within the ALDH1A2 gene and with rare variants at 1p31

The Genetics and Functional Genomics of Osteoarthritis

Osteoarthritis is the most prevalent whole-joint degenerative disorder, and is characterized by the degradation of articular cartilage and the underlying bone structures. Almost 600 million people are affected by osteoarthritis worldwide. No curative treatments are available, and management strategies focus mostly on pain relief. Here, we provide a comprehensive overview of the available human genetic and functional genomics studies for osteoarthritis to date and delineate how these studies have helped shed light on disease etiopathology. We highlight genetic discoveries from genome-wide association studies and provide a detailed overview of molecular-level investigations in osteoarthritis tissues, including methylation-, transcriptomics-, and proteomics-level analyses. We review how functional genomics data from different molecular levels have helped to prioritize effector genes that can be used as drug targets or drug-repurposing opportunities. Finally, we discuss future directions with the potential to drive a step change in osteoarthritis research.

Hand osteoarthritis: A fresh look

Hand osteoarthritis (OA) has been the subject of numerous publications in recent years, particularly in the fields of imaging and therapeutics. The imaging studies revealed a good correlation between the presence of synovitis and/or subchondral edema and arthritic joint pain. Several randomized controlled trials (RCTs) have assessed the efficacy of biologics and conventional DMARDs in patients with symptomatic hand OA. No less than six RCTs have evaluated the symptomatic and, in some cases, structural efficacy of anti-IL-1, anti-TNF or anti-IL-6 drugs. Overall, the results of these trials were disappointing - none of them demonstrated superiority over placebo. There were also two negative trials with hydroxychloroquine. In the end, the only trial that was positive evaluated 10mg oral prednisone versus placebo for 6 weeks in patients with flares of hand OA and synovitis visible on ultrasound. While that trial confirms the role of inflammation in hand OA, it should obviously not encourage the long-term use of corticosteroids as a symptomatic treatment.

A functional variant of ALDH1A2 is associated with hand osteoarthritis in the Chinese population

Genome-wide association study identified common variants within the ALDH1A2 gene as the susceptible loci of hand osteoarthritis (HOA) in UK and Iceland populations. Located in chromosome 15, ALDH1A2 encodes aldehyde dehydrogenase family 1 member A2, which is an enzyme that catalyses the synthesis of retinoic acid from retinaldehyde. Our purposes were to replicate the association of functional variant in ALDH1A2 with the development of HOA in the Chinese population. Variant rs12915901 of ALDH1A2 was genotyped in 872 HOA patients and 1223 healthy controls. Subchondral bone samples were collected from 40 patients who had undergone a trapeziectomy, and the tissue expression of ALDH1A2 was analysed. The chi-square analysis was used to compare the frequency of genotype and risk allele between the HOA cases and controls. The Student t test was used to compare the mRNA expression of ALDH1A2 between patients with genotype AA/AG and those with genotype GG. The frequency of genotype AA was significantly higher in HOA patients than in the controls (7.6% vs. 5.1%, p = .01). The frequency of allele A was significantly higher in the patients than in the controls (28.9% vs. 24.6%, p = .005). The mRNA expression of ALDH1A2 was 1.31-folds higher in patients with genotype GG than in the patients with genotype AA/AG (0.000617 ± 0.000231 vs. 0.000471 ± 0.000198, p = .04). Variant rs12915901 of ALDH1A2 contributed to the susceptibility of HOA in the Chinese population. Allele A of rs12915901 can add to the risk of HOA possibly via down-regulation of ALDH1A2 expression.

Aldehyde dehydrogenase 1 family: A potential molecule target for diseases

Aldehyde dehydrogenase 1 (ALDH1), a crucial aldehyde metabolizing enzyme, has six family members. The ALDH1 family is expressed in various tissues, with a significant presence in the liver. It plays a momentous role in several pathophysiological processes, including aldehyde detoxification, oxidative stress, and lipid peroxidation. Acetaldehyde detoxification is the fundamental function of the ALDH1 family in participating in vital pathological mechanisms. The ALDH1 family can catalyze retinal to retinoic acid (RA) that is a hormone-signaling molecule and plays a vital role in the development and adult tissues. Furthermore, there is a need for further and broader research on the role of the ALDH1 family as a signaling molecule. The ALDH1 family is widely recognized as a cancer stem cell (CSC) marker and plays a significant role in the proliferation, invasion, metastasis, prognosis, and drug resistance of cancer. The ALDH1 family also participates in other human diseases, such as neurodegenerative diseases, osteoarthritis, diabetes, and atherosclerosis. It can inhibit disease progression by inhibiting/promoting the expression/activity of the ALDH1 family. In this review, we comprehensively analyze the tissue distribution, and functions of the ALDH1 family. Additionally, we review the involvement of the ALDH1 family in diseases, focusing on the underlying pathological mechanisms and briefly talk about the current status and development of ALDH1 family inhibitors. The ALDH1 family presents new possibilities for treating diseases, with both its upstream and downstream pathways serving as promising targets for therapeutic intervention. This offers fresh perspectives for drug development in the field of disease research.

Epigenetic mechanisms of osteoarthritis risk in human skeletal development

The epigenome, including the methylation of cytosine bases at CG dinucleotides, is intrinsically linked to transcriptional regulation. The tight regulation of gene expression during skeletal development is essential, with ~1/500 individuals born with skeletal abnormalities. Furthermore, increasing evidence is emerging to link age-associated complex genetic musculoskeletal diseases, including osteoarthritis (OA), to developmental factors including joint shape. Multiple studies have shown a functional role for DNA methylation in the genetic mechanisms of OA risk using articular cartilage samples taken from aged patients. Despite this, our knowledge of temporal changes to the methylome during human cartilage development has been limited. We quantified DNA methylation at ~700,000 individual CpGs across the epigenome of developing human articular cartilage in 72 samples ranging from 7-21 post-conception weeks, a time period that includes cavitation of the developing knee joint. We identified significant changes in 8% of all CpGs, and >9400 developmental differentially methylated regions (dDMRs). The largest hypermethylated dDMRs mapped to transcriptional regulators of early skeletal patterning including MEIS1 and IRX1. Conversely, the largest hypomethylated dDMRs mapped to genes encoding extracellular matrix proteins including SPON2 and TNXB and were enriched in chondrocyte enhancers. Significant correlations were identified between the expression of these genes and methylation within the hypomethylated dDMRs. We further identified 811 CpGs at which significant dimorphism was present between the male and female samples, with the majority (68%) being hypermethylated in female samples. Following imputation, we captured the genotype of these samples at >5 million variants and performed epigenome-wide methylation quantitative trait locus (mQTL) analysis. Colocalization analysis identified 26 loci at which genetic variants exhibited shared impacts upon methylation and OA genetic risk. This included loci which have been previously reported to harbour OA-mQTLs (including GDF5 and ALDH1A2), yet the majority (73%) were novel (including those mapping to CHST3, FGF1 and TEAD1). To our knowledge, this is the first extensive study of DNA methylation across human articular cartilage development. We identify considerable methylomic plasticity within the development of knee cartilage and report active epigenomic mediators of OA risk operating in prenatal joint tissues.

Genetics may affect the risk of undergoing surgery for rhizarthrosis

Osteoarthritis is a prevalent and severe disease. Involvement of the trapeziometacarpal joint is common and can lead to both pain and disability. Genetics are known to affect the risk of osteoarthritis, but it remains unclear how genetics affect disease trajectories. In this study, we investigated whether the genetic associations of trapeziometacarpal osteoarthritis (rhizarthrosis) vary with the need for surgical treatment. The study was conducted as a case-control genome-wide association study using individuals from the Copenhagen Hospital Biobank pain and degenerative musculoskeletal disease study and the Danish Blood Donor Study (N = 208,342). We identified patients diagnosed with rhizarthrosis and grouped them by treatment status, resulting in two case groups: surgical (N = 1083) and nonsurgical (N = 1888). The case groups were tested against osteoarthritis-free controls in two genome-wide association studies. We then compared variants suggestive of association (p < 10-6) in either of these analyses directly between the treatment groups (surgical vs. nonsurgical rhizarthrosis). We identified 10 variants suggestive of association with either surgical (seven variants) or nonsurgical (three variants) rhizarthrosis. None of the variants reached nominal significance in the opposite treatment group (p ≥ 0.14), and all 10 variants were significantly different between the treatment groups at a false discovery rate of 5%. These results suggest possible differences in the genetic associations of rhizarthrosis depending on surgical treatment. Clinical significance: Uncovering genetic differences between clinically distinct patient groups can reveal biological determinants of disease trajectories.

The Genetic Markers of Knee Osteoarthritis in Women from Russia

Osteoarthritis is a chronic progressive joint disease that clinically debuts at the stage of pronounced morphologic changes, which makes treatment difficult. In this regard, an important task is the study of genetic markers of the disease, which have not been definitively established, due to the clinical and ethnic heterogeneity of the studied populations. To find the genetic markers for the development of knee osteoarthritis (OA) in women from the Volga-Ural region of Russia, we conducted research in two stages using different genotyping methods, such as the restriction fragment length polymorphism (RFLP) measurement, TaqMan technology and competitive allele-specific PCR-KASPTM. In the first stage, we studied polymorphic variants of candidate genes (ACAN, ADAMTS5, CHST11, SOX9, COL1A1) for OA development. The association of the *27 allele of the VNTR locus of the ACAN gene was identified (OR = 1.6). In the second stage, we replicated the GWAS results (ASTN2, ALDH1A2, DVWA, CHST11, GNL3, NCOA3, FILIP/SENP1, MCF2L, GLT8D, DOT1L) for knee OA studies. The association of the *T allele of the rs7639618 locus of the DVWA gene was detected (OR = 1.54). Thus, the VNTR locus of ACAN and the rs7639618 locus of DVWA are risk factors for knee OA in women from the Volga-Ural region of Russia.

Erosive Hand Osteoarthritis: Recent Advances and Future Treatments

PURPOSE OF THE REVIEW

Erosive hand osteoarthritis (EHOA) is an aggressive form of hand osteoarthritis that leads to significant disability, and recent data suggests that it is increasing in prevalence. This review provides an update of our current understanding of epidemiology, genetic associations, biomarkers, pathogenesis, and treatment of EHOA, with particular focus on studies published within the last 5 years.

RECENT FINDINGS

New studies of EHOA have identified new genetic loci associated with disease, including variants in genes involved in inflammation and bone remodeling. Preclinical studies implicate pathways of innate immunity, including some that may be causal in the condition. Recent novel studies showed that inflammatory features identified by ultrasound and MRI are associated with development of erosive lesions over time on conventional radiography. In the future, these imaging modalities may be useful in identifying patients at risk of adverse outcomes. Promising new findings in genetics, biomarkers, and treatment targets will hopefully allow for future therapeutic options for this debilitating condition.

Estrogen and testosterone supplementation improves tendon healing and functional recovery after rotator cuff repair

Failure of healing after rotator cuff repair (RCR) is common. The purpose of the current study was to evaluate the effect of systemic estrogen or testosterone supplementation on tendon healing after RCR. Seventy-two adult male mice were utilized for all experiments. The supraspinatus tendon was transected and repaired with 6-0 Prolene suture on the left shoulder of 51 animals. Mice were segregated into three groups postoperative: (1) vehicle group (VG; n = 18), (2) estrogen group (EST; n = 17), and (3) testosterone group (TST; n = 16). An unrepaired control group (unrepaired, n = 21) did not have surgery. Utilizing these animals, histological analysis, activity testing, biomechanical testing and RNA sequencing (RNA-seq) was performed. At 8 weeks post-RCR, TST, and EST supplementation improved the overall histologic structure of the repaired enthesis site. No differences in ultimate failure loads or stiffness were detected between VG, EST, and TST groups after biomechanical testing. RCR caused a reduction in wheel activity compared to unrepaired controls and supplementation with TST restored wheel activity. RNA-seq analysis indicated that estrogen and testosterone regulated different pathways associated with enthesis healing, including a suppression of inflammatory signaling. Supplementation with sex hormones improved the structure of the repaired tendon enthesis and significantly regulated expression of diverse pathways regulating multiple biological processes. Testosterone administration following RCR restored wheel activity without having a detrimental impact on biomechanical strength. Future human studies of sex hormone supplementation after RCR are warranted as supplementation in an animal model may improve tendon enthesis healing.

Osteoarthritis year in review 2023: Biology

Great progress continues to be made in our understanding of the multiple facets of osteoarthritis (OA) biology. Here, we review the major advances in this field and progress towards therapy development over the past year, highlighting a selection of relevant published literature from a PubMed search covering the year from the end of April 2022 to the end of April 2023. The selected articles have been arranged in themes. These include 1) molecular regulation of articular cartilage and implications for OA, 2) mechanisms of subchondral bone remodelling, 3) role of synovium and inflammation, 4) role of age-related changes including cartilage matrix stiffening, cellular senescence, mitochondrial dysfunction, metabolic dysfunction, and impaired autophagy, and 5) peripheral mechanisms of OA pain. Progress in the understanding of the cellular and molecular mechanisms responsible for the multiple aspects of OA biology is unravelling novel therapeutic targets for disease modification.

Exploring Talarozole as a Novel Therapeutic Approach for Osteoarthritis: Insights From Experimental Studies

Osteoarthritis (OA) is a prevalent degenerative joint disease characterized by the slow degeneration of joint components that primarily affects the elderly. There is currently no cure for OA; thus, treatment focuses on symptom reduction. This article investigates the potential of talarozole, a retinoic acid metabolism-blocking agent (RAMBA), as a new treatment for hand OA. Talarozole showed promising results by inhibiting retinoic acid degradation and increasing its levels in the body. Six hours after destabilization of the medial meniscus, talarozole significantly reduced inflammation in mice's cartilage. The findings underscore the importance of the protein encoded by the ALDH1A2 gene in retinoic acid metabolism, shedding light on its potential implications for the management of OA. Maintaining adequate retinoic acid levels may help to reduce mechano-inflammatory gene regulation. Furthermore, RAMBAs like talarozole may emerge as disease-modifying OA therapies, promising improved symptom control and slower disease progression. In conclusion, this research provides critical genetic insights into severe hand OA and promotes talarozole as a prospective therapy option. These findings pave the door for additional research that could revolutionize OA treatment by targeting retinoic acid metabolism to reduce symptoms and slow disease progression.

Retinoic acid signaling in development and differentiation commitment and its regulatory topology

Retinoic acid (RA), the derivative of vitamin A/retinol, is a signaling molecule with important implications in health and disease. It is a well-known developmental morphogen that functions mainly through the transcriptional activity of nuclear RA receptors (RARs) and, uncommonly, through other nuclear receptors, including peroxisome proliferator-activated receptors. Intracellular RA is under spatiotemporally fine-tuned regulation by synthesis and degradation processes catalyzed by retinaldehyde dehydrogenases and P450 family enzymes, respectively. In addition to dictating the transcription architecture, RA also impinges on cell functioning through non-genomic mechanisms independent of RAR transcriptional activity. Although RA-based differentiation therapy has achieved impressive success in the treatment of hematologic malignancies, RA also has pro-tumor activity. Here, we highlight the relevance of RA signaling in cell-fate determination, neurogenesis, visual function, inflammatory responses and gametogenesis commitment. Genetic and post-translational modifications of RAR are also discussed. A better understanding of RA signaling will foster the development of precision medicine to improve the defects caused by deregulated RA signaling.

Genome-Wide Association Studies to Drug: Identifying Retinoic Acid Metabolism Blocking Agents to Suppress Mechanoflammation in Osteoarthritis

Osteoarthritis (OA) is a highly prevalent debilitating joint disease for which there are currently no licensed disease-modifying treatments. The pathogenesis of OA is complex, involving genetic, mechanical, biochemical, and environmental factors. Cartilage injury, arguably the most important driving factor in OA development, is able to activate both protective and inflammatory pathways within the tissue. Recently, >100 genetic risk variants for OA have been identified through Genome Wide Association Studies, which provide a powerful tool to validate existing putative disease pathways and discover new ones. Using such an approach, hypomorphic variants within the aldehyde dehydrogenase 1 family member A2 (ALDH1A2) gene were shown to be associated with increased risk of severe hand OA. ALDH1A2 encodes the enzyme that synthesizes all-trans retinoic acid (atRA), an intracellular signaling molecule. This review summarizes the influence of the genetic variants on expression and function of ALDH1A2 in OA cartilage, its role in the mechanical injury response of cartilage, and its potent anti-inflammatory effect after cartilage injury. In doing so it identifies atRA metabolism-blocking agents as potential treatments for suppressing mechanoflammation in OA.

Bone-cartilage crosstalk informed by aging mouse bone transcriptomics and human osteoarthritis genome-wide association studies

Subchondral bone participates in crosstalk with articular cartilage to maintain joint homeostasis, and disruption of either tissue results in overall joint degeneration. Among the subchondral bone changes observed in osteoarthritis (OA), subchondral bone plate (SBP) thickening has a time-dependent relationship with cartilage degeneration and has recently been shown to be regulated by osteocytes. Here, we evaluate the effect of age on SBP thickness and cartilage degeneration in aging mice. We find that SBP thickness significantly increases by 18-months of age, corresponding temporally with increased cartilage degeneration. To identify factors in subchondral bone that may participate in bone cartilage crosstalk or OA, we leveraged mouse transcriptomic data from one joint tissue compartment - osteocyte-enriched bone - to search for enrichment with human OA in UK Biobank and Arthritis Research UK Osteoarthritis Genetics (arcOGEN) GWAS using the mouse2human (M2H, www.mouse2human.org) strategy. Genes differentially expressed in aging mouse bone are significantly enriched for human OA, showing joint site-specific (knee vs. hip) relationships, exhibit temporal associations with age, and unique gene clusters are implicated in each type of OA. Application of M2H identifies genes with known and unknown functions in osteocytes and OA development that are clinically associated with human OA. Altogether, this work prioritizes genes with a potential role in bone/cartilage crosstalk for further mechanistic study based on their association with human OA in GWAS.

RNA-Seq Reveals the mRNAs, miRNAs, and lncRNAs Expression Profile of Knee Joint Synovial Tissue in Osteoarthritis Patients

Osteoarthritis (OA) is a chronic disease common in the elderly population and imposes significant health and economic burden. Total joint replacement is the only currently available treatment but does not prevent cartilage degeneration. The molecular mechanism of OA, especially the role of inflammation in disease progression, is incompletely understood. We collected knee joint synovial tissue samples of eight OA patients and two patients with popliteal cysts (controls), measured the expression levels of lncRNAs, miRNAs, and mRNAs in these tissues by RNA-seq, and identified differentially expressed genes (DEGs) and key pathways. In the OA group, 343 mRNAs, 270 lncRNAs, and 247 miRNAs were significantly upregulated, and 232 mRNAs, 109 lncRNAs, and 157 miRNAs were significantly downregulated. mRNAs potentially targeted by lncRNAs were predicted. Nineteen overlapped miRNAs were screened based on our sample data and GSE 143514 data. Pathway enrichment and functional annotation analyses showed that the inflammation-related transcripts CHST11, ALDH1A2, TREM1, IL-1β, IL-8, CCL5, LIF, miR-146a-5p, miR-335-5p, lncRNA GAS5, LINC02288, and LOC101928134 were differentially expressed. In this study, inflammation-related DEGs and non-coding RNAs were identified in synovial samples, suggesting that competing endogenous RNAs have a role in OA. TREM1, LIF, miR146-5a, and GAS5 were identified to be OA-related genes and potential regulatory pathways. This research helps elucidate the pathogenesis of OA and identify novel therapeutic targets for this disorder.

A Functional Polymorphism Downstream of Vitamin A Regulator Gene CYP26B1 Is Associated with Hand Osteoarthritis

While genetic analyses have revealed ~100 risk loci associated with osteoarthritis (OA), only eight have been linked to hand OA. Besides, these studies were performed in predominantly European and Caucasian ancestries. Here, we conducted a genome-wide association study in the Han Chinese population to identify genetic variations associated with the disease. We recruited a total of 1136 individuals (n = 420 hand OA-affected; n = 716 unaffected control subjects) of Han Chinese ancestry. We carried out genotyping using Axiom Asia Precisi on Medicine Research Array, and we employed the RegulomeDB database and RoadMap DNase I Hypersensitivity Sites annotations to further narrow down our potential candidate variants. Genetic variants identified were tested in the Geisinger's hand OA cohort selected from the Geisinger MyCode community health initiative (MyCode^®). We also performed a luciferase reporter assay to confirm the potential impact of top candidate single-nucleotide polymorphisms (SNPs) on hand OA. We identified six associated SNPs (p-value = 6.76 × 10^-7-7.31 × 10^-6) clustered at 2p13.2 downstream of the CYP26B1 gene. The strongest association signal identified was rs883313 (p-value = 6.76 × 10^-7, odds ratio (OR) = 1.76), followed by rs12713768 (p-value = 1.36 × 10^-6, OR = 1.74), near or within the enhancer region closest to the CYP26B1 gene. Our findings showed that the major risk-conferring CC haplotype of SNPs rs12713768 and rs10208040 [strong linkage disequilibrium (LD); D' = 1, r² = 0.651] drives 18.9% of enhancer expression activity. Our findings highlight that the SNP rs12713768 is associated with susceptibility to and severity of hand OA in the Han Chinese population and that the suggested retinoic acid signaling pathway may play an important role in its pathogenesis.

Variants in ALDH1A2 reveal an anti-inflammatory role for retinoic acid and a new class of disease-modifying drugs in osteoarthritis

More than 40% of individuals will develop osteoarthritis (OA) during their lifetime, yet there are currently no licensed disease-modifying treatments for this disabling condition. Common polymorphic variants in ALDH1A2, which encodes the key enzyme for synthesis of all-trans retinoic acid (atRA), are associated with severe hand OA. Here, we sought to elucidate the biological significance of this association. We first confirmed that ALDH1A2 risk variants were associated with hand OA in the U.K. Biobank. Articular cartilage was acquired from 33 individuals with hand OA at the time of routine hand OA surgery. After stratification by genotype, RNA sequencing was performed. A reciprocal relationship between ALDH1A2 mRNA and inflammatory genes was observed. Articular cartilage injury up-regulated similar inflammatory genes by a process that we have previously termed mechanoflammation, which we believe is a primary driver of OA. Cartilage injury was also associated with a concomitant drop in atRA-inducible genes, which were used as a surrogate measure of cellular atRA concentration. Both responses to injury were reversed using talarozole, a retinoic acid metabolism blocking agent (RAMBA). Suppression of mechanoflammation by talarozole was mediated by a peroxisome proliferator-activated receptor gamma (PPARγ)-dependent mechanism. Talarozole was able to suppress mechano-inflammatory genes in articular cartilage in vivo 6 hours after mouse knee joint destabilization and reduced cartilage degradation and osteophyte formation after 26 days. These data show that boosting atRA suppresses mechanoflammation in the articular cartilage in vitro and in vivo and identifies RAMBAs as potential disease-modifying drugs for OA.

Novel genetic loci associated with osteoarthritis in multi-ancestry analyses in the Million Veteran Program and UK Biobank

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain.

Familial Clustering and Genetic Analysis of Severe Thumb Carpometacarpal Joint Osteoarthritis in a Large Statewide Cohort

PURPOSE

Our goals were to identify individuals who required surgery for thumb carpometacarpal (CMC) joint osteoarthritis (OA), determine if CMC joint OA clusters in families, define the magnitude of familial risk of CMC joint OA, identify risk factors associated with CMC joint OA, and identify rare genetic variants that segregate with familial CMC joint OA.

METHODS

We searched the Utah Population Database to identify a cohort of CMC joint OA patients who required surgery. Affected individuals were mapped to pedigrees to identify high-risk families with excess clustering of CMC joint OA. Cox regression models were used to calculate familial risk of CMC joint OA in related individuals. Risk factors were evaluated using logistic regression models. Whole exome sequencing was used to identify rare coding variants associated with familial CMC joint OA.

RESULTS

We identified 550 pedigrees with excess clustering of severe CMC joint OA. The relative risk of CMC joint OA requiring surgical treatment was elevated significantly in first- and third-degree relatives of affected individuals, and significant associations with advanced age, female sex, obesity, and tobacco use were observed. We discovered candidate genes that dominantly segregate with severe CMC joint OA in 4 independent families, including a rare variant in Chondroitin Sulfate Synthase 3 (CHSY3).

CONCLUSIONS

Familial clustering of severe CMC joint OA was observed in a statewide population. Our data indicate that genetic and environmental factors contribute to the disease process, further highlighting the multifactorial nature of the disease. Genomic analyses suggest distinct biological processes are involved in CMC joint OA pathogenesis.

CLINICAL RELEVANCE

Awareness of associated comorbidities may guide the diagnosis of CMC joint OA in at-risk populations and help identify individuals who may not do well with nonoperative treatment. Further pursuit of the genes associated with severe CMC joint OA may lead to assays for detection of early stages of disease and have therapeutic potential.

A neutrophil cell membrane-biomimetic nanoplatform based on L-arginine nanoparticles for early osteoarthritis diagnosis and nitric oxide therapy

Osteoarthritis (OA) is a common debilitating disease affecting articular joints for which no effective disease-modifying early diagnosis or medical therapy tools are currently available. The inefficient delivery of drugs into inflamed chondrocytes has restricted the development of anti-OA medication. Evidence has shown that inflammatory neutrophils possess the property of targeting inflammation via inflammatory tissue recruiting. Herein, we report neutrophil-cell-membrane-based biomimetic nanoparticles (NM-LANPs@Ru) as an OA theranostic nanoplatform; they act as a NO delivery system, coating neutrophil cell membrane onto the surface of self-assembled PEGylated L-arginine nanoparticles (LANPs) to act as a NO donor and loading a Ru complex to act as a ROS inducer. NM-LANPs@Ru demonstrated the specific targeting of inflamed OA with low toxicity, good NO release, and excellent fluorescence/photoacoustic (FL/PA) imaging properties. We showed that NM-LANPs@Ru exhibited enhanced cellular association in inflamed chondrocyte cells (C28/I2), much higher than NO release from ROS oxidized LA, and it improved the inhibition of the apoptosis of inflamed C28/I2 cells compared with control treatments. In vivo studies demonstrated that NM-LANPs@Ru effectively targeted inflamed OA, based on real-time dual-modal FL/PA imaging, eventually exhibiting its excellent anti-inflammatory activity. Our study may provide a new approach for the early diagnosis and treatment of osteoarthritis using a neutrophil-cell-membrane-based biomimetic nanoplatform for NO or drug delivery.

Identification of candidate enhancers controlling the transcriptome during the formation of interphalangeal joints

The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

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.

An eQTL variant of ALDH1A2 is associated with Kashin-Beck disease in Chinese population

INTRODUCTION

The aims of the study were to investigate the relationship between aldehyde dehydrogenase 1 family member A2 (ALDH1A2) and Kashin-Beck disease (KBD), explore the effects of the rs3204689 polymorphism and methylation status on the expression levels of ALDH1A2, and further clarify the pathogenesis of KBD.

MATERIALS AND METHODS

The genotype of ALDH1A2 rs3204689 was detected by PCR-RFLP in 103 KBD patients and 109 healthy controls in the whole blood. The mRNA level of ALDH1A2 was measured by qRT-PCR, and the protein expression was detected using IHC staining and Western blotting. The MSP-PCR was used to identify the ALDH1A2 methylation level.

RESULTS

There were significant differences in G/G, G/C, and C/C frequencies of ALDH1A2 rs3204689 between the KBD and control groups (χ² = 7.113, P = 0.029); the minor allele G of ALDH1A2 was associated with the risk of KBD (χ² = 5.984, P = 0.014). The mRNA and protein levels of ALDH1A2 were increased in the whole blood and cartilage of KBD patients compared with the controls (P = 0.049, P < 0.0001, P = 0.019). Meanwhile, a statistically significant difference was observed between G/G, G/C and C/C genotype on mRNA expression (P = 0.039). The methylation level of the ALDH1A2 gene promoter region showed no significant difference between the KBD and control groups (χ² = 0.317, P = 0.573).

CONCLUSION

Our case-control study indicates that the common variant rs3204689 near ALDH1A2 is associated with KBD in Chinese population. The risk allele G of rs3204689 is statistically linked to the high expression of ALDH1A2, which is up-regulated in the cartilage and whole blood of KBD patients. Our findings suggest a potential role of ALDH1A2 in the pathogenesis of KBD.

Genetic variants associated with platelet count are predictive of human disease and physiological markers

Platelets play an important role in hemostasis and other aspects of vascular biology. We conducted a meta-analysis of platelet count GWAS using data on 536,974 Europeans and identified 577 independent associations. To search for mechanisms through which these variants affect platelets, we applied cis-expression quantitative trait locus, DEPICT and IPA analyses and assessed genetic sharing between platelet count and various traits using polygenic risk scoring. We found genetic sharing between platelet count and counts of other blood cells (except red blood cells), in addition to several other quantitative traits, including markers of cardiovascular, liver and kidney functions, height, and weight. Platelet count polygenic risk score was predictive of myeloproliferative neoplasms, rheumatoid arthritis, ankylosing spondylitis, hypertension, and benign prostate hyperplasia. Taken together, these results advance understanding of diverse aspects of platelet biology and how they affect biological processes in health and disease.

Evgenia Mikaelsdottir et al. report a study of variants associated with platelet count among European individuals where they identify 577 associations. They also report a genetic overlap between platelet count and human diseases, including myeloproliferative neoplasms, rheumatoid arthritis, and hypertension, as well as a genetic overlap between platelet count and various physiological markers.

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.

The Genetic Epidemiology of Joint Shape and the Development of Osteoarthritis

Congruent, low-friction relative movement between the articulating elements of a synovial joint is an essential pre-requisite for sustained, efficient, function. Where disorders of joint formation or maintenance exist, mechanical overloading and osteoarthritis (OA) follow. The heritable component of OA accounts for ~ 50% of susceptible risk. Although almost 100 genetic risk loci for OA have now been identified, and the epidemiological relationship between joint development, joint shape and osteoarthritis is well established, we still have only a limited understanding of the contribution that genetic variation makes to joint shape and how this modulates OA risk. In this article, a brief overview of synovial joint development and its genetic regulation is followed by a review of current knowledge on the genetic epidemiology of established joint shape disorders and common shape variation. A summary of current genetic epidemiology of OA is also given, together with current evidence on the genetic overlap between shape variation and OA. Finally, the established genetic risk loci for both joint shape and osteoarthritis are discussed.

Linking chondrocyte and synovial transcriptional profile to clinical phenotype in osteoarthritis

OBJECTIVES

To determine how gene expression profiles in osteoarthritis joint tissues relate to patient phenotypes and whether molecular subtypes can be reproducibly captured by a molecular classification algorithm.

METHODS

We analysed RNA sequencing data from cartilage and synovium in 113 osteoarthritis patients, applying unsupervised clustering and Multi-Omics Factor Analysis to characterise transcriptional profiles. We tested the association of the molecularly defined patient subgroups with clinical characteristics from electronic health records.

RESULTS

We detected two patient subgroups in low-grade cartilage (showing no/minimal degeneration, cartilage normal/softening only), with differences associated with inflammation, extracellular matrix-related and cell adhesion pathways. The high-inflammation subgroup was associated with female sex (OR 4.12, p=0.0024) and prescription of proton pump inhibitors (OR 4.21, p=0.0040). We identified two independent patient subgroupings in osteoarthritis synovium: one related to inflammation and the other to extracellular matrix and cell adhesion processes. A seven-gene classifier including MMP13, APOD, MMP2, MMP1, CYTL1, IL6 and C15orf48 recapitulated the main axis of molecular heterogeneity in low-grade knee osteoarthritis cartilage (correlation ρ=-0.88, p<10^-10) and was reproducible in an independent patient cohort (ρ=-0.85, p<10^-10).

CONCLUSIONS

These data support the reproducible stratification of osteoarthritis patients by molecular subtype and the exploration of new avenues for tailored treatments.

Sources of Individual Differences in Pain

Pain is an immense clinical and societal challenge, and the key to understanding and treating it is variability. Robust interindividual differences are consistently observed in pain sensitivity, susceptibility to developing painful disorders, and response to analgesic manipulations. This review examines the causes of this variability, including both organismic and environmental sources. Chronic pain development is a textbook example of a gene-environment interaction, requiring both chance initiating events (e.g., trauma, infection) and more immutable risk factors. The focus is on genetic factors, since twin studies have determined that a plurality of the variance likely derives from inherited genetic variants, but sex, age, ethnicity, personality variables, and environmental factors are also considered.

Current Epidemiology and Risk Factors for the Development of Hand Osteoarthritis

PURPOSE OF REVIEW

Hand osteoarthritis (hand OA), the most common peripheral arthritis in the world, is less studied than osteoarthritis (OA) of the knee and hip. However, it is uniquely situated to offer novel insight into OA as a disease process by removing weight-bearing as a confounder of systemic disease mechanisms. Here we review the epidemiology of hand OA and key risk factors for its development.

RECENT FINDINGS

Mounting evidence points to obesity as an important risk factor for hand OA development, with new evidence implicating a role for leptin and serum fatty acids. Disease progression in hand OA and specifically the erosive OA subtype may be associated with diabetes. New evidence supports an association between cardiovascular disease progression and symptomatic hand OA. Alcohol use may be associated with increased synovitis and erosive hand OA. Differences in ethnical distributions of hand OA have become more apparent, with a lower prevalence in Black patients compared to White patients. Novel genetic insights implicating the WNT gene pathway and IL-1β have led to novel potential targets in hand OA pathogenesis. Hand OA is a heterogeneous disease with many modifiable and non-modifiable risk factors that can determine disease severity and shed light on disease pathogenesis.

Genetic and Epigenetic Interplay Within a COLGALT2 Enhancer Associated With Osteoarthritis

OBJECTIVE

The osteoarthritis (OA)-associated single-nucleotide polymorphism (SNP) rs11583641 is located in COLGALT2, encoding a posttranslational modifier of collagen. In cartilage, the SNP genotype correlates with DNA methylation in a putative enhancer. This study was undertaken to characterize the mechanistic relationship between rs11583641, the putative enhancer, and COLGALT2 expression using cartilage samples from human patients and a chondrocyte cell model.

METHODS

Nucleic acids were extracted from articular cartilage samples obtained from patients with OA (n = 137). Samples were genotyped, and DNA methylation was quantified at 12 CpGs using pyrosequencing. The putative enhancer was deleted in Tc28a2 chondrocytes using clustered regularly interspaced short palindromic repeat/Cas9, and the impact on nearby gene expression was determined using real-time quantitative polymerase chain reaction. Targeted modulation of the epigenome using catalytically dead Cas9 (dCas9) constructs fused to DNA methyltransferase 3a or ten-eleven translocase 1 allowed for the investigation of a causal relationship between DNA methylation and enhancer activity.

RESULTS

The genotype at rs11583641 correlated with DNA methylation at 3 CpGs, and the presence of the OA risk allele, C, corresponded to reduced levels of methylation. Deletion of the enhancer resulted in a 2.7-fold reduction in COLGALT2 expression. Targeted methylation and demethylation of the CpGs had antagonistic effects on COLGALT2 expression. An allelic imbalance in the expression of COLGALT2 was identified in the cartilage from patients with OA, with relative overexpression of the OA risk allele. Allelic expression ratios correlated with DNA methylation at 4 CpGs.

CONCLUSION

COLGALT2 is a target of OA genetic risk at this locus. The genotype at rs11583641 impacts DNA methylation in a gene enhancer, which, in turn, modulates COLGALT2 expression. COLGALT2 encodes an enzyme that initiates posttranslational glycosylation of collagens and is therefore a compelling OA susceptibility target.

Familial Clustering of Erosive Hand Osteoarthritis in a Large Statewide Cohort

OBJECTIVE

Erosive hand osteoarthritis (OA) is a severe and rapidly progressing subset of hand OA. Its etiology remains largely unknown, which has hindered development of successful treatments. This study was undertaken to test the hypothesis that erosive hand OA demonstrates familial clustering in a large statewide population linked to genealogical records, and to determine the association of potential risk factors with erosive hand OA.

METHODS

Patients diagnosed as having erosive hand OA were identified by searching 4,741,840 unique medical records from a comprehensive statewide database, the Utah Population Database (UPDB). Affected individuals were mapped to pedigrees to identify high-risk families with excess clustering of erosive hand OA as defined by a familial standardized incidence ratio (FSIR) of ≥2.0. The magnitude of familial risk of erosive hand OA in related individuals was calculated using Cox regression models. Association of potential erosive hand OA risk factors was analyzed using multivariate conditional logistic regression and logistic regression models.

RESULTS

We identified 703 affected individuals linked to 240 unrelated high-risk pedigrees with excess clustering of erosive hand OA (FSIR ≥2.0, P < 0.05). The relative risk of developing erosive hand OA was significantly elevated in first-degree relatives (P < 0.001). There were significant associations between a diagnosis of erosive hand OA and age, sex, diabetes, and obesity (all P < 0.05).

CONCLUSION

Familial clustering of erosive hand OA observed in a statewide database indicates a potential genetic contribution to the etiology of the disease. Age, sex, diabetes, and obesity are risk factors for erosive hand OA. Identification of causal gene variants in these high-risk families may provide insight into the genes and pathways that contribute to erosive hand OA onset and progression.

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10-17), arthritis (GDF5 p = 4 × 10-13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

Understanding osteoarthritis pathogenesis: a multiomics system-based approach

PURPOSE OF REVIEW

Osteoarthritis is a heterogeneous, multifactorial condition regulated by complex biological interactions at multiple levels. Comprehensive understanding of these regulatory interactions is required to develop feasible advances to improve patient outcomes. Improvements in technology have made extensive genomic, transcriptomic, epigenomic, proteomic, and metabolomic profiling possible. This review summarizes findings over the past 20 months related to omics technologies in osteoarthritis and examines how using a multiomics approach is necessary for advancing our understanding of osteoarthritis as a disease to improve precision osteoarthritis treatments.

RECENT FINDINGS

Using the search terms 'genomics' or 'transcriptomics' or 'epigenomics' or 'proteomics' or 'metabolomics' and 'osteoarthritis' from January 1, 2018 to August 31, 2019, we identified advances in omics approaches applied to osteoarthritis. Trends include untargeted whole genome, transcriptome, proteome, and metabolome analyses leading to identification of novel molecular signatures, cell subpopulations and multiomics validation approaches.

SUMMARY

To address the complexity of osteoarthritis, integration of multitissue analyses by multiomics approaches with the inclusion of longitudinal clinical data is necessary for a comprehensive understanding of the disease process, and for appropriate development of efficacious diagnostics, prognostics, and biotherapeutics.

Mesenchyme-specific loss of Dot1L histone methyltransferase leads to skeletal dysplasia phenotype in mice

Chromatin modifying enzymes play essential roles in skeletal development and bone maintenance, and deregulation of epigenetic mechanisms can lead to skeletal growth and malformation disorders. Here, we report a novel skeletal dysplasia phenotype in mice with conditional loss of Disruptor of telomeric silencing 1-like (Dot1L) histone methyltransferase in limb mesenchymal progenitors and downstream descendants. Phenotypic characterizations of mice with Dot1L inactivation by Prrx1-Cre (Dot1L-cKO^Prrx1) revealed limb shortening, abnormal bone morphologies, and forelimb dislocations. Our in vivo and in vitro data support a crucial role for Dot1L in regulating growth plate chondrocyte proliferation and differentiation, extracellular matrix production, and secondary ossification center formation. Micro-computed tomography analysis of femurs revealed that partial loss of Dot1L expression is sufficient to impair trabecular bone formation and microarchitecture in young mice. Moreover, RNAseq analysis of Dot1L deficient chondrocytes implicated Dot1L in the regulation of key genes and pathways necessary to promote cell cycle regulation and skeletal growth. Collectively, our data show that early expression of Dot1L in limb mesenchyme provides essential regulatory control of endochondral bone morphology, growth, and stability.

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.

Assessing the genetic relationships between osteoarthritis and human plasma proteins: a large scale genetic correlation scan

Background

Osteoarthritis (OA) is a multifactorial complex disease. The impact of plasma proteins on OA remains elusive now.

Methods

The UK Biobank genome-wide association study data of OA was used here. Genome-wide SNP genotyping was performed using the Affymetrix UK BiLEVE Axiom or UK Biobank Axiom array. Equally, the GWAS summary data of 3,622 plasma proteins was derived from a recently published study. Consequently, linkage disequilibrium score regression (LD score regression) analysis was performed to evaluate the genetic correlation between each plasma protein and different sites of OA.

Results

Several suggestive plasma proteins were identified for OA. For hand OA, evidence of genetic correlation was observed for inter-alpha-trypsin inhibitor heavy chain H1 (coefficient =−0.3854, P value =0.0198), multiple inositol polyphosphate phosphatase 1 (coefficient =−1.1721, P value =0.0303). For hip OA, 7 suggestive genetic correlation signals were observed, such as Transmembrane glycoprotein NMB (coefficient =0.6944, P value =0.0098), Endothelial cell-specific molecule 1 (coefficient =0.6337, P value =0.03). For Knee OA, 12 suggestive genetic correlation signals were identified, including Elafin (coefficient =−0.5562, P value =0.0092), Interleukin-16 (coefficient =0.3949, P value =0.0435).

Conclusions

We investigated the genetic correlations between plasma proteins and different sites of OA in a systematic way. Our results provide novel evidence that OA is a heterogeneous disease.

Development and validation of a prediction model for incident hand osteoarthritis in the HUNT study

OBJECTIVE

To develop and externally validate prediction models for incident hand osteoarthritis (OA) in a large population-based cohort of middle aged and older men and women.

DESIGN

We included 17,153 men and 18,682 women from a population-based cohort, aged 35-70 years at baseline (1995-1997). Incident hand OA were obtained from diagnostic codes in the Norwegian National Patient Register (1995-2018). We studied whether a range of self-reported and clinically measured predictors could predict hand OA, using the Area Under the receiver-operating Curve (AUC) from logistic regression. External validation of an existing prediction model for male hand OA was tested on discrimination in a sample of men. Bootstrapping was used to avoid overfitting.

RESULTS

The model for men showed modest discriminatory ability (AUC = 0.67, 95% CI 0.62-0.71). Adding a genetic risk score did not improve prediction. Similar discrimination was observed in the model for women (AUC = 0.62, 95% CI 0.59-0.64). Prediction was not improved by adding a genetic risk score or hormonal and reproductive factors. Applying external validation, similar results were observed among men in HUNT (The Nord-Trøndelag Health Study) as in the developmental sample (AUC = 0.62, 95% CI 0.57-0.65).

CONCLUSION

We developed prediction models for incident hand OA in men and women. For women, the model included body mass index (BMI), heavy physical work, high physical activity and perceived poor health. The model showed moderate discrimination. For men, we have shown that a prediction model including BMI, education and information on sleep can predict incident hand OA in several populations with moderate discriminative ability.

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.

Generation and characterization of human induced pluripotent stem cells (iPSCs) from hand osteoarthritis patient-derived fibroblasts

Knowledge and research results about hand osteoarthritis (hOA) are limited due to the lack of samples and animal models of the disease. Here, we report the generation of two induced pluripotent stem cell (iPSC)-lines from patients with radiographic hOA. Furthermore, we wondered whether these iPSC-lines carried single nucleotide polymorphisms (SNPs) within genes that have been associated with hOA. Finally, we performed chondrogenic differentiation of the iPSCs in order to prove their usefulness as cellular models of the disease. We performed a non-integrative reprogramming of dermal fibroblasts obtained from two patients with radiographic rhizarthrosis and non-erosive hOA by introducing the transcriptional factors Oct4, Sox2, Klf4 and c-Myc using Sendai virus. After reprogramming, embryonic stem cell-like colonies emerged in culture, which fulfilled all the criteria to be considered iPSCs. Both iPSC-lines carried variants associated with hOA in the four studied genes and showed differences in their chondrogenic capacity when compared with a healthy control iPSC-line. To our knowledge this is the first time that the generation of iPSC-lines from patients with rhizarthrosis and non-erosive hOA is reported. The obtained iPSC-lines might enable us to model the disease in vitro, and to deeper study both the molecular and cellular mechanisms underlying hOA.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

Interleukin 1 receptor antagonist (IL1RN) gene variants predict radiographic severity of knee osteoarthritis and risk of incident disease

OBJECTIVE

In these studies, we examined the association of single nucleotide polymorphisms (SNPs) of the IL1RN gene with radiographic severity of symptomatic knee osteoarthritis (SKOA) and the risk of incident OA. We also explored these genetic polymorphisms in patients with new onset rheumatoid arthritis (RA).

METHODS

Over 1000 subjects who met American College of Rheumatology criteria for tibiofemoral OA were selected from three independent, National Institute of Health (NIH)-funded cohorts. CTA and TTG haplotypes formed from three SNPs of the IL1RN gene (rs419598, rs315952, rs9005) were assessed for association with radiographic severity, and risk for incident radiographic OA (rOA) in a nested case-control cohort. These IL1RN haplotypes were also assessed for association with disease activity (DAS28) and plasma inflammatory markers in patients with RA.

RESULTS

Carriage of the IL1RN TTG haplotype was associated with increased odds of more severe rOA compared with age-matched, sex-matched and body mass index-matched individuals. Examination of the osteoarthritis initiative Incidence Subcohort demonstrated that carriage of the TTG haplotype was associated with 4.1-fold (p=0.001) increased odds of incident rOA. Plasma IL-1Ra levels were lower in TTG carriers, while chondrocytes from TTG carriers exhibited decreased secretion of IL-1Ra. In patients with RA, the TTG haplotype was associated with increased DAS28, decreased plasma IL-1Ra and elevations of plasma inflammatory markers (hsCRP, interleukin 6 (IL-6)).

CONCLUSION

Carriage of the IL1RN TTG haplotype is associated with more severe rOA, increased risk for incident OA, and increased evidence of inflammation in RA. These data suggest that the IL1RN TTG risk haplotype, associated with decreased IL-1Ra plasma levels, impairs endogenous 'anti-inflammatory' mechanisms.

Integrating Genome-Wide Association Studies With Pathway Analysis and Gene Expression Analysis Highlights Novel Osteoarthritis Risk Pathways and Genes

Osteoarthritis (OA) is the most common degenerative joint disorder worldwide. To identify more genetic signals, genome-wide association study (GWAS) has been widely used and elucidated some OA susceptibility genes. However, these susceptibility genes could only explain only a small part of heritability of OA. It is suggested that the identification of disease-related pathways may contribute to understand the genomic etiology of OA. Here, we integrated the GWAS into pathway analysis to identify novel OA risk pathways. In this study, we first selected 187 independent genetic variants identified by GWAS (P < 1.00E−05) and found that most of these genetic variants are noncoding mutations. We then conducted an expression quantitative trait loci analysis and found that 165 of these 187 genetic variants could significantly regulate the expression of nearby genes. Third, we identified OA susceptibility genes corresponding to these genetic variants, conducted a pathway analysis, and identified novel OA-related KEGG pathways, GO biological processes, GO molecular functions, and GO cellular components. In KEGG database, transforming growth factor β signaling pathway is the most significant signal (P = 5.98E−05) and is the only pathway after the BH multiple-test adjustment with false discovery rate (FDR) = 0.02. In GO database, we identified 24 statistically significant GO biological processes, one statistically significant GO molecular function, and five statistically significant GO cellular components (FDR < 0.05). These signals are related with chondrocyte differentiation and development, which are all known biological pathways associated with OA. Finally, we conducted an OA case–control gene expression analysis to evaluate the differential expression of these OA risk genes. Using an OA case–control gene expression analysis, we showed that 44 risk genes were suggestively differentially expressed in OA cases compared with controls (P < 0.05). Three genes, WWP2, COG5, and MAPT, were statistically differentially expressed in OA cases compared with controls (P < 0.05/122 = 4.10E−04). Hence, our findings may contribute to understanding the genomic etiology of OA.

Hand osteoarthritis: clinical phenotypes, molecular mechanisms and disease management

Osteoarthritis (OA) is a highly prevalent condition, and the hand is the most commonly affected site. Patients with hand OA frequently report symptoms of pain, functional limitations and frustration in undertaking everyday activities. The condition presents clinically with changes to the bone, ligaments, cartilage and synovial tissue, which can be observed using radiography, ultrasonography or MRI. Hand OA is a heterogeneous disorder and is considered to be multifactorial in aetiology. This Review provides an overview of the epidemiology, presentation and burden of hand OA, including an update on hand OA imaging (including the development of novel techniques), disease mechanisms and management. In particular, areas for which new evidence has substantially changed the way we understand, consider and treat hand OA are highlighted. For example, genetic studies, clinical trials and careful prospective imaging studies from the past 5 years are beginning to provide insights into the pathogenesis of hand OA that might uncover new therapeutic targets in the disease.

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.

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.

IL-1 in osteoarthritis: time for a critical review of the literature

The concept of interleukin-1 (IL-1) as a target in osteoarthritis (OA) has been an attractive one for many years. It is a highly potent inducer of cartilage degradation, causing the induction of mRNA and controlling the bioavailability of disease-relevant proteases such as ADAMTS5 and MMP13. It drives synovitis and can induce other disease-relevant genes such as nerve growth factor, a key pain sensitiser in OA. However, the quality of evidence for its involvement in disease is modest. Descriptive studies have demonstrated expression of IL-1α and β in OA cartilage and elevated levels in the synovial fluid of some patients. Agnostic transcriptomic and genomic analyses do not identify IL-1 as a key pathway. In vivo models show a conflicting role for this molecule; early studies using therapeutic approaches in large animal models show a benefit, but most murine studies fail to demonstrate protection where the ligands (IL-1α/β), the cytokine activator (IL-1–converting enzyme), or the receptor (IL-1R) have been knocked out. Recently, a number of large double-blind randomised controlled clinical studies targeting IL-1 have failed. Enthusiasm for IL-1 as a target in OA is rapidly dwindling.

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.

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.