A molecular quantitative trait locus map for osteoarthritis

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Vitamin K2 ameliorates osteoarthritis by suppressing ferroptosis and extracellular matrix degradation through activation GPX4's dual functions

Vitamin K2 (VK2) is an effective compound for anti-ferroptosis and anti-osteoporosis, and Semen sojae praeparatum (Dandouchi in Chinese) is the main source of VK2. Chondrocyte ferroptosis and extracellular matrix (ECM) degradation playing a role in the pathogenesis of osteoarthritis (OA). Glutathione peroxidase 4 (GPX4) is the intersection of two mechanisms in regulating OA progression. But no studies have elucidated the therapeutic effects and mechanisms of VK2 on OA. This study utilized an in vivo rat OA model created via anterior cruciate ligament transection (ACLT) and an in vitro chondrocyte oxidative damage model induced by TBHP to investigate the protective effects and mechanisms of action of VK2 in OA. Knee joint pain in mice was evaluated using the Von Frey test. Micro-CT and Safranin O-Fast Green staining were employed to observe the extent of damage to the tibial cartilage and subchondral bone, while immunohistochemistry and PCR were used to examine GPX4 levels in joint cartilage. The effects of VK2 on rat chondrocyte viability were assessed using CCK-8 and flow cytometry assays, and chondrocyte morphology was observed with toluidine blue and alcian blue staining. The impact of VK2 on intracellular ferroptosis-related markers was observed using fluorescent staining and flow cytometry. Protein expression changes were detected by immunofluorescence and Western blot analysis. Furthermore, specific protein inhibitors were applied to confirm the dual-regulatory effects of VK2 on GPX4. VK2 can increase bone mass and cartilage thickness in the subchondral bone of the tibia, and reduce pain and the OARSI score induced by OA. Immunohistochemistry results indicate that VK2 exerts its anti-OA effects by regulating GPX4 to delay ECM degradation. VK2 can inhibit the activation of the MAPK/NFκB signaling pathway caused by reduced expression of intracellular GPX4, thereby decreasing ECM degradation. Additionally, VK2 can reverse the inhibitory effect of RSL3 on GPX4, increase intracellular GSH content and the GSH/GSSG ratio, reduce MDA content, and rescue chondrocyte ferroptosis. The protective mechanism of VK2 may involve its dual-target regulation of GPX4, reducing chondrocyte ferroptosis and inhibiting the MAPK/NFκB signaling pathway to decelerate the degradation of the chondrocyte extracellular matrix.

Study on the treatment of osteoarthritis by acupuncture combined with traditional Chinese medicine based on pathophysiological mechanism: A review

Osteoarthritis (OA) is a major contributor to disability and social costs in the elderly. As the population ages and becomes increasingly obese, the incidence of the disease is higher than in previous decades. In recent years, important progress has been made in the causes and pathogenesis of OA pain. Modern medical treatment modalities mainly include the specific situation of the patient and focus on the core treatment, including self-management and education, exercise, and related weight loss. As an important part of complementary and alternative medicine, TCM has remarkable curative effect, clinical safety, and diversity of treatment methods in the treatment of OA. Traditional Chinese Medicine treatment of OA has attracted worldwide attention. Therefore, this article will study the pathophysiological mechanism of OA based on modern medicine, and explore the treatment of OA by acupuncture combined with Chinese Medicine.

Epigenomic profiling of the infrapatellar fat pad in osteoarthritis

Osteoarthritis is a prevalent, complex disease of the joints, and affects multiple intra-articular tissues. Here, we have examined genome-wide DNA methylation profiles of primary infrapatellar fat pad and matched blood samples from 70 osteoarthritis patients undergoing total knee replacement surgery. Comparing the DNA methylation profiles between these tissues reveal widespread epigenetic differences. We produce the first genome-wide methylation quantitative trait locus (mQTL) map of fat pad, and make the resource available to the wider community. Using two-sample Mendelian randomization and colocalization analyses, we resolve osteoarthritis GWAS signals and provide insights into the molecular mechanisms underpinning disease aetiopathology. Our findings provide the first view of the epigenetic landscape of infrapatellar fat pad primary tissue in osteoarthritis.

Evolution and advancements in genomics and epigenomics in OA research: How far we have come

OBJECTIVE

Osteoarthritis (OA) is the most prevalent musculoskeletal disease affecting articulating joint tissues, resulting in local and systemic changes that contribute to increased pain and reduced function. Diverse technological advancements have culminated in the advent of high throughput "omic" technologies, enabling identification of comprehensive changes in molecular mediators associated with the disease. Amongst these technologies, genomics and epigenomics - including methylomics and miRNomics, have emerged as important tools to aid our biological understanding of disease.

DESIGN

In this narrative review, we selected articles discussing advancements and applications of these technologies to OA biology and pathology. We discuss how genomics, deoxyribonucleic acid (DNA) methylomics, and miRNomics have uncovered disease-related molecular markers in the local and systemic tissues or fluids of OA patients.

RESULTS

