Genome-wide DNA methylation study of hip and knee cartilage reveals embryonic organ and skeletal system morphogenesis as major pathways involved in osteoarthritis

BMC Musculoskeletal Disorders

Guangju Zhai²  Proton Rahman⁴  Guang Sun⁴  Roger Green¹  Andrew Furey³  Glynn Martin³  Patricia E. Harper¹  Ming Liu¹  Yuhua Zhang¹  Erfan Aref-Eshghi¹ 

[1] Discipline of Genetics, Faculty of Medicine, Memorial University of Newfoundland, St. John’s, NL, Canada;Department of Twin Research & Genetic Epidemiology, King’s College London, London, UK;Division of Orthopedics, Faculty of Medicine, Memorial University of Newfoundland, St. John’s, NL, Canada;Disicpline of Medicine, Faculty of Medicine, Memorial University of Newfoundland, St. John’s, NL, Canada

关键词: Knee;    Hip;    Skeletal morphogenesis;    DNA methylation;    Epigenome;    Osteoarthritis;   

Others  :  1230312
DOI  :  10.1186/s12891-015-0745-5

 received in 2015-06-15, accepted in 2015-10-02,  发布年份 2015

【 摘 要 】

Background

Evidence suggests that epigenetics plays a role in osteoarthrits (OA). The aim of the study was to describethe genome wide DNA methylation changes in hip and knee OA and identify novel genes and pathwaysinvolved in OA by comparing the DNA methylome of the hip and knee osteoarthritic cartilage tissues withthose of OA-free individuals.

Methods

Cartilage samples were collected from hip or knee joint replacement patients either due to primary OA or hip fractures as controls. DNA was extracted from the collected cartilage and assayed by Illumina Infinium HumanMethylation450 BeadChip array, which allows for the analysis of >480,000 CpG sites. Student T-test was conducted for each CpG site and those sites with at least 10 % methylation difference and a p value <0.0005 were defined as differentially methylated regions (DMRs) for OA. A sub-analysis was also done for hip and knee OA separately. DAVID v6.7 was used for the functional annotation clustering of the DMR genes. Clustering analysis was done using multiple dimensional scaling and hierarchical clustering methods.

Results

The study included 5 patients with hip OA, 6 patients with knee OA and 7 hip cartilage samples from OA-free individuals. The comparisons of hip, knee and combined hip/knee OA patients with controls resulted in 26, 72, and 103 DMRs, respectively. The comparison between hip and knee OA revealed 67 DMRs. The overall number of the sites after considering the overlaps was 239, among which 151 sites were annotated to 145 genes. One-fifth of these genes were reported in previous studies. The functional annotation clustering of the identified genes revealed clusters significantly enriched in skeletal system morphogenesis and development. The analysis revealed significant difference among OA and OA-free cartilage, but less different between hip OA and knee OA.

Conclusions

We found that a number of CpG sites and genes across the genome were differentially methylated in OA patients, a remarkable portion of which seem to be involved in potential etiologic mechanisms of OA. Genes involved in skeletal developmental pathways and embryonic organ morphogenesis may be a potential area for further OA studies.

【 授权许可】

   
2015 Aref-Eshghi et al.

【 图 表 】

Fig. 1.

Fig. 2.

Figure 2.

【 参考文献 】

