- Research article
- Open Access
A genetic association study between growth differentiation factor 5 (GDF 5) polymorphism and knee osteoarthritis in Thai population
Journal of Orthopaedic Surgery and Researchvolume 6, Article number: 47 (2011)
Osteoarthritis (OA) is a multi-factorial disease and genetic factor is one of the important etiologic risk factors. Various genetic polymorphisms have been elucidated that they might be associated with OA. Recently, several studies have shown an association between Growth Differentiation Factor 5(GDF5) polymorphism and knee OA. However, the role of genetic predisposing factor in each ethnic group cannot be replicated to all, with conflicting data in the literatures. Therefore, the aim of this study was to investigate the association between GDF5 polymorphism and knee OA in Thai population.
Materials and Methods
One hundred and ninety three patients aged 54-88 years who attended Ramathibodi Hospital were enrolled. Ninety cases with knee OA according to American College of Rheumatology criteria and one hundred and three cases in control group gave informed consent. Blood sample (5 ml) were collected for identification of GDF5 (rs143383) single nucleotide polymorphism by PCR/RFLP according to a standard protocol. This study protocol was approved by the Ethics Committee on human experimentation of Ramathibodi Hospital Faculty of Medicine, Mahidol University. Odds ratios (OR) and 95% confidence intervals were calculated for the risk of knee OA by genotype (TT, TC and CC) and allele (T/C) analyses.
The baseline characteristics between two groups including job, smoking and activity were not different, except age and BMI. The entire cases and controls were in Hardy-Weinberg equilibrium (p > 0.05). The OA knee group (n = 90) had genotypic figure which has shown by TT 42.2% (n = 38), TC 45.6% (n = 41) and CC 12% (n = 11), whereas the control group (n = 103) revealed TT 32% (n = 33), TC 45.6% (n = 47), and CC 22.3% (n = 23), respectively. Genotypic TT increased risk of knee OA as compared to CC [OR = 2.41 (P = 0.04, 95%CI = 1.02-5.67)]. In the allele analysis, the T allele was found to be significantly associated with knee OA [OR = 1.53 (P = 0.043, 95%CI = 1.01-2.30)].
These data suggested that GDF5 polymorphism has an association with knee OA in Thai ethnic. This finding also supports the hypothesis that OA has an important genetic component in its etiology, and GDF5 protein might play important role in the pathophysiology of the disease.
It is widely believed that osteoarthritis develops from an imbalance between anabolic and catabolic processes or homeostasis of cartilage metabolism [1, 2]. The etiology of this disease is related to genetic association . Recently, several studies have demonstrated the polymorphism in many genes which might be related to the pathogenesis of osteoarthritis [4–14]. Growth Differentiation Factor 5 or GDF5 gene regulates the expression of the GDF5 protein which is closely related to BMP and is a member of TGF-beta superfamily . It has a role in regulation of the chondrogenesis and the defect in this gene might be correlated to the abnormal joint development . It has been reported in animal study that GDF5 knockout mice develops knee joint anomaly . Moreover, it has been shown that he polymorphism in GDF5 gene is related with low expression of the GDF5 protein in knee joint .
The large scale analysis has shown that the association of this polymorphism (rs143383) in promoter area might be a risk factor of the osteoarthritis . In addition, studies of this genetic variant in China and Japan have shown the association of T allele and knee OA . However, there is a report from Greece which had an inconsistent result and found no an association . The genetic susceptibility of the disease in different ethnic cannot be applied to others because each ethnic has a different genetic background. There is a gap of information about this polymorphism in Thai population, therefore the objective of this study is to determine the association of the SNP rs143383 in GDF5 gene and knee OA in Thai population.
Patients and Methods
Our study was approved by the Ethics Committee of the Faculty of Medicine, Ramathibodi Hospital, Mahidol University, Bangkok, Thailand. All patients recruited in the study were Thai by nationality and had ancestors settled in Thailand for at least three generations. In total, 90 patients with knee OA who underwent total knee arthroplasty (TKA) and 103 patients without knee OA were enrolled from Department of Orthopedics, Ramathibodi Hospital. Informed consent was performed after the purpose of the research project and it had been clearly explained to the patients. The diagnosis of knee OA was based on the American College of Rheumatology criteria . Both OA and control groups were interviewed to obtain demographic data and all of established risk factors. Thereafter, standard weight-bearing antero-posterior and lateral view of knee radiographs were taken to confirm the diagnosis of OA by Kellgren and Lawrence scores (KL scores) .
PCR-RFLP for BsiEI restriction site was used for SNP in GDF-5 identification
Firstly, 5 ml peripheral blood sample was collected from patient using ethylenediamine tetraacetic acid as an anticoagulant and processed for SNP analysis. Genomic DNA was extracted from buffy coat leukocytes using the standard phenol-chloroform method. PCR primers to amplify the promoter area of GDF5 gene were designed by the Primer-3 web-based tool : GATTTTTTCTGAGCACCTGCAGG (forward) and GTGTGTGTTTGTATCCAG (reverse). 50 μl PCR mixture contained 100 ng of genomic DNA, 20 pmol of each primer, 0.2 μ M of each dNTP, 1 unit of Taq DNA polymerase (AmpliTaq®, Applied Biosystem, Foster City, CA), 3.0 mM MgCl2 in 10 × PCR buffer containing 10 mmol of Tri-HCl pH 9.0, 10 mmol KCl and 0.1% Triton X-100 (Invitrogen, Carlsbad, CA). PCR reaction was started with an initial denaturation at 95°C for 5 min, followed by 35 cycles of amplification in a thermocycler (PCR Sprint, Thermofisher, Waltham, MA) with denaturation at 94°C for 1 min, annealing at 58°C for 1 min, extension at 72°C for 1 min, and final extension at 72°C for 10 min. Then, 10 μL of PCR product was incubated at 37°C with 3 units of BsiEI for 4 hours under the manufacturer's recommended conditions (New England Biolabs, Ipswich, MA). The digested product was electrophoresed on 2% agarose gel with ethidium bromide staining before being visualized on a UV transilluminator (Figure 1). The expected fragment length was 104 and 230 bp in CC, 104, 230, and 344 bp in TC, and 344 bp in TT genotypes, respectively.
Analysis of demographic data was performed by Excel 2007 (Microsoft®Excel®). The unpaired t-test was used for continuous data and Chi-square was used for categorical data. Allele frequencies, Odds ratio and the probability for Hardy-Weinberg Equilibrium were estimated and analyzed as explained by following website; http://ihg.gsf.de/cgi-bin/hw/hwa1.pl.
In our study, the average age of patients in the knee OA group was significant older than in the control. BMI which is one of the risk factor of knee OA was also slightly higher in the OA patient group. However, other possible risk factors which might be related to OA as well as patient life style were comparable between these two groups. The baseline characteristics of the patients are shown in Table 1. The prevalence of allele T and C in our sample was normally distributed according to the Hardy-Weinberg Equilibrium (p-value < 0.05). The comparison of the number and odd ratio of genotype TT, TC, TT+TC with CC of GDF5 polymorphism (rs143383) in the promoter area between case and control is demonstrated in Table 2. When analyzing by allele (shown in table 3), the odds ratio of T allele was 1.53 (95%C.I. = 1.01-2.31), which delineated the increase in risk to develop knee OA. Our result shown that this polymorphism, T/C, is inherited by autosomal recessive manner as the TT genotype increases risk of the disease significantly (OR = 2.41, 95%C.I. = 1.02-5.67), whereas the TC has no significant difference.
Discussion and Conclusion
Our study has shown the statistically significant association between polymorphism in promoter of GDF5 gene (rs143383) and knee osteoarthritis in Thai population. We found that T allele in GDF5 polymorphism was a significant risk factor for knee osteoarthritis with Odds ratio 1.53 (95%C.I. = 1.01-2.31). When analysed by genotype, it was found that the TT genotype increased risk of knee OA, whereas TC has not shown any correlation. Regarding these finding, it suggests that this polymorphic gene might be expressed as the autosomal recessive type. The GDF5 gene is located on the chromosome 20q11.2 and it regulates the expression of GDF5 protein. It is categorised in the class of bone morphogenetic protein (BMP) . It involves in development of bone and cartilage, particularly in endochondral ossification process [24, 25]. The mutation of GDF5 associates with generalised osteoarthritis and skeletal-related congenital diseases. Acromesomelic chondrodysplasia; Grebe type and Hunter-Thompson type is autosomal recessive form which is rare hereditary skeletal disorders. They are presented by short status, abnormal limbs development. Brachydactyly type A1, A2 and C and symphalangism proximal syndrome are inherited malformation presented with abnormal morphology of hand and finger. Multiple synostoses syndrome type 2 is an autosomal dominant condition characterised by progressive joint fusions and progressive conductive deafness. Du Pan syndrome is a rare autosomal recessive; the patient presents with absence of the fibulae and severe acromesomelic limb shortening with small, non-functional toes. These abnormal conditions have been reported that they have a mutation in GDF5 gene, thus they support that GDF5 plays a role in skeletal development [26, 27]. Recently, it has been reported that GDF5 deficiency mice had the delay fracture healing process that impaired the cartilaginous matrix deposition in the callus and reduced callus cross-sectional area .
The functional study of this polymorphism has been demonstrated; T allele in rs143383 was associated with the decrease of GFD5 molecule expression and might increase susceptibility to osteoarthritis . Moreover, it has been report that the differential binding of deformed epidermal autoregulatory factor 1 (DEAF-1) can modulate the expression of GDF5 via this polymorphism . It is believed that GDF5 plays role in the regulation of the chondrogenic cell growth and differentiation. These evidences support the function of GDF5 gene which might indicate the importance of this polymorphism in osteoarthritis etiology. As the result, many scientists have been working for study the association of this polymorphism in GDF5 gene. Recently, a large scale association study between GDF5 gene and osteoarthritis has revealed the genetic susceptibility of polymorphism in the GDF5; rs143383 with odds ratio 1.12 (95%CI 0.99-1.31) . This study included the data from diverse ethnics which have been reported in the recent literatures. However, the magnitude of this association is smaller than the Asian ethnic, it has been reported from Japan that odds ratio is 1.79 (P = 1.18 × 10-13) for GDF5 polymorphism with per-risk allele (T) odds ratio by Miyamoto et al . The magnitude of the association of GDF5 polymorphism is different between Caucasian ethnic and Asian ethnic. Furthermore, the association of this polymorphism from Greek and Korean study cannot demonstrate the association of GDF5 polymorphism [20, 30]. We have also reported the polymorphism in ESR1, which is associated with knee OA in Korean population. However, our result did not show any statistically significant difference . Thus, these evidences can imply that the association of ethic and genetic susceptibility in osteoarthritis might not be consistent among the different populations and cannot be applied to others.
In addition, there is a report that GDF5 polymorphism (rs143383) also predisposes to Lumbar disc degeneration in women. Lumbar disc degeneration which is defined by disc space narrowing and the presence of osteophytes significantly associates with GDF5 polymorphism from cohorts from Northern Europe in women, with an odds ratio (OR) of 1.72 (95% CI 1.15-2.57) . The genome-wide association study from Finland and Sardinia shows the common variants in GDF5 contribute to height difference . The GDF5 gene might be involved skeletal growth and development and GDF5 variants might play role in pathogenesis of bone and cartilage diseases.
Although the sample size in our study was small when compared to previous studies, these limited samples were sufficient for demonstrating the statistically significant association of rs143383 in GDF5 core promotor area. Furthermore, the population in our study was in homogeneity and the distribution between the case and control groups was in the Hardy Weinberg's equilibrium. Therefore, our study can represent the significant susceptibility of GDF5 in knee osteoarthritis in Thai ethnic. Our findings have emphasized this association that the susceptibility of the polymorphism in GDF5 among the Asian population because the magnitude of association is close to Japan and Chinese population. It is possible that Thai population ethnic is believed to be close to Chinese ancestry. Finally, we decide to leave some suggestions in genetic susceptibility study. Firstly, it is important to study the genetic susceptibility in common disease as this might be valuable for further investigation in order to understand the disease at molecular level or even apply to use for disease screening. More important, the association between polymorphism and disease might not similar in different ethics; therefore the genetic susceptibility from different ethics is required and the meta-analysis should be conducted in the adjacent area or similar genetic background. Lastly, other modern technology such as microarray technique should be employed for investigation of genetic susceptibility in osteoarthritis.
Felson DT: Osteoarthritis of the knee. N Engl J Med. 2006, 354: 841-8. 10.1056/NEJMcp051726.
Sandell LJ, Aigner T: Articular cartilage and changes in arthritis. An introduction: cell biology of osteoarthritis. Arthritis Res. 2001, 3: 107-13. 10.1186/ar148.
Lanyon P, Muir K, Doherty S, Doherty M: Assessment of a genetic contribution to osteoarthritis of the hip: sibling study. BMJ. 2000, 321: 1179-83. 10.1136/bmj.321.7270.1179.
Loughlin J, Sinsheimer JS, Mustafa Z, Carr AJ, Clipsham K, Bloomfield VA, Chitnavis J, Bailey A, Sykes B, Chapman : Association analysis of the vitamin D receptor gene, the type I collagen gene COL1A1, and the estrogen receptor gene in idiopathic osteoarthritis. J Rheumatol. 2000, 27: 779-84.
Spector TD, MacGregor AJ: Risk factors for osteoarthritis: genetics. Osteoarthritis Cartilage. 2004, 12 (Suppl A): 39-44.
Smith AJ, Keen LJ, Billingham MJ, Perry MJ, Elson CJ, Kirwan JR, Sims JE, Doherty M, Spector TD, Bidwell JL: Extended haplotypes and linkage disequilibrium in the IL1R1-IL1A-IL1B-IL1RN gene cluster: association with knee osteoarthritis. Genes Immun. 2004, 5: 451-60. 10.1038/sj.gene.6364107.
Aigner T, Dudhia J: Genomics of osteoarthritis. Curr Opin Rheumatol. 2003, 15: 634-40. 10.1097/00002281-200309000-00019.
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. 10.1038/ng1496.
Mabuchi A, Ikeda T, Fukuda A, Koshizuka Y, Hiraoka H, Miyoshi K, Haga N, Kawaguchi H, Kawakami A, Yamamoto S, Takatori Y, Nakamura K, Ikegawa S: Identification of sequence polymorphisms of the COMP (cartilage oligomeric matrix protein) gene and association study in osteoarthrosis of the knee and hip joints. J Hum Genet. 2001, 46: 456-462. 10.1007/s100380170045.
Mototani H, Mabuchi A, Saito S, Fujioka M, Iida A, Takatori Y, Kotani A, Kubo T, Nakamura K, Sekine A, Murakami Y, Tsunoda T, Notoya K, Nakamura Y, Ikegawa S: A functional single nucleotide polymorphism in the core promoter region of CALM1 is associated with hip osteoarthritis in Japanese. Hum Mol Genet. 2005, 14: 1009-1017. 10.1093/hmg/ddi093.
Ikeda T, Mabuchi A, Fukuda A, Hiraoka H, Kawakami A, Yamamoto S, Machida H, Takatori Y, Kawaguchi H, Nakamura K, Ikegawa S: Identification of sequence polymorphisms in two sulfation-relatedgenes, PAPSS2 and SLC26A2, and an association analysis with knee osteoarthritis. J Hum Genet. 2001, 46: 538-543. 10.1007/s100380170036.
Jiang Q, Shi D, Yi L, Ikegawa S, Wang Y, Nakamura T, Qiao D, Liu C, Dai J: Replication of the association of the aspartic acid repeat polymorphism in the asporin gene with knee-osteoarthritis susceptibility in Han Chinese. J Hum Genet. 2006, 51: 1068-1072. 10.1007/s10038-006-0065-6.
Valdes AM, Spector TD: The contribution of genes to osteoarthritis. Rheum Dis Clin North Am. 2008, 34: 581-603. 10.1016/j.rdc.2008.04.008.
Valdes AM, Loughlin J, Timms KM, van Meurs JJ, Southam L, Wilson SG, Doherty S, Lories RJ, Luyten FP, Gutin A, Abkevich V, Ge D, Hofman A, Uitterlinden AG, Hart DJ, Zhang F, Zhai G, Egli RJ, Doherty M, Lanchbury J, Spector TD: Genome-wide association scan identifies aprostaglandin-endoperoxide synthase 2 variant involved in risk of knee osteoarthritis. Am J Hum Genet. 2008, 82: 1231-4. 10.1016/j.ajhg.2008.04.006.
Buxton P, Edwards C, Archer CW, Francis-West P: Growth/differentiation factor-5 (GDF-5) and skeletal development. J Bone Joint Surg Am. 2001, 83-A (Suppl 1): 23-30.
Harada M, Takahara M, Zhe P, Otsuji M, Iuchi Y, Takagi M, Ogino T: Developmental failure of the intra-articular ligaments in mice with absence of growth differentiation factor 5. Osteoarthritis Cartilage. 2007, 15: 468-74. 10.1016/j.joca.2006.09.003.
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-33. 10.1038/2005.
Evangelou E, Chapman K, Meulenbelt I: Large-scale analysis of association between GDF5 (rs143383) and FRZB (rs7775 and rs288326) variants and hip, knee and hand osteoarthritis. Arthritis Rheum. 2009, 60: 1710-1721. 10.1002/art.24524.
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-504. 10.1093/hmg/ddn038.
Tsezou A, Satra M, Oikonomou P, Bargiotas K, Malizos KN: The growth differentiation factor 5 (GDF5) core promoter polymorphism is not associated with knee osteoarthritis in the Greek population. J Orthop Res. 2008, 26: 136-40. 10.1002/jor.20464.
Altman R, Asch E, Bloch D, Bole G, Borenstein D, Brandt K, Christry W, Cooke TD, Greenwald R, Hochberg M, Howell D, Kaplan W, Koopman T, Meenan R, Mqskowitz R, Murphy W, Rothschild B, Segal M, Sokoloff L, Wolfe F: Development of criteria for the classification and reporting of osteoarthritis: classification of osteoarthritis of the knee. Diagnostic and Therapeutic Criteria Committee of the American Rheumatism Association. Arthritis Rheum. 1986, 29: 1039-49. 10.1002/art.1780290816.
Kellgren JK, Lawrence JS: Radiological assessment of osteoarthritis. Ann Rheum Dis. 1957, 16: 494-501. 10.1136/ard.16.4.494.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology 2000. Edited by: Krawetz S, Misener S Totowa. 2000, NJ: Humana Press, 365-86.
Francis-West PH, Abdelfattah A, Chen P, Allen C, Parish J, Ladher R, Allen S, MacPherson S, Luyten FP, Archer CW: Mechanisms of GDF-5 action during skeletal development. Development. 1999, 126 (6): 1305-15.
Hatakeyama Y, Tuan RS, Shum L: Distinct functions of BMP4 and GDF5 in the regulation of chondrogenesis. J Cell Biochem. 2004, 91: 1204-17. 10.1002/jcb.20019.
Thomas JT, Lin K, Nandedkar M, Camargo M, Cervenka J, Luyten FP: A human chondrodysplasia due to a mutation in a TGF-b superfamily member. Nat Genet. 1996, 12 (3): 315-7. 10.1038/ng0396-315.
GDF5 growth differentiation factor 5 [Homo sapiens ]. http://www.ncbi.nlm.nih.gov/sites/entrez?Db=gene%26Cmd=ShowDetailView%26TermToSearch=8200
Coleman CM, Scheremeta BH, Boyce AT, Mauck RL, Tuan RS: Delayed fracture Healing in growth differentiation factor 5-deficient mice: A pilot study. Clin Orthop Relat Res.
Egli RJ, Southam L, Wilkins JM, Lorenzen I, Pombo-Suarez M, Gonzalez A, Carr A, Chapman K, Loughlin J: Functional analysis of the osteoarthritis susceptibility-associated GDF5 regulatory polymorphism. Arthritis Rheum. 2009, 60: 2055-64. 10.1002/art.24616.
Cao Z, Lee HS, Song JH, Yoon JW, Park YK, Nam SW, Lee JY, Park WS: Growth Differentiation Factor 5 (GDF5) core promoter polymorphism is not associated with susceptibility to osteoarthritis of the knee in the Korean population. Korean J Pathol. 2010, 44: 404-9. 10.4132/KoreanJPathol.2010.44.4.404.
Tawonsawatruk T, Trachoo O, Channoom T, Sura T, Eu-ahsunthornwattana J, Woratanarat P, Wajanavisit W: Association of Estrogen Receptor-alpha single-nucleotide polymorphism (codon 594 G/A) and Thai patients affected by knee osteoarthritis. J Med Assoc Thai. 2009, 92 (Suppl 6): 45-50.
Williams FM, Popham M, Hart DJ, de Schepper E, Bierma-Zeinstra S, Hofman A, Uitterlinden AG, Arden NK, Cooper C, Spector TD, Valdes AM, van Meurs J: GDF5 single-nucleotide polymorphism rs143383 is associated with lumbar disc degeneration in Northern European women. Arthritis Rheum. 2011, 63 (3): 708-12. 10.1002/art.30169.
Sanna S, Jackson AU, Nagaraja R, Willer CJ, Chen WM, Bonnycastle LL: Common variants in the GDF5-UQCC region are associated with variation in human height. Nat Genet. 2008, 40 (2): 198-203. 10.1038/ng.74.
Acknowledgements and Funding
This work was supported by grants from the biomedical research fund (#50039) of Faculty of Medicine, Ramathibodi Hospital, Bangkok, Thailand.
The authors declare that they have no competing interests.
TT conducted the acquisition, analysis and interpretation of data, carried out the molecular genetic studies and drafted the manuscript. TC collected data. SP carried out the genetic study; RFLP. OT participated in the design of the study and performed the statistical analysis. TS and WW participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.