Genomics investigations into the genetic links of OA, including using genome-wide association studies, have evolved to identify 100+ genetic susceptibility markers of OA. Epigenomic investigations of gene methylation status have identified the importance of methylation to OA-related catabolic gene expression. Furthermore, miRNomic studies have identified key microRNA signatures in various tissues and fluids related to OA disease.

CONCLUSIONS

Sharing of standardized, well-annotated omic datasets in curated repositories will be key to enhancing statistical power to detect smaller and targetable changes in the biological signatures underlying OA pathogenesis. Additionally, continued technological developments and analysis methods, including using computational molecular and regulatory networks, are likely to facilitate improved detection of disease-relevant targets, in-turn, supporting precision medicine approaches and new treatment strategies for OA.

A landscape of gene expression regulation for synovium in arthritis

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

Calycosin ameliorates osteoarthritis by regulating the imbalance between chondrocyte synthesis and catabolism

Osteoarthritis (OA) is a severe chronic inflammatory disease. As the main active component of Astragalus mongholicus Bunge, a classic traditional ethnic herb, calycosin exhibits anti-inflammatory action and its mechanism of exact targets for OA have yet to be determined. In this study, we established an anterior cruciate ligament transection (ACLT) mouse model. Mice were randomized to sham, OA, and calycosin groups. Cartilage synthesis markers type II collagen (Col-2) and SRY-Box Transcription Factor 9 (Sox-9) increased significantly after calycosin gavage. While cartilage matrix degradation index cyclooxygenase-2 (COX-2), phosphor-epidermal growth factor receptor (p-EGFR), and matrix metalloproteinase-9 (MMP9) expression were decreased. With the help of network pharmacology and molecular docking, these results were confirmed in chondrocyte ADTC5 cells. Our results indicated that the calycosin treatment significantly improved cartilage damage, this was probably attributed to reversing the imbalance between chondrocyte synthesis and catabolism.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes

OBJECTIVES

Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome.

METHODS

We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci.

RESULTS

We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation.

CONCLUSIONS

The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis

There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

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.

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis

Multimorbidity is a rising public health challenge with important implications for health management and policy. The most common multimorbidity pattern is the combination of cardiometabolic and osteoarticular diseases. Here, we study the genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis. We find genome-wide genetic correlation between the two diseases and robust evidence for association-signal colocalization at 18 genomic regions. We integrate multi-omics and functional information to resolve the colocalizing signals and identify high-confidence effector genes, including FTO and IRX3, which provide proof-of-concept insights into the epidemiologic link between obesity and both diseases. We find enrichment for lipid metabolism and skeletal formation pathways for signals underpinning the knee and hip osteoarthritis comorbidities with type 2 diabetes, respectively. Causal inference analysis identifies complex effects of tissue-specific gene expression on comorbidity outcomes. Our findings provide insights into the biological basis for the type 2 diabetes-osteoarthritis disease co-occurrence.

GPRC5B protects osteoarthritis by regulation of autophagy signaling

Osteoarthritis (OA) is one of the most common chronic diseases in the world. However, current treatment modalities mainly relieve pain and inhibit cartilage degradation, but do not promote cartilage regeneration. In this study, we show that G protein-coupled receptor class C group 5 member B (GPRC5B), an orphan G-protein-couple receptor, not only inhibits cartilage degradation, but also increases cartilage regeneration and thereby is protective against OA. We observed that Gprc5b deficient chondrocytes had an upregulation of cartilage catabolic gene expression, along with downregulation of anabolic genes in vitro. Furthermore, mice deficient in Gprc5b displayed a more severe OA phenotype in the destabilization of the medial meniscus (DMM) induced OA mouse model, with upregulation of cartilage catabolic factors and downregulation of anabolic factors, consistent with our in vitro findings. Overexpression of Gprc5b by lentiviral vectors alleviated the cartilage degeneration in DMM-induced OA mouse model by inhibiting cartilage degradation and promoting regeneration. We also assessed the molecular mechanisms downstream of Gprc5b that may mediate these observed effects and identify the role of protein kinase B (AKT)-mammalian target of rapamycin (mTOR)-autophagy signaling pathway. Thus, we demonstrate an integral role of GPRC5B in OA pathogenesis, and activation of GPRC5B has the potential in preventing the progression of OA.

Mediation of the Same Epigenetic and Transcriptional Effect by Independent Osteoarthritis Risk-Conferring Alleles on a Shared Target Gene, COLGALT2

OBJECTIVE

Over 100 DNA variants have been associated with osteoarthritis (OA), including rs1046934, located within a linkage disequilibrium block encompassing part of COLGALT2 and TSEN15. The present study was undertaken to determine the target gene(s) and the mechanism of action of the OA locus using human fetal cartilage, cartilage from OA and femoral neck fracture arthroplasty patients, and a chondrocyte cell model.

METHODS

Genotyping and methylation array data of DNA from human OA cartilage samples (n = 87) were used to determine whether the rs1046934 genotype is associated with differential DNA methylation at proximal CpGs. Results were replicated in DNA from human arthroplasty (n = 132) and fetal (n = 77) cartilage samples using pyrosequencing. Allelic expression imbalance (AEI) measured the effects of genotype on COLGALT2 and TSEN15 expression. Reporter gene assays and epigenetic editing determined the functional role of regions harboring differentially methylated CpGs. In silico analyses complemented these experiments.

