A genome-wide association study identifies an osteoarthritis susceptibility locus on chromosome 7q22

The Role of Vitamin D in the Development and Progression of Osteoarthritis

"Low levels of vitamin D in patients with osteoarthritis (ΟΑ) rather adversely affect the structure and function of articular cartilage. Low levels of vitamin D are also associated with joint pain, limited physical activity, quality of life, as well as decreased muscle strength, primarily in the lower extremities with adverse OA progression. A high percentage of low vitamin D levels was found in patients with OA and joint arthroplasty as well as a less satisfactory postoperative follow-up in patients with low vitamin D levels. The administration of vitamin D supplementation in patients with knee osteoarthritis has been associated with improved articular cartilage architecture, and reduced joint pain, combined with improved functionality and quality of life in patients with OA. "

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.

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.

Erosive hand osteoarthritis: latest findings and outlook

Osteoarthritis (OA) most commonly affects knee joints, and the next most commonly affected sites are the hands and hips. Three distinct hand OA phenotypes have been described: erosive hand OA (EHOA), nodal hand OA - also known as non-erosive hand OA (non-EHOA) - and first carpometacarpal joint OA. EHOA predominantly affects women and is the most aggressive form of hand OA, characterized by a severe clinical onset and progression, leading to joint damage, disability and reduction of quality of life. Clinical signs of inflammation associated with EHOA include the acute onset of pain, swelling and redness. Moreover, EHOA is characterized by radiographic features such as central erosion, saw-tooth and gull-wing lesions and, rarely, ankylosis. The aim of this Review is to report the latest findings on epidemiology, clinical features, pathology and aetiopathogenesis, biomarkers, imaging modalities and treatments for EHOA. The ongoing development of new hand OA classification criteria should facilitate standardization between studies.

Nutraceutical Approach to Chronic Osteoarthritis: From Molecular Research to Clinical Evidence

Osteoarthritis (OA) is a degenerative inflammatory condition of the joint cartilage that currently affects approximately 58 million adults in the world. It is characterized by pain, stiffness, and a reduced range of motion with regard to the arthritic joints. These symptoms can cause in the long term a greater risk of overweight/obesity, diabetes mellitus, and falls and fractures. Although the current guidelines for the treatment of OA suggest, as the gold standard for this condition, pharmacological treatment characterized by non-steroidal anti-inflammatory drugs (NSAID), opioids, and cyclooxygenase (COX)-2-specific drugs, a great interest has been applied to nutraceutical supplements, which include a heterogeneous class of molecules with great potential to reduce inflammation, oxidative stress, pain, and joint stiffness and improve cartilage formation. The purpose of this review is to describe the potential application of nutraceuticals in OA, highlighting its molecular mechanisms of actions and data of efficacy and safety (when available).

G Protein-Coupled Receptors in Osteoarthritis: A Novel Perspective on Pathogenesis and Treatment

G protein-coupled receptors (GPCRs) are transmembrane receptor proteins that trigger numerous intracellular signaling pathways in response to the extracellular stimuli. The GPCRs superfamily contains enormous structural and functional diversity and mediates extensive biological processes. Until now, critical roles have been established in many diseases, including osteoarthritis (OA). Existing studies have shown that GPCRs play an important role in some OA-related pathogenesis, such as cartilage matrix degradation, synovitis, subchondral bone remodeling, and osteophyte formation. However, current pharmacological treatments are mostly symptomatic and there is a paucity of disease-modifying OA drugs so far. Targeting GPCRs is capable of inhibiting cartilage matrix degradation and synovitis and up-regulating cartilage matrix synthesis, providing a new therapeutic strategy for OA. In this review, we have comprehensively summarized the structures, biofunctions, and the novel roles of GPCRs in the pathogenesis and treatment of OA, which is expected to lay the foundation for the development of novel therapeutics against OA. Even though targeting GPCRs may ameliorate OA progression, many GPCRs-related therapeutic strategies are still in the pre-clinical stage and require further investigation.

Models of Osteoarthritis: Relevance and New Insights

Osteoarthritis (OA) is a progressive and disabling musculoskeletal disease affecting millions of people and resulting in major healthcare costs worldwide. It is the most common form of arthritis, characterised by degradation of the articular cartilage, formation of osteophytes, subchondral sclerosis, synovial inflammation and ultimate loss of joint function. Understanding the pathogenesis of OA and its multifactorial aetiology will lead to the development of effective treatments, which are currently lacking. Two-dimensional (2D) in vitro tissue models of OA allow affordable, high-throughput analysis and stringent control over specific variables. However, they are linear in fashion and are not representative of physiological conditions. Recent in vitro studies have adopted three-dimensional (3D) tissue models of OA, which retain the advantages of 2D models and are able to mimic physiological conditions, thereby allowing investigation of additional variables including interactions between the cells and their surrounding extracellular matrix. Numerous spontaneous and induced animal models are used to reproduce the onset and monitor the progression of OA based on the aetiology under investigation. This therefore allows elucidation of the pathogenesis of OA and will ultimately enable the development of novel and specific therapeutic interventions. This review summarises the current understanding of in vitro and in vivo OA models in the context of disease pathophysiology, classification and relevance, thus providing new insights and directions for OA research.

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.

Transplantation of human autologous synovial mesenchymal stem cells with trisomy 7 into the knee joint and 5 years of follow-up

Mesenchymal stem cells (MSCs) can show trisomy 7; however, the safety of these cells has not been fully investigated. The purposes of this study were to determine the ratio of patients whose synovial MSCs were transplanted clinically, to intensively investigate MSCs with trisomy 7 from a safety perspective, and to follow up the patients for 5 years after transplantation. Synovial MSCs at passage 0 were transplanted into a knee for degenerative meniscus tears in 10 patients, and the patients were checked at 5 years. The synovial MSCs were evaluated at passages 0 to 15 by G‐bands and digital karyotyping, and trisomy 7 was found in 3 of 10 patients. In those three patients, 5% to 10% of the synovial MSCs showed trisomy 7. The mRNA expressions of representative oncogenes and genes on chromosome 7 did not differ between MSCs with and without trisomy 7. Whole‐genome sequencing and DNA methylation analysis showed similar results for MSCs with and without trisomy 7. Transplantation of human synovial MSCs with trisomy 7 into eight mouse knees did not result in tumor formation under the skin or in the knees after 8 weeks in any mouse, whereas transplanted HT1080 cells formed tumors. In vitro chondrogenic potentials were similar between MSCs with and without trisomy 7. Five‐year follow‐ups revealed no serious adverse events in all 10 human patients, including 3 who had received MSCs with trisomy 7. Overall, our findings indicated that synovial MSCs with trisomy 7 were comparable with MSCs without trisomy 7 from a safety perspective.

