The GDF5 rs143383 polymorphism is associated with osteoarthritis of the knee with genome-wide statistical significance

Transcriptome Sequencing Analysis of Genes Associated with Different Developmental Periods of the Ovarian Follicle in the Duolang Sheep

BACKGROUND

The ovaries are crucial reproductive organs in female mammals, directly influencing the reproductive efficiency and productivity of these animals. The Duolang sheep, native to Xinjiang, is known for its rapid growth and high fertility. However, the mechanisms underlying ovarian follicle development and regulation in sheep remain unclear.

METHODS

Employing transcriptome sequencing technology, this study methodically analyzed ovaries from sheep across various estrous cycles to uncover key genes and signaling pathways that play a role in the development of ovarian follicles.

RESULTS

The results indicated that a total of 130, 183, and 175 differentially expressed genes were identified in the DTA/DTB, DTB/DTC, and DTA/DTC groups, respectively. Key genes like BAG3, GDF5, RHOB, RUNX2, LGALS3, and CDH1, along with pathways such as endoplasmic reticulum protein processing, the NOTCH signaling pathway, and the MAPK signaling pathway, were found to be involved. RT-qPCR confirmed the differential expression of BAG3, RHOB, and RUNX2.

CONCLUSIONS

This research provides insights into the molecular mechanisms of ovarian follicle development and a basis for enhancing the reproductive performance of Duolang sheep.

The Genetics and Functional Genomics of Osteoarthritis

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

Correlation of growth differentiation factor-5 + 104T>C polymorphism with the risk of knee, hand, and hip osteoarthritis: a case-control study and meta-analysis based on 47 case-control studies

Osteoarthritis (OA) arises from a intricate interplay of genetic and environmental factors. Numerous studies have explored the link between the growth differentiation factor 5 (GDF-5) +104T>C polymorphism and OA risk, but the findings have been inconclusive. We carried out a case-control study with 704 OA cases and 418 healthy controls. Furthermore, we conducted a meta-analysis by thoroughly searching the literature for relevant studies published until 1 September, 2023. The combined odds ratio and 95% confidence intervals were used to assess the correlation's strength. A total of 47 independent case-control studies, including 17,602 OA cases and 30,947 controls, were analyzed. Of these, 31 studies (11,176 cases, 16,724 controls) focused on knee OA, 8 studies (3,973 cases, 8,055 controls) examined hip OA, and 6 studies (2244 cases, 5965 controls) investigated hand OA. Overall, our findings suggest that the GDF-5 + 104T>C polymorphism has a protectibe role in development of OA in global scale. Subgroup analyses by ethnicity indicated that this genetic variation provides protection against OA in Caucasian, Asian, and African populations. Further subgroup analysis based on the type of OA showed a decreased risk of knee and hand OA associated with this variation, but not for hip OA. Our combined data indicates that the GDF-5 + 104T>C polymorphism offers protection against the development of OA in general, as well as knee and hand OA. Nevertheless, there was no correlation found between this polymorphism and the development of hip OA.

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.

Transforming growth factor beta signaling and craniofacial development: modeling human diseases in zebrafish

Humans and other jawed vertebrates rely heavily on their craniofacial skeleton for eating, breathing, and communicating. As such, it is vital that the elements of the craniofacial skeleton develop properly during embryogenesis to ensure a high quality of life and evolutionary fitness. Indeed, craniofacial abnormalities, including cleft palate and craniosynostosis, represent some of the most common congenital abnormalities in newborns. Like many other organ systems, the development of the craniofacial skeleton is complex, relying on specification and migration of the neural crest, patterning of the pharyngeal arches, and morphogenesis of each skeletal element into its final form. These processes must be carefully coordinated and integrated. One way this is achieved is through the spatial and temporal deployment of cell signaling pathways. Recent studies conducted using the zebrafish model underscore the importance of the Transforming Growth Factor Beta (TGF-β) and Bone Morphogenetic Protein (BMP) pathways in craniofacial development. Although both pathways contain similar components, each pathway results in unique outcomes on a cellular level. In this review, we will cover studies conducted using zebrafish that show the necessity of these pathways in each stage of craniofacial development, starting with the induction of the neural crest, and ending with the morphogenesis of craniofacial elements. We will also cover human skeletal and craniofacial diseases and malformations caused by mutations in the components of these pathways (e.g., cleft palate, craniosynostosis, etc.) and the potential utility of zebrafish in studying the etiology of these diseases. We will also briefly cover the utility of the zebrafish model in joint development and biology and discuss the role of TGF-β/BMP signaling in these processes and the diseases that result from aberrancies in these pathways, including osteoarthritis and multiple synostoses syndrome. Overall, this review will demonstrate the critical roles of TGF-β/BMP signaling in craniofacial development and show the utility of the zebrafish model in development and disease.

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

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

Decreased Risk of Knee Osteoarthritis with Taller Height in an East Asian Population: A Nationwide Cohort Study

BACKGROUND

Numerous studies have explored factors impacting osteoarthritis (OA), but its relationship with height remains uncertain. This study investigates the relationship between height and osteoarthritis risk in South Korea.

METHODS

Participants aged 50 or older who underwent health screenings in 2009 were selected from the National Health Insurance System database. A total of 1,138,904 subjects were divided into height quartiles (Q1-Q4) based on age and gender. Cox proportional hazard models were used to assess knee osteoarthritis incidence risk, with the shortest quartile (Q1) as the reference.

RESULTS

After adjusting for age, sex, income, smoking, drinking, exercise, hypertension, diabetes mellitus, dyslipidemia, and body mass index (BMI), no significant difference in OA incidence risk based on height was observed. However, when adjusted for weight instead of BMI, we observed a gradual decrease in hazard ratio with increasing height. The hazard ratio for the tallest group was 0.787 (95% CI, 0.781~0.795). Similar results were obtained in all subgroups.

CONCLUSIONS

Compared to previous studies, our findings present a clear distinction. Therefore, there may be racial differences in the association between height and knee OA risk, and our study provides evidence that, in East Asian populations, taller individuals have a reduced risk of knee OA.

A systematic review of genome-wide association studies for pain, nociception, neuropathy, and pain treatment responses

ABSTRACT