RESULTS

Three differentially methylated CpGs residing within regulatory regions were detected in the human OA cartilage array data, and 2 of these were replicated in human arthroplasty and fetal cartilage. AEI was detected for COLGALT2 and TSEN15, with associations between expression and methylation for COLGALT2. Reporter gene assays confirmed that the CpGs are in chondrocyte enhancers, with epigenetic editing results directly linking methylation with COLGALT2 expression.

CONCLUSION

COLGALT2 is a target of this OA locus. We previously characterized another OA locus, marked by rs11583641, that independently targets COLGALT2. The genotype of rs1046934, like rs11583641, mediates its effect by modulating expression of COLGALT2 via methylation changes to CpGs located in enhancers. Although the single-nucleotide polymorphisms, CpGs, and enhancers are distinct between the 2 independent OA risk loci, their effect on COLGALT2 is the same. COLGALT2 is the target of independent OA risk loci sharing a common mechanism of action.

A Genome-Wide Association Study Meta-Analysis of Alpha Angle Suggests Cam-Type Morphology May Be a Specific Feature of Hip Osteoarthritis in Older Adults

OBJECTIVE

To examine the genetic architecture of cam morphology using alpha angle (AA) as a proxy measure and conduct an AA genome-wide association study (GWAS) followed by Mendelian randomization (MR) to evaluate its causal relationship with hip osteoarthritis (OA).

METHODS

Observational analyses examined associations between AA measurements derived from hip dual x-ray absorptiometry (DXA) scans from the UK Biobank study and radiographic hip OA outcomes and subsequent total hip replacement. Following these analyses, an AA GWAS meta-analysis was performed (N = 44,214) using AA measurements previously derived in the Rotterdam Study. Linkage disequilibrium score regression assessed the genetic correlation between AA and hip OA. Genetic associations considered significant (P < 5 × 10-8 ) were used as AA genetic instrument for 2-sample MR analysis.

RESULTS

DXA-derived AA showed expected associations between AA and radiographic hip OA (adjusted odds ratio [OR] 1.63 [95% confidence interval (95% CI) 1.58, 1.67]) and between AA and total hip replacement (adjusted hazard ratio 1.45 [95% CI 1.33, 1.59]) in the UK Biobank study cohort. The heritability of AA was 10%, and AA had a moderate genetic correlation with hip OA (rg  = 0.26 [95% CI 0.10, 0.43]). Eight independent genetic signals were associated with AA. Two-sample MR provided weak evidence of causal effects of AA on hip OA risk (inverse variance weighted OR 1.84 [95% CI 1.14, 2.96], P = 0.01). In contrast, genetic predisposition for hip OA had stronger evidence of a causal effect on increased AA (inverse variance weighted β = 0.09 [95% CI 0.04, 0.13], P = 4.58 × 10-5 ).

CONCLUSION

Expected observational associations between AA and related clinical outcomes provided face validity for the DXA-derived AA measurements. Evidence of bidirectional associations between AA and hip OA, particularly for risk of hip OA on AA, suggests that hip shape modeling secondary to a genetic predisposition to hip OA contributes to the well-established relationship between hip OA and cam morphology in older adults.

Integrated Analysis of Transcriptome Changes in Osteoarthritis: Gene Expression, Pathways and Alternative Splicing

OBJECTIVE

Osteoarthritis (OA) is the most prevalent joint disease characterized by the degeneration of articular cartilage and the remodeling of its underlying bones, resulting in pain and loss of function in the knees and hips. As far as we know, no curative treatments are available except for the joint replacement. The precise molecular mechanisms which are involved in the degradation of cartilage matrix and development of osteoarthritis are still unclear.

DESIGN

By analyzing RNA-seq data, we found the molecular changes at the transcriptome level such as alternative splicing, gene expression, and molecular pathways in OA knees cartilage.

RESULTS

Expression analysis have identified 457 differential expressed genes including 266 up-regulated genes such as TNFSF15, ST6GALNAC5, TGFBI, ASPM, and TYM, and 191 down-regulated genes such as ADM, JUN, IRE2, PIGA, and MAFF. Gene set enrichment analysis (GSEA) analysis identified down-regulated pathways related to translation, transcription, immunity, PI3K/AKT, and circadian as well as disturbed pathways related to extracellular matrix and collagen. Splicing analysis identified 442 differential alternative splicing events within 284 genes in osteoarthritis, including genes involved in extracellular matrix (ECM) and alternative splicing, and TIA1 was identified as a key regulator of these splicing events.

CONCLUSIONS

These findings provide insights into disease etiology, and offer favorable information to support the development of more effective interventions in response to the global clinical challenge of osteoarthritis.

Osteoarthritis genetic risk acting on the galactosyltransferase gene COLGALT2 has opposing functional effects in articulating joint tissues

BACKGROUND