Familial Clustering of Erosive Hand Osteoarthritis in a Large Statewide Cohort

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Accelerating functional gene discovery in osteoarthritis

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease.

Association between the MMP-1-1607 1G/2G Polymorphism and Osteoarthritis Risk: A Systematic Review and Meta-Analysis

Background

Osteoarthritis (OA) is a common disease characterized by articular cartilage degeneration and secondary hyperosteogenesis. Genetic factors are associated with the occurrence of OA. While several studies have shown that the matrix metalloproteinase-1- (MMP-1-) 1607 1G/2G (rs1799750) polymorphism may be related to the occurrence and development of OA, there is inconsistency in the literature. To better estimate the relationship between the MMP-1 gene polymorphism and OA, a comprehensive meta-analysis of relevant literature was carried out.

Results

In total, seven studies comprising 1245 OA patients and 1230 controls were included in this meta-analysis. The combined results revealed no significant association between the MMP-1-1607 1G/2G polymorphism and risk of OA in the five genetic models. However, after Bonferroni correction, the results of subgroup analysis revealed a significant correlation between the MMP-1-1607 1G/2G polymorphism and OA susceptibility in the temporomandibular joint (TMJ) OA subgroup (allelic: 2G vs. 1G: OR = 1.575, 95%CI = 1.259–1.972, P < 0.01; recessive: 2G2G vs. 1G1G+1G2G: OR = 2.411, 95%CI = 1.658–3.504, P < 0.01; and homozygote: 2G2G vs. 1G1G: OR = 2.313, 95%CI = 1.341, 3.991, P = 0.003), the younger subgroup (aged less than 60 years; allelic: 2G vs. 1G: OR = 1.635, 95%CI = 1.354, 1.974, P < 0.01; dominant: 2G1G+2G2G vs. 1G1G: OR = 1.622, 95%CI = 1.158, 2.271, P = 0.005; recessive: 2G2G vs. 1G1G+1G2G: OR = 2.209, 95%CI = 1.718, 2.840, P < 0.01; and homozygote: 2G2G vs. 1G1G: OR = 2.578, 95%CI = 1.798, 3.696, P < 0.01), the larger subgroup (N > 300), and the hospital-based case-control study (HCC) subgroup. The sensitivity analysis suggested that the results of the meta-analysis were stable and reliable. Begg's funnel plot and Egger's test indicated that there was no publication bias in this study.

Conclusion

Our meta-analysis indicated that although the MMP-1-1607 1G/2G polymorphism was not significantly associated with OA susceptibility among the whole sample, it played a key role in the etiology and development of TMJ OA and OA in people aged less than 60 years.

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.

Genetic Susceptibility Loci for Cardiovascular Disease and Their Impact on Atherosclerotic Plaques

Supplemental Digital Content is available in the text.

Background:

Atherosclerosis is a chronic inflammatory disease in part caused by lipid uptake in the vascular wall, but the exact underlying mechanisms leading to acute myocardial infarction and stroke remain poorly understood. Large consortia identified genetic susceptibility loci that associate with large artery ischemic stroke and coronary artery disease. However, deciphering their underlying mechanisms are challenging. Histological studies identified destabilizing characteristics in human atherosclerotic plaques that associate with clinical outcome. To what extent established susceptibility loci for large artery ischemic stroke and coronary artery disease relate to plaque characteristics is thus far unknown but may point to novel mechanisms.

Methods:

We studied the associations of 61 established cardiovascular risk loci with 7 histological plaque characteristics assessed in 1443 carotid plaque specimens from the Athero-Express Biobank Study. We also assessed if the genotyped cardiovascular risk loci impact the tissue-specific gene expression in 2 independent biobanks, Biobank of Karolinska Endarterectomy and Stockholm Atherosclerosis Gene Expression.

Results:

A total of 21 established risk variants (out of 61) nominally associated to a plaque characteristic. One variant (rs12539895, risk allele A) at 7q22 associated to a reduction of intraplaque fat, P=5.09×10⁻⁶ after correction for multiple testing. We further characterized this 7q22 Locus and show tissue-specific effects of rs12539895 on HBP1 expression in plaques and COG5 expression in whole blood and provide data from public resources showing an association with decreased LDL (low-density lipoprotein) and increase HDL (high-density lipoprotein) in the blood.

Conclusions:

Our study supports the view that cardiovascular susceptibility loci may exert their effect by influencing the atherosclerotic plaque characteristics.

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.

The role of GPCRs in bone diseases and dysfunctions. Bone Res. 2019;7:19. [PMID: 31646011 PMCID: PMC6804689 DOI: 10.1038/s41413-019-0059-6]

The superfamily of G protein-coupled receptors (GPCRs) contains immense structural and functional diversity and mediates a myriad of biological processes upon activation by various extracellular signals. Critical roles of GPCRs have been established in bone development, remodeling, and disease. Multiple human GPCR mutations impair bone development or metabolism, resulting in osteopathologies. Here we summarize the disease phenotypes and dysfunctions caused by GPCR gene mutations in humans as well as by deletion in animals. To date, 92 receptors (5 glutamate family, 67 rhodopsin family, 5 adhesion, 4 frizzled/taste2 family, 5 secretin family, and 6 other 7TM receptors) have been associated with bone diseases and dysfunctions (36 in humans and 72 in animals). By analyzing data from these 92 GPCRs, we found that mutation or deletion of different individual GPCRs could induce similar bone diseases or dysfunctions, and the same individual GPCR mutation or deletion could induce different bone diseases or dysfunctions in different populations or animal models. Data from human diseases or dysfunctions identified 19 genes whose mutation was associated with human BMD: 9 genes each for human height and osteoporosis; 4 genes each for human osteoarthritis (OA) and fracture risk; and 2 genes each for adolescent idiopathic scoliosis (AIS), periodontitis, osteosarcoma growth, and tooth development. Reports from gene knockout animals found 40 GPCRs whose deficiency reduced bone mass, while deficiency of 22 GPCRs increased bone mass and BMD; deficiency of 8 GPCRs reduced body length, while 5 mice had reduced femur size upon GPCR deletion. Furthermore, deficiency in 6 GPCRs induced osteoporosis; 4 induced osteoarthritis; 3 delayed fracture healing; 3 reduced arthritis severity; and reduced bone strength, increased bone strength, and increased cortical thickness were each observed in 2 GPCR-deficiency models. The ever-expanding number of GPCR mutation-associated diseases warrants accelerated molecular analysis, population studies, and investigation of phenotype correlation with SNPs to elucidate GPCR function in human diseases.

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

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

