Email updates

Keep up to date with the latest news and content from Arthritis Research & Therapy and BioMed Central.

Open Access Research article

Genetic, clinical and radiographic signs in knee osteoarthritis susceptibility

Luigi Minafra1*, Valentina Bravatà1, Michele Saporito2, Francesco P Cammarata1, Giusi I Forte1, Salvatore Caldarella1, Michele D’Arienzo2, Maria C Gilardi134, Cristina Messa145 and Filippo Boniforti6

Author Affiliations

1 IBFM CNR-LATO, Contrada Pietrapollastra Pisciotto, 90015 Cefalù, PA, Italy

2 Clinica Ortopedica e Traumatologica, Policlinico Universitario P. Giaccone-Università degli Studi di Palermo, via del Vespro 143, 90127 Palermo, Italy

3 Nuclear Medicine, San Raffaele Scientific Institute, via Olgettina 60, 20132 Milan, Italy

4 Department of Health Sciences, Tecnomed Foundation, University of Milano-Bicocca, via Pergolesi 33, 20900 Monza (MB), Italy

5 Nuclear Medicine Center, San Gerardo Hospital, via Pergolesi 33, 20900 Monza (MB), Italy

6 Unità Operativa di Ortopedia, San Raffaele Hospital ‘G. Giglio’, Contrada Pietrapollastra Pisciotto, 90015 Cefalù (PA), Italy

For all author emails, please log on.

Arthritis Research & Therapy 2014, 16:R91  doi:10.1186/ar4535


The electronic version of this article is the complete one and can be found online at: http://arthritis-research.com/content/16/2/R91


Received:22 October 2013
Accepted:25 March 2014
Published:9 April 2014

© 2014 Minafra et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Introduction

Osteoarthritis (OA) is considered to be a multifactorial and polygenic disease and diagnosis is mainly clinical and radiological. Correlation between radiographic data and clinical status has been reported. However, very few studies, especially in Caucasian people, describe the association between the Kellgren and Lawrence OA grading scale (KL) and genetic alterations to better understand OA etiopathogenesis and susceptibility. In order to update the knee OA grading, in this study we assessed the associations between KL grade, clinical features such as American Knee Society Score (AKSS), age, and polymorphisms in the principal osteoarthritis susceptibility (OS) genes in Sicilian individuals.

Methods

In 66 Sicilian individuals affected by primary knee OA, the clinical and radiographic evaluation was performed using 2 sub-scores of AKSS (knee score (KS) and function score (FS)) and KL. The patients were also classified according to age. Online Mendelian Inheritance in Man (OMIM) and Database of Single Nucleotide Polymorphisms (dbSNP) Short Genetic Variations databases were used to select gene regions containing the following polymorphisms to analyze: FRZB rs288326 and rs7775, MATN3 rs77245812, ASPN D14 repeats, PTHR2 rs76758470, GDF5 rs143383 and DVWA rs11718863. Patient genotypes were obtained using Sanger DNA sequencing analysis.

Results

In our cohort of patients a statistical association between the variables analyzed was reported in all associations tested (KL versus KS, FS and age). We observed that a mild to severe OA radiographic grade is related to severe clinical conditions and loss of articular function and that the severity of symptoms increases with age. Concerning the genotyping analysis, our results revealed a significant statistical association between KL grading and GDF5 rs143383 and DVWA rs11718863 genetic alterations. The latter was also associated with a more severe radiographic grade, displaying its predictive role as OA marker progression. Statistically significant association between clinical, radiographic and genetic signs observed, suggests extending the actual grading of knee OA based mainly on X-ray features.

Conclusions

This work represents a multidisciplinary and translational medicine approach to study OA where clinical, radiological, and OS5 and OS6 SNPs evaluation could contribute to better define grading and progression of OA and to the development of new therapies.

Introduction

Osteoarthritis (OA) is a multifactorial, inflammatory and degenerative disorder of the joints [1-5]. OA involves the synovial tissues and articular cartilage, resulting in symptoms that cause a decrease in the quality of life and disability. Among the risk factors, age contributes to a substantially increased risk of knee OA onset and progression [6-9], even if the association of age with the progression of knee OA is sometimes conflicting [10].

OA diagnosis is mainly clinical and radiological, and X-ray images are the gold standard to confirm the clinical diagnosis and to grade the disease [11-13]. A critical point for OA diagnosis is to identify an early onset and an early progression of this disease. Many studies analyzed the correlation between knee OA radiographic data and clinical status of the affected joint by using specific clinical scores and radiographic grading scales. A strong statistical association between pain and the Kellgren and Lawrence osteoarthritis grading scale (KL) of knee OA was observed by Muraki and colleagues in an older Japanese population [14]. Moreover, a study conducted by Cho and colleagues in a Korean population showed that knee OA symptoms became more severe with high KL levels and also that KL values were related to sex: women had worse symptoms than men at the same radiographic grade [15]. In addition, Cubukcu and colleagues investigated the possible relation between KL and the Western Ontario and McMaster Universities subscore of function, pain and stiffness, but the association was not often verified [16]. These findings were in agreement with data reported by recent observations that knee radiographic OA results cannot always be related to knee pain, due to radiographic view extension, definition of pain and variability of the cohort selected for the studies [17]. The Ahlbäck OA radiographic grading scale was used by Hernàndez-Vaquero and Fernández-Carreira in a Spanish study to analyze the relation of OA with clinical status expressed by the Hospital for Special Surgery Knee Score. However, a small variation of the clinical status score compared with the radiographic grades was observed [18].

