Martha C. Castaño-Betancourt

High Bone Mineral Density and Fracture Risk in Type 2 Diabetes as Skeletal Complications of Inadequate Glucose Control: The Rotterdam Study

OBJECTIVE Individuals with type 2 diabetes have increased fracture risk despite higher bone mineral density (BMD). Our aim was to examine the influence of glucose control on skeletal complications. RESEARCH DESIGN AND METHODS Data of 4,135 participants of the Rotterdam Study, a prospective population-based cohort, were available (mean follow-up 12.2 years). At baseline, 420 participants with type 2 diabetes were classified by glucose control (according to HbA1c calculated from fructosamine), resulting in three comparison groups: adequately controlled diabetes (ACD; n = 203; HbA1c <7.5%), inadequately controlled diabetes (ICD; n = 217; HbA1c ≥7.5%), and no diabetes (n = 3,715). Models adjusted for sex, age, height, and weight (and femoral neck BMD) were used to test for differences in bone parameters and fracture risk (hazard ratio [HR] [95% CI]). RESULTS The ICD group had 1.1–5.6% higher BMD, 4.6–5.6% thicker cortices, and −1.2 to −1.8% narrower femoral necks than ACD and ND, respectively. Participants with ICD had 47–62% higher fracture risk than individuals without diabetes (HR 1.47 [1.12–1.92]) and ACD (1.62 [1.09–2.40]), whereas those with ACD had a risk similar to those without diabetes (0.91 [0.67–1.23]). CONCLUSIONS Poor glycemic control in type 2 diabetes is associated with fracture risk, high BMD, and thicker femoral cortices in narrower bones. We postulate that fragility in apparently “strong” bones in ICD can result from microcrack accumulation and/or cortical porosity, reflecting impaired bone repair.

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.

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.

The DOT1L rs12982744 polymorphism is associated with osteoarthritis of the hip with genome-wide statistical significance in males

Osteoarthritis (OA) of the hip is a major cause of pain, disability and use of healthcare resources.1 Although OA is multifactorial, it is known to have significant genetic contribution and a number of studies have attempted to dissect such contribution (see2 for review). The DOT1L gene encodes the DOT1-like histone H3 methyltransferaseis, a potentially dedicated enzyme for Wnt target gene activation in leukaemia recently shown to be associated with endochondral bone formation.3 A polymorphism (rs12982744) in DOT1L has been found to be strongly associated with minimum joint space width (minJSW) at the hip. This exact same single nucleotide polymorphism (SNP) was previously identified to be associated with increased height.4 ,5 The C allele associated with lower minJSW and lower height was associated with hip OA, although this association did not reach genome-wide significance (GWS) (OR 1.14, CI 1.06 to 1.22; p=1.5  ×  10−4).3 The GWS level of p<5×10−8 is the threshold at which genetic associations are considered credible.6 The aim of our study was to prove that common genetic variation in the DOT1L gene is important …

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.