Role of SREBP2 gene polymorphism on knee osteoarthritis in the South Indian Hyderabad Population: A hospital based study with G595C variant

Introduction

Osteoarthritis (OA) is a multifactorial disease with genetic factors playing a crucial role, and it has been associated with a family history of obesity. G595C polymorphism in the sterol regulatory element-binding protein 2 (SREBP2) gene has demonstrated an association with knee osteoarthritis (KOA) patients. However, this polymorphism has been never explored in an Indian population. Hence, the current study aimed to examine whether G595C (rs2228314) polymorphism in SREBP2 gene was associated with KOA susceptibility in the South Indian Hyderabad population.

Methods

G595C polymorphism was genotyped with 200 KOA cases and 200 healthy controls using polymerase chain reaction-restriction fragment length polymorphism analysis.

Results

A significant association was observed between age, body mass index (BMI), and family histories in KOA cases and controls (p < 0.05). The current allele (C vs G; OR-2.8 [95%CI = 2.1-3.7]; p < 0.0001) and genotype analysis confirms the significant association with (GC + CC vs GG; OR-3.5 [95%CI = 2.3-5.3]; p < 0.0001 & GC vs GG + CC; OR-1.7 [95%CI = 1.0-2.9]; p = 0.02) KOA vs. control subjects. On stratification analysis, genotype CC and C allele were associated with KOA. Gender association failed to demonstrate positive genotype frequencies (p > 0.05). Multifactor-dimensionality reduction (MDR) analysis showed a positive association with BMI and G595C genotypes (p < 0.05); 51% of the homozygous variant CC genotypes were present in obesity subjects.

Conclusion

In conclusion, our findings suggest that G595C polymorphism in SREBP2 gene is associated with KOA in the South Indian Hyderabad population and presents scope for further investigation of the gene's function in KOA.

Association between ACE gene I/D polymorphism and knee osteoarthritis in a Chinese population

Osteoarthritis (OA) is a degenerative joint disease characterized by joint destruction with cartilage loss and occasional gross derangement of joint integrity. In recent years, several studies have reported the association between angiotensin-converting enzyme (ACE) gene insertion/deletion (I/D) polymorphism and knee OA. However, the results were conflicting. To determine the association between ACE gene I/D polymorphism and knee OA, we conducted a hospital-based case–control study with 282 knee OA cases and 316 controls to investigate the association between ACE gene I/D polymorphism and knee OA susceptibility in a Chinese Han population. The present study found that DD genotype or D allele carriers of ACE gene I/D polymorphism increased the risk of knee OA. Stratification analyses of sex, age, and body mass index (BMI) showed significant associations amongst the groups of females, ≥55 years, and abnormal BMI. In addition, the present study made analysis between ACE I/D polymorphism and some clinical features of OA, and found DD genotype of I/D polymorphism was associated with arthralgia. Furthermore, we undertook a meta-analysis together with the present study between this single nucleotide polymorphism (SNP) and knee OA risk. This meta-analysis found that ACE gene I/D polymorphism was associated with increased risk for OA. Stratification analysis of ethnicity in this meta-analysis indicated that I/D polymorphism increased the risk of knee OA amongst the Asians and Caucasians. In conclusion, this case–control study and meta-analysis suggest that ACE gene I/D polymorphism is associated with increased risk for knee OA.

Association between the polymorphisms of CALM1 gene and osteoarthritis risk: a meta-analysis based on observational studies

The existing studies on the association between polymorphisms of Calmodulin 1 (CALM1) gene and the risk of osteoarthritis (OA, a complex multifactorial disease and a major degenerative form of arthritis) in different populations have yielded conflicting findings. Therefore, we conducted a meta-analysis by systematically searching PubMed, Embase, Medline, Cochrane Library and Google Scholar, and assessing this association by calculating pooled odds ratios with 95% confidence intervals. Subgroup analyses stratified by ethnicity, OA type, and genotype were also conducted. Six studies (2752 cases and 3259 controls) involving six single nucleotide polymorphisms were included. Our data suggested that the T allele and genotype TT of the rs12885713 polymorphism, and the C allele of the rs2300496 polymorphism in the CALM1 gene all increased the risk of OA. The pooled results revealed no significant association between the CALM1 rs3213718 polymorphism and the risk of OA. Stratification analyses by ethnicity and OA type showed that the rs12885713 polymorphism increased the risk of OA among Asians and in knee OA, respectively. In conclusion, the rs12885713 and rs2300496 polymorphisms of the CALM1 gene may both increase the risk of OA. Owing to the limitations of the present study, this finding should be further confirmed in future well-designed studies.

The association between MMP-1 gene rs1799750 polymorphism and knee osteoarthritis risk

Matrix metalloproteinase 1 (MMP-1) degrades cartilage, which may result in osteoarthritis (OA) development. Several studies have explored the association between MMP-1 gene rs1799750 polymorphism and OA in different populations. However, the results are inconsistent. The aim of this case–control study was to investigate the association between MMP-1 gene rs1799750 polymorphism and knee OA in a Chinese population. The present study included 308 cases and 404 controls. Genotyping was performed using standard polymerase chain reaction and restriction fragment length polymorphism. The present study found that 2G2G genotype (2G2G vs 1G1G: OR & 95% CI, 2.28 (1.47–3.53), P<0.001; 2G2G + 1G2G vs 1G1G: OR & 95% CI, 1.61 (1.15–2.24), P=0.005; 2G2G vs 1G2G + 1G1G: OR & 95% CI, 1.84 (1.26–2.68), P=0.002) or 2G allele carriers (2G vs 1G: OR & 95% CI, 1.48 (1.20–1.83), P<0.001) of MMP-1 gene rs1799750 polymorphism increased the risk of OA. In conclusion, this case–control study confirms that MMP-1 gene rs1799750 polymorphism increases the risk of knee OA in Chinese Han population.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Identification and Analysis of Blood Gene Expression Signature for Osteoarthritis With Advanced Feature Selection Methods