Investigation of cartilage and chondrocytes has revealed that the osteoarthritis risk marked by the independent DNA variants rs11583641 and rs1046934 mediate  their effects by decreasing the methylation status of CpG dinucleotides in enhancers and increasing the expression of shared target gene COLGALT2. We set out to investigate if these functional effects operate in a non-cartilaginous joint tissue.

METHODS

Nucleic acids were extracted from the synovium of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at CpGs within the COLGALT2 enhancers. CpGs were tested for enhancer effects using a synovial cell line and a reporter gene assay. DNA methylation was altered using epigenetic editing, with the impact on gene expression determined using quantitative polymerase chain reaction. In silico analysis complemented laboratory experiments.

RESULTS

The rs1046934 genotype did not associate with DNA methylation or COLGALT2 expression in the synovium, whereas the rs11583641 genotype did. Surprisingly, the effects for rs11583641 were opposite to those previously observed in cartilage. Epigenetic editing in synovial cells revealed that enhancer methylation is causally linked to COLGALT2 expression.

CONCLUSIONS

This is the first direct demonstration for osteoarthritis genetic risk of a functional link between DNA methylation and gene expression operating in opposite directions between articular joint tissues. It highlights pleiotropy in the action of osteoarthritis risk and provides a cautionary note in the application of future genetically based osteoarthritis therapies: an intervention that decreases the detrimental effect of a risk allele in one joint tissue may inadvertently increase its detrimental effect in another joint tissue.

Therapeutic Perspectives for Inflammation and Senescence in Osteoarthritis Using Mesenchymal Stem Cells, Mesenchymal Stem Cell-Derived Extracellular Vesicles and Senolytic Agents

Osteoarthritis (OA) is the most common cause of disability worldwide among the elderly. Alarmingly, the incidence of OA in individuals less than 40 years of age is rising, likely due to the increase in obesity and post-traumatic osteoarthritis (PTOA). In recent years, due to a better understanding of the underlying pathophysiology of OA, several potential therapeutic approaches targeting specific molecular pathways have been identified. In particular, the role of inflammation and the immune system has been increasingly recognized as important in a variety of musculoskeletal diseases, including OA. Similarly, higher levels of host cellular senescence, characterized by cessation of cell division and the secretion of a senescence-associated secretory phenotype (SASP) within the local tissue microenvironments, have also been linked to OA and its progression. New advances in the field, including stem cell therapies and senolytics, are emerging with the goal of slowing disease progression. Mesenchymal stem/stromal cells (MSCs) are a subset of multipotent adult stem cells that have demonstrated the potential to modulate unchecked inflammation, reverse fibrosis, attenuate pain, and potentially treat patients with OA. Numerous studies have demonstrated the potential of MSC extracellular vesicles (EVs) as cell-free treatments that comply with FDA regulations. EVs, including exosomes and microvesicles, are released by numerous cell types and are increasingly recognized as playing a critical role in cell-cell communication in age-related diseases, including OA. Treatment strategies for OA are being developed that target senescent cells and the paracrine and autocrine secretions of SASP. This article highlights the encouraging potential for MSC or MSC-derived products alone or in combination with senolytics to control patient symptoms and potentially mitigate the progression of OA. We will also explore the application of genomic principles to the study of OA and the potential for the discovery of OA phenotypes that can motivate more precise patient-driven treatments.

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.

Responding to ACL Injury and its Treatments: Comparative Gene Expression between Articular Cartilage and Synovium

The relationship between cartilage and synovium is a rapidly growing area of osteoarthritis research. However, to the best of our knowledge, the relationships in gene expression between these two tissues have not been explored in mid-stage disease development. The current study compared the transcriptomes of these two tissues in a large animal model one year following posttraumatic osteoarthritis induction and multiple surgical treatment modalities. Thirty-six Yucatan minipigs underwent transection of the anterior cruciate ligament. Subjects were randomized to no further intervention, ligament reconstruction, or ligament repair augmented with an extracellular matrix (ECM) scaffold, followed by RNA sequencing of the articular cartilage and synovium at 52 weeks after harvest. Twelve intact contralateral knees served as controls. Across all treatment modalities, the primary difference in the transcriptomes was that the articular cartilage had greater upregulation of genes related to immune activation compared to the synovium-once baseline differences between cartilage and synovium were adjusted for. Oppositely, synovium featured greater upregulation of genes related to Wnt signaling compared to articular cartilage. After adjusting for expression differences between cartilage and synovium seen following ligament reconstruction, ligament repair with an ECM scaffold upregulated pathways related to ion homeostasis, tissue remodeling, and collagen catabolism in cartilage relative to synovium. These findings implicate inflammatory pathways within cartilage in the mid-stage development of posttraumatic osteoarthritis, independent of surgical treatment. Moreover, use of an ECM scaffold may exert a chondroprotective effect over gold-standard reconstruction through preferentially activating ion homeostatic and tissue remodeling pathways within cartilage.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

Osteoarthritis year in review: genetics, genomics and epigenetics