Pain is the leading cause of disability worldwide, imposing an enormous burden on personal health and society. Pain is a multifactorial and multidimensional problem. Currently, there is (some) evidence that genetic factors could partially explain individual susceptibility to pain and interpersonal differences in pain treatment response. To better understand the underlying genetic mechanisms of pain, we systematically reviewed and summarized genome-wide association studies (GWASes) investigating the associations between genetic variants and pain/pain-related phenotypes in humans. We reviewed 57 full-text articles and identified 30 loci reported in more than 1 study. To check whether genes described in this review are associated with (other) pain phenotypes, we searched 2 pain genetic databases, Human Pain Genetics Database and Mouse Pain Genetics Database. Six GWAS-identified genes/loci were also reported in those databases, mainly involved in neurological functions and inflammation. These findings demonstrate an important contribution of genetic factors to the risk of pain and pain-related phenotypes. However, replication studies with consistent phenotype definitions and sufficient statistical power are required to validate these pain-associated genes further. Our review also highlights the need for bioinformatic tools to elucidate the function of identified genes/loci. We believe that a better understanding of the genetic background of pain will shed light on the underlying biological mechanisms of pain and benefit patients by improving the clinical management of pain.

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.

R399E, A Mutated Form of Growth and Differentiation Factor 5, for Disease Modification of Osteoarthritis

OBJECTIVE

To preclinically characterize a mutant form of growth and differentiation factor 5, R399E, with reduced osteogenic properties as a potential disease-modifying osteoarthritis (OA) drug.

METHODS

Cartilage, synovium, and meniscus samples from patients with OA were used to evaluate anabolic and antiinflammatory properties of R399E. In the rabbit joint instability model, 65 rabbits underwent transection of the anterior cruciate ligament plus partial meniscectomy. Three intraarticular (IA) R399E doses were administered biweekly 6 times, and static incapacitance was determined to assess joint pain. OA was evaluated 13 weeks after surgery. In sheep, medial meniscus transection was performed to induce OA, dynamic weight bearing was measured in-life, and OA was assessed after 13 weeks.

RESULTS

Intermittent exposure to R399E (1 week per month) was sufficient to induce cell proliferation and release of anabolic markers in 3-dimensional chondrocyte cultures. R399E also inhibited the release of interleukin-1β (IL-1β), IL-6, and prostaglandin E₂ from cartilage with synovium, meniscal cell, and synoviocyte cultures. In rabbits, the mean difference (95% confidence interval [95% CI]) in weight bearing for R399E compared to vehicle was -5.8 (95% confidence interval [95% CI] -9.54, -2.15), -7.2 (95% CI -10.93, -3.54), and -7.7 (95% CI -11.49, -3.84) for the 0.6, 6, and 60 μg doses, respectively, 6 hours after the first IA injection, and was statistically significant through the entire study for all doses. Cartilage surface structure improved with the 6-μg dose. Structural and symptomatic improvement with the same dose was confirmed in the sheep model of OA.

CONCLUSION

R399E influences several pathologic processes contributing to OA, highlighting its potential as a disease-modifying therapy.

Phenotype-Genotype analysis of caucasian patients with high risk of osteoarthritis

Background: Osteoarthritis (OA) is a common cause of disability and pain around the world. Epidemiologic studies of family history have revealed evidence of genetic influence on OA. Although many efforts have been devoted to exploring genetic biomarkers, the mechanism behind this complex disease remains unclear. The identified genetic risk variants only explain a small proportion of the disease phenotype. Traditional genome-wide association study (GWAS) focuses on radiographic evidence of OA and excludes sex chromosome information in the analysis. However, gender differences in OA are multifactorial, with a higher frequency in women, indicating that the chromosome X plays an essential role in OA pathology. Furthermore, the prevalence of comorbidities among patients with OA is high, indicating multiple diseases share a similar genetic susceptibility to OA. Methods: In this study, we performed GWAS of OA and OA-associated key comorbidities on 3366 OA patient data obtained from the Osteoarthritis Initiative (OAI). We performed Mendelian randomization to identify the possible causal relationship between OA and OA-related clinical features. Results: One significant OA-associated locus rs2305570 was identified through sex-specific genome-wide association. By calculating the LD score, we found OA is positively correlated with heart disease and stroke. A strong genetic correlation was observed between knee OA and inflammatory disease, including eczema, multiple sclerosis, and Crohn's disease. Our study also found that knee alignment is one of the major risk factors in OA development, and we surprisingly found knee pain is not a causative factor of OA, although it was the most common symptom of OA. Conclusion: We investigated several significant positive/negative genetic correlations between OA and common chronic diseases, suggesting substantial genetic overlaps between OA and these traits. The sex-specific association analysis supports the critical role of chromosome X in OA development in females.

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.

Epigenetic Regulation of Chondrocytes and Subchondral Bone in Osteoarthritis

The aim of this review is to provide an updated review of the epigenetic factors involved in the onset and development of osteoarthritis (OA). OA is a prevalent degenerative joint disease characterized by chronic inflammation, ectopic bone formation within the joint, and physical and proteolytic cartilage degradation which result in chronic pain and loss of mobility. At present, no disease-modifying therapeutics exist for the prevention or treatment of the disease. Research has identified several OA risk factors including mechanical stressors, physical activity, obesity, traumatic joint injury, genetic predisposition, and age. Recently, there has been increased interest in identifying epigenetic factors involved in the pathogenesis of OA. In this review, we detail several of these epigenetic modifications with known functions in the onset and progression of the disease. We also review current therapeutics targeting aberrant epigenetic regulation as potential options for preventive or therapeutic treatment.

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.

The effect of common variants in GDF5 gene on the susceptibility to chronic postsurgical pain

Background

The growth differentiation factor 5 (GDF5) gene regulates the growth of neuronal axons and dendrites and plays a role in the inflammatory response and tissue damage. The gene may also be associated with chronic postsurgical pain. This study aimed to reveal the relationship between SNPs in the GDF5 gene and orthopedic chronic postsurgical pain in Han Chinese population based on a case-control study.

Methods

We genotyped 8 SNPs within GDF5 gene in 1048 surgical patients with chronic postsurgical pain as the case group and 2062 surgical patients who were pain free as the control group. SNP and haplotypic analyses were performed, and stratified analyses were conducted to determine the correlations between significant SNPs and clinical characteristics.

Results

Only rs143384 in the 5′UTR of GDF5 was identified as significantly associated with increased susceptibility to chronic postsurgical pain, and the risk of A allele carriers was increased approximately 1.35-fold compared with that of G allele carriers. Haplotypes AGG and GGG in the LD block rs143384-rs224335-rs739329 also showed similar association patterns. Furthermore, we found that rs143384 was significantly correlated with chronic postsurgical pain in the subgroup aged ≤ 61 years, subgroup with a BMI ≤ 26, subgroup with no-smoking or no pain history, and subgroup with a drinking history.

