Genome-wide association and functional studies identify the DOT1L gene to be involved in cartilage thickness and hip osteoarthritis.

Hip osteoarthritis (HOA) is one of the most disabling and common joint disorders with a large genetic component that is, however, still ill-defined. To date, genome-wide association studies (GWAS) in osteoarthritis (OA) and specifically in HOA have yielded only few loci, which is partly explained by heterogeneity in the OA definition. Therefore, we here focused on radiographically measured joint-space width (JSW), a proxy for cartilage thickness and an important underlying intermediate trait for HOA. In a GWAS of 6,523 individuals on hip-JSW, we identified the G allele of rs12982744 on chromosome 19p13.3 to be associated with a 5% larger JSW (P = 4.8 × 10(-10)). The association was replicated in 4,442 individuals from three United Kingdom cohorts with an overall meta-analysis P value of 1.1 × 10(-11). The SNP was also strongly associated with a 12% reduced risk for HOA (P = 1 × 10(-4)). The SNP is located in the DOT1L gene, which is an evolutionarily conserved histone methyltransferase, recently identified as a potentially dedicated enzyme for Wnt target-gene activation in leukemia. Immunohistochemical staining of the DOT1L protein in mouse limbs supports a role for DOT1L in chondrogenic differentiation and adult articular cartilage. DOT1L is also expressed in OA articular chondrocytes. Silencing of Dot1l inhibited chondrogenesis in vitro. Dot1l knockdown reduces proteoglycan and collagen content, and mineralization during chondrogenesis. In the ATDC5 chondrogenesis model system, DOT1L interacts with TCF and Wnt signaling. These data are a further step to better understand the role of Wnt-signaling during chondrogenesis and cartilage homeostasis. DOT1L may represent a therapeutic target for OA.

A variant in MCF2L is associated with osteoarthritis.

Osteoarthritis (OA) is a prevalent, heritable degenerative joint disease with a substantial public health impact. We used a 1000-Genomes-Project-based imputation in a genome-wide association scan for osteoarthritis (3177 OA cases and 4894 controls) to detect a previously unidentified risk locus. We discovered a small disease-associated set of variants on chromosome 13. Through large-scale replication, we establish a robust association with SNPs in MCF2L (rs11842874, combined odds ratio [95% confidence interval] 1.17 [1.11-1.23], p = 2.1 × 10(-8)) across a total of 19,041 OA cases and 24,504 controls of European descent. This risk locus represents the third established signal for OA overall. MCF2L regulates a nerve growth factor (NGF), and treatment with a humanized monoclonal antibody against NGF is associated with reduction in pain and improvement in function for knee OA patients.

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

OBJECTIVES: Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS: We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS: We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS: Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

Bone parameters across different types of hip osteoarthritis and their relationship to osteoporotic fracture risk.

OBJECTIVE: The atrophic type of hip osteoarthritis (OA) is characterized by cartilage degradation without the formation of osteophytes. Individuals with atrophic OA have been less well studied, and it is unknown whether this OA type differs from the osteophytic types with regard to bone tissue. The purpose of this study was to examine bone mineral density (BMD), hip structural properties, and fracture risk in individuals with the atrophic type of OA as compared to those with the osteophytic types (normotrophic/hypertrophic) as well as individuals without OA. METHODS: This study is part of the Rotterdam Study, a large prospective population-based cohort study. We examined 5,006 participants who had been assessed for OA, BMD, and geometric measures at baseline and for incident nonvertebral osteoporotic fractures (mean followup 9.6 years). We estimated the differences in bone characteristics between the OA groups and the controls (no joint space narrowing or osteophytes). Cox proportional hazards regression was used to calculate osteoporotic fracture risk. RESULTS: Participants with atrophic OA had systemically lower BMD as compared to those with normotrophic OA and as compared to the controls (6.5% and 9% for total body BMD; 4% and 5% for skull BMD, respectively). Participants with osteophytic OA had ∼4% and ∼5% higher total body and skull BMD, respectively, a wider femoral neck, and greater bone strength (12% and 5% higher section modulus, respectively) as compared to the controls or to those with atrophic OA. The risk of osteoporotic fractures was almost 50% higher in those with atrophic OA as compared to the controls (hazard risk 1.48, P = 0.008). This difference was not explained by differences in the BMD, number of falls, degree of disability, or use of corticosteroids. CONCLUSION: Individuals with atrophic hip OA have an increased risk of osteoporotic fractures that is not fully explained by systemically lower BMD as compared to controls.

Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study.

