Whole-genome sequencing identifies rare genotypes in COMP and CHADL associated with high risk of hip osteoarthritis

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.

Elevated secretion of pro-collagen I-alpha and vascular endothelial growth factor as biomarkers of acetabular labrum degeneration and calcification in hip osteoarthritis: An explant study

Background

Hip osteoarthritis (OA) involves structural degeneration of different joint compartments, including femoral head cartilage, periarticular ligaments and the acetabular labrum. However, the molecular mechanisms underlying labrum degeneration in hip OA remain poorly understood.

Aim

To assess secretion of putative biomarkers for OA from explanted human labrum tissues under basal and inflammatory conditions and to determine whether these could differentiate between OA and calcification status compared to fracture controls.

Methods

Intact labrum specimens were collected from patients undergoing joint arthroplasty for primary hip OA (n ​= ​15, mean age 70) or non-OA femoral neck fracture (n ​= ​5, mean age 64). Tissues were dissected in equal-sized samples and explanted for one week. To mimic activation of inflammatory signaling by endogenous damage-associated molecular patterns (DAMP) tissue were stimulated with a toll-like receptor 4 (TLR4) agonist (1 ​μg/mL LPS). The involvement of transforming growth factor-beta (TGF-beta) signaling was evaluated by treatment with a TGF-beta type 1 receptor inhibitor (10 ​μM SB-505124). Secretion of aggrecan (ACAN), pro-collagen-I alpha (Pro-Col-Iα), cartilage oligomeric matrix protein (COMP), interleukin-6 (IL-6) and vascular endothelial growth factor (VEGF) was assessed by enzyme-linked immunosorbent assay (ELISA). Labrum calcification was evaluated by 3D whole mount fluorescent microscopy of ethyl cinnamate-based optically cleared tissues stained with Alcian blue/Alizarin red.

Results

Whole mount microscopy revealed non-OA fracture controls were non-calcified, whereas six OA labra (40%) were partially calcified or ossified. Basal secretion of Pro-Col-Iα and VEGF was increased four-fold in OA versus non-OA labra. Pro-Col-Iα levels were correlated with those of VEGF (r ​= ​0.65) and COMP (r ​= ​0.54). Stimulation of DAMP signaling through TLR4 affected secretion of IL-6, VEGF, COMP and Pro-Col-Iα, with distinct responses between non-OA and OA tissues. Inhibition of TGF-beta signaling specifically reduced elevated secretion of Pro-Col- Iα and VEGF in calcified OA labrum.

Conclusions

Secretion of the putative OA biomarkers Pro-Col-Iα and VEGF is elevated in degenerated human acetabular labrum and may serve as indicators of OA and calcification status. Secretion of both factors was partially regulated by TGF-beta signaling in calcified OA labrum tissues.The Translational potential of this article:Our findings suggest that a biomarker panel consisting of Pro-Col-Iα/VEGF/COMP may be valuable for assessing subradiographic labrum degeneration and calcification in hip OA. Targeting TGF-beta signaling may offer a means to reduce vascular invasion and fibrosis in acetabular labrum tissue.

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.

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.

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.

Lack of Association between (AAT)n Polymorphism of the CNR1 Gene Encoding the Cannabinoid Receptor (CB1) and Patient’s Quality of Life

Genetic factors may predispose persons to decreased pain excitability. One of the interesting modulators affecting pain perception may be polymorphisms of the cannabinoid receptor type 1 (CNR1) gene. In this study, we examined the association between three-nucleotide repeats (AAT) polymorphism located in the 3′UTR non-translational region of CNR1 and the patient’s quality of life after total hip arthroplasty. Our study examined the degree of pain sensation, hip function, and the patient’s performance at defined intervals after elective hip replacement due to degenerative changes. The study included 198 patients (128 women and 70 men). The average age was 67 years. PCR genotyping assay was used to identify the (AAT)n triplet repeat polymorphism in the CNR1 gene. The (AAT)n repeat number was determined by sequencing using a standard sequencing protocol. Our study found no statistically significant association between the degree of pain, hip function, and the change in the degree of disability and the (AAT)n polymorphism in the CNR1 gene, no statistically significant correlations between clinical symptoms, the patient’s age, and the number of AAT repeats, no association between the length of the allele and the degree of pain, hip function, and the change in disability.

Phosphoproteomics reveals the BRAF-ERK1/2 axis as an important pathogenic signaling node in cartilage degeneration

OBJECTIVE

Osteoarthritis (OA) causes gradual cellular alterations, structural anomalies and joint dysfunction. Progressive decline of chondrocyte function plays a vital role on OA pathogenesis. Although protein phosphorylation controls cartilage metabolism, its regulation mechanism in OA remains unclear. Thus, proteomic methods were used to investigate phosphorylation changes in preserved and OA articular cartilage samples, and to explore the intervention targets of phosphorylated kinase.

METHODS

Preserved (control) and lesioned (OA) cartilage samples from OA cases were assessed by phosphoproteomics. Immobilized metal affinity chromatography was performed for phosphopeptide enrichment. Quantitated phosphosites were comparatively assessed in the cartilage sample pair. Kinase-substrate enrichment analyses were carried out for identifying OA-related kinases. BRAF expression in cartilage tissues was assessed by immunohistochemical staining. The effects of BRAF inhibitor on cartilage degeneration were examined in mouse chondrocytes and OA mouse model.

RESULTS