• [1]Vos T, Flaxman AD, Naghavi M, Lozano R, Michaud C, Ezzati M et al.. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012; 380:2163-96.
• [2]Kean WF, Kean R, Buchanan WW. Osteoarthritis: symptoms, signs and source of pain. Inflammopharmacology. 2004; 12:3-31.
• [3]Guccione AA, Felson DT, Anderson JJ, Anthony JM, Zhang Y, Wilson PW et al.. The effects of specific medical conditions on the functional limitations of elders in the Framingham Study. Am J Public Health. 1994; 84:351-8.
• [4]Felson DT. An update on the pathogenesis and epidemiology of osteoarthritis. Radiol Clin North Am. 2004; 42:1-9.
• [5]Hunter DJ, March L, Sambrook PN. Knee osteoarthritis: the influence of environmental factors. Clin Exp Rheumatol. 2002; 20:93-100.
• [6]Martel-Pelletier J, Boileau C, Pelletier JP, Roughley PJ. Cartilage in normal and osteoarthritis conditions. Best Pract Res Clin Rheumatol. 2008; 22:351-84.
• [7]Zhai G, Aref-Eshghi E. Biomarkers for osteoarthritis: investigation, identification, and prognosis. Curr Biomark Find. 2012; 2:19-28.
• [8]Reynard LN, Loughlin J. Genetics and epigenetics of osteoarthritis. Maturitas. 2012; 71:200-4.
• [9]Rodriguez-Paredes M, Esteller M. Cancer epigenetics reaches mainstream oncology. Nat Med. 2011; 17:330-9.
• [10]Imagawa K, de Andrés MC, Hashimoto K, Itoi E, Otero M, Roach HI et al.. Association of reduced type IX collagen gene expression in human osteoarthritic chondrocytes with epigenetic silencing by DNA hypermethylation. Arthritis Rheumatol. 2014; 66:3040-51.
• [11]De Andres MC, Imagawa K, Hashimoto K, Gonzalez A, Roach HI, Goldring MB et al.. Loss of methylation in CpG sites in the NF-κB enhancer elements of inducible nitric oxide synthase is responsible for gene induction in human articular chondrocytes. Arthritis Rheum. 2013; 65:732-42.
• [12]Reynard LN, Bui C, Canty-Laird EG, Young DA, Loughlin J. Expression of the osteoarthritis-associated gene GDF5 is modulated epigenetically by DNA methylation. Hum Mol Genet. 2011; 20:3450-60.
• [13]Jeffries MA, Donica M, Baker LW, Stevenson ME, Annan AC, Humphrey MB et al.. Genome-wide DNA methylation study identifies significant epigenomic changes in osteoarthritic cartilage. Arthritis Rheumatol. 2014; 66:2804-15.
• [14]Rushton MD, Reynard LN, Barter MJ, Refaie R, Rankin KS, Young DA et al.. Characterization of the cartilage DNA methylome in knee and hip osteoarthritis. Arthritis Rheumatol. 2014; 66:2450-2460.
• [15]Aref-Eshghi E, Rahman P, Zhang H, Martin G, Furey A, Green R et al.. Attempt to replicate the published osteoarthritis-associated genetic variants in the Newfoundland& Labrador Population. J Orthopedics Rheumatol. 2014; 1:5.
• [16]Zhang W, Likhodii S, Aref-Eshghi E, Zhang Y, Harper PE, Randell E et al.. Relationship between blood plasma and synovial fluid metabolite concentrations in patients with osteoarthritis. J Rheumatol. 2015; 42:859-865.
• [17]Zhang W, Likhodii S, Zhang Y, Aref-Eshghi E, Harper PE, Randell E et al.. Classification of osteoarthritis phenotypes by metabolomics analysis. BMJ Open. 2014; 4:e006286.
• [18]Altman R, Alarcón G, Appelrouth D, Bloch D, Borenstein D, Brandt K et al.. The American College of Rheumatology criteria for the classification and reporting of osteoarthritis of the hip. Arthritis Rheum. 1991; 34:505-514.
• [19]Altman R, Asch E, Bloch D, Bole G, Borenstein D, Brandt K et al.. Development of the criteria for the classification and reporting of osteoarthritis, classification of osteoarthritis of the knee. Arthritis Rheum. 1986; 29:1039-1049.
• [20]Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. Evaluation of the Infinium Methylation 450 K technology. Epigenomics. 2011; 3:771-784.
• [21]Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM et al.. High density DNA methylation array with single CpG site resolution. Genomics. 2011; 98(4):288-295.
• [22]Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014; 30(10):1363-1369.
• [23]Hansen KD, Aryee M, Timp W. minfiData: Example data for the Illumina Methylation 450 k array. R package version 0.7, 1.
• [24]Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, Lin SM. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC bioinformatics. 2010; 11(1):587. BioMed Central Full Text
• [25]Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D et al.. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013; 29:189-96.
• [26]Laird PW. Principles and challenges of genome-wide DNA methylation analysis. Nat Rev Genet. 2010; 11:191-203.
• [27]Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009; 4:44-57.
• [28]den Hollander W, Ramos YF, Bos SD, Bomer N, van der Breggen R, Lakenberg N et al.. Knee and hip articular cartilage have distinct epigenomic landscapes: implications for future cartilage regeneration approaches. Ann Rheum Dis. 2014; 73:2208-12.
• [29]Fernández-Tajes J, Soto-Hermida A, Vázquez-Mosquera ME, Cortés-Pereira E, Mosquera A, Fernández-Moreno M et al.. Genome-wide DNA methylation analysis of articular chondrocytes reveals a cluster of osteoarthritic patients. Ann Rheum Dis. 2014; 73:668-77.
• [30]Wang X, Manner PA, Horner A, Shum L, Tuan RS, Nuckolls GH. Regulation of MMP-13 expression by RUNX2 and FGF2 in osteoarthritic cartilage. Osteoarthritis Cartilage. 2004; 12:963-73.
• [31]Ma XJ, Guo HH, Hao SW, Sun SX, Yang XC, Yu B et al.. Association of single nucleotide polymorphisms (SNPs) in leptin receptor gene with knee osteoarthritis in the Ningxia Hui population. Yi Chuan. 2013; 35:359-64.
• [32]Ye XL, Lu CF. Association of polymorphisms in the leptin and leptin receptor genes with inflammatory mediators in patients with osteoporosis. Endocrine. 2013; 44:481-8.
• [33]Pettersson U, Albagha OM, Mirolo M, Taranta A, Frattini A, McGuigan FE et al.. Polymorphisms of the CLCN7 gene are associated with BMD in women. J Bone Miner Res. 2005; 20:1960-7.
• [34]Elfassihi L, Giroux S, Bureau A, Laflamme N, Cole DE, Rousseau F. Association with replication between estrogen-related receptor gamma (ESRRgamma) polymorphisms and bone phenotypes in women of European ancestry. J Bone Miner Res. 2010; 25:901-11.
• [35]Ettinger R, Kuchen S, Lipsky PE. Interleukin 21 as a target of intervention in autoimmune disease. Ann Rheum Dis. 2008; 67 Suppl 3:iii83-6.
• [36]Windsor RE, Strauss SJ, Kallis C, Wood NE, Whelan JS. Germline genetic polymorphisms may influence chemotherapy response and disease outcome in osteosarcoma: a pilot study. Cancer. 2012; 118:1856-67.
• [37]Fei Q, Wu ZH, Yuan SM, Wang H, Zhou X, Liu Z et al.. Association of PAX1 gene polymorphisms with susceptibility to congenital scoliosis in Chinese Han population. Zhonghua Yi Xue Za Zhi. 2008; 88:2597-602.
• [38]Moazedi-Fuerst FC, Hofner M, Gruber G, Weinhaeusel A, Stradner MH, Angerer H et al.. Epigenetic differences in human cartilage between mild and severe OA. J Orthop Res. 2014; 32:1636-45.
• [39]Delgado-Calle J, Fernández AF, Sainz J, Zarrabeitia MT, Sañudo C, García-Renedo R et al.. Genome-wide profiling of bone reveals differentially methylated regions in osteoporosis and osteoarthritis. Arthritis Rheum. 2013; 65:197-205.
• [40]Wagner EF, Karsenty G. Genetic control of skeletal development. Curr Opin Genet Dev. 2001; 11:527-532.
• [41]Lane NE, Lin P, Christiansen L, Gore LR, Williams EN, Hochberg MC et al.. Association of mild acetabular dysplasia with an increased risk of incident hip osteoarthritis in elderly white women: the Study of Osteoporotic Fractures. Arthritis Rheum. 2000; 43:400-404.
• [42]Shepstone L, Rogers J, Kirwan JR, Silverman BW. Shape of the intercondylar notch of the human femur: a comparison of osteoarthritic and non-osteoarthritic bones from a skeletal sample. Ann Rheum Dis. 2001; 60:968-973.
• [43]Valdes AM, Spector TD. The contribution of genes to osteoarthritis. Rheum Dis Clin North Am. 2008; 34:581-603.