Conclusion

Our study provided supportive evidence that genetic variations in the GDF5 gene are potential genetic factors that can increase the risk of chronic postsurgical pain in the Han Chinese population, but further research is necessary to elucidate the underlying mechanism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-021-02549-5.

Frequency of Growth Differentiation Factor 5 rs143383 and asporin D-repeat polymorphisms in patients with hand and knee osteoarthritis in Kurdistan province, Iran

AIM

Osteoarthritis (OA) is the most common chronic joint disorder, resulting from the breakdown of joint cartilage. It occurs in the knees, hands, and hips, leading to pain, stiffness, inflammation, and swelling.

METHODS

In this study, 100 hand and knee OA patients, meeting the American College of Rheumatology criteria were included in the case group, and 100 healthy individuals were allocated to the control group. Blood samples were collected from the participants. After DNA extraction, genotyping was carried out for GDF5 rs143383 C/T polymorphism by allele-specific polymerase chain reaction (ASPCR) and for D-repeat alleles of asporin (ASPN) by conventional PCR assay.

RESULTS

The results showed that the frequency of the D14 allele of ASPN was significantly higher than other alleles in the case group (P = .0001). Also, the frequency of the D14 allele among women was significantly higher than in men (P = .004). Moreover, the frequency of the TT allele in GDF5 rs143383 C/T polymorphism was significantly higher than the CC and CT alleles in the case group, compared with the control group (P = .001). A significant difference was found between the TT allele and other alleles in female and male patients compared with the control group (P = .02).

CONCLUSIONS

The D14 allele of the ASPN gene and TT allele of the GDF5 gene (rs143383 + 104T/C) are associated with hand and knee OA in the Kurdish population, indicating that these alleles could be risk factors for OA, at least in our populations.

Constitutional morphological features and risk of hip osteoarthritis: a case-control study using standard radiographs

OBJECTIVES

To evaluate the risk of association with hip osteoarthritis (OA) of 14 morphological features measured on standard antero-posterior pelvis radiographs.

METHODS

A case-control study of 566 symptomatic unilateral hip OA cases and 1108 controls without hip OA, using the Genetics of OA and Lifestyle database. Unaffected hips of cases were assumed to reflect pre-OA morphology of the contralateral affected hip. ORs with 95% CI adjusted for confounding factors were calculated using logistic regression. Hierarchical clustering on principal component method was used to identify clusters of morphological features. Proportional risk contribution (PRC) of these morphological features in the context of other risk factors of hip OA was estimated using receiver operating characteristic analysis.

RESULTS

All morphological features showed right-left symmetry in controls. Each feature was associated with hip OA after adjusting for age, gender and body mass index. Increased sourcil angle had the strongest association (OR: 6.93, 95% CI 5.16 to 9.32). Three clusters were identified. The PRC varied between individual features, as well as between clusters. It was 35% (95% CI 31% to 40%) for all 14 morphological features, compared to 21% (95% CI 19% to 24%) for all other well-established risk factors.

CONCLUSIONS

Constitutional morphological variation strongly associates with hip OA development and may explain much of its heritability. Relevant morphological measures can be assessed readily on standard radiographs to help predict risk of hip OA. Prospective studies are required to provide further support for causality.

Exogenous Application of Proteoglycan to the Cell Surface Microenvironment Facilitates to Chondrogenic Differentiation and Maintenance

Osteoarthritis (OA), a disease that greatly impacts quality of life, has increasing worldwide prevalence as the population ages. However, its pathogenic mechanisms have not been fully elucidated and current therapeutic treatment strategies are inadequate. In recent years, abnormal endochondral ossification in articular cartilage has received attention as a pathophysiological mechanism in OA. Cartilage is composed of abundant extracellular matrix components, which are involved in tissue maintenance and regeneration, but how these factors affect endochondral ossification is not clear. Here, we show that the application of aggrecan-type proteoglycan from salmon nasal cartilage (sPG) exhibited marked proliferative capacity through receptor tyrosine kinases in chondroprogenitor cells, and also exhibited differentiation and three-dimensional structure formation via phosphorylation of Insulin-like Growth Factor-1 Receptor and Growth Differentiation Factor 5 expression. Furthermore, sPG inhibited calcification via expression of Runx2 and Col10 (factors related to induction of calcification), while increasing Mgp, a mineralization inhibitory factor. As a result of analyzing the localization of sPG applied to the cells, it was localized on the surface of the cell membrane. In this study, we found that sPG, as a biomaterial, could regulate cell proliferation, differentiation and calcification inhibition by acting on the cell surface microenvironment. Therefore, sPG may be the foundation for a novel therapeutic approach for cartilage maintenance and for improved symptoms in OA.

Identification of TGFβ signatures in six murine models mimicking different osteoarthritis clinical phenotypes

OBJECTIVE

TGFβ is a key player in cartilage homeostasis and OA pathology. However, few data are available on the role of TGFβ signalling in the different OA phenotypes. Here, we analysed the TGFβ pathway by transcriptomic analysis in six mouse models of OA.

METHOD

We have brought together seven expert laboratories in OA pathophysiology and, used inter-laboratories standard operating procedures and quality controls to increase experimental reproducibility and decrease bias. As none of the available OA models covers the complexity and heterogeneity of the human disease, we used six different murine models of knee OA: from post-traumatic/mechanical models (meniscectomy (MNX), MNX and hypergravity (HG-MNX), MNX and high fat diet (HF-MNX), MNX and seipin knock-out (SP-MNX)) to aging-related OA and inflammatory OA (collagenase-induced OA (CIOA)). Four controls (MNX-sham, young, SP-sham, CIOA-sham) were added. OsteoArthritis Research Society International (OARSI)-based scoring of femoral condyles and ribonucleic acid (RNA) extraction from tibial plateau samples were done by single operators as well as the transcriptomic analysis of the TGFβ family pathway by Custom TaqMan® Array Microfluidic Cards.

RESULTS

The transcriptomic analysis revealed specific gene signatures in each of the six models; however, no gene was deregulated in all six OA models. Of interest, we found that the combinatorial Gdf5-Cd36-Ltbp4 signature might discriminate distinct subgroups of OA: Cd36 upregulation is a hallmark of MNX-related OA while Gdf5 and Ltbp4 upregulation is related to MNX-induced OA and CIOA.

CONCLUSION

These findings stress the OA animal model heterogeneity and the need of caution when extrapolating results from one model to another.

Association between growth differentiation factor 5 rs143383 genetic polymorphism and the risk of knee osteoarthritis among Caucasian but not Asian: a meta-analysis