Nowadays, molecular genetic investigations have gained an increasingly significant role in the knowledge of OA etiology and have provided evidence for a genetic component to OA [19-21]. The completion of large genome-wide association studies introduced single nucleotide polymorphisms (SNPs) as risk factors for human disease. Several gene association analyses, either genome wide or a gene candidate approach, identified different genes related to the OA susceptibility, onset and progression. OA may thus be placed into the category of polygenic diseases [22-24], but the functional influence of specific SNPs on OA needs further research in order to contribute to the etiopathogenesis. Several association studies between SNPs and OA disease remain unconfirmed or controversial, due to bias in patient enrolling criteria, differences in OA-affected joint sites, in classification and staging modes, in the radiographic evaluation scales used and in subjective differences in patient’s pain evaluation scoring. Furthermore, it becomes of interest to explore the geographical and ethnic allele distribution, which is extremely important in fully understanding the SNP variant effects. Precisely, Sicilian individuals have a specific genetic background and different allele distribution compared with the rest of Europe and with the rest of Italy (north–south genetic trend), due to distinct gene–environment interactions and, certainly, due to deep human migration movements, which have occurred in Sicily over the centuries as described by several authors [25-28].

For all these reasons and the lack of comparison with studies conducted on Italian cohorts, we have chosen to initiate this study on our cohort of Sicilian OA patients starting from loci listed and described by the Online Mendelian Inheritance in Man database, widely recognized as a comprehensive and authoritative compendium of human genes and genetic phenotypes.

Particularly, in the Online Mendelian Inheritance in Man database – which collects known genetic lesions responsible for human inherited diseases – the following principal loci of osteoarthritis susceptibility (OS) and the associated polymorphisms, SNPs and aspartic acid (D) repeats, are reported: frizzled-related protein (FRZB) rs288326 (OS1A) and rs7775 (OS1B), MATN3 rs77245812 (OS2), ASPN D14 repeats (OS3), parathyroid hormone 2 (PTHR2) rs76758470 (OS4), growth and differentiation factor 5 (GDF5) rs143383 (OS5) and DVWA rs11718863 (OS6).

In particular, the FRZB gene is a member of a family of the soluble Wingless (Wnt) antagonist. Recent evidence has demonstrated that products of Wnt and Frizzled play a key role in the development and maintenance of bones and joints [29,30]. The rs7775 and rs288326 FRZB SNPs variants showed an increased frequency in subjects with generalized radiographic OA, as confirmed by other studies in Caucasian individuals [31,32].

Articular cartilage is composed of specialized cells, the chondrocytes, that produce a large amount of extracellular matrix composed of collagen fibers. MATN3 encodes a noncollagenous extracellular matrix protein expressed during the development of the skeletal system and in the cartilage [33]. Another extracellular matrix component deregulated in the articular cartilage of OA patients is asporin protein, encoded by the ASPN gene and expressed at high levels in knee and hip cartilage of individuals with ASPN D14 repeats [34]. The PTHR2 receptor is a member of the G-protein coupled receptor family 2 and its functional role in OA is based on the observation that PTHR2 is expressed in a number of endocrine cell types and regulates pituitary hormone secretion and specifically growth hormone [30]. The role of GDF5 in the development and maintenance of bone and cartilage has also been recognized. Mutations in GDF5 are involved in several disorders of skeletal development and also in hip and knee OA progression [35-37]. Finally, the DVWA gene, which encodes for a protein containing two von Willebrand A domains, was found to harbor the rs11718863 SNP, showing a consistent association with knee OA in Japanese and Chinese OA cohorts [38].

Although several studies described the association between these specific polymorphisms and susceptibility to OA, to our knowledge no studies have examined their simultaneous presence in OA patients, especially in the European people groups. The purpose of this pilot study is to highlight possible associations between KL grade, clinical features (American Knee Society Score (AKSS), age) and the abovementioned genetic polymorphisms in order to update the knee OA grading and to improve a personalized treatment program in the future.

Methods

Patients

On admission to hospital, 66 Sicilian patients affected by primary OA, aged 54 to 86 years and candidates for knee surgery of arthroscopy or arthroplasty, were enrolled in this project. The patients were grouped, according to age, into two groups: young (from 54 to 65 years old) and old (from 66 to 86 years old). Among these, 61 patients were selected for genotyping analysis due to availability of a blood sample. This study (named OA_BIOMOL_1) was approved by the Ethical Committee of the San Raffaele G. Giglio Hospital, Cefalù, Italy (number of protocol: CE 2011/63) and the patients gave their written informed consent according to the Helsinki Declaration.

Clinical evaluation

The clinical evaluation was performed for each patient before surgery using the AKSS, which includes two subscores: knee score (KS) and function score (FS). Each subscore ranges from 0 to 100 points. For KS evaluation, pain, range of motion, anterior–posterior and medio-lateral stability, flexion contracture, leg extension and varus–valgus alignment were investigated. FS evaluates knee function from a patient’s point of view, describing walking ability, climbing stairs ability and the use of walking aids. The AKSS was classified into three levels for both KS and FS: high (HKS, HFS), medium (MKS, MFS) and low (LKS, LFS).