Osteoarthritis (OA) is a complex disease that affects articular joints and may cause disability. The incidence of OA is extremely high. Most elderly people have the symptoms of osteoarthritis. The physiotherapy of OA is time consuming, and the chances of full recovery from OA are very minimal. The most effective way of fighting OA is early diagnosis and early intervention. Liquid biopsy has become a popular noninvasive test. To find the blood gene expression signature for OA, we reanalyzed the publicly available blood gene expression profiles of 106 patients with OA and 33 control samples using an automatic computational pipeline based on advanced feature selection methods. Finally, a compact 23-gene set was identified. On the basis of these 23 genes, we constructed a Support Vector Machine (SVM) classifier and evaluated it with leave-one-out cross-validation. Its sensitivity (Sn), specificity (Sp), accuracy (ACC), and Mathew's correlation coefficient (MCC) were 0.991, 0.909, 0.971, and 0.920, respectively. Obviously, the performance needed to be validated in an independent large dataset, but the in-depth biological analysis of the 23 biomarkers showed great promise and suggested that mRNA surveillance pathway and multicellular organism growth played important roles in OA. Our results shed light on OA diagnosis through liquid biopsy.

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.

Osteoarthritis: Trauma vs Disease

Osteoarthritis (OA) is the most prevalent joint disease characterized by pain and degenerative lesions of the cartilage, subchondral bone, and other joint tissues. The causes of OA remain incompletely understood. Over the years, it has become recognized that OA is a multifactorial disease. In particular, aging and trauma are the main risk factors identified for the development of OA; however, other factors such as genetic predisposition, obesity, inflammation, gender and hormones, or metabolic syndrome contribute to OA development and lead to a more severe outcome. While this disease mainly affects people older than 60 years, OA developed after joint trauma affects all range ages and has a particular impact on young individuals and people who have highest levels of physical activity such as athletes. Traumatic injury to the joint often results in joint instability or intra-articular fractures which lead to posttraumatic osteoarthritis (PTOA). In response to injury, several molecular mechanisms are activated, increasing the production and activation of different factors that contribute to the progression of OA.In this chapter, we have focused on the interactions and contribution of the multiple factors involved in joint destruction and progression of OA. In addition, we overview the main changes and molecular mechanisms related to OA pathogenesis.

Genome-wide association and functional studies identify a role for matrix Gla protein in osteoarthritis of the hand

OBJECTIVE

Osteoarthritis (OA) is the most common form of arthritis and the leading cause of disability in the elderly. Of all the joints, genetic predisposition is strongest for OA of the hand; however, only few genetic risk loci for hand OA have been identified. Our aim was to identify novel genes associated with hand OA and examine the underlying mechanism.

METHODS

We performed a genome-wide association study of a quantitative measure of hand OA in 12 784 individuals (discovery: 8743, replication: 4011). Genome-wide significant signals were followed up by analysing gene and allele-specific expression in a RNA sequencing dataset (n=96) of human articular cartilage.

RESULTS

We found two significantly associated loci in the discovery set: at chr12 (p=3.5 × 10-10) near the matrix Gla protein (MGP) gene and at chr12 (p=6.1×10-9) near the CCDC91 gene. The DNA variant near the MGP gene was validated in three additional studies, which resulted in a highly significant association between the MGP variant and hand OA (rs4764133, Betameta=0.83, Pmeta=1.8*10-15). This variant is high linkage disequilibrium with a coding variant in MGP, a vitamin K-dependent inhibitor of cartilage calcification. Using RNA sequencing data from human primary cartilage tissue (n=96), we observed that the MGP RNA expression of the hand OA risk allele was significantly lowercompared with the MGP RNA expression of the reference allele (40.7%, p<5*10-16).

CONCLUSIONS

Our results indicate that the association between the MGP variant and increased risk for hand OA is caused by a lower expression of MGP, which may increase the burden of hand OA by decreased inhibition of cartilage calcification.

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.

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.

Association study of candidate genes for susceptibility to Kashin-Beck disease in a Tibetan population

Background

Many osteoarthritis (OA) susceptibility genes have been identified in recent years. Given the overlap in the phenotype of joint inflammation between OA and Kashin-Beck disease (KBD), the aim of this study is to explore whether the reported OA susceptibility genes and two genes that may link to OA pathophysiology are associated with KBD in the Tibetan population.

Method

Fifteen single-nucleotide polymorphisms (SNPs) in 12 candidate genes previously reported as OA susceptibility loci were selected for investigation. Genotyping was performed using the SNaPshot method for these SNPs in a Tibetan population composed of 849 KBD patients and 565 normal controls. Meanwhile, the coding regions of two genes, COL10A1 and HABP2, which may involve in the pathological mechanism of OA/KBD, were sequenced by Sanger sequencing to identify susceptibility coding variants for KBD in the Tibetan population.

Results

The two arthritis-susceptible candidate SNPs, rs7775 (p.Arg324Gly) in the FRZB gene and rs7033979 in the ASPN gene, showed associations with KBD (OR = 1.568, P = 4 × 10⁻³ and OR = 0.744, P = 8 × 10⁻³, respectively). The coding variants rs142463796 (p.Asp128Asn) and rs2228547 (p.Gly545Arg) in the COL10A1 gene (OR = 9.832 and P = 6 × 10⁻³ and OR = 1.242, P = 0.043, respectively) and rs548354451 (p.Asp272Glu) in the HABP2 gene (OR = 2.813, P = 0.010) were associated with KBD patients.

Conclusion

These finding suggested that rs7775 in the FRZB gene may increase susceptibility to KBD, while rs7033979 in the ASPN gene may play a protective role in susceptibility to KBD in Tibetans. Moreover, genetic variants in chondrogenesis-related genes COL10A1 and HABP2 may play a role in the risk of developing KBD in the Tibetan population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-017-0423-6) contains supplementary material, which is available to authorized users.

TIMP-2 SNPs rs7342880 and rs4789936 are linked to risk of knee osteoarthritis in the Chinese Han Population