BACKGROUND: Osteoarthritis is the most common form of arthritis worldwide and is a major cause of pain and disability in elderly people. The health economic burden of osteoarthritis is increasing commensurate with obesity prevalence and longevity. Osteoarthritis has a strong genetic component but the success of previous genetic studies has been restricted due to insufficient sample sizes and phenotype heterogeneity. METHODS: We undertook a large genome-wide association study (GWAS) in 7410 unrelated and retrospectively and prospectively selected patients with severe osteoarthritis in the arcOGEN study, 80% of whom had undergone total joint replacement, and 11,009 unrelated controls from the UK. We replicated the most promising signals in an independent set of up to 7473 cases and 42,938 controls, from studies in Iceland, Estonia, the Netherlands, and the UK. All patients and controls were of European descent. FINDINGS: We identified five genome-wide significant loci (binomial test p≤5·0×10(-8)) for association with osteoarthritis and three loci just below this threshold. The strongest association was on chromosome 3 with rs6976 (odds ratio 1·12 [95% CI 1·08-1·16]; p=7·24×10(-11)), which is in perfect linkage disequilibrium with rs11177. This SNP encodes a missense polymorphism within the nucleostemin-encoding gene GNL3. Levels of nucleostemin were raised in chondrocytes from patients with osteoarthritis in functional studies. Other significant loci were on chromosome 9 close to ASTN2, chromosome 6 between FILIP1 and SENP6, chromosome 12 close to KLHDC5 and PTHLH, and in another region of chromosome 12 close to CHST11. One of the signals close to genome-wide significance was within the FTO gene, which is involved in regulation of bodyweight-a strong risk factor for osteoarthritis. All risk variants were common in frequency and exerted small effects. INTERPRETATION: Our findings provide insight into the genetics of arthritis and identify new pathways that might be amenable to future therapeutic intervention. FUNDING: arcOGEN was funded by a special purpose grant from Arthritis Research UK.

Association of atherosclerosis with presence and progression of osteoarthritis: the Rotterdam Study.

OBJECTIVE: This study examined whether vascular alterations are associated with the presence and progression of osteoarthritis of the knee, the hip and the different hand joints in a large prospective cohort study. METHODS: In this population-based study involving participants aged 55 years and older (Rotterdam Study I), men (n=2372) and women (n=3278) were analysed separately. x-Rays of the knee, hip and hand were scored using the Kellgren and Lawrence score for osteoarthritis at baseline, after 6.6 years and 10 years. Measures of atherosclerosis (carotid intima media thickness (IMT) and carotid plaque) and data on covariates (age, body mass index, hypertension, cholesterol ratio, diabetes mellitus and smoking) were collected at baseline. Multivariate logistic regression models with generalised estimated equations were used to calculate OR and corresponding 95% CI. Secondary multiple comparison adjustment resulted in a significance level of p<0.0021. RESULTS: In women, IMT showed an independent association with the prevalence of knee osteoarthritis (adjusted OR (aOR) 1.7, 1.1 to 2.7), and carotid plaque with distal interphalangeal (DIP) osteoarthritis (aOR 1.4, 1.2 to 1.7) and with metacarpophalangeal osteoarthritis (aOR 1.5, 1.1 to 2.2). An independent association for IMT with progression of metacarpophalangeal osteoarthritis was found in women (aOR 2.9, 1.18 to 6.93). Additional adjustment for multiple testing yielded a significant association between carotid plaque and DIP osteoarthritis in women (p<0.001). CONCLUSIONS: This study showed independent associations of atherosclerosis with osteoarthritis of the knee and hand joints in women. The evidence was most solid for a relation with DIP osteoarthritis. More research is needed to confirm the associations and examine the differential association with various joints.

Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region.