The patients with LKS and LFS had scores between 0 and 49 points. The patients with MKS and MFS had scores between 50 and 69 points. The patients with HKS and HFS had scores from 70 to 100.

Radiographic evaluation

The radiographic evaluation was performed on antero-posterior and lateral X-ray views of the knee by a single investigator using the KL, which includes four grades: grade 1, possible narrowing of joint space and possible presence of osteophytes; grade 2, definite narrowing of joint space and definite osteophytes; grade 3, definite narrowing of joint space, multiple osteophytes, sclerosis, cysts and possible deformity of bone contour; and grade 4, marked narrowing of joint space, large osteophytes, severe sclerosis, cysts and definite deformity of bone contour [39,40]. The evaluation was undertaken on an X-ray performed no more than 4 months before surgery. In our study we grouped grade 1 and grade 2 into a single grade because the radiographic differences in our cohort were considered not relevant compared with those between KL grade 3 and grade 4. The KL classification was therefore summarized into three groups: group A (grades 1 and 2), group B (grade 3), and group C (grade 4).

Genetic analysis

The patients were genotyped by sequencing analysis, for the following genetic polymorphisms associated with OS, SNPs and D repeats: FRZB rs288326 (OS1A) and rs7775 (OS1B), MATN3 rs77245812 (OS2), PTHR2 rs76758470 (OS3), ASPN D14 repeats (OS4), GDF5 rs143383 (OS5) and DVWA rs11718863 (OS6). The Human Gene Mutation Database [41] and the dbSNP Short Genetic Variations database [42] were used to analyze gene regions containing the selected SNPs. Genomic DNA was extracted from peripheral blood using the QIAamp DNA blood mini kit, according to the manufacturer’s specifications (Qiagen Inc., Valencia, CA, USA).

After quality and quantity analysis, DNA was polymerase chain reaction amplified using primers designed by the Primer3 software [43] and listed in Table 1. Polymerase chain reaction reactions were performed with 50 ng genomic DNA in a total volume of 50 μl containing 1× PCR Gold Buffer, 1.5 mM di-MgCl2, 200 μM dNTPs, 200 nM forward and reverse primer mix and 1.25 U AmpliTaq Gold DNA Polymerase (Life Technologies Monza, MB, Italy). The thermal cycle profile employed a 5-minute denaturing step at 94°C, followed by 35 cycles at 94°C for 45 seconds, 59°C for 45 seconds and 72°C for 45 seconds, and a final extension step of 5 minutes at 72°C. The quality and quantity of polymerase chain reaction products were assessed on the Bioanalyzer instrument (Agilent Technologies, Santa Clara, CA, USA) and were purified using the QIAquick PCR purification kit, according to the manufacturer’s specifications (Qiagen Inc., Valencia, CA, USA).

Table 1. Primers sequence used for genotyping analysis

To perform DNA sequencing, purified amplicons were labeled with the BigDye Terminator v3.1 Cycle Sequencing Kit following the manufacturer’s standard protocol (Applied Biosystems). The thermal cycle profile employed a 1-minute denaturing step at 96°C, followed by 25 cycles at 96°C for 10 seconds, 54°C for 5 seconds and 60°C for 3 minutes. Labeled samples were purified with the X-terminator purification kit according to the manufacturer’s standard protocol and loaded in a 3500-Dx Genetic Analyzer (Applied Biosystems) for separation by capillary electrophoresis. Electropherograms and sequence files were analyzed using Sequencing Analysis and SeqScape software (Applied Biosystems).

Statistical analysis

The association between the clinical data (KS, FS, age) and the radiographic data (KL) and the association between genotypes and KL groups (A, B, C) were analyzed using GraphPad InStat software version 3.05 [44]. The Mann–Whitney U test, the chi-square test and Fisher’s exact test were performed. Differences in groups were considered significant when P ≤ 0.05. Hardy–Weinberg equilibrium was evaluated as described previously [28].

Results

Clinical and radiographic evaluation

We recruited 66 cases (37 females and 29 males), of which 24 were young (54 to 65 years) and 42 were old (66–86 years), and they were divided into three groups (A, B, C) depending on the degree of radiographic knee OA.

According to the clinical scores we classified the patients as follows. Group A consisted of 24 patients (11 females and 13 males, 14 young and 10 old); KS was poor in 13 cases and fair in 11 cases, and the average FS score was 51 points. Group B consisted of 21 patients (15 females and six males, eight young and 13 old); KS was poor in 19 cases and fair in two cases, and the average FS score was 41 points. Group C consisted of 21 patients (11 females and 10 males, two young and 19 old); KS was low in all cases, and the average FS was 35 points.

Regarding the treatment, 22 patients of group A underwent arthroscopy and two patients arthroplasty, two patients of group B underwent arthroscopy and 19 patients arthroplasty, and 21 patients of group C underwent arthroplasty (Table 2).

Table 2. Clinical features and treatment for each radiographic group of patients

According to the KS and FS scores, the patients in the LKS group were in the majority (n = 46), while the 36 patients in the LFS group generally had severe symptoms and high disability. There were 19 patients in the MKS group and 23 patients in the MFS group, and thus one-third of patients had moderate to severe symptoms and disability. One patient was in the HKS group and seven patients were in the HFS group, with zero to mild symptoms and disability (Table 3).

Table 3. Patient classification according to knee score and function score