Assessment of Gene-by-Sex Interaction Effect on Bone Mineral Density

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene‐by‐sex autosomal interactions genome‐wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome‐wide meta‐analysis of gene‐by‐sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single‐nucleotide polymorphisms (SNPs; p < 1 × 10−5) in an additional set of 24,763 individuals. Gene‐by‐sex interaction and sex‐specific effects were examined in these 12 SNPs. We detected one novel genome‐wide significant interaction associated with LS‐BMD at the Chr3p26.1‐p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10−5; female effect = −0.007 and p = 3.3 × 10−2), and 11 suggestive loci associated with either FN‐ or LS‐BMD in discovery cohorts. However, there was no evidence for genome‐wide significant (p < 5 × 10−8) gene‐by‐sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome‐wide significant evidence for gene‐by‐sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene‐by‐sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2,666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at P<5x10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6 x 10-8. However, the association was not significant across 5,720 cases and 21,791 controls from 14 studies. Fixed-effects meta analyses summary estimate was 1.06 (95% CI: 0.98-1.14; P=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p= 0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (P=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size > 1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions are needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis.

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10−8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

Genome-wide association and functional studies identify a role for IGFBP3 in hip osteoarthritis

Objectives To identify genetic associations with hip osteoarthritis (HOA), we performed a meta-analysis of genome-wide association studies (GWAS) of HOA. Methods The GWAS meta-analysis included approximately 2.5 million imputed HapMap single nucleotide polymorphisms (SNPs). HOA cases and controls defined radiographically and by total hip replacement were selected from the Osteoporotic Fractures in Men (MrOS) Study and the Study of Osteoporotic Fractures (SOF) (654 cases and 4697 controls, combined). Replication of genome-wide significant SNP associations (p ≤5×10−8) was examined in five studies (3243 cases and 6891 controls, combined). Functional studies were performed using in vitro models of chondrogenesis and osteogenesis. Results The A allele of rs788748, located 65 kb upstream of the IGFBP3 gene, was associated with lower HOA odds at the genome-wide significance level in the discovery stage (OR 0.71, p=2×10−8). The association replicated in five studies (OR 0.92, p=0.020), but the joint analysis of discovery and replication results was not genome-wide significant (p=1×10−6). In separate study populations, the rs788748 A allele was also associated with lower circulating IGFBP3 protein levels (p=4×10−13), suggesting that this SNP or a variant in linkage disequilibrium could be an IGFBP3 regulatory variant. Results from functional studies were consistent with association results. Chondrocyte hypertrophy, a deleterious event in OA pathogenesis, was largely prevented upon IGFBP3 knockdown in chondrocytes. Furthermore, IGFBP3 overexpression induced cartilage catabolism and osteogenic differentiation. Conclusions Results from GWAS and functional studies provided suggestive links between IGFBP3 and HOA.

Association of lumbar disc degeneration with osteoporotic fractures; the Rotterdam study and meta-analysis from systematic review

OBJECTIVE To investigate the relation between lumbar disc degeneration (LDD) and all type of osteoporotic (OP) fractures including vertebral. METHODS This study is part of the Rotterdam study, a large prospective population-based cohort study among men and women aged 55years and over. In 2819 participants spine radiographs were scored for LDD (osteophytes and disc space narrowing (DSN)) from L1 till S1, using the Lane atlas. Osteoporotic (OP) fracture data were collected and verified by specialists during 12.8years. We considered two types of vertebral fractures (VFx): Clinical VFx (symptomatic fractures recorded by medical practitioners) and Radiographic VFx (using the McCloskey-Kanis method). Meta-analysis of published studies reporting an association of LDD features and VFx was performed. Differences in Bone Mineral Density (BMD) between participants with and without LDD features were analyzed using ANOVA. Risk of OP-fractures was analyzed using Cox regression. RESULTS In a total of 2385 participants, during 12.8years follow-up, 558 suffered an OP-fracture. Subjects with LDD had an increased OP fracture risk compared to subjects without LDD (HR: 1.29, CI: 1.04-1.60). LDD-cases have between 0.3 and 0.72 standard deviations more BMD than non-cases in all analyzed regions including total body BMD and skull BMD (P<0.001). Only males with LDD had increased risk for OP-fractures compared to males without LDD (adjusted-HR: 1.80, 95%CI: 1.20-2.70, P=0.005). The risk was also higher for VFx in males (HR: 1.64, CI: 1.03-2.60, P: 0.04). The association LDD-OP-fractures in females was lower and not significant (adjusted-HR: 1.08, 95%CI: 0.82-1.41). Meta-analyses showed that the risk of VFx in subjects with LDD has been studied only in women and there is not enough evidence to confidently analyze the relationship between LDD-features (DSN or/and OPH) and VFx due to low power and heterogeneity in phenotype definition in the collected studies. CONCLUSIONS Male subjects with LDD have a higher osteoporotic fracture risk, in spite of systemically higher BMD.