BACKGROUND

A few months ago, the Bioscience Reports journal showed that growth differentiation factor 5 (GDF5) rs143383 genetic polymorphism increases the susceptibility of knee osteoarthritis (KOA), but previous studies' results have debates about available data. Considering the availability of more recent data, we focus on clarifying the relationship of KOA and GDF5 rs143383 genetic polymorphism by a meta-analysis of case-control trial data.

METHODS

The eligible studies from the time of database established to Oct. 2019 were collected from PubMed, Springer, Cochrane library, Web of Science, China National Knowledge Infrastructure (CNKI), and Wan Fang library. Odds ratios (OR) and 95% confidence intervals (CI) were used to estimate the association between these polymorphisms and KOA risk. The meta-analysis was completed by STATA 18.0 software.

RESULTS

A total of 196 studies were collected, 16 of them included in final meta-analysis (7997 cases and 12,684 controls). There was significant association between GDF5 rs143383 polymorphism and KOA in all genetic models (for Allele model (C versus T): OR = 0.84 (95% CI = 0.76-0.91); dominate model (CC+CT versus TT): OR = 0.80 (95% CI = 0.72-0.90); recessive model (CC versus CT+TT): OR = 0.79 (95% CI = 0.68-0.92); heterozygote model (CT versus CC+TT): OR = 0.89 (95% CI = 0.80-0.97); homozygous model (CC versus TT): OR = 0.71 (95% CI = 0.60-0.85)). In the subgroup analysis, we obtained the results that there is no significance among Asians.

CONCLUSION

GDF5 rs143383 genetic polymorphism increases the risk of KOA among Caucasians; CC genotype and C allele are protective factors for the susceptibility of KOA among Caucasians.

Donor genetic backgrounds contribute to the functional heterogeneity of stem cells and clinical outcomes

Stable and sustainable stem cell sources for stem cell‐based therapies are scarce and a key bottleneck for clinical applications. The regenerative potential of stem cells is usually attributed to several allogeneic or even autologous donor‐related factors. Genetic background and epigenetic variations in different individuals may significantly affect the functional heterogeneity of stem cells. Particularly, single‐nucleotide polymorphisms (SNPs) have been implicated in diseases with monogenetic or multifactorial and complex genetic etiologies. However, the possible effects of individual SNPs on donor stem cells remain far from fully elucidated. In this Perspective, we will discuss the roles played by donor genetic traits in the functional heterogeneity of induced pluripotent stem cells, mesenchymal stem cells, and hematopoietic stem cells and their implications for regenerative medicine and therapy.

Epigenetic Therapies for Osteoarthritis

Osteoarthritis (OA) is an age-associated disease characterized by chronic joint pain resulting from degradation of articular cartilage, inflammation of the synovial lining, and changes to the subchondral bone. Despite the wide prevalence, no FDA-approved disease-modifying drugs exist. Recent evidence has demonstrated that epigenetic dysregulation of multiple molecular pathways underlies OA pathogenesis, providing a new mechanistic and therapeutic axis with the advantage of targeting multiple deregulated pathways simultaneously. In this review, we focus on the epigenetic regulators that have been implicated in OA, their individual roles, and potential crosstalk. Finally, we discuss the pharmacological molecules that can modulate their activities and discuss the potential advantages and challenges associated with epigenome-based therapeutics for OA.

The GDF-5 mutant M1673 exerts robust anabolic and anti-catabolic effects in chondrocytes

The growth and differentiation factor 5 (GDF‐5) is known to play a key role in cartilage morphogenesis and homeostasis, and a single‐nucleotide polymorphism in its promoter sequence was found to be associated with osteoarthritis (OA). In addition, GDF‐5 was shown to promote extracellular matrix (ECM) production in healthy chondrocytes, to stimulate chondrogenesis of mesenchymal stem cells (MSCs) and to protect against OA progression in vivo. Therefore, GDF‐5 appears to be a promising treatment for osteoarthritis. However, GDF‐5 also promotes osteogenesis and hypertrophy, limiting its therapeutic utility. To circumvent this, a GDF‐5 mutant with lower hypertrophic and osteogenic properties was engineered: M1673. The present study aimed to evaluate and compare the effects of GDF‐5 and M1673 on primary porcine and human OA chondrocytes. We found that both GDF‐5 and M1673 can robustly stimulate ECM accumulation, type II collagen and aggrecan expression in porcine and human OA chondrocytes in 3D culture. In addition, both molecules also down‐regulated MMP13 and ADAMTS5 expression. These results suggest that M1673 retained the anabolic and anti‐catabolic effects of GDF‐5 on chondrocytes and is an alternative to GDF‐5 for osteoarthritis.

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.

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.

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.

Inhibition of 4-aminobutyrate aminotransferase protects against injury-induced osteoarthritis in mice

Recently we demonstrated that ablation of the DNA methyltransferase enzyme, Dnmt3b, resulted in catabolism and progression of osteoarthritis (OA) in murine articular cartilage through a mechanism involving increased mitochondrial respiration. In this study, we identify 4-aminobutyrate aminotransferase (Abat) as a downstream target of Dnmt3b. Abat is an enzyme that metabolizes γ-aminobutyric acid to succinate, a key intermediate in the tricarboxylic acid cycle. We show that Dnmt3b binds to the Abat promoter, increases methylation of a conserved CpG sequence just upstream of the transcriptional start site, and inhibits Abat expression. Dnmt3b deletion in articular chondrocytes results in reduced methylation of the CpG sequence in the Abat promoter, which subsequently increases expression of Abat. Increased Abat expression in chondrocytes leads to enhanced mitochondrial respiration and elevated expression of catabolic genes. Overexpression of Abat in murine knee joints via lentiviral injection results in accelerated cartilage degradation following surgical induction of OA. In contrast, lentiviral-based knockdown of Abat attenuates the expression of IL-1β-induced catabolic genes in primary murine articular chondrocytes in vitro and also protects against murine articular cartilage degradation in vivo. Strikingly, treatment with the FDA-approved small-molecule Abat inhibitor, vigabatrin, significantly prevents the development of injury-induced OA in mice. In summary, these studies establish Abat as an important new target for therapies to prevent OA.

Genome-wide association study of knee pain identifies associations with GDF5 and COL27A1 in UK Biobank