Association between Kellgren and Lawrence osteoarthritis grading and knee score, function score and age

Association analyses were performed to verify the possible association between clinical data (KS, FS, age) and radiographic data (KL). A statistical association between the variables analyzed was observed (Table 4).

Table 4. Association between Kellgren and Lawrence osteoarthritis grading and knee score, function score and age

Kellgren and Lawrence osteoarthritis grading versus knee score

In group A, we observed eight patients (33.3%) with LKS, 15 (62.5%) with MKS and one (4.2%) with HKS. In group B, we observed 17 patients (81%) with LKS, four (19%) with MKS and none with HKS. In group C, we observed all patients (n = 21) with LKS. The highest number of patients with LKS were therefore in groups B and C and the radiographic findings are related to clinical pictures expressed by the KS score (P = 0.0001).

Kellgren and Lawrence osteoarthritis grading versus function score

In group A, we observed eight patients (33.3%) with LKS, 10 (41.7%) with MFS and six (25%) with HFS. In group B, we observed 13 patients (62%) with LFS, seven (33.3%) with MFS and one (4.7%) with HFS. In group C, we observed 15 patients (71.4%) with LFS, six (28.6%) with MFS and none with HFS. These data show that an increase of the OA radiographic severity corresponds to a decrease of the function score (P = 0.022).

Kellgren and Lawrence osteoarthritis grading versus age

In group A, we observed 14 young (58.3%) and 10 old (41.7%) patients. In group B, we observed eight young (38.1%) and 13 old (61.9%) patients. In group C, we observed two young (9.5%) and 19 old (90.5%) patients. So, it is more common to observe a medium to high OA radiographic grade in the population over 65 years old, and a low to medium in adults under the age of 65 years old (P = 0.0011).

Mutational analysis of osteoarthritis susceptibility genes

The OA patients were genotyped for the following polymorphisms associated with OS, such as SNPs and D repeats: FRZB rs288326 and rs7775, MATN3 rs77245812, ASPN D14, PTHR2 rs76758470, GDF5 rs143383, and DVWA rs11718863. Percentages of the wild type, heterozygote and homozygote genotypes for each polymorphism were calculated. We reported genotyping data of the three radiographic groups (A, B, C) and the number of individuals for each genotype (Table 5). In each group, deviations of Hardy–Weinberg equilibrium for all polymorphisms analyzed were not observed.

Table 5. Genetic analysis results

Kellgren and Lawrence osteoarthritis grading and genotype association analysis

To evaluate a potential association between genotypes, wild-type group or mutated (heterozygote and homozygote) group and the KL groups (A, B, C), the Mann–Whitney U test, the chi-square test and Fisher’s exact test were performed (Table 6). Analysis showed a statistically significant association between genotype and KL grade for the GDF5 rs143383 and the DVWA rs11718863 polymorphisms (P = 0.02 and P = 0.03, respectively). These results are in line with the study of Valdes and colleagues where GDF5 rs143383 and DVWA rs11718863 polymorphisms are consistently associated with the risk of knee OA in the Caucasian population [45], but to our knowledge this is the first study that reports the simultaneous presence of these two polymorphisms associated with KL in a European group. Unfortunately, concerning the other four OS SNPs, no genotype showed any significant association with KL data, as revealed by statistical analysis.

Table 6. Association between Kellgren and Lawrence osteoarthritis grading and genotype

Finally, it is possible to note in Table 6 that the DVWA rs11718863 polymorphism (genotype heterozygote + homozygote) is more represented in group C (55%), compared with the other two groups A (25%) and B (19%), suggesting that OS6 can be associated with a more severe OA radiographic grade.

Discussion

OA grading is commonly based on the radiographic classification of Kellgren and Lawrence [39] and can be supported by general or joint specific clinical and functional scores (for example, Western Ontario and McMaster Universities score, AKSS, and so forth) [46].

The aim of our pilot study was to update the knee OA grading with further clinical and genetic associated data, since OA is nowadays considered a polygenic and multifactorial disease [47].

We investigated the possible association between KL grade, clinical features (AKSS, age) and susceptibility polymorphisms to OA, such as FRZB (OS1A and OS1B), MATN3 (OS2), PTHR2 (OS3), ASPN D14 (OS4), GDF5 (OS5) and DVWA (OS6), in order to better define the grading of this disorder. As far as we know, this is the first study that simultaneously evaluates the association between clinical and radiographic data and the presence of the known six OA susceptibility SNPs reported in the Online Mendelian Inheritance in Man database.

In our cohort of 66 patients, a statistical association between the variables analyzed was observed (KL data versus: KS, FS, age). In particular, statistical association between KL grade versus KS and FS showed that KL group A can be associated with a medium clinical score, while KL group B and KL group C are related with low KS and FS. This suggests that a mild to severe OA radiographic grade is linked to severe clinical conditions and loss of articular function. In addition, association analysis between KL grade and age of patients confirms that severity of symptoms increases with age: the majority of our patients with KL grade A were 54 to 65 years old, and most of the patients with grade C were 65 to 86 years old.

Concerning the mutational analysis, we genotyped the patients for the abovementioned OS SNPs and our results revealed a significant statistical association between KL grade and GDF5 rs143383 (OS5) and DVWA rs11718863 (OS6) genetic alterations. The presence of these SNPs could thus contribute to better define OA grading and progression.

