J.B. van Meurs

Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study

Summary Background Osteoporosis is diagnosed by the measurement of bone mineral density, which is a highly heritable and multifactorial trait. We aimed to identify genetic loci that are associated with bone mineral density. Methods In this genome-wide association study, we identified the most promising of 314 075 single nucleotide polymorphisms (SNPs) in 2094 women in a UK study. We then tested these SNPs for replication in 6463 people from three other cohorts in western Europe. We also investigated allelic expression in lymphoblast cell lines. We tested the association between the replicated SNPs and osteoporotic fractures with data from two studies. Findings We identified genome-wide evidence for an association between bone mineral density and two SNPs (p<5×10−8). The SNPs were rs4355801, on chromosome 8, near to the TNFRSF11B (osteoprotegerin) gene, and rs3736228, on chromosome 11 in the LRP5 (lipoprotein-receptor-related protein) gene. A non-synonymous SNP in the LRP5 gene was associated with decreased bone mineral density (rs3736228, p=6·3×10−12 for lumbar spine and p=1·9×10−4 for femoral neck) and an increased risk of both osteoporotic fractures (odds ratio [OR] 1·3, 95% CI 1·09–1·52, p=0·002) and osteoporosis (OR 1·3, 1·08–1·63, p=0·008). Three SNPs near the TNFRSF11B gene were associated with decreased bone mineral density (top SNP, rs4355801: p=7·6×10−10 for lumbar spine and p=3·3×10−8 for femoral neck) and increased risk of osteoporosis (OR 1·2, 95% CI 1·01–1·42, p=0·038). For carriers of the risk allele at rs4355801, expression of TNFRSF11B in lymphoblast cell lines was halved (p=3·0×10−6). 1883 (22%) of 8557 people were at least heterozygous for these risk alleles, and these alleles had a cumulative association with bone mineral density (trend p=2·3×10−17). The presence of both risk alleles increased the risk of osteoporotic fractures (OR 1·3, 1·08–1·63, p=0·006) and this effect was independent of bone mineral density. Interpretation Two gene variants of key biological proteins increase the risk of osteoporosis and osteoporotic fracture. The combined effect of these risk alleles on fractures is similar to that of most well-replicated environmental risk factors, and they are present in more than one in five white people, suggesting a potential role in screening. Funding Wellcome Trust, European Commission, NWO Investments, Arthritis Research Campaign, Chronic Disease Research Foundation, Canadian Institutes of Health Research, European Society for Clinical and Economic Aspects of Osteoporosis, Genome Canada, Genome Quebéc, Canada Research Chairs, National Health and Medical Research Council of Australia, and European Union.

Insights into the genetic architecture of osteoarthritis from stage 1 of the arcOGEN study

Objectives The genetic aetiology of osteoarthritis has not yet been elucidated. To enable a well-powered genome-wide association study (GWAS) for osteoarthritis, the authors have formed the arcOGEN Consortium, a UK-wide collaborative effort aiming to scan genome-wide over 7500 osteoarthritis cases in a two-stage genome-wide association scan. Here the authors report the findings of the stage 1 interim analysis. Methods The authors have performed a genome-wide association scan for knee and hip osteoarthritis in 3177 cases and 4894 population-based controls from the UK. Replication of promising signals was carried out in silico in five further scans (44 449 individuals), and de novo in 14 534 independent samples, all of European descent. Results None of the association signals the authors identified reach genome-wide levels of statistical significance, therefore stressing the need for corroboration in sample sets of a larger size. Application of analytical approaches to examine the allelic architecture of disease to the stage 1 genome-wide association scan data suggests that osteoarthritis is a highly polygenic disease with multiple risk variants conferring small effects. Conclusions Identifying loci conferring susceptibility to osteoarthritis will require large-scale sample sizes and well-defined phenotypes to minimise heterogeneity.

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.

Gremlin 1, frizzled-related protein, and dkk-1 are key regulators of human articular cartilage homeostasis

OBJECTIVE The development of osteoarthritis (OA) may be caused by activation of hypertrophic differentiation of articular chondrocytes. Healthy articular cartilage is highly resistant to hypertrophic differentiation, in contrast to other hyaline cartilage subtypes, such as growth plate cartilage. The purpose of this study was to elucidate the molecular mechanism responsible for the difference in the propensity of human articular cartilage and growth plate cartilage to undergo hypertrophic differentiation. METHODS Whole-genome gene-expression microarray analysis of healthy human growth plate and articular cartilage derived from the same adolescent donors was performed. Candidate genes, which were enriched in the articular cartilage, were validated at the messenger RNA (mRNA) and protein levels and examined for their potential to inhibit hypertrophic differentiation in two models. In addition, we studied a possible genetic association with OA. RESULTS Pathway analysis demonstrated decreased Wnt signaling in articular cartilage as compared to growth plate cartilage. This was at least partly due to increased expression of the bone morphogenetic protein and Wnt antagonists Gremlin 1, Frizzled-related protein (FRP), and Dkk-1 at the mRNA and protein levels in articular cartilage. Supplementation of these proteins diminished terminal hypertrophic differentiation without affecting chondrogenesis in long-bone explant cultures and chondrogenically differentiating human mesenchymal stem cells. Additionally, we found that single-nucleotide polymorphism rs12593365, which is located in a genomic control region of GREM1, was significantly associated with a 20% reduced risk of radiographic hip OA in 2 population-based cohorts. CONCLUSION Taken together, our study identified Gremlin 1, FRP, and Dkk-1 as natural brakes on hypertrophic differentiation in articular cartilage. As hypertrophic differentiation of articular cartilage may contribute to the development of OA, our findings may open new avenues for therapeutic intervention.