Knee pain is one of the most common musculoskeletal complaints that brings people to medical attention. Approximately 50% of individuals over the age of 50 report an experience of knee pain within the past 12 months. We sought to identify the genetic variants associated with knee pain in 171,516 subjects from the UK Biobank cohort and seek supporting evidence in cohorts from 23andMe, the Osteoarthritis Initiative, and the Johnston County Osteoarthritis Project. We identified two loci that reached genome-wide significance in the UK Biobank: rs143384, located in GDF5 (P = 1.32 × 10-12), a gene previously implicated in osteoarthritis; and rs2808772, located near COL27A1 (P = 1.49 × 10-8). These findings were supported in cohorts with self-reported osteoarthritis/radiographic knee osteoarthritis without pain information. In this report on genome-wide association of knee pain, we identified two loci in or near GDF5 and COL27A1 that are associated with knee pain.

Integrating transcriptome-wide study and mRNA expression profiles yields novel insights into the biological mechanism of chondropathies

Background

Chondropathies are a group of cartilage diseases, which share some common pathogenetic features. The etiology of chondropathies is still largely obscure now.

Methods

A transcriptome-wide association study (TWAS) was performed using the UK Biobank genome-wide association study (GWAS) data of chondropathies (including 1314 chondropathy patients and 450,950 controls) with gene expression references of muscle skeleton (MS) and peripheral blood (YBL). The candidate genes identified by TWAS were further compared with three gene expression profiles of osteoarthritis (OA), cartilage tumor (CT), and spinal disc herniation (SDH), to confirm the functional relevance between the chondropathies and the candidate genes identified by TWAS. Functional mapping and annotation (FUMA) was used for the gene ontology enrichment analyses. Immunohistochemistry (IHC) was conducted to validate the accuracy of integrative analysis results.

Results

Integrating TWAS and mRNA expression profiles detected 84 candidate genes for knee OA, such as DDX20 (P_{TWAS YBL} = 1.79 × 10^− 3, fold change (FC) = 2.69), 10 candidate genes for CT, such as SRGN (P_{TWAS YBL} = 1.46 × 10^− 3, FC = 3.36), and 4 candidate genes for SDH, such as SUPV3L1 (P_{TWAS YBL} = 3.59 × 10^− 3, FC = 3.22). Gene set enrichment analysis detected 73 GO terms for knee OA, 3 GO terms for CT, and 1 GO term for SDH, such as mitochondrial protein complex (P = 7.31 × 10^− 5) for knee OA, cytokine for CT (P = 1.13 × 10^− 4), and ion binding for SDH (P = 3.55 × 10^− 4). IHC confirmed that the protein expression level of DDX20 was significantly different between knee OA cartilage and healthy control cartilage (P = 0.0358).

Conclusions

Multiple candidate genes and GO terms were detected for chondropathies. Our findings may provide a novel insight in the molecular mechanisms of chondropathies.

Electronic supplementary material

The online version of this article (10.1186/s13075-019-1978-8) contains supplementary material, which is available to authorized users.

Genetic Association Between Growth Differentiation Factor 5 Single Nucleotide Polymorphism and Primary Knee Osteoarthritis in a Group of Egyptian Patients: A Pilot Study

Aim

This study aimed to determine the genetic association between Growth Differentiation Factor 5 (GDF5) gene (rs143383 T/C) single nucleotide polymorphism (SNP) and primary knee osteoarthritis (OA) in a group of Egyptian patients.

Patients and Methods

The study included 47 patients with primary knee OA and 40 apparently healthy control subjects. The disease was assessed using Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) score and Health Assessment Questionnaire (HAQ). Radiological assessment was done by Kellgren-Laurence (K/L) grading system. The genetic association of the SNP with primary knee OA was assessed by restriction fragment length polymorphism - polymerase chain reaction (RFLP-PCR).

Results

The mean total WOMAC index was significantly higher in patients with TT genotype as compared to patients with CC and CT genotypes (P<0.001). Similarly, the HAQ score was significantly higher among patients with TT genotype when compared to patients with CT and CC genotypes (P<0.001). There was a statistically significant association between different GDF5 genotypes and K/L radiological grading of knee OA among the studied patients (P=0.029). No statistically significant association was detected on comparing the frequency distribution of GDF5 alleles and genotypes frequencies of the SNP in patients and healthy controls.

Conclusion

There is a possible genetic association between GDF5 (rs143383) SNP and severity of primary knee OA, which might facilitate the detection of patients with high risk for disease progression. The present study did not detect an association between the SNP and development of primary knee OA.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Nutritional, metabolic and genetic considerations to optimise regenerative medicine outcome for knee osteoarthritis

Knee osteoarthritis (KOA) is a multifactorial degenerative disorder of joints, affecting the world's population over the age of 65 and with a higher prevalence in females. KOA is responsible for many age associated joint problems such as stiffness and pain. Conventional methods for managing KOA such as nonsteroidal anti-inflammatory drugs (NSAID) may not improve pain or alter the disease progression and may have adverse side effects. Non-pharmacological management of OA is fundamental to management of functional limitations and provides effective symptom relief but has not shown that disease progression can be altered. Regenerative medicine is a relatively new approach which aims to induce cellular regeneration and promote self-healing through minimally invasive methods. The use of regenerative medicine slowed the progression of KOA and revealed significant improvements, yet further investigations are required to optimize the outcomes. Nutritional and metabolic aspects such as supplementations, vitamins and minerals were proven to have an impact on the progression of KOA. Genetic variations are rapidly inspected to identify any potential influence of these variations in the predisposition and diagnosis of KOA. Further supporting evidence suggests the potential influence of metabolic, nutritional and genetic aspects in optimizing the outcomes of regenerative medicine in the management of KOA.

Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis

Osteoarthritis has a highly negative impact on quality of life because of the associated pain and loss of joint function. Here we describe the largest meta-analysis so far of osteoarthritis of the hip and the knee in samples from Iceland and the UK Biobank (including 17,151 hip osteoarthritis patients, 23,877 knee osteoarthritis patients, and more than 562,000 controls). We found 23 independent associations at 22 loci in the additive meta-analyses, of which 16 of the loci were novel: 12 for hip and 4 for knee osteoarthritis. Two associations are between rare or low-frequency missense variants and hip osteoarthritis, affecting the genes SMO (rs143083812, frequency 0.11%, odds ratio (OR) = 2.8, P = 7.9 × 10^-12, p.Arg173Cys) and IL11 (rs4252548, frequency 2.08%, OR = 1.30, P = 2.1 × 10^-11, p.Arg112His). A common missense variant in the COL11A1 gene also associates with hip osteoarthritis (rs3753841, frequency 61%, P = 5.2 × 10^-10, OR = 1.08, p.Pro1284Leu). In addition, using a recessive model, we confirm an association between hip osteoarthritis and a variant of CHADL¹ (rs117018441, P = 1.8 × 10^-25, OR = 5.9). Furthermore, we observe a complex relationship between height and risk of osteoarthritis.