As described by several authors, the OS5 rs143383 SNP is still the most robustly replicated polymorphism associated with OA and the only locus successfully validated across diverse Asian and European populations. The OS5 rs143383 polymorphism, localized in the 5′-untranslated region, causes a decrease in the transcriptional activity of GDF5. This gene, expressed in the adult human articular cartilage, encodes growth differentiation factor 5, a bone morphogenic protein involved in the development, homeostasis and repairing of the bone, cartilage and other articular tissues [48].

To our knowledge, genetic contribution to radiographic severity in knee OA was previously described only by Valdes and colleagues in the UK OA patients and no information is available on individuals from other European countries. For the first time, our work analyzes the association between KL grade, clinical features (AKSS, age) and the rs143383 GDF5 SNP in Sicilian OA patients. In order to support this variant role as a progression marker in knee OA, we meta-analyzed our results with those previously published by Valdes and colleagues [49]. Nonetheless, the UK OA cohorts are much larger than our Sicilian cohort but all the populations are comparable for age and body mass index, whereas a slight reduction of the rs143383 T-allele percentage is observable in Sicilian cases with respect to the other UK populations (data not shown). Overall, the Sicilian OA patients’ odds ratio of 1.53 (confidence interval, 1.11 to 2.11) describes a positive association between rs143383 GDF5 and KL (KL ≤ 2 vs. KL > 2), a trend in line with Valdes and colleagues [48], supporting this variant as an OA progression marker.

The OS6 rs11718863 polymorphism is localized in an exonic region of the DWVA gene and causes a missense mutation with a consequent amino acidic substitution Tyr169Asn [50]. The DWVA protein, involved in cellular adhesion and protein to protein interactions, interacts with β-tubulin of microtubules and has an important role in the regulation of chondrocyte differentiation, protecting articulate joints from OA onset. In particular, the OS6 rs11718863 SNP induces a decreasing interaction between DVWA and β-tubulin [51,52]. The DVWA rs11718863 SNP is reported to be strongly associated with the risk of knee OA (odds ratio = 1.43, P < 7 × 10−11) and able to influence β-tubulin binding in Asian populations. Nevertheless, Valdes and colleagues showed an association between this genetic alteration and knee OA in the UK Nottingham total knee replacement OA group (P < 0.046) [32]. Meulenbelt and colleagues also described a moderately significant association in the UK sample (P = 0.046), not confirmed for other European countries. However, the higher risk allele frequency in the European samples highlighted the ethnic different penetrance of OS genes and, once again, the need to evaluate the alleles’ geographic distribution [38].

In addition, in our cohort of patients this genetic alteration was more represented in the KL group C (55%) compared with the other groups, KL group A (25%) and KL group B (19%), respectively. We therefore suggest that OS6 can be associated with a more severe OA radiographic grade, displaying its predictive role as OA marker progression.

In summary, according to our results we propose a model where we highlight the relationship between all data obtained, clinical and radiographic groups, and genetic factors. Our data pointed out that patients who belong to the KL group A (affected by mild OA) are younger and show a better joint function. In contrast, patients belonging to KL group B have a gradual progression of the disease. Finally, patients of KL group C have a more severe OA, they are older, they have poor function of the joint and most of them carry the rs11718863 OS6 polymorphism. Furthermore, according to our evidence, we show the DNA contribution, represented by rs143383 GDF5 and rs11718863 DVWA SNPs in defining the KL grade in OA prognosis and progression (Figure 1).

thumbnailFigure 1. Model of relationship between Kellgren and Lawrence osteoarthritis grading scale groups and clinical and genetic factors described in our study. FS, function score; KS, knee score; OA, osteoarthritis; OS, osteoarthritis susceptibility; SNP, single nucleotide polymorphism.

Conclusion

The statistically significant association between clinical, radiographic and genetic signs observed suggests the extension of the actual grading of knee OA based mainly on X-ray features.

This work highlights the importance of a multidisciplinary and translational medicine approach to study OA where clinical, radiological, and OS5 and OS6 SNP analysis could contribute to better define grading and progression of OA and to the development of new therapies.

Abbreviations

AKSS: American Knee Society Score; D: aspartic acid; FRZB: frizzled-related protein; FS: function score; GDF5: growth and differentiation factor 5; KL: Kellgren and Lawrence osteoarthritis grading scale; KS: knee score; OA: osteoarthritis; OS: osteoarthritis susceptibility; PTHR2: parathyroid hormone 2; SNP: single nucleotide polymorphism; Wnt: Wingless.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LM and VB were responsible for conception and design, genotyping experiments, data collection and analysis, interpretation, results elaboration, manuscript writing, and final approval of the manuscript. GIF and FPC were responsible for genotyping experiments, manuscript revision, and final approval of the manuscript. MS was responsible for data collection, patient enrolment, clinical and radiographical evaluations, manuscript writing, and final approval of the manuscript. FB was responsible for data collection, patient enrolment and clinical and radiographical evaluations, manuscript revision, financial support, and final approval of the manuscript. SC was responsible for data collection and analysis, manuscript revision, and final approval of the manuscript. MD’A was responsible for clinical evaluations, manuscript revision, and final approval of the manuscript. MCG and CM were responsible for result elaborations, financial support, manuscript revision, and final approval of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by PON01_01059 MRgFUS.