BACKGROUND AND OBJECTIVES: Chronic widespread pain (CWP) is a common disorder affecting ∼10% of the general population and has an estimated heritability of 48-52%. In the first large-scale genome-wide association study (GWAS) meta-analysis, we aimed to identify common genetic variants associated with CWP. METHODS: We conducted a GWAS meta-analysis in 1308 female CWP cases and 5791 controls of European descent, and replicated the effects of the genetic variants with suggestive evidence for association in 1480 CWP cases and 7989 controls. Subsequently, we studied gene expression levels of the nearest genes in two chronic inflammatory pain mouse models, and examined 92 genetic variants previously described associated with pain. RESULTS: The minor C-allele of rs13361160 on chromosome 5p15.2, located upstream of chaperonin-containing-TCP1-complex-5 gene (CCT5) and downstream of FAM173B, was found to be associated with a 30% higher risk of CWP (minor allele frequency=43%; OR=1.30, 95% CI 1.19 to 1.42, p=1.2×10(-8)). Combined with the replication, we observed a slightly attenuated OR of 1.17 (95% CI 1.10 to 1.24, p=4.7×10(-7)) with moderate heterogeneity (I2=28.4%). However, in a sensitivity analysis that only allowed studies with joint-specific pain, the combined association was genome-wide significant (OR=1.23, 95% CI 1.14 to 1.32, p=3.4×10(-8), I2=0%). Expression levels of Cct5 and Fam173b in mice with inflammatory pain were higher in the lumbar spinal cord, not in the lumbar dorsal root ganglions, compared to mice without pain. None of the 92 genetic variants previously described were significantly associated with pain (p>7.7×10(-4)). CONCLUSIONS: We identified a common genetic variant on chromosome 5p15.2 associated with joint-specific CWP in humans. This work suggests that CCT5 and FAM173B are promising targets in the regulation of pain.

Meta-analysis of genome-wide association studies confirms a susceptibility locus for knee osteoarthritis on chromosome 7q22.

OBJECTIVES: Osteoarthritis (OA) is the most prevalent form of arthritis and accounts for substantial morbidity and disability, particularly in older people. It is characterised by changes in joint structure, including degeneration of the articular cartilage, and its aetiology is multifactorial with a strong postulated genetic component. METHODS: A meta-analysis was performed of four genome-wide association (GWA) studies of 2371 cases of knee OA and 35 909 controls in Caucasian populations. Replication of the top hits was attempted with data from 10 additional replication datasets. RESULTS: With a cumulative sample size of 6709 cases and 44 439 controls, one genome-wide significant locus was identified on chromosome 7q22 for knee OA (rs4730250, p=9.2 × 10⁻9), thereby confirming its role as a susceptibility locus for OA. CONCLUSION: The associated signal is located within a large (500 kb) linkage disequilibrium block that contains six genes: PRKAR2B (protein kinase, cAMP-dependent, regulatory, type II, ß), 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 (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.

Serum C reactive protein levels and genetic variation in the CRP gene are not associated with the prevalence, incidence or progression of osteoarthritis independent of body mass index.

OBJECTIVE: To study the relationship between serum C reactive protein (CRP) levels, genetic variation in the CRP gene and the prevalence, incidence and progression of radiographic osteoarthritis (ROA) in the Rotterdam Study-I (RS-I). A systematic review of studies assessing the relationship between osteoarthritis (OA) and CRP levels was also performed. METHODS: The association between CRP levels and genetic variation in the CRP gene and ROA was examined in 861 patients with hand OA, 718 with knee OA, 349 with hip OA and 2806 controls in the RS-I using one-way analysis of covariance and logistic regression, respectively. PubMed was searched for articles published between January 1992 and August 2009 assessing the relationship between CRP levels and OA. RESULTS: In RS-I the prevalence of knee OA, but not hip OA or hand OA, was associated with 14% higher serum CRP levels compared with controls (p=0.001). This association disappeared after adjustment for age and especially body mass index (BMI) (p=0.33). Genetic variation of the CRP gene was not consistently associated with the prevalence, incidence or progression of OA within RS-I. The systematic review included 18 studies (including RS-I) on serum CRP levels and the prevalence, incidence or progression of OA. Consistently higher crude CRP levels were found in cases of prevalent knee OA compared with controls. No association was observed between serum CRP levels and the prevalence of knee OA following adjustment for BMI (n=3 studies, meta-analysis p value=0.61). CONCLUSION: There is no evidence of association between serum CRP levels or genetic variation in the CRP gene with the prevalence, incidence or progression of OA independent of BMI.

Common genetic variation in the Estrogen Receptor Beta (ESR2) gene and osteoarthritis: results of a meta-analysis.

BACKGROUND: The objective of this study was to examine the relationship between common genetic variation of the ESR2 gene and osteoarthritis. METHODS: 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. RESULTS: 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). CONCLUSION: 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.