Statin use is associated with reduced incidence and progression of knee osteoarthritis in the Rotterdam study

Background Osteoarthritis is the most frequent chronic joint disease causing pain and disability. Besides biomechanical mechanisms, the pathogenesis of osteoarthritis may involve inflammation, vascular alterations and dysregulation of lipid metabolism. As statins are able to modulate many of these processes, this study examines whether statin use is associated with a decreased incidence and/or progression of osteoarthritis. Methods Participants in a prospective population-based cohort study aged 55 years and older (n=2921) were included. x-Rays of the knee/hip were obtained at baseline and after on average 6.5 years, and scored using the Kellgren and Lawrence score for osteoarthritis. Any increase in score was defined as overall progression (incidence and progression). Data on covariables were collected at baseline. Information on statin use during follow-up was obtained from computerised pharmacy databases. The overall progression of osteoarthritis was compared between users and non-users of statins. Using a multivariate logistic regression model with generalised estimating equation, OR and 95% CI were calculated after adjusting for confounding variables. Results Overall progression of knee and hip osteoarthritis occurred in 6.9% and 4.7% of cases, respectively. The adjusted OR for overall progression of knee osteoarthritis in statin users was 0.43 (95% CI 0.25 to 0.77, p=0.01). The use of statins was not associated with overall progression of hip osteoarthritis. Conclusions Statin use is associated with more than a 50% reduction in overall progression of osteoarthritis of the knee, but not of the hip.

High prevalence of vertebral fractures despite normal bone mineral density in patients with long-term controlled acromegaly.

OBJECTIVE To establish the prevalence of osteoporosis, vertebral fractures (VFs), and non-VFs in acromegaly patients with long-term controlled disease and factors potentially influencing fracture risk. DESIGN Case-control study. Patients and measurements Eighty-nine patients (46% male, mean age: 58 years) were included. We studied VFs and non-VFs, bone mineral density (BMD), and markers of bone turnover. In 48 patients, BMD assessment was also obtained 7 years prior to the current study. To compare VF prevalence, data from a sample of the Dutch population (n=3469) were used. RESULTS VF prevalence was 59% (men 64% and women 54%), significantly increased when compared with controls (odds ratio up to 6.5), and independent of the duration of disease control, BMD, markers of bone turnover, and acromegalic disease characteristics. Mean number of VFs per patient was 3.4±0.3 (range 1-8). There was no relationship between the number and severity of fractures, parameters of bone turnover, and follow-up BMD measurements. BMD did not change during prolongation of follow-up by 7 years of controlled acromegaly. CONCLUSION There is a very high prevalence of VFs in acromegaly patients with long-term controlled disease, independently of BMD. In view of the significant morbidity and mortality associated with VFs in general and the inability of BMD to predict fracture risk in acromegalic patients, we propose to include VF assessment, for example by lateral conventional radiographs of the spine in the screening of patients with acromegaly, both at diagnosis and during follow-up after establishment of disease control.

Prediction model for knee osteoarthritis incidence, including clinical, genetic and biochemical risk factors

Objective To develop and validate a prognostic model for incident knee osteoarthritis (KOA) in a general population and determine the value of different risk factor groups to prediction. Methods The prognostic model was developed in 2628 individuals from the Rotterdam Study-I (RS-I). Univariate and multivariate analyses were performed for questionnaire/easily obtainable variables, imaging variables, genetic and biochemical markers. The extended multivariate model was tested on discrimination (receiver operating characteristic curve and area under the curve (AUC)) in two other population-based cohorts: Rotterdam Study-II and Chingford Study. Results In RS-I, there was moderate predictive value for incident KOA based on the genetic score alone in subjects aged <65 years (AUC 0.65), while it was only 0.55 for subjects aged ≥65 years. The AUC for gender, age and body mass index (BMI) in prediction for KOA was 0.66. Addition of the questionnaire variables, genetic score or biochemical marker urinary C-terminal cross-linked telopeptide of type II collagen to the model did not change the AUC. However, when adding the knee baseline KL score to the model the AUC increased to 0.79. Applying external validation, similar results were observed in the Rotterdam Study-II and the Chingford Study. Conclusions Easy obtainable ‘Questionnaire’ variables, genetic markers, OA at other joint sites and biochemical markers add only modestly to the prediction of KOA incidence using age, gender and BMI in an elderly population. Doubtful minor radiographic degenerative features in the knee, however, are a very strong predictor of future KOA. This is an important finding, as many radiologists do not report minor degenerative changes in the knee.