This "year in review" provides a summary of the research findings on the topic of genetics, genomics and epigenetics for osteoarthritis (OA) between Mar 2021-Apr 2022. A search routine of the literature in PubMed for the keyword, osteoarthritis, together with topics on genetics, genomics, epigenetics, polymorphism, DNA methylation, noncoding RNA, lncRNA, proteomics, and single cell RNA sequencing, returned key research articles and relevant reviews. Following filtering of duplicates across search routines, 695 unique research articles and 112 reviews were identified. We manually curated these articles and selected 90 as references for this review. However, we were unable to refer to all these articles, and only used selected articles to highlight key outcomes and trends. The trend in genetics is on the meta-analysis of existing cohorts with comparable genetic and phenotype characterisation of OA; in particular, clear definition of endophenotypes to enhance the genetic power. Further, many researchers are realizing the power of big data and multi-omics approaches to gain molecular insights for OA, and this has opened innovative approaches to include transcriptomics and epigenetics data as quantitative trait loci (QTLs). Given that most of the genetic loci for OA are not located within coding regions of genes, implying the impact is likely to be on gene regulation, epigenetics is a hot topic, and there is a surge in studies relating to the role of miRNA and long non-coding RNA on cartilage biology and pathology. The findings are exciting and new insights are provided in this review to summarize a year of research and the road map to capture all new innovations to achieve the desired goal in OA prevention and treatment.

Causal Effects of Plasma Proteome on Osteoporosis and Osteoarthritis

The two-sample Mendelian randomization (MR) study revealed a causal association of plasma proteins with osteoporosis (OP) and osteoarthritis (OA). Bone mineral density (BMD) is the gold standard for the clinical assessment of OP. Recent studies have shown that plasma proteins play an essential role in the regulation of bone development. However, the causal association of plasma proteins with BMD and OA remains unclear. We estimated the effects of 2889 plasma proteins on 2 BMD phenotypes and 6 OA phenotypes using two-sample MR analysis based on the genome-wide association study summary statistics. Then, we performed sensitivity analysis and reverse-direction MR analysis to evaluate the robustness of the MR analysis results, followed by gene ontology (GO) enrichment analysis and KEGG pathway analysis to explore the functional relevance of the identified plasma proteins. Overall, we observed a total of 257 protein-estimated heel BMD associations, 17 protein-total-body BMD associations, 2 protein-all-OA associations, and 2 protein-knee-OA associations at P_FDR < 0.05. Reverse-direction MR analysis demonstrated that there was little evidence of the causal association of BMD and OA with plasma proteins. GO enrichment analysis and KEGG pathway analysis identified multiple pathways, which may be involved in the development of OP and OA. Our findings recognized plasma proteins that could be used to regulate changes in OP and OA, thus, providing new insights into protein-mediated mechanisms of bone development.

Pitfalls and opportunities for applying latent variables in single-cell eQTL analyses