Finger length pattern as a biomarker for osteoarthritis and chronic joint pain: a population-based study and meta-analysis after systematic review.

OBJECTIVE: Type 3 finger length pattern (longer fourth digit than second digit) is influenced by prenatal androgens and has been studied previously as a biomarker for sexually dimorphic traits. Because osteoarthritis (OA) and chronic pain are known to be sexually dimorphic traits, we evaluated the association between finger length pattern and OA and chronic joint pain. METHODS: This study was part of the Rotterdam Study, a prospective population-based cohort study. We examined 4,784 participants. Associations between type 3 finger length and radiologic knee, hip, and hand OA and chronic joint pain were analyzed using a logistic regression model. Our results for OA were combined with previously published data in a meta-analysis. RESULTS: Participants with type 3 finger length pattern had an odds ratio of 1.64 for hand OA (P = 1.06 × 10(-7)). No associations with radiologic knee or hip OA were observed in the Rotterdam Study. The meta-analysis of previously published data and our novel data showed a significant association between type 3 finger length pattern and clinical symptomatic knee OA, but no association was found with radiologic knee OA. In addition, within the Rotterdam Study, we observed an odds ratio of 1.41 for individuals having joint pain at multiple sites (P = 1.4 × 10(-3)). CONCLUSION: Type 3 finger length pattern, as an indicator of prenatal androgen exposure, was associated with having symptomatic knee OA, chronic pain, and hand OA. Therefore, it may be applicable as an easy measurable biomarker to identify susceptible subjects for these traits.

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

OBJECTIVE: To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA. METHODS: A total of 199 candidate genes for association with OA were identified using Human Genome Epidemiology (HuGE) Navigator. All of their single-nucleotide polymorphisms (SNPs) with an allele frequency of >5% were assessed by fixed-effects meta-analysis of 9 genome-wide association studies (GWAS) that included 5,636 patients with knee OA and 16,972 control subjects and 4,349 patients with hip OA and 17,836 control subjects of European ancestry. An additional 5,921 individuals were genotyped for significantly associated SNPs in the meta-analysis. After correction for the number of independent tests, P values less than 1.58 × 10(-5) were considered significant. RESULTS: SNPs at only 2 of the 199 candidate genes (COL11A1 and VEGF) were associated with OA in the meta-analysis. Two SNPs in COL11A1 showed association with hip OA in the combined analysis: rs4907986 (P = 1.29 × 10(-5) , odds ratio [OR] 1.12, 95% confidence interval [95% CI] 1.06-1.17) and rs1241164 (P = 1.47 × 10(-5) , OR 0.82, 95% CI 0.74-0.89). The sex-stratified analysis also showed association of COL11A1 SNP rs4908291 in women (P = 1.29 × 10(-5) , OR 0.87, 95% CI 0.82-0.92); this SNP showed linkage disequilibrium with rs4907986. A single SNP of VEGF, rs833058, showed association with hip OA in men (P = 1.35 × 10(-5) , OR 0.85, 95% CI 0.79-0.91). After additional samples were genotyped, association at one of the COL11A1 signals was reinforced, whereas association at VEGF was slightly weakened. CONCLUSION: Two candidate genes, COL11A1 and VEGF, were significantly associated with OA in this focused meta-analysis. The remaining candidate genes were not associated.

Factors for pain in patients with different grades of knee osteoarthritis.

OBJECTIVE: Discordance between having pain and radiologic osteoarthritis (OA) is a well-established fact. It is suggested that this particularly applies to the less severe grades of OA. However, some people with a Kellgren/Lawrence (K/L) grade of 3 or 4 for OA are without pain. This study aimed to identify factors and differences in the factors associated with pain in persons with different grades of knee OA. METHODS: We stratified the knees of more than 5,000 participants of a population-based cohort study, the Rotterdam Study, based on the grade of knee OA. Multivariate generalized estimating equation analysis was used to analyze the association with knee pain. We tested several factors not directly related to structural damage of the knee. RESULTS: As expected, an increasing percentage of participants did not report pain with decreasing severity of knee OA: 25.8% for K/L grade 3 or 4 and 84.5% for no knee OA. Being a woman, having widespread pain, reporting general health symptoms, familial OA, and morning stiffness are factors for knee pain, but not specific for a particular grade of radiographic knee OA. Depression and hip OA showed significant interactions with the grade of OA being a factor for knee pain in knees without OA (K/L grade 0), but not in knees with OA. In addition, increasing age is protective for reporting pain in general. CONCLUSION: Several factors are associated with knee pain, but are not specific for a grade of radiographic knee OA. Two factors were associated with knee pain in the knee without signs of OA.

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