This study aimed to investigate whether functional polymorphisms in the tissue inhibitors of metalloproteinase-2 (TIMP-2) gene are associated with susceptibility to knee osteoarthritis (OA) in the Chinese Han population. Six TIMP-2 single nucleotide polymorphisms (SNPs) were assayed using MassARRAY in 300 patients clinically and radiographically diagnosed with knee OA and in 428 controls. Allelic and genotypic frequencies were compared between groups. Logistic regression adjusting for age and gender was used to estimate risk associations between specific genotypes and knee OA by computing odds ratios (ORs) and 95% confidence intervals (95% CIs). We found that allele "A" in rs7342880 was significantly associated with increased risk of knee OA (OR = 1.44, 95%CI = 1.09-1.91, p = 0.035). In addition, in the over-dominant model, rs4789936 correlated with reduced risk of knee OA, adjusting for age and gender (OR = 0.69, 95%CI = 0.49-0.98, p = 0.036). Finally, rs7342880 correlated with increased risk of knee OA in females. This study provides evidence that TIMP-2 is a knee OA susceptibility gene in the Chinese population and a potential diagnostic and preventive marker for the disease.

Association between Knee Osteoarthritis, Cardiovascular Risk Factors, and the Framingham Risk Score in South Koreans: A Cross-Sectional Study

Background

Osteoarthritis is a significant burden on personal health and for social cost, and its prevalence is rising. Recent research has revealed an association between osteoarthritis and cardiovascular disease, and this study uses the Framingham risk score (FRS), which is widely used as a composite index of cardiovascular risk factors, to investigate the association between osteoarthritis and various cardiovascular risk factors.

Methods

A total 9,514 participants aged 50 years or older who received knee X-ray diagnosis of the 5th Korean National Health and Nutrition Examination Survey (total surveyees = 24,173) released by the Korean Centers for Disease Control and Prevention was included for analysis. Knee osteoarthritis patients were defined as participants with K-L grade ≥2 on knee X-ray regardless of knee pain. The association between major cardiovascular risk factors (blood pressure, diabetes, cholesterol, and smoking habits), FRS, and knee osteoarthritis was analyzed, adjusting for various covariates.

Results

Prevalence of knee osteoarthritis in Koreans aged ≥50 years was 36.6%, and higher in women (men: 24.9%, women: 45.4%). Prevalence of knee osteoarthritis in participants with hypertension was significantly higher than those without hypertension (fully adjusted odds ratio (OR) 1.26; 95% confidence interval (CI) 1.08–1.48). Knee osteoarthritis prevalence was also higher in participants with impaired fasting glucose or diabetes than those without (age, sex adjusted OR 1.19; 95% CI 1.00–1.41). Also, OR values increased statistically significantly with FRS as a continuous variable (fully adjusted OR 1.007; 95% CI 1.00–1.01).

Conclusions

Prevalence of knee osteoarthritis was associated with hypertension and diabetes, which are major cardiovascular risk factors, and the FRS. Further studies on FRS pertaining to its relationship with osteoarthritis are warranted.

Identification of novel genetic variations affecting osteoarthritis patients

Background

Osteoarthritis (OA) is a progressive joint disease characterized by gradual degradation of extracellular matrix (ECM) components in the cartilage and bone. The ECM of cartilage is a highly specified structure that is mainly composed of type II collagen and provides tensile strength to the tissue via aggrecan and proteoglycans. However, changes in the ECM composition and structure can lead to loss of collagen type II and network integrity. Several risk factors have been correlated with OA including age, genetic predisposition, hereditary factors, obesity, mechanical injuries, and joint trauma. Certain genetic association studies have identified several genes associated with OA using genome-wide association studies (GWASs).

Results

We identified several novel genetic variants affecting genes that function in several candidate causative pathways including immune responses, inflammatory and cartilage degradation such as SELP, SPN, and COL6A6.

Conclusions

The approach of whole-exome sequencing can be a promising method to identify genetic mutations that can influence the OA disease.

Angiotensin-Converting Enzyme Insertion/Deletion Polymorphism and Susceptibility to Osteoarthritis of the Knee: A Case-Control Study and Meta-Analysis

Background

Studies of angiotensin-converting enzyme insertion/deletion (ACE I/D) polymorphisms and the risks of knee osteoarthritis (OA) have yielded conflicting results.

Objective

To determine the association between ACE I/D and knee OA, we conducted a combined case-control study and meta-analysis.

Methods

For the case-control study, 447 knee OA cases and 423 healthy controls were recruited between March 2010 and July 2011. Knee OA cases were defined using the Kellgren-Lawrence grading system, and the ACE I/D genotype was determined using a standard polymerase chain reaction. The association between ACE I/D and knee OA was detected using allele, genotype, dominant, and recessive models. For the meta-analysis, PubMed and Embase databases were systematically searched for prospective observational studies published up until August 2015. Studies of ACE I/D and knee OA with sufficient data were selected. Pooled results were expressed as odds ratios (ORs) with corresponding 95% confidence intervals (CI) for the D versus I allele with regard to knee OA risk.

Results

We found no significant association between the D allele and knee OA [OR: 1.09 (95% CI: 0.76–1.89)] in the present case-control study, and the results of other genetic models were also nonsignificant. Five current studies were included, and there were a total of six study populations after including our case-control study (1165 cases and 1029 controls). In the meta-analysis, the allele model also yielded nonsignificant results [OR: 1.37 (95% CI: 0.95–1.99)] and a high heterogeneity (I²: 87.2%).

Conclusions

The association between ACE I/D and knee OA tended to yield negative results. High heterogeneity suggests a complex, multifactorial mechanism, and an epistasis analysis of ACE I/D and knee OA should therefore be conducted.

Genetic Contribution to the Development of Radiographic Knee Osteoarthritis in a Population Presenting with Nonacute Knee Symptoms a Decade Earlier

This study examined the contribution of the osteoarthritis (OA) susceptibility genes ASPN, GDF5, DIO2, and the 7q22 region to the development of radiographic knee OA in patients with a mean age of 40.6 ± 7.9 years (standard deviation) and who suffered from nonacute knee complaints a decade earlier. Dose–response associations of four single nucleotide polymorphisms(SNPs) in the susceptibility genes were determined by comparing 36 patients who showed the development of OA on radiographs (Kellgren and Lawrence score ≥1) with 88 patients having normal cartilage with no development of OA on radiographs. Multivariate logistic regression analysis including the variables such as age, gender, body mass index, and reported knee trauma was performed. A dose–response association of DIO2 SNP rs225014: odds ratio (OR) 2.3, 95% confidence interval (CI) 1.1–4.5 (P = 0.019) and GDF5 SNP rs143383: OR 2.0, 95% CI 1.1–3.8 (P = 0.031) was observed with knee OA development. The ASPN and 7q22 SNPs were not associated with OA development.