Association of GDF-5 rs143383 polymorphism with radiographic defined knee osteoarthritis: A systematic review and meta-analysis

Objective

To assess the association of GDF-5 rs143383 polymorphism with radiographic defined knee osteoarthritis (OA), a systematic review and meta-analysis was conducted.

Methods

A total of 17 relevant case-control studies with 7424 cases and 11,310 controls was collected from several electronic databases up to June 2018.

Results

The pooled results showed that GDF-5 rs143383 polymorphism was significantly associated with radiographic defined knee OA in overall and stratified analysis by ethnicity, source of controls and genotyping techniques.

Conclusions

The GDF-5 rs143383 polymorphism might be used as a relevant risk estimate for radiographic defined Knee OA.

Lin28a overexpression reveals the role of Erk signaling in articular cartilage development

Adult articular cartilage shows limited tissue turnover, and therefore development of the proper structure of articular cartilage is crucial for life-long joint function. However, the mechanism by which the articular cartilage structure is developmentally regulated is poorly understood. In this study, we show evidence that activation of extracellular signal-regulated kinases (Erk1/2) in articular chondrocyte progenitors during developmental stages control articular cartilage thickness. We found that overexpression of Lin28a, an RNA-binding protein that regulates organismal growth and metabolism, in articular chondrocyte progenitor cells upregulated Erk signaling and increased articular cartilage thickness. Overexpression of a constitutively active Kras mimicked Lin28a overexpression, and inhibition of Erk signaling during embryonic stages normalized the cartilage phenotype of both Kras- and Lin28a-overexpressing mice. These results suggest that articular cartilage thickness is mainly determined during the process of embryonic synovial joint development, which is positively regulated by Erk signaling.

A machine-learned analysis of human gene polymorphisms modulating persisting pain points to major roles of neuroimmune processes

Background

Human genetic research has implicated functional variants of more than one hundred genes in the modulation of persisting pain. Artificial intelligence and machine‐learning techniques may combine this knowledge with results of genetic research gathered in any context, which permits the identification of the key biological processes involved in chronic sensitization to pain.

Methods

Based on published evidence, a set of 110 genes carrying variants reported to be associated with modulation of the clinical phenotype of persisting pain in eight different clinical settings was submitted to unsupervised machine‐learning aimed at functional clustering. Subsequently, a mathematically supported subset of genes, comprising those most consistently involved in persisting pain, was analysed by means of computational functional genomics in the Gene Ontology knowledgebase.

Results

Clustering of genes with evidence for a modulation of persisting pain elucidated a functionally heterogeneous set. The situation cleared when the focus was narrowed to a genetic modulation consistently observed throughout several clinical settings. On this basis, two groups of biological processes, the immune system and nitric oxide signalling, emerged as major players in sensitization to persisting pain, which is biologically highly plausible and in agreement with other lines of pain research.

Conclusions

The present computational functional genomics‐based approach provided a computational systems‐biology perspective on chronic sensitization to pain. Human genetic control of persisting pain points to the immune system as a source of potential future targets for drugs directed against persisting pain. Contemporary machine‐learned methods provide innovative approaches to knowledge discovery from previous evidence.

Significance

We show that knowledge discovery in genetic databases and contemporary machine‐learned techniques can identify relevant biological processes involved in Persitent pain.

The Involvement of MicroRNAs in Osteoarthritis and Recent Developments: A Narrative Review

Background:

Osteoarthritis (OA) is the most common chronic joint disease and it may progressively cause disability and compromise quality of life. Lately, the role of miRNAs in the pathogenesis of OA has drawn a lot of attention. miRNAs are small, single-stranded, non-coding molecules of RNA which regulate gene expression at post-transcriptional level. The dysregulation of the expression of several miRNAs affects pathways involved in OA pathogenesis.

Objective:

The purpose of this article is to review the literature on the involvement of miRNAs in the pathogenesis of OA and the implications on its diagnosis and treatment.

Materials and Methods:

An extensive electronic literature search was conducted by two researchers from January 2008 to August 2017. Titles and abstracts of papers were screened by the authors for further inclusion in the present work. Finally, full texts of the selected articles were retrieved.

Results:

Abnormally expressed miRNAs enhance the production of cartilage degrading enzymes, inhibit the expression of cartilage matrix components, increase the production of proinflammatory cytokines, facilitate chondrocyte apoptosis, suppress autophagy in chondrocytes and are involved in pain-related pathways. miRNAs are also incorporated in extra-cellular membranous vesicles such as exosomes and participate in the intercellular communication in osteoarthritic joints.

Conclusion:

Ongoing research on miRNAs has potential implications in the diagnosis and treatment of OA. Their different levels in peripheral blood and synovial fluid between OA patients and healthy population makes them candidates for being used as biomarkers of the disease, while targeting miRNAs may be a novel therapeutic strategy in OA.

Association between EN1 rs4144782 and susceptibility of knee osteoarthritis: A case-control study

Osteoarthritis (OA) is a complex disease that affects the whole joint, resulting from the combined influence of biomechanical factors and genetic factors. The heritable component for primary OA accounts for about 60% of variation in population liability to the disease. So far, genome-wide association studies (GWAS) and candidate gene studies have established many OA-related loci. However, these findings account for only a rather small fraction of the genetic component. To further reveal the genetic architecture of OA, we conducted this case-control study to explore the association of locus EN1 rs4144782 and knee OA susceptibility in a Chinese population. EN1 rs4144782 was significantly associated with increased risk of knee OA (OR=1.26; 95% CI: 1.05-1.50, P value=0.012). In dominant model, compared with carriers of GG genotype, those with AG or AA genotype have an 1.44-fold increased risk of OA (OR: 1.44; 95% CI: 1.10-1.88; P value=0.008). Subgroup analyses didn't change the results materially. This should be the first association study of EN1 locus on risk of OA, and our finding suggested that the EN1 rs4144782 might contribute to the susceptibility of knee OA.