High-sensitivity mass spectrometry-based proteomics revealed 7,471 peptides and 4,375 phosphorylated peptides differing between preserved and lesioned cartilage samples, which represented the significant alteration of kinase hubs and transduction pathways. Phosphoproteomics identified BRAF may be involved in developing OA. BRAF regulated the downstream ERK signaling pathway. In addition, BRAF was upregulated in human OA cartilage as shown by immunohistochemistry. Remarkably, BRAF inhibition alleviated cartilage degradation in a mouse model of OA through its downstream of ERK pathway activation.

CONCLUSIONS

Jointly, these findings provide an overview of phosphoproteomic alterations occurring during cartilage degeneration, identifying the BRAF-ERK1/2 Axis signaling as a potential signaling pathway involved in OA.

Physical therapy as a promising treatment for osteoarthritis: A narrative review

Osteoarthritis (OA) is the most prevalent joint disease and a leading cause of disability in older adults. With an increasing population ageing and obesity, OA is becoming even more prevalent than it was in previous decades. Evidence indicates that OA is caused by the breakdown of joint tissues from mechanical loading and inflammation, but the deeper underlying mechanism of OA pathogenesis remains unclear, hindering efforts to prevent and treat this disease. Pharmacological treatments are mostly related to relieving symptoms, and there is no drug for radical cure. However, compelling evidence suggests that regular practice of resistance exercise may prevent and control the development of several musculoskeletal chronic diseases including OA, which may result in improved quality of life of the patients. In this review, we introduced the current understanding of the mechanism and clinical treatments of OA pathogenesis. We also reviewed the recent study of physical therapy in the treatment of skeletal system disorders, especially in OA. Finally, we discuss the present challenges and promising advantages of physical therapy in OA treatment.