Correlation of the osteoarthritis susceptibility variants that map to chromosome 20q13 with an expression quantitative trait locus operating on NCOA3 and with functional variation at the polymorphism rs116855380

Objective

To functionally characterize the osteoarthritis (OA) susceptibility variants that map to a region of high linkage disequilibrium (LD) on chromosome 20q13 marked by the single‐nucleotide polymorphism (SNP) rs6094710 and encompassing NCOA3 and SULF2.

Methods

Nucleic acids were extracted from the cartilage of OA patients. Overall and allelic expression of NCOA3 and SULF2 were measured by quantitative reverse transcription–polymerase chain reaction and pyrosequencing, respectively. The functional effect of SNPs within the 20q13 locus was assessed in vitro using luciferase reporter constructs and electrophoretic mobility shift assays (EMSAs). The in vivo effect of nuclear receptor coactivator 3 (NCOA3) protein depletion on primary human OA articular cartilage chondrocytes was assessed using RNA interference.

Results

Expression of NCOA3 correlated with the genotype at rs6094710 (P = 0.006), and the gene demonstrated allelic expression imbalance (AEI) in individuals heterozygous for the SNP (mean AEI 1.21; P < 0.0001). In both instances, expression of the OA‐associated allele was reduced. In addition, there was reduced enhancer activity of the OA‐associated allele of rs116855380, a SNP in perfect LD with rs6094710 in luciferase assays (P < 0.001). EMSAs demonstrated a protein complex binding with reduced affinity to this allele. Depletion of NCOA3 led to significant changes (all P < 0.05) in the expression of genes involved in cartilage homeostasis.

Conclusion

NCOA3 is subject to a cis‐acting expression quantitative trait locus in articular cartilage, which correlates with the OA association signal and with the OA‐associated allele of the functional SNP rs116855380, a SNP that is located only 10.3 kb upstream of NCOA3. These findings elucidate the effect of the association of the 20q13 region on OA cartilage and provide compelling evidence of a potentially causal candidate SNP.

Single Nucleotide Polymorphisms and Osteoarthritis: An Overview and a Meta-Analysis

Osteoarthritis (OA) is a complex disorder characterized by degenerative articular cartilage and is largely attributed to genetic risk factors. Single nucleotide polymorphisms (SNPs) are common DNA variants that have shown promising and efficiency, compared with positional cloning, to map candidate genes of complex diseases, including OA. In this study, we aim to provide an overview of multiple SNPs from a number of genes that have recently been linked to OA susceptibility. We also performed a comprehensive meta-analysis to evaluate the association of SNP rs7639618 of double von Willebrand factor A domains (DVWA) gene with OA susceptibility. A systematic search of studies on the association of SNPs with susceptibility to OA was conducted in PubMed and Google scholar. Studies subjected to meta-analysis include human and case-control studies that met the Hardy-Weinberg equilibrium model and provide sufficient data to calculate an odds ratio (OR). A total of 9500 OA cases and 9365 controls in 7 case-control studies relating to SNP rs7639618 were included in this study and the ORs with 95% confidence intervals (CIs) were calculated. Over 50 SNPs from different genes have been shown to be associated with either hip (23), or knee (20), or both (13) OA. The ORs of these SNPs for OA and the subtypes are not consistent. As to SNP rs7639618 of DVWA, increased knee OA risk was observed in all genetic models analyzed. Specifically, people from Asian with G-allele showed significantly increased risk of knee OA (A versus G: OR = 1.28, 95% CI 1.13-1.46; AA versus GG: OR = 1.60, 95% CI 1.25-2.05; GA versus GG: OR = 1.31, 95% CI 1.18-1.44; AA versus GA+GG: OR = 1.34, 95% CI 1.12-1.61; AA+GA versus GG: OR = 1.40, 95% CI 1.19-1.64), but not in Caucasians or with hip OA. Our results suggest that multiple SNPs play different roles in the pathogenesis of OA and its subtypes; SNP rs7639618 of DVWA gene is associated with a significantly increased risk of knee OA in Asians. Given the limited sample size, further studies are needed to evaluate this observation.

Methylation quantitative trait locus analysis of osteoarthritis links epigenetics with genetic risk

Osteoarthritis (OA) is a common, painful and debilitating disease of articulating joints resulting from the age-associated loss of cartilage. Well-powered genetic studies have identified a number of DNA polymorphisms that are associated with OA susceptibility. Like most complex trait loci, these OA loci are thought to influence disease susceptibility through the regulation of gene expression, so-called expression quantitative loci, or eQTLs. One mechanism through which eQTLs act is epigenetic, by modulating DNA methylation. In such cases, there are quantitative differences in DNA methylation between the two alleles of the causal polymorphism, with the association signal referred to as a methylation quantitative trait locus, or meQTL. In this study, we aimed to investigate whether the OA susceptibility loci identified to date are functioning as meQTLs by integrating genotype data with whole genome methylation data of cartilage DNA. We investigated potential genotype-methylation correlations within a 1.0-1.5 Mb region surrounding each of 16 OA-associated single-nucleotide polymorphisms (SNPs) in 99 cartilage samples and identified four that function as meQTLs. Three of these replicated in an additional cohort of up to 62 OA patients. These observations suggest that OA susceptibility loci regulate the level of DNA methylation in cis and provide a mechanistic explanation as to how these loci impact upon OA susceptibility, further increasing our understanding of the role of genetics and epigenetics in this common disease.

Targets, models and challenges in osteoarthritis research