GDF5 single-nucleotide polymorphism rs143383 is associated with lumbar disc degeneration in Northern European women

Objective Lumbar disc degeneration (LDD) is a serious social and medical problem which has been shown to be highly heritable. It has similarities with peripheral joint osteoarthritis (OA) in terms of both epidemiology and pathologic processes. A few known genetic variants have been identified using a candidate gene approach, but many more are thought to exist. GDF5 is a gene whose variants have been shown to play a role in skeletal height as well as predisposing to peripheral joint OA. In vitro, the gene product growth differentiation factor 5 has been shown to promote growth and repair of animal disc. This study was undertaken to investigate whether the GDF5 gene plays a role in LDD. Methods We investigated whether the 5′ upstream single-nucleotide polymorphism (SNP) variant rs143383 was associated with LDD, using plain radiography and magnetic resonance imaging to identify disc space narrowing and osteophytes, in 5 population cohorts from Northern Europe. Results An association between LDD and the SNP rs143383 was identified in women, with the same risk allele as in knee and hip OA (odds ratio 1.72 [95% confidence interval 1.15–2.57], P = 0.008). Conclusion Our findings in 5 population cohorts from Northern Europe indicate that a variant in the GDF5 gene is a risk factor for LDD in women. Many more such variants are predicted to exist, but this result highlights the growth and differentiation cellular pathway as a possible route to a better understanding of the process behind lumbar disc degeneration.

Genetic variation in the GDF5 region is associated with osteoarthritis, height, hip axis length and fracture risk: the Rotterdam study.

Background: A polymorphism (rs143383; T to C) near the GDF5 gene has been associated with height and osteoarthritis (OA), but debate exists about whether its primary biological action is directed to cartilage or bone. Objective: To study the association between genetic variation in the GDF5 region and radiographic osteoarthritis (ROA) susceptibility, height, bone size parameters and fracture risk in a large population-based cohort of Caucasian elderly subjects. Methods: 6365 men and women had genotype data available. ROA was defined as a Kellgren/Lawrence (K/L) score ⩾2 for hand, knee and hip joints. CTX-II levels, height, bone mineral density (BMD), bone size and fracture risk were also assessed. Results: rs143383 and three highly correlated single nucleotide polymorphisms (SNPs) in the GDF5 region were found to be independently associated with OA, height, bone size and fracture risk in women. Women with homozygotes for the rs143383 C allele had a 37% lower risk for hand OA (p = 8×10−6) and a 28% lower risk for knee OA (p = 0.003). In addition, they were 1.1 cm taller (p = 0.001), had a larger hip axis length (HAL) (p = 4×10−4) and had a 29% increased risk of incident non-vertebral fractures (p = 0.02). No associations with hip OA or BMD were detected. No associations were found in men. Conclusion: This population-based study shows that GDF5 gene variants are associated with hand OA, knee OA and fracture risk in elderly women. It also replicates previous association between GDF5 variation and height. Furthermore, our findings for HAL suggest that GDF5 action is primarily directed to the long bones, rather than the axial skeleton.

Radiographic osteoarthritis at three joint sites and FRZB, LRP5, and LRP6 polymorphisms in two population-based cohorts

OBJECTIVE To examine the association of genetic variation in key players in the Wnt signaling pathway with aspects of osteoarthritis (OA) in two population-based cohort studies: the Rotterdam Study and the Chingford Study. METHODS Radiographic OA (ROA) was defined as a Kellgren/Lawrence score (K/L) score > or = 2 for the knee and hip. Total hip replacement (THR) was scored. Hand OA was defined as presence of ROA (K/L > or = 2) in two out of three hand joint groups [distal interphalangeal (DIPs), proximal interphalangeal (PIPs), first carpometacarpal (CMC1)/trapezio-scaphoid joint (TS)] of each hand. The concentration of urinary C-terminal cross-linked telopeptide of type II collagen (CTX-II) was standardized to the total urine creatinine. Genotypes for the amino acid variants, Arg200Trp and Arg324Gly of Frizzled-Related protein gene (FRZB), Ala1330Val of Low-density lipoprotein receptor-related protein 5 (LRP5) and Ile1062Val of Low-density lipoprotein receptor-related protein 6 (LRP6), were obtained using the Taqman allelic discrimination assay. A meta-analysis was performed for the FRZB Arg324Gly polymorphism and hip- and knee-OA using RevMan version 4.3. RESULTS No consistent associations were observed between the FRZB, LRP5 and LRP6 amino acid variants and radiographic hip-, knee-, or hand-OA or THR, in either study population. While power was limited for most studies to date, a meta-analysis of all published studies regarding the FRZB Arg324Gly polymorphism was performed for hip- and knee-OA separately. This showed no significant associations between the Gly324 allele and risk for hip- or knee OA, although there was large heterogeneity between studies for hip OA in females. CONCLUSION No association was seen between FRZB, LRP5 and LRP6 variants with radiographic osteoarthritic outcomes in two population-based cohorts. In future studies, increased power and standardization of OA-phenotypes are highly recommended for replication studies and to allow meta-analysis.