References

  1. Goldring MB, Goldring SR: Osteoarthritis.

    J Cell Physiol 2007, 213:626-634. PubMed Abstract | Publisher Full Text OpenURL

  2. Berenbaum F: Osteoarthritis as an inflammatory disease (osteoarthritis is not osteoarthrosis!).

    Osteoarthritis Cartilage 2013, 21:16-21. PubMed Abstract | Publisher Full Text OpenURL

  3. Blagojevic M, Jinks C, Jeffery A, Jordan KP: Risk factors for onset of osteoarthritis of the knee in older adults: a systematic review and meta-analysis.

    Osteoarthritis Cartilage 2010, 18:24-33. PubMed Abstract | Publisher Full Text OpenURL

  4. Arden N, Nevitt M: Osteoarthritis: epidemiology.

    Best Pract Res Clin Rheumatol 2006, 20:3-25. PubMed Abstract | Publisher Full Text OpenURL

  5. Suri P, Morgenroth DC, Hunter DJ: Epidemiology of osteoarthritis and associated comorbidities.

    PM R 2012, 4:S10-S19. PubMed Abstract | Publisher Full Text OpenURL

  6. Lane NE, Brandt K, Hawker G, Peeva E, Schreyer E, Tsuji W, Hochberg MC: OARSI-FDA initiative: defining the disease state of osteoarthritis.

    Osteoarthritis Cartilage 2011, 19:478-482. PubMed Abstract | Publisher Full Text OpenURL

  7. Pereira D, Peleteiro B, Araújo J, Branco J, Santos RA, Ramos E: The effect of osteoarthritis definition on prevalence and incidence estimates: a systematic review.

    Osteoarthritis Cartilage 2011, 19:1270-1285. PubMed Abstract | Publisher Full Text OpenURL

  8. Srikanth VK, Fryer JL, Zhai G, Winzenberg TM, Hosmer D, Jones G: A meta-analysis of sex differences prevalence, incidence and severity of osteoarthritis.

    Osteoarthritis Cartilage 2005, 13:769-781. PubMed Abstract | Publisher Full Text OpenURL

  9. Oliveria SA, Felson DT, Reed JI, Cirillo PA, Walker AM: Incidence of symptomatic hand, hip, and knee osteoarthritis among patients in a health maintenance organization.

    Arthritis Rheum 1995, 38:1134-1141. PubMed Abstract | Publisher Full Text OpenURL

  10. Belo JN, Berger MY, Reijman M, Koes BW, Bierma-Zeinstra SM: Prognostic factors of progression of osteoarthritis of the knee: a systematic review of observational studies.

    Arthritis Rheum 2007, 57:13-26. PubMed Abstract | Publisher Full Text OpenURL

  11. Braun HJ, Gold GE: Diagnosis of osteoarthritis: imaging.

    Bone 2012, 51:278-288. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Bijlsma JW, Berenbaum F, Lafeber FP: Osteoarthritis: an update with relevance for clinical practice.

    Lancet 2011, 377:2115-2126. PubMed Abstract | Publisher Full Text OpenURL

  13. Busija L, Bridgett L, Williams SR, Osborne RH, Buchbinder R, March L, Fransen M: Osteoarthritis.

    Best Pract Res Clin Rheumatol 2010, 24:757-768. PubMed Abstract | Publisher Full Text OpenURL

  14. Muraki S, Oka H, Akune T, Mabuchi A, En-yo Y, Yoshida M, Saika A, Suzuki T, Yoshida H, Ishibashi H, Yamamoto S, Nakamura K, Kawaguchi H, Yoshimura N: Prevalence of radiographic knee osteoarthritis and its association with knee pain in the elderly of Japanese population-based cohorts: the ROAD study.

    Osteoarthritis Cartilage 2009, 17:1137-1143. PubMed Abstract | Publisher Full Text OpenURL

  15. Cho HJ, Chang CB, Yoo JH, Kim SJ, Kim TK: Gender differences in the correlation between symptom and radiographic severity in patients with knee osteoarthritis.

    Clin Orthop Relat Res 2010, 468:1749-1758. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Cubukcu D, Sarsan A, Alkan H: Relationships between pain, function and radiographic findings in osteoarthritis of the knee: a cross-sectional study.

    Arthritis 2012, 2012:984060. PubMed Abstract | PubMed Central Full Text OpenURL

  17. Bedson J, Croft PR: The discordance between clinical and radiographic knee osteoarthritis: a systematic search and summary of the literature.

    BMC Musculoskelet Disord 2008, 9:116. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Hernández-Vaquero D, Fernández-Carreira JM: Relationship between radiological grading and clinical status in knee osteoarthritis.

    A multicentric study. BMC Musculoskelet Disord 2012, 13:194. BioMed Central Full Text OpenURL

  19. Van der Kraan PM: Osteoarthritis year 2012 in review: biology.

    Osteoarthritis Cartilage 2012, 20:1447-1450. PubMed Abstract | Publisher Full Text OpenURL

  20. Van Meurs JB, Uitterlinden AG: Osteoarthritis year 2012 in review: genetics and genomics.

    Osteoarthritis Cartilage 2012, 20:1470-1476. PubMed Abstract | Publisher Full Text OpenURL

  21. Loughlin J: Genome studies and linkage in primary osteoarthritis.

    Rheum Dis Clin N Am 2002, 28:95-109. Publisher Full Text OpenURL

  22. Karlsson C, Dehne T, Lindahl A, Brittberg M, Pruss A, Sittinger M, Ringe J: Genome-wide expression profiling reveals new candidate genes associated with osteoarthritis.

    Osteoarthr Cartil 2010, 18:581-592. PubMed Abstract | Publisher Full Text OpenURL

  23. Loughlin J: Genetics of osteoarthritis.

    Curr Opin Rheumatol 2011, 23:479-483. PubMed Abstract | Publisher Full Text OpenURL

  24. Rodriguez-Fontenla C, López-Golán Y, Calaza M, Pombo-Suarez M, Gómez-Reino JJ, González A: Genetic risk load and age at symptom onset of knee osteoarthritis.

    J Orthop Res 2012, 30:905-909. PubMed Abstract | Publisher Full Text OpenURL

  25. Forte GI, Piccione M, Scola L, Crivello A, Galfano C, Corsi MM, Chiappelli M, Candore G, Giuffrè M, Verna R, Licastro F, Corsello G, Caruso C, Lio D: Apolipoprotein E genotypic frequencies among down syndrome patients imply early unsuccessful aging for ApoE4 carriers.

    Rejuvenation Res 2007, 10:293-299. PubMed Abstract | Publisher Full Text OpenURL

  26. Bosco P, Guéant-Rodríguez RM, Anello G, Spada RS, Romano A, Caraci F, Ferri R, Guéant JL: Allele epsilon 4 of APOE is a stronger predictor of Alzheimer risk in Sicily than in continental South Italy.

    Neurosci Lett 2005, 388:168-172. PubMed Abstract | Publisher Full Text OpenURL

  27. Wilcken B, Bamforth F, Li Z, Zhu H, Ritvanen A, Ritvanen A, Renlund M, Stoll C, Alembik Y, Dott B, Czeizel AE, Gelman-Kohan Z, Scarano G, Bianca S, Ettore G, Tenconi R, Bellato S, Scala I, Mutchinick OM, López MA, de Walle H, Hofstra R, Joutchenko L, Kavteladze L, Bermejo E, Martínez-Frías ML, Gallagher M, Erickson JD, Vollset SE, Mastroiacovo P, et al.: Geographical and ethnic variation of the 677C > T allele of 5,10 methylenetetrahydrofolate reductase (MTHFR): findings from over 7000 newborns from 16 areas worldwide.

    J Med Genet 2003, 40:619-625. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Bravatà V, Stefano A, Cammarata FP, Minafra L, Russo G, Nicolosi S, Pulizzi S, Gelfi C, Gilardi MC, Messa C: Genotyping analysis and 18FDG uptake in breast cancer patients: a preliminary research.

    J Exp Clin Cancer Res 2013, 32:23. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Lane NE, Lian K, Nevitt MC, Zmuda JM, Lui L, Li J, Wang J, Fontecha M, Umblas N, Rosenbach M, de Leon P, Corr M: Frizzled-related protein variants are risk factors for hip osteoarthritis.

    Arthritis Rheum 2006, 54:1246-1254. PubMed Abstract | Publisher Full Text OpenURL

  30. Evangelou E, Chapman K, Meulenbelt I, Karassa FB, Loughlin J, Carr A, Doherty M, Doherty S, Gómez-Reino JJ, Gonzalez A, Halldorsson BV, Hauksson VB, Hofman A, Hart DJ, Ikegawa S, Ingvarsson T, Jiang Q, Jonsdottir I, Jonsson H, Kerkhof HJ, Kloppenburg M, Lane NE, Li J, Lories RJ, van Meurs JB, Näkki A, Nevitt MC, Rodriguez-Lopez J, Shi D, Slagboom PE, et al.: Large-scale analysis of association between GDF5 and FRZB variants and osteoarthritis of the hip, knee, and hand.

    Arthritis Rheum 2009, 60:1710-1721. PubMed Abstract | Publisher Full Text OpenURL

  31. Min JL, Meulenbelt I, Riyazi N, Kloppenburg M, Houwing-Duistermaat JJ, Seymour AB, Pols HA, van Duijn CM, Slagboom PE: Association of the frizzled-related protein gene with symptomatic osteoarthritis at multiple sites.

    Arthritis Rheum 2005, 52:1077-1080. PubMed Abstract | Publisher Full Text OpenURL

  32. Valdes AM, Loughlin J, Oene MV, Chapman K, Surdulescu GL, Doherty M, Spector TD: Sex and ethnic differences in the association of ASPN, CALM1, COL2A1, COMP, and FRZB with genetic susceptibility to osteoarthritis of the knee.

    Arthritis Rheum 2007, 56:137-146. PubMed Abstract | Publisher Full Text OpenURL

  33. Stefánsson SE, Jónsson H, Ingvarsson T, Manolescu I, Jónsson HH, Olafsdóttir G, Pálsdóttir E, Stefánsdóttir G, Sveinbjörnsdóttir G, Frigge ML, Kong A, Gulcher JR, Stefánsson K: Genomewide scan for hand osteoarthritis: a novel mutation in matrilin-3.

    Am J Hum Genet 2003, 72:1448-1459. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Kizawa H, Kou I, Iida A, Sudo A, Miyamoto Y, Fukuda A, Mabuchi A, Kotani A, Kawakami A, Yamamoto S, Uchida A, Nakamura K, Notoya K, Nakamura Y, Ikegawa S: An aspartic acid repeat polymorphism in asporin inhibits chondrogenesis and increases susceptibility to osteoarthritis.

    Nat Genet 2005, 37:138-144. PubMed Abstract | Publisher Full Text OpenURL

  35. Min JL, Meulenbelt I, Kloppenburg M, van Duijn CM, Slagboom PE: Mutation analysis of candidate genes within the 2q33.3 linkage area for familial early-onset generalised osteoarthritis.

    Eur J Hum Genet 2007, 15:791-799. PubMed Abstract | Publisher Full Text OpenURL

  36. Southam L, Rodriguez-Lopez J, Wilkins JM, Pombo-Suarez M, Sarah S, Gomez-Reino JJ, Chapman K, Gonzalez A, Loughlin J: An SNP in the 5′-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage.

    Hum Mol Genet 2007, 16:2226-2232. PubMed Abstract | Publisher Full Text OpenURL

  37. Chapman K, Takahashi A, Meulenbelt I, Watson C, Rodriguez-Lopez J, Egli R, Tsezou A, Malizos KN, Kloppenburg M, Shi D, Southam L, van der Breggen R, Donn R, Qin J, Doherty M, Slagboom PE, Wallis G, Kamatani N, Jiang Q, Gonzalez A, Loughlin J, Ikegawa S: A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5′ UTR of GDF5 with osteoarthritis susceptibility.

    Hum Mol Genet 2008, 17:1497-1504. PubMed Abstract | Publisher Full Text OpenURL

  38. Meulenbelt I, Chapman K, Dieguez-Gonzalez R, Shi D, Tsezou A, Dai J, Malizos KN, Kloppenburg M, Carr A, Nakajima M, van der Breggen R, Lakenberg N, Gomez-Reino JJ, Jiang Q, Ikegawa S, Gonzalez A, Loughlin J, Slagboom EP: Large replication study and meta-analyses of DVWA as an osteoarthritis susceptibility locus in European and Asian populations.

    Hum Mol Genet 2009, 18:1518-1523. PubMed Abstract | Publisher Full Text OpenURL

  39. Kellgren JH, Lawrence JS: Radiological assessment of osteoathrosis.

    Ann Rheum Dis 1957, 16:494-502. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Reijman M, Hazes JM, Pols HA, Bernsen RM, Koes BW, Bierma-Zeinstra SM: Validity and reliability of three definitions of hip osteoarthritis: cross sectional and longitudinal approach.

    Ann Rheum Dis 2004, 63:1427-1433. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Human Gene Mutation Database [http://www.hgmd.cf.ac.uk webcite]

  42. dbSNP Short Genetic Variations database [http://www.ncbi.nlm.nih.gov/snp webcite]

  43. Primer3 [http://fokker.wi.mit.edu/primer3 webcite]

  44. GraphPad InStat software version 3.05 [http://www.graphpad.com webcite]

  45. Valdes AM, Spector TD, Doherty S, Wheeler M, Hart DJ, Doherty M: Association of the DVWA and GDF5 polymorphisms with osteoarthritis in UK populations.

    Ann Rheum Dis 2009, 68:1916-1920. PubMed Abstract | Publisher Full Text OpenURL

  46. Insall J, Dorr LD, Scott RD, Scott DW: Rationale of the Knee Society Clinical Rating System.

    Clin Orthop 1989, 248:13-24. PubMed Abstract OpenURL

  47. Reynard LN, Loughlin J: The genetics and functional analysis of primary osteoarthritis susceptibility.

    Expert Rev Mol Med 2013, 15:e2. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Miyamoto Y, Mabuchi A, Shi D, Kubo T, Takatori Y, Saito S, Fujioka M, Sudo A, Uchida A, Yamamoto S, Ozaki K, Takigawa M, Tanaka T, Nakamura Y, Jiang Q, Ikegawa S: A functional polymorphism in the 5′-UTR of GDF5 is associated with susceptibility to osteoarthritis.

    Nat Genet 2007, 39:529-533. PubMed Abstract | Publisher Full Text OpenURL

  49. Valdes AM, Doherty S, Muir KR, Zhang W, Maciewicz RA, Wheeler M, Arden N, Cooper C, Doherty M: Genetic contribution to radiographic severity in osteoarthritis of the knee.

    Ann Rheum Dis 2012, 71:1537-1540. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Miyamoto Y, Shi D, Nakajima M, Ozaki K, Sudo A, Kotani A, Uchida A, Tanaka T, Fukui N, Tsunoda T, Takahashi A, Nakamura Y, Jiang Q, Ikegawa S, Miyamoto Y: Common variants in DVWA on chromosome 3p24.3 are associated with susceptibility to knee osteoarthritis.

    Nat Genet 2008, 40:994-998. PubMed Abstract | Publisher Full Text OpenURL

  51. Nakajima M, Miyamoto Y, Ikegawa S: Cloning and characterization of the osteoarthritis-associated gene DVWA.

    J Bone Miner Metab 2011, 29:300-308. PubMed Abstract | Publisher Full Text OpenURL

  52. Wagener R, Gara SK, Kobbe B, Paulsson M, Zaucke F: The knee osteoarthritis susceptibility locus DVWA on chromosome 3p24.3 is the 5′ part of the split COL6A4 gene.

    Matrix Biol 2009, 28:307-310. PubMed Abstract | Publisher Full Text OpenURL