Using latent variables in gene expression data can help correct unobserved confounders and increase statistical power for expression quantitative trait Loci (eQTL) detection. The probabilistic estimation of expression residuals (PEER) and principal component analysis (PCA) are widely used methods that can remove unwanted variation and improve eQTL discovery power in bulk RNA-seq analysis. However, their performance has not been evaluated extensively in single-cell eQTL analysis, especially for different cell types. Potential challenges arise due to the structure of single-cell RNA-seq data, including sparsity, skewness, and mean-variance relationship. Here, we show by a series of analyses that PEER and PCA require additional quality control and data transformation steps on the pseudo-bulk matrix to obtain valid latent variables; otherwise, it can result in highly correlated factors (Pearson's correlation r = 0.63 ~ 0.99). Incorporating valid PFs/PCs in the eQTL association model would identify 1.7 ~ 13.3% more eGenes. Sensitivity analysis showed that the pattern of change between the number of eGenes detected and fitted PFs/PCs varied significantly in different cell types. In addition, using highly variable genes to generate latent variables could achieve similar eGenes discovery power as using all genes but save considerable computational resources (~ 6.2-fold faster).

Insights from multi-omics integration in complex disease primary tissues

Genome-wide association studies (GWAS) have provided insights into the genetic basis of complex diseases. In the next step, integrative multi-omics approaches can characterize molecular profiles in relevant primary tissues to reveal the mechanisms that underlie disease development. Here, we highlight recent progress in four examples of complex diseases generated by integrative studies: type 2 diabetes (T2D), osteoarthritis, Alzheimer's disease (AD), and systemic lupus erythematosus (SLE). High-resolution methodologies such as single-cell and spatial omics techniques will become even more important in the future. Furthermore, we emphasize the urgent need to include as yet understudied cell types and increase the diversity of studied populations.

3D chromatin structure in chondrocytes identifies putative osteoarthritis risk genes

Genome-wide association studies have identified over 100 loci associated with osteoarthritis risk, but the majority of osteoarthritis risk variants are noncoding, making it difficult to identify the impacted genes for further study and therapeutic development. To address this need, we used a multiomic approach and genome editing to identify and functionally characterize potential osteoarthritis risk genes. Computational analysis of genome-wide association studies and ChIP-seq data revealed that chondrocyte regulatory loci are enriched for osteoarthritis risk variants. We constructed a chondrocyte-specific regulatory network by mapping 3D chromatin structure and active enhancers in human chondrocytes. We then intersected these data with our previously collected RNA-seq dataset of chondrocytes responding to fibronectin fragment, a known osteoarthritis trigger. Integration of the 3 genomic datasets with recently reported osteoarthritis genome-wide association study variants revealed a refined set of putative causal osteoarthritis variants and their potential target genes. One of the putative target genes identified was SOCS2, which was connected to a putative causal variant by a 170-kb loop and is differentially regulated in response to fibronectin fragment. CRISPR-Cas9-mediated deletion of SOCS2 in primary human chondrocytes from 3 independent donors led to heightened expression of inflammatory markers after fibronectin fragment treatment. These data suggest that SOCS2 plays a role in resolving inflammation in response to cartilage matrix damage and provides a possible mechanistic explanation for its influence on osteoarthritis risk. In total, we identified 56 unique putative osteoarthritis risk genes for further research and potential therapeutic development.

Identification of TMEM129, encoding a ubiquitin-protein ligase, as an effector gene of osteoarthritis genetic risk

BACKGROUND

Osteoarthritis is highly heritable and genome-wide studies have identified single nucleotide polymorphisms (SNPs) associated with the disease. One such locus is marked by SNP rs11732213 (T > C). Genotype at rs11732213 correlates with the methylation levels of nearby CpG dinucleotides (CpGs), forming a methylation quantitative trait locus (mQTL). This study investigated the regulatory activity of the CpGs to identify a target gene of the locus.

METHODS

Nucleic acids were extracted from the articular cartilage of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at 14 CpGs within a 259-bp interval. CpGs were tested for enhancer effects in immortalised chondrocytes using a reporter gene assay. DNA methylation at the locus was altered using targeted epigenome editing, with the impact on gene expression determined using quantitative polymerase chain reaction.

RESULTS

rs11732213 genotype correlated with DNA methylation at nine CpGs, which formed a differentially methylated region (DMR), with the osteoarthritis risk allele T corresponding to reduced levels of methylation. The DMR acted as an enhancer and demethylation of the CpGs altered expression of TMEM129. Allelic imbalance in TMEM129 expression was identified in cartilage, with under-expression of the risk allele.

CONCLUSIONS

TMEM129 is a target of osteoarthritis genetic risk at this locus. Genotype at rs11732213 impacts DNA methylation at the enhancer, which, in turn, modulates TMEM129 expression. TMEM129 encodes an enzyme involved in protein degradation within the endoplasmic reticulum, a process previously implicated in osteoarthritis. TMEM129 is a compelling osteoarthritis susceptibility target.

An epigenome-wide view of osteoarthritis in primary tissues

Osteoarthritis is a complex degenerative joint disease. Here, we investigate matched genotype and methylation profiles of primary chondrocytes from macroscopically intact (low-grade) and degraded (high-grade) osteoarthritis cartilage and from synoviocytes collected from 98 osteoarthritis-affected individuals undergoing knee replacement surgery. We perform an epigenome-wide association study of knee cartilage degeneration and report robustly replicating methylation markers, which reveal an etiologic mechanism linked to the migration of epithelial cells. Using machine learning, we derive methylation models of cartilage degeneration, which we validate with 82% accuracy in independent data. We report a genome-wide methylation quantitative trait locus (mQTL) map of articular cartilage and synovium and identify 18 disease-grade-specific mQTLs in osteoarthritis cartilage. We resolve osteoarthritis GWAS loci through causal inference and colocalization analyses and decipher the epigenetic mechanisms that mediate the effect of genotype on disease risk. Together, our findings provide enhanced insights into epigenetic mechanisms underlying osteoarthritis in primary tissues.

Cell-free fat extract attenuates osteoarthritis via chondrocytes regeneration and macrophages immunomodulation

Background

The prevalence of osteoarthritis (OA) is increasing, yet clinically effective and economical treatments are unavailable. We have previously proposed a cell-free fat extract (CEFFE) containing multiple cytokines, which possessed antiapoptotic, anti-oxidative, and proliferation promotion functions, as a “cell-free” strategy. In this study, we aimed to evaluate the therapeutic effect of CEFFE in vivo and in vitro.

Methods

In vivo study, sodium iodoacetate-induced OA rats were treated with CEFFE by intra-articular injections for 8 weeks. Behavioral experiments were performed every two weeks. Histological analyses, anti-type II collagen, and toluidine staining provided structural evaluation. Macrophage infiltration was assessed by anti-CD68 and anti-CD206 staining. In vitro study, the effect of CEFFE on macrophage polarization and secretory factors was evaluated by flow cytometry, immunofluorescence, and quantitative reverse-transcription polymerase chain reaction (qRT-PCR). The effect of CEFFE on cartilage regeneration was accessed by cell counting kit-8 assay and qRT-PCR. The generation of reactive oxygen species (ROS) and levels of ROS-related enzymes were investigated by qRT-PCR and western blotting.

Results

In rat models with sodium iodoacetate (MIA)-induced OA, CEFFE increased claw retraction pressure while decreasing bipedal pressure in a dose-dependent manner. Moreover, CEFFE promoted cartilage structure restoration and increased the proportion of CD206⁺ macrophages in the synovium. In vitro, CEFFE decreased the proportion of CD86⁺ cells and reduced the expression of pro-inflammatory factors in LPS + IFN-γ induced Raw 264.7. In addition, CEFFE decreased the expression of interleukin-6 and ADAMTs-5 and promoted the expression of SOX-9 in mouse primary chondrocytes. Besides, CEFFE reduced the intracellular levels of reactive oxygen species in both in vitro models through regulating ROS-related enzymes.

Conclusions

CEFFE inhibits the progression of OA by promoting cartilage regeneration and limiting low-grade joint inflammation.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02813-3.

Translating osteoarthritis genetics research: challenging times ahead

The ultimate goal of molecular genetic studies of human diseases is to translate the discoveries for patient benefit. For diseases that lack licensed disease-modifying therapeutics, such as osteoarthritis (OA), the need is acute. OA is polygenic and affects older individuals, with a recent genome-wide study of over 800 000 individuals adding 52 novel association signals to those already reported on for this common arthritis. Many of the predicted effector genes of these signals encode proteins that are targets of drugs for other indications, highlighting repurposing opportunities. Here, the potential for OA genetic data to translate is discussed, including whether the developmental origin of OA will limit the application of genetic risk data for disease-modification purposes.

Osteoarthritis year in review: genetics, genomics, epigenetics

OBJECTIVE

In this review, we have highlighted the advances over the past year in genetics, genomics and epigenetics in the field of osteoarthritis (OA).

METHODS

A literature search of PubMed was performed using the criteria: "osteoarthritis" and one of the following terms "genetic(s), genomic(s), epigenetic(s), polymorphism, noncoding ribonucleic acid (RNA), microRNA, long noncoding RNA, lncRNA, circular RNA, RNA sequencing (RNA-seq), single cell sequencing, transcriptomics, or deoxyribonucleic acid (DNA) methylation between April 01, 2020 and April 30, 2021.

RESULTS

In total we identified 765 unique publications, which eventually reduced to 380 of relevance to the field as judged by two assessors. Many of these studies included multiple search terms. We summarised advances relating to genetics, functional genetics, genomics and epigenetics, focusing on our personal key papers during the year.

CONCLUSIONS

This year few studies have identified new genetic variants contributing to OA susceptibility, but a focus has been on refining risk loci or their functional validation. The use of new technologies together with investigating the cross-talk between multiple tissue types, greater sample sizes and/or better patient classification (OA subtypes) will continue to increase our knowledge of disease mechanisms and progress towards understanding and treating OA.

A molecular map of long non-coding RNA expression, isoform switching and alternative splicing in osteoarthritis

Osteoarthritis is a prevalent joint disease and a major cause of disability worldwide with no curative therapy. Development of disease-modifying therapies requires a better understanding of the molecular mechanisms underpinning disease. A hallmark of osteoarthritis is cartilage degradation. To define molecular events characterizing osteoarthritis at the whole transcriptome level, we performed deep RNA sequencing in paired samples of low- and high-osteoarthritis grade knee cartilage derived from 124 patients undergoing total joint replacement. We detected differential expression between low- and high-osteoarthritis grade articular cartilage for 365 genes and identified a 38-gene signature in osteoarthritis cartilage by replicating our findings in an independent dataset. We also found differential expression for 25 novel long non-coding RNA genes (lncRNAs) and identified potential lncRNA interactions with RNA-binding proteins in osteoarthritis. We assessed alterations in the relative usage of individual gene transcripts and identified differential transcript usage for 82 genes, including ABI3BP, coding for an extracellular matrix protein, AKT1S1, a negative regulator of the mTOR pathway and TPRM4, coding for a transient receptor potential channel. We further assessed genome-wide differential splicing, for the first time in osteoarthritis, and detected differential splicing for 209 genes, which were enriched for extracellular matrix, proteoglycans and integrin surface interactions terms. In the largest study of its kind in osteoarthritis, we find that isoform and splicing changes, in addition to extensive differences in both coding and non-coding sequence expression, are associated with disease and demonstrate a novel layer of genomic complexity to osteoarthritis pathogenesis.

Insights into the molecular landscape of osteoarthritis in human tissues

PURPOSE OF REVIEW

To provide an overview of recent developments in the field of osteoarthritis research with a focus on insights gleaned from the application of different -omic technologies.

RECENT FINDINGS

We searched for osteoarthritis-relevant studies focusing on transcriptomics, epigenomics, proteomics and metabolomics, published since November of 2019. Study designs showed a trend towards characterizing the genomic profile of osteoarthritis-relevant tissues with high resolution, for example either by using single-cell technologies or by considering several -omic levels and disease stages.

SUMMARY

Multitissue interactions (cartilage-subchondral bone; cartilage-synovium) are prevalent in the pathophysiology of osteoarthritis, which is characterized by substantial matrix remodelling in an inflammatory milieu. Subtyping approaches using -omic technologies have contributed to the identification of at least two osteoarthritis endotypes. Studies using data integration approaches have provided molecular maps that are tissue-specific for osteoarthritis and pave the way for expanding these data integration approaches towards a more comprehensive view of disease aetiopathogenesis.

Evaluating Causal Relationship Between Metabolites and Six Cardiovascular Diseases Based on GWAS Summary Statistics

Cardiovascular diseases (CVDs) remain the main cause of morbidity and mortality worldwide. The pathological mechanism and underlying biological processes of these diseases with metabolites remain unclear. In this study, we conducted a two-sample Mendelian randomization (MR) analysis to evaluate the causal effect of metabolites on these diseases by making full use of the latest GWAS summary statistics for 486 metabolites and six major CVDs. Extensive sensitivity analyses were implemented to validate our MR results. We also conducted linkage disequilibrium score regression (LDSC) and colocalization analysis to investigate whether MR findings were driven by genetic similarity or hybridization between LD and disease-associated gene loci. We identified a total of 310 suggestive associations across all metabolites and CVDs, and finally obtained four significant associations, including bradykinin, des-arg(9) (odds ratio [OR] = 1.160, 95% confidence intervals [CIs]: 1.080-1.246, false discovery rate [FDR] = 0.022) on ischemic stroke, N-acetylglycine (OR = 0.946, 95%CIs: 0.920-0.973, FDR = 0.023), X-09026 (OR = 0.845, 95%CIs: 0.779-0.916, FDR = 0.021) and X-14473 (OR = 0.938, 95%CIs = 0.907-0.971, FDR = 0.040) on hypertension. Sensitivity analyses showed that these causal associations were robust, the LDSC and colocalization analyses demonstrated that the identified associations were unlikely confused by LD. Moreover, we identified 15 important metabolic pathways might be involved in the pathogenesis of CVDs. Overall, our work identifies several metabolites that have a causal relationship with CVDs, and improves our understanding of the pathogenesis and treatment strategies for these diseases.

Spatiotemporally controlled calcitonin delivery: Long-term and targeted therapy of skeletal diseases

Bone is a connective tissue that support the entire body and protect the internal organs. However, there are great challenges on curing intractable skeletal diseases such as hypercalcemia, osteoporosis and osteoarthritis. To address these issues, calcitonin (CT) therapy is an effective treatment alternative to regulate calcium metabolism and suppress inflammation response, which are closely related to skeletal diseases. Traditional calcitonin formulation requires frequent administration due to the low bioavailability resulting from the short half-life and abundant calcitonin receptors distributed through the whole body. Therefore, long-term and targeted calcitonin delivery systems (LCDS and TCDS) have been widely explored as the popular strategies to overcome the intrinsic limitations of calcitonin and improve the functions of calcium management and inflammation inhibition in recent years. In this review, we first explain the physiological effects of calcitonin on bone remodeling: (i) inhibitory effects on osteoclasts and (ii) facilitated effects on osteoblasts. Then we summarized four strategies for spatiotemporally controlled delivery of calcitonin: micro-/nanomedicine (e.g. inorganic micro-/nanomedicine, polymeric micro-/nanomedicine and supramolecular assemblies), hydrogels (especially thermosensitive hydrogels), prodrug (PEGylation and targeting design) and hybrid biomaterials. Subsequently, we discussed the application of LCDS and TCDS in treating hypercalcemia, osteoporosis, and arthritis. Understanding and analyzing these advanced calcitonin delivery applications are essential for future development of calcitonin therapies toward skeletal diseases with superior efficacy in clinic.

Distinct Concentration-Dependent Molecular Pathways Regulate Bone Cell Responses to Cobalt and Chromium Exposure from Joint Replacement Prostheses

Systemic cobalt (Co) and chromium (Cr) concentrations may be elevated in patients with metal joint replacement prostheses. Several studies have highlighted the detrimental effects of this exposure on bone cells in vitro, but the underlying mechanisms remain unclear. In this study, we use whole-genome microarrays to comprehensively assess gene expression in primary human osteoblasts, osteoclast precursors and mature resorbing osteoclasts following exposure to clinically relevant circulating versus local periprosthetic tissue concentrations of Co²⁺ and Cr³⁺ ions and CoCr nanoparticles. We also describe the gene expression response in osteoblasts on routinely used prosthesis surfaces in the presence of metal exposure. Our results suggest that systemic levels of metal exposure have no effect on osteoblasts, and primarily inhibit osteoclast differentiation and function via altering the focal adhesion and extracellular matrix interaction pathways. In contrast, periprosthetic levels of metal exposure inhibit both osteoblast and osteoclast activity by altering HIF-1α signaling and endocytic/cytoskeletal genes respectively, as well as increasing inflammatory signaling with mechanistic implications for adverse reactions to metal debris. Furthermore, we identify gene clusters and KEGG pathways for which the expression correlates with increasing Co²⁺:Cr³⁺ concentrations, and has the potential to serve as early markers of metal toxicity. Finally, our study provides a molecular basis for the improved clinical outcomes for hydroxyapatite-coated prostheses that elicit a pro-survival osteogenic gene signature compared to grit-blasted and plasma-sprayed titanium-coated surfaces in the presence of metal exposure.