Whole genome sequencing of pairwise human subjects reveals DNA mutations specific to developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is a common congenital malformation characterized by mismatch in shape between the femoral head and acetabulum, and leads to hip dysplasia. To date, the pathogenesis of DDH is poorly understood and may involve multiple factors, including genetic predisposition. However, comprehensive genetic analysis has not been applied to investigate a genetic component of DDH. In the present study, 10 pairs of healthy fathers and DDH daughters were enrolled to identify genetic hallmarks of DDH using high throughput whole genome sequencing. The DDH-specific DNA mutations were found in each patient. Overall 1344 genes contained DDH-specific mutations. Functional enrichment analysis showed that these genes played important roles in the cytoskeleton, microtubule cytoskeleton, sarcoplasm and microtubule associated complex. These functions affected osteoblast and osteoclast development. Therefore, we proposed that the DDH-specific mutations might affect bone development, and caused DDH. Our pairwise high throughput sequencing results comprehensively delineated genetic hallmarks of DDH. Further research into the biological impact of these mutations may inform the development of DDH diagnostic tools and allow neonatal gene screening.

Functional polymorphisms in asporin and CILP together with joint loading predispose to hand osteoarthritis

Background

Osteoarthritis (OA) is the most common degenerative joint disease afflicting people in the Western world and has a strong genetic influence. The aim of this study was to examine the association of two known functional polymorphisms in the TGF-β inhibiting genes, asporin (ASPN) and cartilage intermediate layer protein (CILP), with hand OA and potential gene-occupational hand loading interaction.

Results

Statistically significant interaction of the CILP rs2073711 T and ASPN D15 alleles with hand OA was observed (OR = 2.48, 95% CI 1.27–4.85, p = 0.008) in a Finnish hand OA cohort of 543 women (aged 45–63). When stratified by variation in working tasks, low variation of working tasks increased the risk further (OR = 3.00, 95% CI 1.35–6.66, p = 0.007). Based on the analysis of ASPN and CILP protein-coding regions, functional studies were performed with one observed variant, rs41278695 in the ASPN gene. Analyses showed that bone morphogenetic protein 2 (BMP2) mediated expression of aggrecan (Agc1) and type II collagen (Col2a1) was significantly suppressed (p = 0.011 and p = 0.023, respectively) in a murine chondrocytic cell line (ATDC5) with cells stably expressing ASPN rs41278695.

Conclusions

The carriage of either ASPN D15 or CILP rs2073711 TT is associated with increased risk of symmetrical hand OA, particularly in individuals with low variation in work tasks. ASPN rs41278695 SNP had an effect on Agc1 and Col2a1 gene expression when induced with BMP-2 suggesting an effect on the cartilage extracellular matrix composition.

Genetic Determinants of Radiographic Knee Osteoarthritis in African Americans

OBJECTIVE

The etiology of knee osteoarthritis (OA), the most common form of arthritis, is complex and may differ by race or ethnicity. In recent years, genetic studies have identified many genetic variants associated with OA, but nearly all the studies were conducted in European whites and Asian Americans. Few studies have focused on the genetics of knee OA in African Americans.

METHODS

We performed a genome-wide association study of radiographic knee OA in 1217 African Americans from 2 North American cohort studies: 590 subjects from the Johnston County Osteoarthritis Project and 627 subjects from the Osteoarthritis Initiative. Analyses were conducted in each cohort separately and combined in an inverse variance fixed effects metaanalysis, which were then included in pathway analyses. We additionally tested 12 single-nucleotide polymorphisms robustly associated with OA in European white populations for association in African Americans.

RESULTS

We identified a genome-wide significant variant in LINC01006 (minor allele frequency 12%; p = 4.11 × 10-9) that is less common in European white populations (minor allele frequency < 3%). Five other independent loci reached suggestive significance (p < 1 × 10-6). In pathway analyses, dorsal/ventral neural tube patterning and iron ion transport pathways were significantly associated with knee OA in African Americans (false discovery rate < 0.05). We found no evidence that previously reported OA susceptibility variants in European whites were associated with knee OA in African Americans.

CONCLUSION

These results highlight differences in the genetic architecture of knee OA between African American and European whites. This finding underscores the need to include more diverse populations in OA genetics studies.

Cushioning the cartilage: a canonical Wnt restricting matter

Wnt signalling pathways have key roles in joint development, homeostasis and disease, particularly in osteoarthritis. New data is starting to reveal the importance of tightly regulating canonical Wnt signalling pathway activation to maintain homeostasis and health in articular cartilage. In addition to the presence of different Wnt antagonists that limit pathway activation in articular cartilage, the reciprocal crosstalk between the canonical and non-canonical cascades and competitive antagonism between different Wnt ligands seem to be critical in restraining excessive Wnt pathway activation. Changes in transcriptional complex assembly upon Wnt pathway activation, epigenetic modulation of target gene transcription, in particular through histone modifications, and complex interactions between the Wnt signalling pathway and other signalling pathways, are also instrumental in adjusting Wnt signalling. In this Review, the cellular and molecular mechanisms involved in fine-tuning canonical Wnt signalling in the joint are updated, with a focus on the articular cartilage. The interventions for preventing or treating osteoarthritis are also discussed, which should aim to limit disease-associated excessive canonical Wnt activity to avoid joint damage.

Association of BMP-14 rs143383 ploymorphism with its susceptibility to osteoarthritis: A meta-analysis and systematic review according to PRISMA guideline

BACKGROUND

Osteoarthritis (OA) is a complex disease which can be caused by both environmental and genetic factors. A functional locus rs143383 of bone morphogenetic protein-14 (BMP-14) has been pointed out to be associated with OA etiology, but conflicting conclusions have been reached. To provide a more comprehensive conclusion about this issue, we performed this meta-analysis.

METHODS

Relevant studies were searched from electronic databases including PubMed, Chinese National Knowledge Infrastructure, Embase, and Wanfang. The strength of correlations was examined with pooled odds ratios (ORs) and 95% confidence intervals (95% CIs). Subgroup analyses stratified by ethnicity and source of control were carried out. All statistical analyses were performed with STATA software (version 12.0).

RESULTS

Overall, BMP-14 rs143383 polymorphism was negatively correlated with the susceptibility to knee OA and hand OA under genetic contrasts of CC versus TT, CC + TC versus TT, CC versus TT + TC, C versus T, TC versus TT (OR = 0.71, 95% CI = 0.65-0.79; OR = 0.81, 95% CI = 0.73-0.89; OR = 0.79, 95% CI = 0.71-0.86; OR = 0.85, 95% CI = 0.81-0.90; OR = 0.84, 95% CI = 0.75-0.93), and TC versus TT, CC versus TT + TC, C versus T (OR = 0.76, 95% CI = 0.65-0.89; OR = 0.79, 95% CI = 0.68-0.92; OR = 0.90, 95% CI = 0.85-0.95), respectively; similar results were observed in subgroups after stratification analyses. Additionally, the polymorphism also reduced hip OA risk in Asian group after stratified analysis by ethnicity.