Osteoarthritis is a chronic degenerative disorder of the joint and represents one of the most common diseases worldwide. Its prevalence and severity are increasing owing to aging of the population, but treatment options remain largely limited to painkillers and anti-inflammatory drugs, which only provide symptomatic relief. In the late stages of the disease, surgical interventions are often necessary to partially restore joint function. Although the focus of osteoarthritis research has been originally on the articular cartilage, novel findings are now pointing to osteoarthritis as a disease of the whole joint, in which failure of different joint components can occur. In this Review, we summarize recent progress in the field, including data from novel ‘omics’ technologies and from a number of preclinical and clinical trials. We describe different in vitro and in vivo systems that can be used to study molecules, pathways and cells that are involved in osteoarthritis. We illustrate that a comprehensive and multisystem approach is necessary to understand the complexity and heterogeneity of the disease and to better guide the development of novel therapeutic strategies for osteoarthritis.

Association between genetic variants of DVWA and osteoarthritis of the knee and hip: a comprehensive meta-analysis

Recently, double von Willebrand factor domain A (DVWA) gene, a previously unknown gene, was revealed to contain several single nucleotide polymorphisms (SNPs) that showed consistent association with knee osteoarthritis (OA) in Japanese and Chinese cohorts. However, subsequent studies failed to confirm this result in several different populations. To deal with the issues raised by inconsistent results among those studies, we investigated the association between DVWA and OA using meta-analytic techniques, combining all published data up to December 2014. 10 independent samples from 4 teams contributed data for a possible association between SNP rs7639618 and knee or hip OA. The total number of cases and controls of this SNP was respectively 4,142 versus 6,575 for knee OA, and 2,325 versus 2,914 for hip OA. A trend of significant association was observed in the combined population with knee OA (P=0.06), and a significant difference was identified between patients with knee OA and controls for the G-allele of rs7639618 (P=0.02). Together with the reported functional studies, our results indicate that DVWA may have a small but strong effect on the susceptibility to knee OA, at least in Asian population. Further functional studies are needed to determine the underlying variation of DVWA and to relate this to the pathophysiology of OA.

Genetic polymorphisms of interleukin-16 and risk of knee osteoarthritis

BACKGROUND

Interleukin-16 (IL-16), a pleiotropic cytokine, plays a fundamental role in inflammatory diseases. This study investigates the association between IL-16 polymorphisms and the risk of knee osteoarthritis (OA) in a Chinese population.

METHODS

The IL-16 rs11556218, rs4072111, and rs4778889 polymorphisms were determined in 150 knee OA cases and 147 healthy controls through polymerase chain reaction-restriction fragment length polymorphism.

RESULTS

The results suggested that the variants in IL-16 gene rs11556218 site were associated with a decreased knee OA risk after adjusting for age, sex, BMI, and smoking and drinking status (TG vs. TT: OR, 0.69; 95% CI, 0.53-0.89; P = 0.006; GG vs. TT: OR, 0.64; 95% CI, 0.45-0.90; P = 0.042; dominant model: OR, 0.68; 95% CI, 0.29-0.87; P = 0.002; G vs. T allele: OR, 0.77; 95% CI, 0.66-0.90; P = 0.003). Similarly, subjects bearing the rs4072111 variant genotypes and alleles also had a lower susceptibility to knee OA compared with those bearing the wild-type (CT vs. CC: OR, 0.66; 95% CI, 0.53-0.83; P = 0.002; TT vs. CC: OR, 0.57; 95% CI, 0.40-0.82; P = 0.027; dominant model: OR, 0.65; 95%, CI 0.52-0.80; P <0.001; T vs. C allele: OR, 0.69; 95% CI, 0.58-0.81; P <0.001). Further, the C allele and the combined genotype (CC+CT) of rs4778889 were associated with a slightly decreased risk of knee OA. In addition, we found two high-risk haplotypes: TTT (OR, 3.70) and GCC (OR, 6.22). Finally, serum IL-16 levels of knee OA patients were significantly higher than those of controls (P = 0.001).

CONCLUSIONS

Despite the small sample size, this is the first study suggesting IL-16 gene polymorphisms to be associated with the risk of knee OA.

A genome-wide association study of copy-number variation identifies putative loci associated with osteoarthritis in Koreans

Background

OA is a complex disease caused by environmental and genetic risk factors. The purpose of this study is to identify candidate copy number variations (CNVs) associated with OA.

Methods

We performed a genome-wide association study of CNV to identify potential loci that confer susceptibility to or protection from OA. CNV genotyping was conducted using NimbleGen HD2 3 × 720K comparative hybridization array and included samples from 371 OA patients and 467 healthy controls. The putative CNV regions identified were confirmed with a TaqMan assay.

Results

We identified six genomic regions associated with OA encompassing CNV loci. None of six loci had previously been reported in genome-wide association studies with OA, although a genetic analysis suggested that they have functional effects. The protein product of a candidate risk gene for obesity, TNKS, targets Wnt inhibition, and this gene was significantly associated with hand and knee OA. Copy number deletion on TNKS was associated with a 1.37-fold decreased risk for OA. In addition, CA10, which shows a strong association with osteoporosis, was also significant in our study. Copy number deletion on this gene was associated with a 1.69-fold decreased risk for OA.

Conclusion

We identified several CNV loci that may contribute to OA susceptibility in Koreans. Further functional investigations of these genes are warranted to fully characterize their putative association.

Electronic supplementary material

The online version of this article (doi:10.1186/s12891-015-0531-4) contains supplementary material, which is available to authorized users.

Adult mice lacking the type 2 iodothyronine deiodinase have increased subchondral bone but normal articular cartilage

BACKGROUND

Although osteoarthritis (OA) is the commonest joint disorder and has a rising prevalence as the population ages, no drugs are available that prevent or delay the onset and progression of disease. Recent studies identified the DIO2 gene encoding type 2 deiodinase (D2) as a susceptibility locus for OA, and further data suggest deiodinase-regulated local availability of triiodothyronine (T3) in the joint plays an important role in cartilage maintenance and repair. To investigate the hypothesis that reduced tissue T3 availability protects joints from development of OA, the joint phenotypes of adult mice lacking D2 (D2KO) or lacking both D1 and D2 (D1D2KO), the only enzymes that catalyze conversion of the prohormone thyroxine to active T3, were determined.

METHODS

Knee joints were prepared from male 16-week-old adult wild type (WT; n=9), D2KO (n=5), and D1D2KO (n=3) mice. Articular cartilage pathology was scored using the Osteoarthritis Research Society International (OARSI) histopathology scale for murine OA to determine the severity and extent of disease. Digital X-ray microradiography was used to determine the area and mineral content of subchondral bone immediately beneath the articular cartilage surface.

RESULTS

