GDF5 encodes an extracellular signalling molecule that is essential for normal skeletal development. The rs144383 C to T SNP located in the 5ʹUTR of this gene is functional and has a pleiotropic effect on the musculoskeletal system, being a risk factor for knee-osteoarthritis (OA), congenital hip dysplasia, lumbar disc degeneration and Achilles tendon pathology. rs143383 exerts a joint-wide effect on GDF5 expression, with expression of the OA-associated T allele being significantly reduced relative to the C allele, termed allelic expression imbalance. We have previously reported that the GDF5 locus is subject to DNA methylation and that allelic imbalance of rs143383 is mediated by SP1, SP3 and DEAF1 transcriptional repressors. In this study, we have assayed GDF5 methylation in normal and osteoarthritic cartilage, and investigated the effect of methylation on the allelic imbalance of rs143383. We observed demethylation of the GDF5 5ʹUTR in OA knee cartilage relative to both OA (p = 0.009) and non-OA (p = 0.001) hip cartilage, with the most significant demethylation observed at the highly conserved +37 CpG site located 4 bp upstream of rs143383. Methylation modulates the level and direction of allelic imbalance of rs143383, with methylation of the +37 CpG dinucleotide within the SP1/SP3 binding site having an allele-specific effect on SP1 and SP3 binding. Furthermore, methylation attenuated the repressive effects of SP1, SP3 and DEAF1 on GDF5 promoter activity. This data suggest that the differential methylation of the +37 CpG site between osteoarthritic hip and knee cartilage may be responsible for the knee-specific effect of rs143383 on OA susceptibility.
Electronic supplementary material
The online version of this article (doi:10.1007/s00439-014-1447-z) contains supplementary material, which is available to authorized users.
An osteoarthritis (OA) susceptibility locus has been mapped to chromosome 3p21, to a region of high linkage disequilibrium encompassing twelve genes. Six of these genes are expressed in joint tissues and we therefore assessed whether any of the six were subject to cis-acting regulatory polymorphisms active in these tissues and which could therefore account for the association signal.
We measured allelic expression using pyrosequencing assays that can distinguish mRNA output from each allele of a transcript single nucleotide polymorphism. We assessed RNA extracted from the cartilage and other joint tissues of OA patients who had undergone elective joint replacement surgery. A two-tailed Mann–Whitney exact test was used to test the significance of any allelic differences.
GNL3 and SPCS1 demonstrated significant allelic expression imbalance (AEI) in OA cartilage (GNL3, mean AEI = 1.04, p = 0.0002; SPCS1, mean AEI = 1.07, p < 0.0001). Similar results were observed in other tissues. Expression of the OA-associated allele was lower than that of the non-associated allele for both genes.
cis-acting regulatory polymorphisms acting on GNL3 and SPCS1 contribute to the OA association signal at chromosome 3p21, and these genes therefore merit further investigation.
Osteoarthritis; Genetic risk; Single nucleotide polymorphism; Allelic imbalance; Linkage disequilibrium
A genetic deficit mediated by SNP rs143383 that leads to reduced expression of GDF5 is strongly associated with large-joint osteoarthritis. We speculated that this deficit could be attenuated by the application of exogenous GDF5 protein and as a first step we have assessed what effect such application has on primary osteoarthritis chondrocyte gene expression. Chondrocytes harvested from cartilage of osteoarthritic patients who had undergone joint replacement were cultured with wildtype recombinant mouse and human GDF5 protein. We also studied variants of GDF5, one that has a higher affinity for the receptor BMPR-IA and one that is insensitive to the GDF5 antagonist noggin. As a positive control, chondrocytes were treated with TGF-β1. Chondrocytes were cultured in monolayer and micromass and the expression of genes coding for catabolic and anabolic proteins of cartilage were measured by quantitative PCR. The expression of the GDF5 receptor genes and the presence of their protein products was confirmed and the ability of GDF5 signal to translocate to the nucleus was demonstrated by the activation of a luciferase reporter construct. The capacity of GDF5 to elicit an intracellular signal in chondrocytes was demonstrated by the phosphorylation of intracellular Smads. Chondrocytes cultured with TGF-β1 demonstrated a consistent down regulation of MMP1, MMP13 and a consistent upregulation of TIMP1 and COL2A1 with both culture techniques. In contrast, chondrocytes cultured with wildtype GDF5, or its variants, did not show any consistent response, irrespective of the culture technique used. Our results show that osteoarthritis chondrocytes do not respond in a predictable manner to culture with exogenous GDF5. This may be a cause or a consequence of the osteoarthritis disease process and will need to be surmounted if treatment with exogenous GDF5 is to be advanced as a potential means to overcome the genetic deficit conferring osteoarthritis susceptibility at this gene.
rs143383 is a C to T transition SNP located in the 5′untranslated region (5′UTR) of the growth differentiation factor 5 gene GDF5. The T allele of the SNP is associated with increased risk of osteoarthritis (OA) in Europeans and in Asians. This susceptibility is mediated by the T allele producing less GDF5 transcript relative to the C allele, a phenomenon known as differential allelic expression (DAE). The aim of this study was to identify trans-acting factors that bind to rs143383 and which regulate this GDF5 DAE. Protein binding to the gene was investigated by two experimental approaches: 1) competition and supershift electrophoretic mobility shift assays (EMSAs) and 2) an oligonucleotide pull down assay followed by quantitative mass spectrometry. Binding was then confirmed in vivo by chromatin immunoprecipitation (ChIP), and the functional effects of candidate proteins investigated by RNA interference (RNAi) and over expression. Using these approaches the trans-acting factors Sp1, Sp3, P15, and DEAF-1 were identified as interacting with the GDF5 5′UTR. Knockdown and over expression of the factors demonstrated that Sp1, Sp3, and DEAF-1 are repressors of GDF5 expression. Depletion of DEAF-1 modulated the DAE of GDF5 and this differential allelic effect was confirmed following over expression, with the rs143383 T allele being repressed to a significantly greater extent than the rs143383 C allele. In combination, Sp1 and DEAF-1 had the greatest repressive activity. In conclusion, we have identified four trans-acting factors that are binding to GDF5, three of which are modulating GDF5 expression via the OA susceptibility locus rs143383.
GDF5 is an important growth factor that plays a vital role in the development and repair of articulating joints. rs143383 is a polymorphism within the regulatory region of the GDF5 gene and has two allelic forms, C and T. Genetic studies have demonstrated that the T allele is associated with an increased risk of osteoarthritis in a range of ethnic populations whilst previous functional studies revealed that this allele mediates its effect by producing less GDF5 transcript than the C allele. In this study, we sought to identify transcription factors that are binding to rs143383 and that are responsible for mediating this differential level of expression. Using two different approaches we have identified four factors and our functional studies have revealed that three of these factors repress GDF5 expression and that DEAF-1 modulates the differential expression of the two rs143383 alleles. The factors that we have identified could serve as novel therapeutic targets, with their depletion restoring the expression levels of GDF5 in patients with the osteoarthritis susceptibility T allele. The relevance of our results extends beyond osteoarthritis, since the T allele of rs143383 is also a risk factor for a number of other musculoskeletal diseases.
Osteoarthritis (OA) is the most prevalent form of arthritis and accounts for substantial morbidity and disability, particularly in the elderly. It is characterized by changes in joint structure including degeneration of the articular cartilage and its etiology is multifactorial with a strong postulated genetic component. We performed a meta-analysis of four genome-wide association (GWA) studies of 2,371 knee OA cases and 35,909 controls in Caucasian populations. Replication of the top hits was attempted with data from additional ten replication datasets. With a cumulative sample size of 6,709 cases and 44,439 controls, we identified one genome-wide significant locus on chromosome 7q22 for knee OA (rs4730250, p-value=9.2×10−9), thereby confirming its role as a susceptibility locus for OA. The associated signal is located within a large (500kb) linkage disequilibrium (LD) block that contains six genes; PRKAR2B (protein kinase, cAMP-dependent, regulatory, type II, beta), HPB1 (HMG-box transcription factor 1), COG5 (component of oligomeric golgi complex 5), GPR22 (G protein-coupled receptor 22), DUS4L (dihydrouridine synthase 4-like), and BCAP29 (the B-cell receptor-associated protein 29). Gene expression analyses of the (six) genes in primary cells derived from different joint tissues confirmed expression of all the genes in the joint environment.
The single nucleotide polymorphism (SNP) rs2615977 is associated with osteoarthritis (OA) and is located in intron 31 of COL11A1, a strong candidate gene for this degenerative musculoskeletal disease. Furthermore, the common non-synonymous COL11A1 SNP rs1676486 is associated with another degenerative musculoskeletal disease, lumbar disc herniation (LDH). rs1676486 is a C-T transition mediating its affect on LDH susceptibility by modulating COL11A1 expression. The risk T-allele of rs1676486 leads to reduced expression of the COL11A1 transcript, a phenomenon known as allelic expression imbalance (AEI). We were keen therefore to assess whether the effect that rs1676486 has on COL11A1 expression in LDH is also observed in OA and whether the rs2615977 association to OA also marked AEI.
Using RNA from OA cartilage, we assessed whether either SNP correlated with COL11A1 AEI by 1) measuring COL11A1 expression and stratifying the data by genotype at each SNP; and 2) quantifying the mRNA transcribed from each allele of the two SNPs. We also assessed whether rs1676486 was associated with OA susceptibility using a case–control cohort of over 18,000 individuals.
We observed significant AEI at rs1676486 (p < 0.0001) with the T-allele correlating with reduced COL11A1 expression. This corresponded with observations in LDH but the SNP was not associated with OA. We did not observe AEI at rs2615977.
COL11A1 is subject to AEI in OA cartilage. AEI at rs1676486 is a risk factor for LDH, but not for OA. These two diseases therefore share a common functional phenotype, namely AEI of COL11A1, but this appears to be a disease risk only in LDH. Other functional effects on COL11A1 presumably account for the OA susceptibility that maps to this gene.
Osteoarthritis; Lumbar disc herniation; Genetics; Susceptibility; COL11A1; Allelic expression
Recent genome-wide association scans (GWASs) along with several adequately powered
candidate gene studies have yielded a number of risk alleles for osteoarthritis (OA). This
number is now sufficiently large to allow conclusions to be drawn regarding the nature of
genetic susceptibility, including the fact that the risk alleles have variable effects
depending on sex, ethnicity and on the skeletal site of the disease. Several of the
alleles that have emerged from the GWASs are within or close to highly plausible candidate
genes, including RUNX2 and CHST11. However, the majority
of risk alleles do not map to genes previously reported to play a role in musculoskeletal
biology, indicating that the GWAS datasets are telling us something new about the OA
disease process. Functional studies have so far revealed that effects on gene expression
are likely to be one of the main mechanisms through which OA susceptibility is acting.
Epigenetic mechanisms such as DNA methylation also influence OA risk, and integration of
genetic, transcriptomic and epigenetic data will allow us to use the genetic discoveries
for informed development of new OA biological treatments.
We aimed to explore the involvement of a multiallelic functional polymorphism in knee osteoarthritis (OA) susceptibility as a prototype of possible genetic factors escaping GWAS detection.
OA patients and controls from three European populations (Greece, Spain and the UK) adding up to 1003 patients (716 women, 287 men) that had undergone total knee joint replacement (TKR) due to severe primary OA and 1543 controls (758 women, 785 men) lacking clinical signs or symptoms of OA were genotyped for the D6S1276 microsatellite in intron 1 of BMP5. Genotype and mutiallelic trend tests were used to compare cases and controls.
Significant association was found between the microsatellite and knee OA in women (P from 3.1 x10-4 to 4.1 x10-4 depending on the test), but not in men. Three of the alleles showed significant differences between patients and controls, one of them of increased risk and two of protection. The gender association and the allele direction of change were very concordant with those previously reported for hip OA.
We have found association of knee OA in women with the D6S1276 functional microsatellite that modifies in cis the expression of BMP5 making this a sounder OA genetic factor and extending its involvement to other joints. This result also shows the interest of analysing other multiallelic polymorphisms.
To identify genes involved in osteoarthritis (OA), the most prevalent form of joint disease, we performed a genome-wide association study (GWAS) in which we tested 500,510 Single Nucelotide Polymorphisms (SNPs) in 1341 OA cases and 3496 Dutch Caucasian controls. SNPs associated with at least two OA-phenotypes were analysed in 14,938 OA cases and approximately 39,000 controls. The C-allele of rs3815148 on chromosome 7q22 (MAF 23%, 172 kb upstream of the GPR22 gene) was consistently associated with a 1.14-fold increased risk (95%CI: 1.09–1.19) for knee- and/or hand-OA (p=8×10−8), and also with a 30% increased risk for knee-OA progression (95%CI: 1.03–1.64, p=0.03). This SNP is in almost complete linkage disequilibrium with rs3757713 (located 68 kb upstream of GPR22) which is associated with GPR22 expression levels in lymphoblast cell lines (p=4×10−12). GPR22 encodes an G-protein coupled receptor with unkown ligand (orphan receptor). Immunohistochemistry experiments showed absence of GPR22 in normal mouse articular cartilage or synovium. However, GPR22 positive chondrocytes were found in the upper layers of the articular cartilage of mouse knee joints that were challenged by in vivo papain treatment or in the presence of interleukin-1 driven inflammation. GRP22 positive chondrocyte-like cells were also found in osteophytes in instability-induced OA. In addition, GPR22 is also present in areas of the brain involved in locomotor function. Our findings reveal a novel common variant on chromosome 7q22 to influence susceptibility for prevalence and progression of OA.
A genome-wide association scan with subsequent replication study that involved over 67,000 individuals of European ancestry has produced evidence of association of single nucleotide polymorphism rs2277831 to primary osteoarthritis (OA) with a P-value of 2.9 × 10-5. rs2277831, an A/G transition, is located in an intron of MICAL3. This gene is located on chromosome 22q11.21 and the association signal encompasses two additional genes, BCL2L13 and BID. It is becoming increasingly apparent that many common complex traits are mediated by cis-acting regulatory polymorphisms that influence, in a tissue-specific manner, gene expression or transcript stability.
We used total and allelic expression analysis to assess whether the OA association to rs2277831 is mediated by an influence on MICAL3, BCL2L13 or BID expression. Using RNA extracted from joint tissues of 60 patients who had undergone elective joint replacement surgery, we assessed whether rs2277831 correlated with allelic expression of either of the three genes by: 1) measuring the expression of each gene by quantitative PCR and then stratifying the data by genotype at rs2277831 and 2) accurately discriminating and quantifying the mRNA synthesised from the alleles of OA patients using allelic-quantitative PCR.
We found no evidence for a correlation between gene expression and genotype at rs2277831, with P-values of 0.09 for BCL2L13, 0.07 for BID and 0.33 for MICAL3. In the allelic expression analysis we observed several examples of significant (p < 0.05) allelic imbalances, with an allelic expression ratio of 2.82 observed in BCL2L13 (P = 0.004), 2.09 at BID (P = 0.001) and the most extreme case being at MICAL3, with an allelic expression ratio of 5.47 (P = 0.001). However, there was no correlation observed between the pattern of allelic expression and the genotype at rs2277831.
In the tissues that we have studied, our data do not support our hypothesis that the association between rs2277831 and OA is due to the effect this SNP has on MICAL3, BCL2L13 or BID gene expression. Instead, our data point towards other functional effects accounting for the OA associated signal.
Genetics highlights relationships between biological systems, and as the number of defined osteoarthritis susceptibility alleles increases, there is the natural tendency to assess whether the alleles influence other musculoskeletal phenotypes. That has proven to be the case for the GDF5 polymorphism rs143383, a risk factor for knee osteoarthritis and several other common conditions, including lumbar-disc degeneration and developmental dysplasia of the hip. Another interesting example has recently emerged in the repeat polymorphism of the asporin gene, ASPN, which is also associated with these three phenotypes. Such discoveries increase our understanding of shared disease etiology but also emphasize the complexity of common genetic risk.
Imputation is an extremely valuable tool in conducting and synthesising genome-wide association studies (GWASs). Directly typed SNP quality control (QC) is thought to affect imputation quality. It is, therefore, common practise to use quality-controlled (QCed) data as an input for imputing genotypes. This study aims to determine the effect of commonly applied QC steps on imputation outcomes. We performed several iterations of imputing SNPs across chromosome 22 in a dataset consisting of 3177 samples with Illumina 610k (Illumina, San Diego, CA, USA) GWAS data, applying different QC steps each time. The imputed genotypes were compared with the directly typed genotypes. In addition, we investigated the correlation between alternatively QCed data. We also applied a series of post-imputation QC steps balancing elimination of poorly imputed SNPs and information loss. We found that the difference between the unQCed data and the fully QCed data on imputation outcome was minimal. Our study shows that imputation of common variants is generally very accurate and robust to GWAS QC, which is not a major factor affecting imputation outcome. A minority of common-frequency SNPs with particular properties cannot be accurately imputed regardless of QC stringency. These findings may not generalise to the imputation of low frequency and rare variants.
genome-wide association study; imputation; quality control; single nucleotide polymorphism
The objective of this study was to examine the relationship between common genetic variation of the ESR2 gene and osteoarthritis.
In the discovery study, the Rotterdam Study-I, 7 single nucleotide polymorphisms (SNPs) were genotyped and tested for association with hip (284 cases, 2772 controls), knee (665 cases, 2075 controls), and hand OA (874 cases, 2184 controls) using an additive model. In the replication stage one SNP (rs1256031) was tested in an additional 2080 hip, 1318 knee and 557 hand OA cases and 4001, 2631 and 1699 controls respectively. Fixed- and random-effects meta-analyses were performed over the complete dataset including 2364 hip, 1983 knee and 1431 hand OA cases and approximately 6000 controls.
The C allele of rs1256031 was associated with a 36% increased odds of hip OA in women of the Rotterdam Study-I (OR 1.36, 95% CI 1.08-1.70, p = 0.009). Haplotype analysis and analysis of knee- and hand OA did not give additional information. With the replication studies, the meta-analysis did not show a significant effect of this SNP on hip OA in the total population (OR 1.06, 95% CI 0.99-1.15, p = 0.10). Stratification according to gender did not change the results. In this study, we had 80% power to detect an odds ratio of at least 1.14 for hip OA (α = 0.05).
This study showed that common genetic variation in the ESR2 gene is not likely to influence the risk of osteoarthritis with effects smaller than a 13% increase.
In a previous study carried out by our group, the genotyping of 36 microsatellite markers from within a narrow interval of chromosome 6p12.3-q13 generated evidence for linkage and for association to female hip osteoarthritis (OA), with the most compelling association found for a marker within intron 1 of the bone morphogenetic protein 5 gene (BMP5). In this study, we aimed to further categorize the association of variants within intron 1 of BMP5 with OA through an expanded genetic association study of the intron and subsequent functional analysis of associated polymorphisms.
We genotyped 18 common polymorphisms including 8 microsatellites and 9 single nucleotide polymorphisms (SNPs) and 1 insertion/deletion (INDEL) from within highly conserved regions between human and mouse within intron 1 of BMP5. These markers were then tested for association to OA by a two-stage approach in which the polymorphisms were initially genotyped in a case-control cohort comprising 361 individuals with associated polymorphisms (P ≤ 0.05) then genotyped in a second case-control cohort comprising 1185 individuals.
Two BMP5 intron 1 polymorphisms demonstrated association in the combined case-control cohort of 1546 individuals (765 cases and 781 controls): microsatellite D6S1276 (P = 0.018) and SNP rs921126 (P = 0.013). Functional analyses in osteoblastic, chondrocytic, and adipocytic cell lines indicated that allelic variants of D6S1276 have significant effects on the transcriptional activity of the BMP5 promoter in vitro.
Variability in gene expression of BMP5 may be an important contributor to OA genetic susceptibility.
A genetic association of the ENPP1 gene with primary hand osteoarthritis was recently reported in this journal. ENPP1 encodes an enzyme that regulates soft tissue calcification. The study as it stands is far from complete because the actual causal variant(s) within ENPP1 has not been identified and no functional study on the activity of the enzyme in hand osteoarthritis was presented. Nevertheless, the study stimulates interest and will encourage others in the field to test ENPP1 as a possible osteoarthritis susceptibility gene in their cohorts. The genetic basis of osteoarthritis is slowly being uncovered, and this report constitutes another interesting find.
Osteoarthritis (OA) has a complex aetiology with a strong genetic component. Genome-wide association studies implicate several nuclear genes in the aetiology, but a major component of the heritability has yet to be defined at the molecular level. Initial studies implicate maternally inherited variants of mitochondrial DNA (mtDNA) in subgroups of patients with OA based on gender and specific joint involvement, but these findings have not been replicated.
The authors studied 138 maternally inherited mtDNA variants genotyped in a two cohort genetic association study across a total of 7393 OA cases from the arcOGEN consortium and 5122 controls genotyped in the Wellcome Trust Case Control consortium 2 study.
Following data quality control we examined 48 mtDNA variants that were common in cohort 1 and cohort 2, and found no association with OA. None of the phenotypic subgroups previously associated with mtDNA haplogroups were associated in this study.
We were not able to replicate previously published findings in the largest mtDNA association study to date. The evidence linking OA to mtDNA is not compelling at present.
Gene Polymorphism; Osteoarthritis; Pharmacogenetics
Obesity as measured by body mass index (BMI) is one of the major risk factors for osteoarthritis. In addition, genetic overlap has been reported between osteoarthritis and normal adult height variation. We investigated whether this relationship is due to a shared genetic aetiology on a genome-wide scale.
We compared genetic association summary statistics (effect size, p value) for BMI and height from the GIANT consortium genome-wide association study (GWAS) with genetic association summary statistics from the arcOGEN consortium osteoarthritis GWAS. Significance was evaluated by permutation. Replication of osteoarthritis association of the highlighted signals was investigated in an independent dataset. Phenotypic information of height and BMI was accounted for in a separate analysis using osteoarthritis-free controls.
We found significant overlap between osteoarthritis and height (p=3.3×10−5 for signals with p≤0.05) when the GIANT and arcOGEN GWAS were compared. For signals with p≤0.001 we found 17 shared signals between osteoarthritis and height and four between osteoarthritis and BMI. However, only one of the height or BMI signals that had shown evidence of association with osteoarthritis in the arcOGEN GWAS was also associated with osteoarthritis in the independent dataset: rs12149832, within the FTO gene (combined p=2.3×10−5). As expected, this signal was attenuated when we adjusted for BMI.
We found a significant excess of shared signals between both osteoarthritis and height and osteoarthritis and BMI, suggestive of a common genetic aetiology. However, only one signal showed association with osteoarthritis when followed up in a new dataset.
Osteoarthritis; Gene Polymorphism; Epidemiology
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.
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.
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).
Novel genetic loci for hip OA were found in this meta-analysis of GWAS.
Epidemiology; Gene Polymorphism; Osteoarthritis
The aim of this study was to characterize the genome-wide DNA methylation profile of chondrocytes from knee and hip cartilage obtained from patients with osteoarthritis (OA) and hip cartilage obtained from patients with femoral neck fracture, providing the first comparison of DNA methylation between OA and non-OA hip cartilage, and between OA hip and OA knee cartilage.
The study was performed using the Illumina Infinium HumanMethylation450 BeadChip array, which allows the annotation of ∼480,000 CpG sites. Genome-wide methylation was assessed in chondrocyte DNA extracted from 23 hip OA patients, 73 knee OA patients, and 21 healthy hip control patients with femoral neck fracture.
Analysis revealed that chondrocytes from the hip cartilage of OA patients and healthy controls have unique methylation profiles, with 5,322 differentially methylated loci (DMLs) identified between the 2 groups. In addition, a comparison between hip and knee OA chondrocytes revealed 5,547 DMLs between the 2 groups, including DMLs in several genes known to be involved in the pathogenesis of OA. Hip OA samples were found to cluster into 2 groups. A total of 15,239 DMLs were identified between the 2 clusters, with an enrichment of genes involved in inflammation and immunity. Similarly, we confirmed a previous report of knee OA samples that also clustered into 2 groups.
We demonstrated that global DNA methylation using a high-density array can be a powerful tool in the characterization of OA at the molecular level. Identification of pathways enriched in DMLs between OA and OA-free cartilage highlight potential etiologic mechanisms that are involved in the initiation and/or progression of the disease and that could be therapeutically targeted.