Multi-omics molecular biomarkers and database of osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis in the adult population and is a leading cause of disability. OA-related genetic loci may play an important role in clinical diagnosis and disease progression. With the rapid development of diverse technologies and omics methods, many OA-related public data sets have been accumulated. Here, we retrieved a diverse set of omics experimental results from 159 publications, including genome-wide association study, differentially expressed genes and differential methylation regions, and 2405 classified OA-related gene markers. Meanwhile, based on recent single-cell RNA-seq data from different joints, 5459 cell-type gene markers of joints were collected. The information has been integrated into an online database named OAomics and molecular biomarkers (OAOB). The database (http://ibi.zju.edu.cn/oaobdb/) provides a web server for OA marker genes, omics features and so on. To our knowledge, this is the first database of molecular biomarkers for OA.

Obesity, Inflammation, and Immune System in Osteoarthritis

Obesity remains the most important risk factor for the incidence and progression of osteoarthritis (OA). The leading cause of OA was believed to be overloading the joints due to excess weight which in turn leads to the destruction of articular cartilage. However, recent studies have proved otherwise, various other factors like adipose deposition, insulin resistance, and especially the improper coordination of innate and adaptive immune responses may lead to the initiation and progression of obesity-associated OA. It is becoming increasingly evident that multiple inflammatory cells are recruited into the synovial joint that serves an important role in pathological changes in the synovial joint. Polarization of macrophages and macrophage-produced mediators are extensively studied and linked to the inflammatory and destructive responses in the OA synovium and cartilage. However, the role of other major innate immune cells such as neutrophils, eosinophils, and dendritic cells in the pathogenesis of OA has not been fully evaluated. Although cells of the adaptive immune system contribute to the pathogenesis of obesity-induced OA is still under exploration, a quantity of literature indicates OA synovium has an enriched population of T cells and B cells compared with healthy control. The interplay between a variety of immune cells and other cells that reside in the articular joints may constitute a vicious cycle, leading to pathological changes of the articular joint in obese individuals. This review addresses obesity and the role of all the immune cells that are involved in OA and summarised animal studies and human trials and knowledge gaps between the studies have been highlighted. The review also touches base on the interventions currently in clinical trials, different stages of the testing, and their shortcomings are also discussed to understand the future direction which could help in understanding the multifactorial aspects of OA where inflammation has a significant function.

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.

Engineering Closed-Loop, Autoregulatory Gene Circuits for Osteoarthritis Cell-Based Therapies

PURPOSE OF REVIEW

Genetic engineering offers the possibility to simultaneously target multiple cellular pathways in the joints affected by osteoarthritis (OA). The purpose of this review is to summarize the ongoing efforts to develop disease-modifying osteoarthritis drugs (DMOADs) using genetic engineering, including targeting approaches, genome editing techniques, and delivery methods.

RECENT FINDINGS

Several gene circuits have been developed that reprogram cells to autonomously target inflammation, and their efficacy has been demonstrated in chondrocytes and stem cells. Gene circuits developed for metabolic disorders, such as those targeting insulin resistance and obesity, also have the potential to mitigate the impact of these conditions on OA onset and/or progression. Despite the strides made in characterizing the inflammatory environment of the OA joint, our incomplete understanding of how the multiple regulators interact to control signal transduction, gene transcription, and translation to protein limits the development of targeted disease-modifying therapeutics. Continuous advances in targeted genome editing, combined with online toolkits that simplify the design and production of gene circuits, have the potential to accelerate the discovery and clinical application of multi-target gene circuits with disease-modifying properties for the treatment of OA.

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

Background

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

Methods

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

Results

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

Conclusions

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

Prospero registration

CRD42020187664.

Supplementary Information

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

Tetra-primers ARMS-PCR Based Association Analyses of Synonymous and Intronic Variants in the ADAM12 Gene with Susceptibility to Knee Osteoarthritis: A Case-Control Study

Genetic variations in a disintegrin and metalloprotease 12 (ADAM12) gene may contribute to develop Osteoarthritis (OA) that is characterized by cartilage matrix degradation and osteophytes formation. Therefore, the aim of present study was to analyze the association between the ADAM12 gene variants and knee OA predisposition. Tetra-primers ARMS-PCR was employed, to genotype the ADAM12 gene polymorphisms (rs1044122 and rs1871054) in 400 knee OA patients and equal number of age-matched controls. The association between ADAM12 gene variants and OA susceptibility was estimated using the Chi-square, logistic regression, haplotypes and linkage analyses. A significant association of rs1044122 (genotype: χ² = 18.94; P < 0.001, allele: χ² = 19.10; P < 0.001) and rs1871054 (genotype: χ² = 10.04; P = 0.007, allele: χ² = 10.57; P = 0.001) was observed with increased OA susceptibility. The variant genotype of rs1044122 increased OA risk more than twice [odds ratio (OR) 2.20; P = 0.001] and the risk was higher in females (OR 2.43; P = 0.001). The variant genotype of rs1871054 was perceived to almost double the risk in females (OR 1.97; P = 0.003). Moreover, a significant association of rs1044122 and rs1871054 under the additive genetic model (P < 0.001 and P = 0.002, respectively) was observed. The targeted ADAM12 gene polymorphisms, showed significant association with knee OA susceptibility. Females harboring the polymorphisms might be at risk. Besides, the haplotype CC of rs1044122 and rs1871054 in the ADAM12 gene may double knee OA risk. These findings may help in determining the etiology of OA and recognizing the people at risk of developing knee OA.

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.

Description of Joint Alterations Observed in a Family Carrying p.Asn453Ser COMP Variant: Clinical Phenotypes, In Silico Prediction of Functional Impact on COMP Protein and Stability, and Review of the Literature

The role of genetics in the development of osteoarthritis is well established but the molecular bases are not fully understood. Here, we describe a family carrying a germline mutation in COMP (Cartilage Oligomeric Matrix Protein) associated with three distinct phenotypes. The index case was enrolled for a familial form of idiopathic early-onset osteoarthritis. By screening potential causal genes for osteoarthritis, we identified a heterozygous missense mutation of COMP (c.1358C>T, p.Asn453Ser), absent from genome databases, located on a highly conserved residue and predicted to be deleterious. Molecular dynamics simulation suggests that the mutation destabilizes the overall COMP protein structure and consequently the calcium releases from neighboring calcium binding sites. This mutation was once reported in the literature as causal for severe multiple epiphyseal dysplasia (MED). However, no sign of dysplasia was present in the index case. The mutation was also identified in one of her brothers diagnosed with MED and secondary osteoarthritis, and in her sister affected by an atypical syndrome including peripheral inflammatory arthritis of unknown cause, without osteoarthritis nor dysplasia. This article suggests that this mutation of COMP is not only causal for idiopathic early-onset osteoarthritis or severe MED, but can also be associated to a broad phenotypic variability with always joint alterations.

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.

"Lessons from Rare Forms of Osteoarthritis"

Osteoarthritis (OA) is one of the most prevalent conditions in the world, particularly in the developed world with a significant increase in cases and their predicted impact as we move through the twenty-first century and this will be exacerbated by the covid pandemic. The degeneration of cartilage and bone as part of this condition is becoming better understood but there are still significant challenges in painting a complete picture to recognise all aspects of the condition and what treatment(s) are most appropriate in individual causes. OA encompasses many different types and this causes some of the challenges in fully understanding the condition. There have been examples through history where much has been learnt about common disease(s) from the study of rare or extreme phenotypes, particularly where Mendelian disorders are involved. The often early onset of symptoms combined with the rapid and aggressive pathogenesis of these diseases and their predictable outcomes give an often-under-explored resource. It is these "rarer forms of disease" that William Harvey referred to that offer novel insights into more common conditions through their more extreme presentations. In the case of OA, GWAS analyses demonstrate the multiple genes that are implicated in OA in the general population. In some of these rarer forms, single defective genes are responsible. The extreme phenotypes seen in conditions such as Camptodactyly Arthropathy-Coxa Vara-pericarditis Syndrome, Chondrodysplasias and Alkaptonuria all present potential opportunities for greater understanding of disease pathogenesis, novel therapeutic interventions and diagnostic imaging. This review examines some of the rarer presenting forms of OA and linked conditions, some of the novel discoveries made whilst studying them, and findings on imaging and treatment strategies.

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.

Feingold syndrome type 2 in a patient from China

Feingold syndrome type 2 (FGLDS2, MIM614326) is a genetic congenital malformation syndrome, caused by germline heterozygous deletion of MIR17HG on chromosome 13q31, which is extremely rare worldwide. To date, less than 25 patients have been described in the literature. Here, we report on a 3-year-old girl presented with hip dysplasia, polysyndactyly of the left thumb, brachymesophalangy of the fifth digit, microcephaly, intellectual disability, and growth delay. This is likely to be the first case of Feingold syndrome type 2 ever discovered among Chinese population. Through genetic testing and pedigree analysis, she was identified to have a de novo 4.8-Mb microdeletion at chromosome 13q31.3-q32.1, encompassing MIR17HG, GPC5, and GPC6. Additionally, we detected two common compound heterozygous variants (c.919-2A>G and c.147C>G) in SLC26A4 encoding pendrin protein, as well as a novel heterozygous variant c.985_988del in COMP encoding cartilage oligomeric matrix protein. This case report aims to analyze the microdeletion and the three types of variant detected in the patient, and to explore the association between the genotype and phenotype in patients with Feingold syndrome type 2, which may contribute to further understanding and future diagnosis of this disorder.

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.

RNA Sequencing Reveals Interacting Key Determinants of Osteoarthritis Acting in Subchondral Bone and Articular Cartilage: Identification of IL11 and CHADL as Attractive Treatment Targets

OBJECTIVE

To identify key determinants of the interactive pathophysiologic processes in subchondral bone and cartilage in osteoarthritis (OA).

METHODS

We performed RNA sequencing on macroscopically preserved and lesional OA subchondral bone from patients in the Research Arthritis and Articular Cartilage study who underwent joint replacement surgery due to OA (n = 24 sample pairs: 6 hips and 18 knees). Unsupervised hierarchical clustering and differential expression analyses were conducted. Results were combined with data on previously identified differentially expressed genes in cartilage (partly overlapping samples) as well as data on recently identified OA risk genes.

RESULTS

We identified 1,569 genes that were significantly differentially expressed between lesional and preserved subchondral bone, including CNTNAP2 (fold change [FC] 2.4, false discovery rate [FDR] 3.36 × 10^-5 ) and STMN2 (FC 9.6, FDR 2.36 × 10^-3 ). Among these 1,569 genes, 305 were also differentially expressed, and with the same direction of effect, in cartilage, including the recently recognized OA susceptibility genes IL11 and CHADL. Upon differential expression analysis with stratification for joint site, we identified 509 genes that were exclusively differentially expressed in subchondral bone of the knee, including KLF11 and WNT4. These genes that were differentially expressed exclusively in the knee were enriched for involvement in epigenetic processes, characterized by, e.g., HIST1H3J and HIST1H3H.

CONCLUSION

IL11 and CHADL were among the most consistently differentially expressed genes OA pathophysiology-related genes in both bone and cartilage. As these genes were recently also identified as robust OA risk genes, they classify as attractive therapeutic targets acting on 2 OA-relevant tissues.

Recent Developments in Clinical Applications of Mesenchymal Stem Cells in the Treatment of Rheumatoid Arthritis and Osteoarthritis

Mesenchymal stem cell (MSC) therapies have been used as cell-based treatments for decades, owing to their anti-inflammatory, immunomodulatory, and regenerative properties. With high expectations, many ongoing clinical trials are investigating the safety and efficacy of MSC therapies to treat arthritic diseases. Studies on osteoarthritis (OA) have shown positive clinical outcomes, with improved joint function, pain level, and quality of life. In addition, few clinical MSC trials conducted on rheumatoid arthritis (RA) patients have also displayed some optimistic outlook. The largely positive outcomes in clinical trials without severe side effects establish MSCs as promising tools for arthritis treatment. However, further research is required to investigate its applicability in clinical settings. This review discusses the most recent advances in clinical studies on MSC therapies for OA and RA.

Fine mapping genetic variants associated with age at puberty and sow fertility using SowPro90 genotyping array

Sow fertility traits, such as litter size and the number of lifetime parities produced (reproductive longevity), are economically important. Selection for these traits is difficult because they are lowly heritable and expressed late in life. Age at puberty (AP) is an early indicator of reproductive longevity. Here, we utilized a custom Affymetrix single-nucleotide polymorphisms (SNPs) array (SowPro90) enriched with positional candidate genetic variants for AP and a haplotype-based genome-wide association study to fine map the genetic sources associated with AP and other fertility traits in research (University of Nebraska-Lincoln [UNL]) and commercial sow populations. Five major quantitative trait loci (QTL) located on four Sus scrofa chromosomes (SSC2, SSC7, SSC14, and SSC18) were discovered for AP in the UNL population. Negative correlations (r = -0.96 to -0.10; P < 0.0001) were observed at each QTL between genomic estimated breeding values for AP and reproductive longevity measured as lifetime number of parities (LTNP). Some of the SNPs discovered in the major QTL regions for AP were located in candidate genes with fertility-associated gene ontologies (e.g., P2RX3, NR2F2, OAS1, and PTPN11). These SNPs showed significant (P < 0.05) or suggestive (P < 0.15) associations with AP, reproductive longevity, and litter size traits in the UNL population and litter size traits in the commercial sows. For example, in the UNL population, when the number of favorable alleles of an SNP located in the 3' untranslated region of PTPN11 (SSC14) increased, AP decreased (P < 0.0001), while LTNP increased (P < 0.10). Additionally, a suggestive difference in the observed NR2F2 isoforms usage was hypothesized to be the source of the QTL for puberty onset mapped on SSC7. It will be beneficial to further characterize these candidate SNPs and genes to understand their impact on protein sequence and function, gene expression, splicing process, and how these changes affect the phenotypic variation of fertility traits.

Multi-tissue epigenetic analysis of the osteoarthritis susceptibility locus mapping to the plectin gene PLEC

OBJECTIVE

In cartilage, the osteoarthritis (OA) associated single nucleotide polymorphism (SNP) rs11780978 correlates with differential expression of PLEC, and with differential methylation of PLEC CpG dinucleotides, forming eQTLs and mQTLs respectively. This implies that methylation links chondrocyte genotype and phenotype, thus driving the functional effect of this genetic risk signal. PLEC encodes plectin, a cytoskeletal protein that enables tissues to respond to mechanical forces. We sought to assess whether these PLEC functional effects were cartilage specific.

METHOD

Cartilage, fat pad, synovium and peripheral blood were collected from patients undergoing arthroplasty. PLEC CpGs were analysed for mQTLs and allelic expression imbalance (AEI) was performed to test for eQTLs. Plectin was knocked down in a mesenchymal stem cell (MSC) line using CRISPR/Cas9 and cells phenotyped by RNA-sequencing.

RESULTS

mQTLs were discovered in fat pad, synovium and blood. Their effects were however stronger in the joint tissues and of comparable effect between these tissues. We observed AEI in synovium in the same direction as for cartilage and correlations between methylation and PLEC expression. Knocking-down plectin impacted on pathways reported to have a role in OA, including Wnt signalling, glycosaminoglycan biosynthesis and immune regulation.

CONCLUSIONS

Synovium is also a target of the rs11780978 OA association functionally operating on PLEC. In fat pad, mQTLs were identified but these did not correlate with PLEC expression, suggesting the functional effect is not joint-wide. Our study highlights interplay between genetic risk, DNA methylation and gene expression in OA, and reveals clear differences between tissues from the same diseased joint.

Evaluation of relationship between DNA methyltransferase 3 β gene and the risk of hip osteoarthritis: A case-control study based on a Han Chinese population

AIM

Osteoarthritis (OA) is a chronic degenerative joint disease. Early studies have indicated that genetic and environmental factors contribute to the risk of OA. However, the etiology of OA remains unknown. Our study aimed to evaluate the association of DNMT3B gene with the risk of hip OA in Han Chinese individuals.

METHODS

A total of 2070 subjects were recruited into the study, including 658 patients with hip OA and 1412 healthy controls. Twelve tag single nucleotide polymorphisms (SNPs) were selected and genotyped in our samples. Genetic associations between DNMT3B gene and the risk of hip OA were examined at both the single marker and haplotype levels. Cis-expression quantitative trait loci signals that achieve genome-wide significance of targeted SNPs from multiple types of human tissues were extracted from the GTEx database.

RESULTS

Significant signals were identified for SNP rs2424905 in 4 genetic models. The T allele was significantly associated with an increased risk of hip OA (odds ratio = 1.53; 95% CI = 1.28-1.83). The T allele was also significantly associated with higher Kellgren-Lawrence grade in the patients with hip OA (χ2  = 32.70, P = 1.37 × 10-6 ). Moreover, SNP rs2424905 was significantly associated with the gene expression level of multiple genes, including DNMT3B, C20orf203, COMMD7, EFCAB8, MAPRE1, and RP5-1085F17.3, from several types of human tissues.

CONCLUSION

Our results indicated that rs2424905 of DNMT3B gene contributed to the risk of hip OA and its clinical severity in a Han Chinese population. These findings suggested that rs2424905 of DNMT3B could be a promising genetic marker to assess susceptibility to hip OA in Han Chinese populations.

Multi-level consistent changes of the ECM pathway identified in a typical keratoconus twin's family by multi-omics analysis

BACKGROUND

Keratoconus (KC) is a common, degenerative disorder of the cornea, and genetic factors play a key role in its development. However, the genetic etiology of KC is still unclear. This study used the family of twins as material, using, for the first time, multi-omics analysis, to systematically display the changes in KC candidate factors in patients at the DNA, RNA, and protein levels. These can evaluate candidate pathogenic factors in depth and lock onto pathogenic targets.

RESULTS

The twins in this study presented classic phenotypes, clear diagnoses, complete case data, and clinical samples, which are excellent materials for genetically studying KC. Whole-exome sequencing was conducted on both the twins and their parents. Transcriptome sequencing was conducted on proband's and health individual's primary human corneal fibroblast cells. Quantitative Real-time PCR and western blot were used to validate the differential gene expressions between the proband and controls. By integrating genomics, transcriptome, and protein level data, multiple consecutive events of KC were systematically analyzed to help better understand the molecular mechanism and genetic basis of KC. The results showed that the accumulation of rare, micro-effect risk variants was the pathogenic factor in this Chinese KC family. Consistent changes in extracellular matrices (ECMs) at the DNA and RNA levels suggested that ECM related changes play a key role in KC pathogenesis. The major gene variants (WNT16, CD248, COL6A2, COL4A3 and ADAMTS3) may affect the expression of related collagens or ECM proteins, thus reducing the amount of ECM in corneas and resulting in KC.

CONCLUSIONS

This study, the first to explore the genetic etiology of KC via multi-omics analysis under the polygenetic model, has provided new insights into the genetic mechanisms underlying KC and an effective strategy for studying KC pathogenesis in the future.

Multiple epiphyseal dysplasia and related disorders: Molecular genetics, disease mechanisms, and therapeutic avenues

For the vast majority of the 6000 known rare disease the pathogenic mechanisms are poorly defined and there is little treatment, leading to poor quality of life and high healthcare costs. Genetic skeletal diseases (skeletal dysplasias) are archetypal examples of rare diseases that are chronically debilitating, often life-threatening and for which no treatments are currently available. There are more than 450 unique phenotypes that, although individually rare, have an overall prevalence of at least 1 per 4000 children. Multiple epiphyseal dysplasia (MED) is a clinically and genetically heterogeneous disorder characterized by disproportionate short stature, joint pain, and early-onset osteoarthritis. MED is caused by mutations in the genes encoding important cartilage extracellular matrix proteins, enzymes, and transporter proteins. Recently, through the use of various cell and mouse models, disease mechanisms underlying this diverse phenotypic spectrum are starting to be elucidated. For example, ER stress induced as a consequence of retained misfolded mutant proteins has emerged as a unifying disease mechanisms for several forms of MED in particular and skeletal dysplasia in general. Moreover, targeting ER stress through drug repurposing has become an attractive therapeutic avenue.

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.

Common variants in IL17F gene contributed to the risk of hip osteoarthritis susceptibility in Han Chinese population

AIM

The prevalence of hip and knee osteoarthritis (OA) varies by ethnicity, suggesting genetic heterogeneity in populations and predilection sites. Given the unknown mechanism of IL17F gene in the etiology of OA, it is necessary to examine the potential shared susceptibility of IL17F gene between knee OA and hip OA (HOA). This study aimed to evaluate the association of the IL17F gene and susceptibility to HOA in a Han Chinese population.

METHODS

A total of 2650 study subjects, comprising 796 HOA patients and 1854 controls, were recruited into the present study. Seven tag single nucleotide polymorphisms (SNPs) were selected for genotyping. Single marker-based genetic association analyses were conducted at both the genotypic and allelic levels. χ² statistics were calculated for statistical testing, and odds ratios were obtained to estimate the effects of genotypes and alleles for each SNP.

RESULTS

The SNP rs763780 was identified to be significantly associated with the risk of HOA at both genotypic (χ²  = 12.45, P = .002) and allelic levels (χ²  = 11.83, P = .0006). A linkage disequilibrium (LD) block comprised of 3 SNPs (rs12201582-rs12203736-rs722323) was also significantly associated with the risk of HOA. In addition, rs2294835 was identified to be associated with HOA severity (χ²  = 12.10, P = .02).

CONCLUSION

Our results showed that IL17F gene contributed to the risk of HOA in a Han Chinese population, which would help to elucidate the pathogenesis of OA and facilitate the development of novel medicines and treatments for OA.

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

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

Evaluation of Common Variants in the AKNA Gene and Susceptibility to Knee Osteoarthritis Among the Han Chinese

Background: Osteoarthritis (OA) is a complex degenerative joint disease that is associated with both genetic and environmental factors. The AKNA gene, located at 9q32, has recently been identified as being associated with knee osteoarthritis (KOA) in the Mexican population. Our aim was to investigate the relationship of common variants in this gene with the risk of KOA in a large Han Chinese population. Methods: A total of 2,500 Han Chinese subjects were recruited, consisting of 824 KOA patients and 1,676 controls. Eight tag single nucleotide polymorphisms (SNPs) located within the ANKA gene were selected for genotyping. Single marker-based association analyses were conducted using multiple modes of inheritance, including genotypic, allelic, dominant, and recessive. Haplotype-based association analyses were also performed. Plink was utilized for genetic association analyses. In addition, we examined the GTEx database to test the expression quantitative loci effects of the significant SNPs within the AKNA gene. Results: Among these eight SNPs evaluated we identified one, rs10817595, as being significantly associated with the risk of KOA. Compared to the CC genotype at this locus, the odds ratio (95% confidence interval) for KOA with the AA genotype was 1.58 (1.23-2.01)-fold greater. A linkage disequilibrium block that included this SNP was also determined to be significantly associated with the risk of KOA (χ² = 25.08, p = 3.58 × 10^-6). In general, the minor allele A of SNP rs10817595 was associated with an increased risk of KOA. Conclusion: This study is the first to present evidence for a potential link between the risk of KOA and an AKNA gene polymorphism among persons with a Han Chinese ancestry. Future functional analyses based on animal models and sequencing-based population studies are needed to elucidate the biological plausibility and genetic architecture of AKNA for KOA susceptibility.

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.

[Biochemical predictors of cartilaginous tissue metabolic disorders for early osteoarthrosis evidence diagnostics.]

The complete laboratory and clinical instrumental examination was conducted, it included serum COMP test, circadian excretion of type II collagen C-terminal telopeptides Urine CartiLaps (СТХ II) and Т2 relaxometry in 29 patients of both sexes of the main group with early (0-I) X-ray osteoarthrosis stages, 30 subjects of comparison group with no X-ray osteoarthrosis evidences aged 44.7±5.9 years and 25 healthy subjects aged 26.3±2.6 years of the control group. The increase (р<0,05) of COMP and Urine CartiLaps levels as well as the increase of Т2 relaxation signal was found at early osteoarthrosis evidences. It was proven that there was (р<0.01) a connection (R=0.8) between COMP and Urine CTX II levels as well as (р<0.05) results of Т2 relaxometry (R=0.8). It was proven that collagen anisotropy and formation of chondromalacia areas as Т2 relaxometry showed in patients with early OA evidences were connected with accumulation of serum COMP and increase of type II collagen circadian renal excretion. The combination of laboratory and radiological methods of articular hyaline cartilage assessment may be used for finding early osteoarthrosis stages.

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.

The Genetic Architecture of High Bone Mass

The phenotypic trait of high bone mass (HBM) is an excellent example of the nexus between common and rare disease genetics. HBM may arise from carriage of many 'high bone mineral density [BMD]'-associated alleles, and certainly the genetic architecture of individuals with HBM is enriched with high BMD variants identified through genome-wide association studies of BMD. HBM may also arise as a monogenic skeletal disorder, due to abnormalities in bone formation, bone resorption, and/or bone turnover. Individuals with monogenic disorders of HBM usually, though not invariably, have other skeletal abnormalities (such as mandible enlargement) and thus are best regarded as having a skeletal dysplasia rather than just isolated high BMD. A binary etiological division of HBM into polygenic vs. monogenic, however, would be excessively simplistic: the phenotype of individuals carrying rare variants of large effect can still be modified by their common variant polygenic background, and by the environment. HBM disorders-whether predominantly polygenic or monogenic in origin-are not only interesting clinically and genetically: they provide insights into bone processes that can be exploited therapeutically, with benefits both for individuals with these rare bone disorders and importantly for the many people affected by the commonest bone disease worldwide-i.e., osteoporosis. In this review we detail the genetic architecture of HBM; we provide a conceptual framework for considering HBM in the clinical context; and we discuss monogenic and polygenic causes of HBM with particular emphasis on anabolic causes of HBM.

Engineered Cartilage from Human Chondrocytes with Homozygous Knockout of Cell Cycle Inhibitor p21

Osteoarthritis (OA) is a highly prevalent disease with limited treatment options. The search for disease-modifying OA therapies would benefit from a more comprehensive knowledge of the genetic variants that contribute to chondrocyte dysfunction and the barriers to cartilage regeneration. One goal of this study was to establish a system for producing engineered cartilage tissue from genetically defined primary human chondrocytes through genome editing and single-cell expansion. This process was utilized to investigate the functional effect of biallelic knockout of the cell cycle inhibitor p21. The use of ribonucleoprotein (RNP) CRISPR/Cas9 complexes targeting two sites in the coding region of p21 resulted in a high frequency (16%) of colonies with homozygous p21 knockout. Chondrogenic pellet cultures from expanded chondrocytes with complete loss of p21 produced more glycosaminoglycans (GAG) and maintained a higher cell number. Single-cell-derived colonies retained the potential for robust matrix production after expansion, allowing for analysis of colony variability from the same population of targeted cells. The effect of enhanced cartilage matrix production in p21 knockout chondrocytes persisted when matrix production from individual colonies was analyzed. Chondrocytes had lower levels of p21 protein with further expansion, and the difference in GAG production with p21 knockout was strongest at early passages. These results support previous findings that implicate p21 as a barrier to cartilage matrix production and regenerative capacity. Furthermore, this work establishes the use of genome-edited human chondrocytes as a promising approach for engineered tissue models containing user-defined gene knockouts and other genetic variants for investigation of OA pathogenesis. Impact Statement This work provides two important advances to the field of tissue engineering. One is the demonstration that engineered cartilage tissue can be produced from genetically defined populations of primary human chondrocytes. While CRISPR/Cas-9 genome editing has been extensively used in cell lines that divide indefinitely, this work extends the technique to an engineered tissue model system to support investigation of genetic changes that affect cartilage production. A second contribution is the finding that chondrocytes with p21 knockout synthesized more cartilage matrix tissue than unedited controls. This supports the continued investigation of p21 as a potential barrier to effective cartilage regeneration.

Association of IL-17A-197G/A and IL-17F-7488T/C polymorphisms and osteoarthritis susceptibility: A meta-analysis

AIM

This meta-analysis was conducted with the aim of collecting and synthesizing the existing evidence on the association of interleukin-17A (IL-17A)-197G/A and IL-17F-7488T/C gene polymorphisms and osteoarthritis (OA) susceptibility.

METHODS

Six electronic databases including PubMed, EMBASE, Web of Science, Cochrane Library, CNKI, and Wanfang were systematically searched for potentially relevant studies previous to June 2019. The strengths of association were estimated by summary odds ratios (ORs) and 95% confidence intervals (95%CIs) in a model-free approach. Heterogeneity test and sensitivity analysis were also conducted to guarantee the reliability of this study.

RESULTS

Six eligible case-control studies comprising 1989 OA patients and 2116 healthy controls were obtained for the meta-analysis. Dominant model was confirmed to be the optimal genetic model (MM + Mm vs mm). The pooled estimate supported that IL-17A-197G/A and IL-17F-7488T/C polymorphisms were significantly associated with OA susceptibility in the overall population (IL-17A-197G/A: GG + GA vs AA, OR = 0.69, 95%CI 0.57-0.80; P < .001; IL-17F-7488T/C, TT + TC vs CC, OR = 0.46, 95%CI 0.29-0.71, P < .001). However, subgroup analyses suggested the association only existed in Asians and knee OA.

CONCLUSION

The findings of the present study indicate that IL-17A-197G/A and IL-17F-7488T/C polymorphisms are positively associated to reduced risk of knee OA, especially in Asian populations.

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

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

SUMMARY

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

Discovery and analysis of methylation quantitative trait loci (mQTLs) mapping to novel osteoarthritis genetic risk signals

OBJECTIVE

Osteoarthritis (OA) is polygenic with over 90 independent genome-wide association loci so far reported. A key next step is the identification of target genes and the molecular mechanisms through which this genetic risk operates. The majority of OA risk-conferring alleles are predicted to act by modulating gene expression. DNA methylation at CpG dinucleotides may be a functional conduit through which this occurs and is detectable by mapping methylation quantitative trait loci, or mQTLs. This approach can therefore provide functional insight into OA risk and will prioritize genes for subsequent investigation. That was our goal, with a focus on the largest set of OA loci yet to be reported.

METHOD

We investigated DNA methylation, genotype and RNA sequencing data derived from the cartilage of patients who had undergone arthroplasty and combined this with in silico analyses of expression quantitative trait loci, epigenomes and chromatin interactions.

RESULTS

We investigated 42 OA risk loci and in ten of these we identified 24 CpGs in which methylation correlated with genotype (false discovery rate (FDR) P-values ranging from 0.049 to 1.73x10^-25). In silico analyses of these mQTLs prioritised genes and regulatory elements at the majority of the ten loci, with COLGALT2 (encoding a collagen galactosyltransferase), COL11A2 (encoding a polypeptide chain of type XI collagen) and WWP2 (encoding a ubiquitin ligase active during chondrogenesis) emerging as particularly compelling target genes.

CONCLUSION

We have highlighted the pivotal role of DNA methylation as a link between genetic risk and OA and prioritized genes for further investigation.

Translational genomics and precision medicine: Moving from the lab to the clinic

Translational genomics aims to improve human health by building on discoveries made through genetics research and applying them in the clinical setting. This progress has been made possible by technological advances in genomics and analytics and by the digital revolution. Such advances should enable the development of prognostic markers, tailored interventions, and the design of prophylactic preventive approaches. We are at the cusp of predicting disease risk for some disorders by means of polygenic risk scores integrated with classical epidemiological risk factors. This should lead to better risk stratification and clinical decision-making. A deeper understanding of the link between genome-wide sequence and association with well-characterized phenotypes will empower the development of biomarkers to aid diagnosis, inform disease progression trajectories, and allow better targeting of treatments to those patients most likely to respond.

Functional rare variants influence the clinical response to anti-TNF therapy in Crohn's disease

BACKGROUND

The effect of low-frequency functional variation on anti-tumor necrosis factor alpha (TNF) response in Crohn's disease (CD) patients remains unexplored. The objective of this study was to investigate the impact of functional rare variants in clinical response to anti-TNF therapy in CD.

METHODS

CD anti-TNF naïve patients starting anti-TNF treatment due to active disease [Crohn's Disease Activity Index (CDAI > 150)] were included. The whole genome was sequenced using the Illumina Hiseq4000 platform. Clinical response was defined as a CDAI score <150 at week 14 of anti-TNF treatment. Low-frequency variants were annotated and classified according to their damaging potential. The whole genome of CD patients was screened to identify homozygous loss-of-function (LoF) variants. The TNF signaling pathway was tested for overabundance of damaging variants using the SKAT-O method. Functional implication of the associated rare variation was evaluated using cell-type epigenetic enrichment analyses.

RESULTS

A total of 41 consecutive CD patients were included; 3250 functional rare variants were identified (2682 damaging and 568 LoF variants). Two homozygous LoF mutations were found in HLA-B and HLA-DRB1 genes associated with lack of response and remission, respectively. Genome-wide LoF variants were enriched in epigenetic marks specific for the gastrointestinal tissue (colon, p = 4.11e-4; duodenum, p = 0.011). The burden of damaging variation in the TNF signaling pathway was associated with response to anti-TNF therapy (p = 0.016); damaging variants were enriched in epigenetic marks from CD8⁺ (p = 6.01e-4) and CD4+ (p = 0.032) T cells.

CONCLUSIONS

Functional rare variants are involved in the response to anti-TNF therapy in CD. Cell-type enrichment analysis suggests that the gut mucosa and CD8⁺ T cells are the main mediators of this response.

Inhibition of cyclooxygenase-2 activity in subchondral bone modifies a subtype of osteoarthritis

Osteoarthritis (OA) causes the destruction of joints. Its pathogenesis is still under investigation, and there is no effective disease-modifying therapy. Here, we report that elevated cyclooxygenase-2 (COX-2) expression in the osteocytes of subchondral bone causes both spontaneous OA and rheumatoid arthritis (RA). The knockout of COX-2 in osteocytes or treatment with a COX-2 inhibitor effectively rescues the structure of subchondral bone and attenuates cartilage degeneration in spontaneous OA (STR/Ort) mice and tumor necrosis factor-α transgenic RA mice. Thus, elevated COX-2 expression in subchondral bone induces both OA-associated and RA-associated joint cartilage degeneration. The inhibition of COX-2 expression can potentially modify joint destruction in patients with arthritis.

Alternative Interleukin 17A/F Locus Haplotypes are Associated With Increased Risk to Hip and Knee Osteoarthritis

We studied the genetic epidemiology of primary large-joint (hip and knee) osteoarthritis (OA), in order to find disease risk factors by a candidate-gene approach. We used case-control study in the Croatian Caucasian population. We genotyped 500 OA patients (260 hip, 240 knee; both with total joint replacements) and 597 healthy individuals for single-nucleotide polymorphisms (SNPs) in interleukin 17A (IL17A) (rs2275913) and IL17F (rs763780 and rs1889570) genes. On the basis of our population and allelic and genotypic frequencies haplotypes were predicted by PHASE software and compared between patients and controls. The three-SNP haplotype (rs2275913-rs763780-rs1889570) G-C-A confers predisposition to hip (p < 0.005) but not knee OA. The three-SNP haplotype having opposed nucleotides A-T-G was found significantly associated with 2.6 times higher risk for developing knee (p < 0.02) but not hip OA. The haplotype G-T (IL17A-IL17F; rs2275913-rs763780) is associated with protection to the disease in hip OA (p < 0.01). Our analyses show that two disparate haplotypes within the IL17A-F gene locus are associated with higher risk to developing hip and knee OA in the Croatian population. The data might suggest a difference in the etiology of hip OA from that of the knee OA, perhaps due to an unknown dissimilarity in vulnerability of these joints to the actions of IL17. Alternatively, other differences in genetic factors like the long non-protein coding region LINCMD1 and/or microRNA species like miR133b and miR206 found in the vicinity of the IL17 locus might be involved in the observed risk. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1972-1978, 2019.

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.