There were no differences in maximum and standardized OA scores, cartilage erosion indices, or articular cartilage cellularity among WT, D2KO, and D1D2KO mice. Subchondral bone area did not differ among genotypes, but mineral content was markedly increased in both D2KO and D1D2KO mice compared to WT.

CONCLUSIONS

Although adult D2KO mice have normal articular cartilage and no other features of spontaneous joint damage, they exhibit increased subchondral bone mineral content.

Role of Vitamin D in Osteoarthritis: Molecular, Cellular, and Clinical Perspectives

Osteoarthritis is a debilitating and degenerative disease which affects millions of people worldwide. The causes and mechanisms of osteoarthritis remain to be fully understood. Vitamin D has been hypothesised to play essential roles in a number of diseases including osteoarthritis. Many cell types within osteoarthritic joints appear to experience negative effects often at increased sensitivity to vitamin D. These findings contrast clinical research which has identified vitamin D deficiency to have a worryingly high prevalence among osteoarthritis patients. Randomised-controlled trial is considered to be the most rigorous way of determining the effects of vitamin D supplementation on the development of osteoarthritis. Studies into the effects of low vitamin D levels on pain and joint function have to date yielded controversial results. Due to the apparent conflicting effects of vitamin D in knee osteoarthritis, further research is required to fully elucidate its role in the development and progression of the disease as well as assess the efficacy and safety of vitamin D supplementation as a therapeutic strategy.

A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip

OBJECTIVES

Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects.

METHODS

We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used.

RESULTS

We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis).

CONCLUSIONS

Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

Orphan G-protein coupled receptor 22 (Gpr22) regulates cilia length and structure in the zebrafish Kupffer's vesicle

GPR22 is an orphan G protein-coupled receptor (GPCR). Since the ligand of the receptor is currently unknown, its biological function has not been investigated in depth. Many GPCRs and their intracellular effectors are targeted to cilia. Cilia are highly conserved eukaryotic microtubule-based organelles that protrude from the membrane of most mammalian cells. They are involved in a large variety of physiological processes and diseases. However, the details of the downstream pathways and mechanisms that maintain cilia length and structure are poorly understood. We show that morpholino knock down or overexpression of gpr22 led to defective left-right (LR) axis formation in the zebrafish embryo. Specifically, defective LR patterning included randomization of the left-specific lateral plate mesodermal genes (LPM) (lefty1, lefty2, southpaw and pitx2a), resulting in randomized cardiac looping. Furthermore, gpr22 inactivation in the Kupffer’s vesicle (KV) alone was still able to generate the phenotype, indicating that Gpr22 mainly regulates LR asymmetry through the KV. Analysis of the KV cilia by immunofluorescence and transmission electron microscopy (TEM), revealed that gpr22 knock down or overexpression resulted in changes of cilia length and structure. Further, we found that Gpr22 does not act upstream of the two cilia master regulators, Foxj1a and Rfx2. To conclude, our study characterized a novel player in the field of ciliogenesis.

Genes involved in the osteoarthritis process identified through genome wide expression analysis in articular cartilage; the RAAK study

Objective

Identify gene expression profiles associated with OA processes in articular cartilage and determine pathways changing during the disease process.

Methods

Genome wide gene expression was determined in paired samples of OA affected and preserved cartilage of the same joint using microarray analysis for 33 patients of the RAAK study. Results were replicated in independent samples by RT-qPCR and immunohistochemistry. Profiles were analyzed with the online analysis tools DAVID and STRING to identify enrichment for specific pathways and protein-protein interactions.

Results

Among the 1717 genes that were significantly differently expressed between OA affected and preserved cartilage we found significant enrichment for genes involved in skeletal development (e.g. TNFRSF11B and FRZB). Also several inflammatory genes such as CD55, PTGES and TNFAIP6, previously identified in within-joint analyses as well as in analyses comparing preserved cartilage from OA affected joints versus healthy cartilage were among the top genes. Of note was the high up-regulation of NGF in OA cartilage. RT-qPCR confirmed differential expression for 18 out of 19 genes with expression changes of 2-fold or higher, and immunohistochemistry of selected genes showed a concordant change in protein expression. Most of these changes associated with OA severity (Mankin score) but were independent of joint-site or sex.

Conclusion

We provide further insights into the ongoing OA pathophysiological processes in cartilage, in particular into differences in macroscopically intact cartilage compared to OA affected cartilage, which seem relatively consistent and independent of sex or joint. We advocate that development of treatment could benefit by focusing on these similarities in gene expression changes and/or pathways.

Assessment of osteoarthritis candidate genes in a meta-analysis of nine genome-wide association studies

Objective

To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA.

Methods

A total of 199 candidate genes for association with OA were identified using Human Genome Epidemiology (HuGE) Navigator. All of their single-nucleotide polymorphisms (SNPs) with an allele frequency of >5% were assessed by fixed-effects meta-analysis of 9 genome-wide association studies (GWAS) that included 5,636 patients with knee OA and 16,972 control subjects and 4,349 patients with hip OA and 17,836 control subjects of European ancestry. An additional 5,921 individuals were genotyped for significantly associated SNPs in the meta-analysis. After correction for the number of independent tests, P values less than 1.58 × 10⁻⁵ were considered significant.

Results

SNPs at only 2 of the 199 candidate genes (COL11A1 and VEGF) were associated with OA in the meta-analysis. Two SNPs in COL11A1 showed association with hip OA in the combined analysis: rs4907986 (P = 1.29 × 10⁻⁵, odds ratio [OR] 1.12, 95% confidence interval [95% CI] 1.06−1.17) and rs1241164 (P = 1.47 × 10⁻⁵, OR 0.82, 95% CI 0.74−0.89). The sex-stratified analysis also showed association of COL11A1 SNP rs4908291 in women (P = 1.29 × 10⁻⁵, OR 0.87, 95% CI 0.82−0.92); this SNP showed linkage disequilibrium with rs4907986. A single SNP of VEGF, rs833058, showed association with hip OA in men (P = 1.35 × 10⁻⁵, OR 0.85, 95% CI 0.79−0.91). After additional samples were genotyped, association at one of the COL11A1 signals was reinforced, whereas association at VEGF was slightly weakened.

Conclusion

Two candidate genes, COL11A1 and VEGF, were significantly associated with OA in this focused meta-analysis. The remaining candidate genes were not associated.