OBJECTIVE: To identify novel genes involved in osteoarthritis (OA), by means of a genome-wide association study. METHODS: We tested 500,510 single-nucleotide polymorphisms (SNPs) in 1,341 Dutch Caucasian OA cases and 3,496 Dutch Caucasian controls. SNPs associated with at least 2 OA phenotypes were analyzed in 14,938 OA cases and approximately 39,000 controls. Meta-analyses were performed using the program Comprehensive Meta-analysis, with P values <1 x 10(-7) considered genome-wide significant. RESULTS: The C allele of rs3815148 on chromosome 7q22 (minor allele frequency 23%; intron 12 of the COG5 gene) was associated with a 1.14-fold increased risk (95% confidence interval 1.09-1.19) of knee and/or hand OA (P = 8 x 10(-8)) and also with a 30% increased risk of knee OA progression (95% confidence interval 1.03-1.64) (P = 0.03). This SNP is in almost complete linkage disequilibrium with rs3757713 (68 kb upstream of GPR22), which is associated with GPR22 expression levels in lymphoblast cell lines (P = 4 x 10(-12)). Immunohistochemistry experiments revealed that G protein-coupled receptor protein 22 (GPR22) was absent 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 with in vivo papain treatment or methylated bovine serum albumin treatment. GPR22-positive chondrocyte-like cells were also found in osteophytes in instability-induced OA. CONCLUSION: Our findings identify a novel common variant on chromosome 7q22 that influences susceptibility to prevalence and progression of OA. Since the GPR22 gene encodes a G protein-coupled receptor, this is potentially an interesting therapeutic target.

Large-scale analysis of association between GDF5 and FRZB variants and osteoarthritis of the hip, knee, and hand.

OBJECTIVE: GDF5 and FRZB have been proposed as genetic loci conferring susceptibility to osteoarthritis (OA); however, the results of several studies investigating the association of OA with the rs143383 polymorphism of the GDF5 gene or the rs7775 and rs288326 polymorphisms of the FRZB gene have been conflicting or inconclusive. To examine these associations, we performed a large-scale meta-analysis of individual-level data. METHODS: Fourteen teams contributed data on polymorphisms and knee, hip, and hand OA. For rs143383, the total number of cases and controls, respectively, was 5,789 and 7,850 for hip OA, 5,085 and 8,135 for knee OA, and 4,040 and 4,792 for hand OA. For rs7775, the respective sample sizes were 4,352 and 10,843 for hip OA, 3,545 and 6,085 for knee OA, and 4,010 and 5,151 for hand OA, and for rs288326, they were 4,346 and 8,034 for hip OA, 3,595 and 6,106 for knee OA, and 3,982 and 5,152 for hand OA. For each individual study, sex-specific odds ratios (ORs) were calculated for each OA phenotype that had been investigated. The ORs for each phenotype were synthesized using both fixed-effects and random-effects models for allele-based effects, and also for haplotype effects for FRZB. RESULTS: A significant random-effects summary OR for knee OA was demonstrated for rs143383 (1.15 [95% confidence interval 1.09-1.22]) (P=9.4x10(-7)), with no significant between-study heterogeneity. Estimates of effect sizes for hip and hand OA were similar, but a large between-study heterogeneity was observed, and statistical significance was borderline (for OA of the hip [P=0.016]) or absent (for OA of the hand [P=0.19]). Analyses for FRZB polymorphisms and haplotypes did not reveal any statistically significant signals, except for a borderline association of rs288326 with hip OA (P=0.019). CONCLUSION: Evidence of an association between the GDF5 rs143383 polymorphism and OA is substantially strong, but the genetic effects are consistent across different populations only for knee OA. Findings of this collaborative analysis do not support the notion that FRZB rs7775 or rs288326 has any sizable genetic effect on OA phenotypes.