CONCLUSION

BMP-14 rs143383 polymorphism may be a protective factor against OA occurrence.

Association of ADAM12 gene polymorphisms with knee osteoarthritis susceptibility

Previous studies that evaluated the association between a disintegrin and metalloprotease 12 (ADAM12) gene polymorphisms and knee osteoarthritis (KOA) have given controversial and indefinite results. Therefore, we performed a meta-analysis to confirm this correlation. We searched the PubMed, Embase, and SinoMed databases for all papers published up to April 11, 2017. Overall, five different studies, totaling 2,353 cases and 3,668 controls, were retrieved on the basis of the search criteria for KOA susceptibility related to four polymorphisms (rs3740199, rs1278279, rs1871054, and rs1044122) in the ADAM12 gene. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the strength of this association. Publication bias was assessed using Egger's and Begg's tests. The rs3740199 G/C polymorphism was found to be associated with increased KOA risk in men (e.g., CG versus GG: OR = 1.44, 95% CI = 1.02–2.04, P = 0.040), but not in the overall analysis and in analyses of other subgroups. Significantly increased associations were also found for the rs1871054 polymorphism (e.g., C versus T allele: OR = 1.85, 95% CI = 1.49–2.30, P < 0.001). However, there were no associations for the rs1278279 and rs1044122 polymorphisms. Furthermore, no obvious evidence of publication bias was detected. Our study indicated that the rs1871054 polymorphism of ADAM12 was significantly associated with increased KOA risk.

Knee osteoarthritis has doubled in prevalence since the mid-20th century

Knee osteoarthritis (OA) is believed to be highly prevalent today because of recent increases in life expectancy and body mass index (BMI), but this assumption has not been tested using long-term historical or evolutionary data. We analyzed long-term trends in knee OA prevalence in the United States using cadaver-derived skeletons of people aged ≥50 y whose BMI at death was documented and who lived during the early industrial era (1800s to early 1900s; n = 1,581) and the modern postindustrial era (late 1900s to early 2000s; n = 819). Knee OA among individuals estimated to be ≥50 y old was also assessed in archeologically derived skeletons of prehistoric hunter-gatherers and early farmers (6000-300 B.P.; n = 176). OA was diagnosed based on the presence of eburnation (polish from bone-on-bone contact). Overall, knee OA prevalence was found to be 16% among the postindustrial sample but only 6% and 8% among the early industrial and prehistoric samples, respectively. After controlling for age, BMI, and other variables, knee OA prevalence was 2.1-fold higher (95% confidence interval, 1.5-3.1) in the postindustrial sample than in the early industrial sample. Our results indicate that increases in longevity and BMI are insufficient to explain the approximate doubling of knee OA prevalence that has occurred in the United States since the mid-20th century. Knee OA is thus more preventable than is commonly assumed, but prevention will require research on additional independent risk factors that either arose or have become amplified in the postindustrial era.

Genome-Wide Association Study of Radiographic Knee Osteoarthritis in North American Caucasians

OBJECTIVE

A major barrier to genetic studies of osteoarthritis (OA) is the need to obtain large numbers of individuals with standardized radiographic evaluations for OA. To address this gap, we performed a genome-wide association study (GWAS) of radiographically defined tibiofemoral knee OA in 3,898 cases and 3,168 controls from 4 well-characterized North American cohorts, and we performed replication analysis of previously reported OA loci.

METHODS

We performed meta-analysis using a 2-stage design. Stage 1 (discovery) consisted of a GWAS meta-analysis of radiographic knee OA carried out in the Osteoarthritis Initiative and the Johnston County Osteoarthritis Project. Knee OA was defined as definite osteophytes and possible joint space narrowing or total joint replacement in one or both knees. Stage 2 (validation) was performed in the Multicenter Osteoarthritis Study and the Genetics of Osteoarthritis study. We genotyped lead meta-analysis variants (P ≤ 1 × 10^-4 ) from stage 1 and tested the association between these variants and knee OA. We then combined results from all cohorts in a meta-analysis.

RESULTS

Lead variants from stage 1, representing 49 unique loci, were analyzed in stage 2; none met genome-wide significance in the combined analysis of stage 1 and stage 2. We validated 1 locus (rs4867568 near LSP1P3) with nominal significance (P < 0.05), which was also our top finding in the combined meta-analysis (odds ratio [OR] 0.84 [95% confidence interval (95% CI) 0.79-0.91], P = 3.02 × 10^-6 ). We observed nominally significant associations (P < 0.05) with 3 previously reported OA loci: rs143383 in GDF5 (OR 1.12 [95% CI 1.04-1.21], P = 2.13 × 10^-3 ), rs835487 in CHST11 (OR 0.93 [95% CI 0.85-0.99], P = 0.03), and rs8044769 in FTO (OR 1.10 [95% CI 1.03-1.19], P = 6.13 × 10^-3 ).

CONCLUSION

These findings provide suggestive evidence of a novel knee OA locus and confirm previously reported associations in GDF5, CHST11, and FTO.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis

OBJECTIVE

Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA.

METHODS

A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage.

RESULTS

We report suggestive evidence (p<5×10-6) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the LRCH1 gene (rs754106, p=1.49×10-7, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between STT3B and GADL1 (rs6766414, p=3.13×10-6, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). STT3B was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10-4).

CONCLUSIONS

Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

DNA methyltransferase 3b regulates articular cartilage homeostasis by altering metabolism

Osteoarthritis (OA) is the most common form of arthritis worldwide. It is a complex disease affecting the whole joint but is generally characterized by progressive degradation of articular cartilage. Recent genome-wide association screens have implicated distinct DNA methylation signatures in OA patients. We show that the de novo DNA methyltransferase (Dnmt) 3b, but not Dnmt3a, is present in healthy murine and human articular chondrocytes and its expression decreases in OA mouse models and in chondrocytes from human OA patients. Targeted deletion of Dnmt3b in murine articular chondrocytes results in an early-onset and progressive postnatal OA-like pathology. RNA-Seq and methylC-Seq analyses of Dnmt3b loss-of-function chondrocytes show that cellular metabolic processes are affected. Specifically, TCA metabolites and mitochondrial respiration are elevated. Importantly, a chondroprotective effect was found following Dnmt3b gain of function in murine articular chondrocytes in vitro and in vivo. This study shows that Dnmt3b plays a significant role in regulating postnatal articular cartilage homeostasis. Cellular pathways regulated by Dnmt3b in chondrocytes may provide novel targets for therapeutic approaches to treat OA.