Objective

To investigate hip shape by active shape modeling (ASM) as a potential predictor of incident radiographic hip osteoarthritis (RHOA) and symptomatic hip osteoarthritis (SRHOA).

Methods

All hips developing RHOA from baseline (Kellgren/Lawrence [K/L] grade 0/1) to mean 6-year followup (K/L grade ≥2, 190 hips) and 1:1 control hips (K/L grade 0/1 at both times, 192 hips) were included. Proximal femur shape was defined on baseline anteroposterior pelvis radiographs and submitted to ASM, producing a mean shape and continuous variables representing independent modes of shape variation. Mode scores (n = 14, explaining 95% of shape variance) were simultaneously included in logistic regression models with incident RHOA and SRHOA as dependent variables, adjusted for intraperson correlations, sex, race, body mass index (BMI), baseline K/L grade, and/or symptoms.

Results

We evaluated 382 hips from 342 individuals: 61% women and 83% white, with mean age 62 years and mean BMI 29 kg/m(2) . Several modes differed by sex and race, but no modes were associated with incident RHOA overall. Among men only, modes 1 and 2 were significantly associated (for a 1-SD decrease in mode 1 score: odds ratio [OR] 1.7 [95% confidence interval (95% CI) 1.1-2.5] and for a 1-SD increase in mode 2 score: OR 1.5 [95% CI 1.0-2.2]) with incident RHOA. A 1-SD decrease in mode 2 or 3 score increased the odds of SRHOA by 50%.

Conclusion

This study confirms other reports that variations in proximal femur shape have a modest association with incident hip OA. The observation of proximal femur shape associations with hip symptoms requires further investigation.

Objective

Hip shape by statistical shape modeling (SSM) is associated with hip radiographic osteoarthritis (rOA). We examined associations between hip shape and knee rOA given the biomechanical interrelationships between these joints.

Methods

Bilateral baseline hip shape assessments [for those with at least 1 hip with a Kellgren-Lawrence arthritis grading scale (KL) 0 or 1] from the Johnston County Osteoarthritis Project were available. Proximal femur shape was defined on baseline pelvis radiographs and evaluated by SSM, producing mean shape and continuous variables representing independent modes of variation (14 modes = 95% of shape variance). Outcomes included prevalent [baseline KL ≥ 2 or total knee replacement (TKR)], incident (baseline KL 0/1 with followup ≥ 2), and progressive knee rOA (KL increase of ≥ 1 or TKR). Limb-based logistic regression models for ipsilateral and contralateral comparisons were adjusted for age, sex, race, body mass index (BMI), and hip rOA, accounting for intraperson correlations.

Results

We evaluated 681 hips and 682 knees from 342 individuals (61% women, 83% white, mean age 62 yrs, BMI 29 kg/m(2)). Ninety-nine knees (15%) had prevalent rOA (4 knees with TKR). Lower modes 2 and 3 scores were associated with ipsilateral prevalent knee rOA, and only lower mode 3 scores were associated with contralateral prevalent knee rOA. No statistically significant associations were seen for incident or progressive knee rOA.

Conclusion

Variations in hip shape were associated with prevalent, but not incident or progressive, knee rOA in this cohort, and may reflect biomechanical differences between limbs, genetic influences, or common factors related to both hip shape and knee rOA.

Osteoarthritis (OA) risk is widely recognized to be heritable but few loci have been identified. Observational studies have identified higher systemic bone mineral density (BMD) to be associated with an increased risk of radiographic knee osteoarthritis. With this in mind, we sought to evaluate whether well-established genetic loci for variance in BMD are associated with risk for radiographic OA in the Osteoarthritis Initiative (OAI) and the Johnston County Osteoarthritis (JoCo) Project. Cases had at least one knee with definite radiographic OA, defined as the presence of definite osteophytes with or without joint space narrowing (Kellgren-Lawrence [KL] grade ≥ 2) and controls were absent for definite radiographic OA in both knees (KL grade ≤ 1 bilaterally). There were 2014 and 658 Caucasian cases, respectively, in the OAI and JoCo Studies, and 953 and 823 controls. Single nucleotide polymorphisms (SNPs) were identified for association analysis from the literature. Genotyping was carried out on Illumina 2.5M and 1M arrays in Genetic Components of Knee OA (GeCKO) and JoCo, respectively and imputation was done. Association analyses were carried out separately in each cohort with adjustments for age, body mass index (BMI), and sex, and then parameter estimates were combined across the two cohorts by meta-analysis. We identified four SNPs significantly associated with prevalent radiographic knee OA. The strongest signal (p = 0.0009; OR = 1.22; 95% CI, 1.08-1.37) maps to 12q3, which contains a gene coding for SP7. Additional loci map to 7p14.1 (TXNDC3), 11q13.2 (LRP5), and 11p14.1 (LIN7C). For all four loci the allele associated with higher BMD was associated with higher odds of OA. A BMD risk allele score was not significantly associated with OA risk. This meta-analysis demonstrates that several genomewide association studies (GWAS)-identified BMD SNPs are nominally associated with prevalent radiographic knee OA and further supports the hypothesis that BMD, or its determinants, may be a risk factor contributing to OA development. © 2014 American Society for Bone and Mineral Research.

Background

Osteoarthritis (OA) is a chronic joint disease, with increasing global burden of disability and healthcare utilisation. Recent meta-analyses have shown a range of effects of OA on mortality, reflecting different OA definitions and study methods. We seek to overcome limitations introduced when using aggregate results by gathering individual participant-level data (IPD) from international observational studies and standardising methods to determine the association of knee OA with mortality in the general population.

Methods

Seven community-based cohorts were identified containing knee OA-related pain, radiographs, and time-to-mortality, six of which were available for analysis. A two-stage IPD meta-analysis framework was applied: (1) Cox proportional hazard models assessed time-to-mortality of participants with radiographic OA (ROA), OA-related pain (POA), and a combination of pain and ROA (PROA) against pain and ROA-free participants; (2) hazard ratios (HR) were then pooled using the Hartung-Knapp modification for random-effects meta-analysis.

Findings

10,723 participants in six cohorts from four countries were included in the analyses. Multivariable models (adjusting for age, sex, race, BMI, smoking, alcohol consumption, cardiovascular disease, and diabetes) showed a pooled HR, compared to pain and ROA-free participants, of 1.03 (0.83, 1.28) for ROA, 1.35 (1.12, 1.63) for POA, and 1.37 (1.22, 1.54) for PROA.

Discussion

Participants with POA or PROA had a 35-37% increased association with reduced time-to-mortality, independent of confounders. ROA showed no association with mortality, suggesting that OA-related knee pain may be driving the association with time-to-mortality.

Funding

Versus Arthritis Centre for Sport, Exercise and Osteoarthritis and Osteoarthritis Research Society International.

Objective

A major barrier to genetic studies of osteoarthritis (OA) is the need to obtain large numbers of individuals with standardized radiographic evaluations for OA. To address this gap, we performed a genome-wide association study (GWAS) of radiographically defined tibiofemoral knee OA in 3,898 cases and 3,168 controls from 4 well-characterized North American cohorts, and we performed replication analysis of previously reported OA loci.

Methods

We performed meta-analysis using a 2-stage design. Stage 1 (discovery) consisted of a GWAS meta-analysis of radiographic knee OA carried out in the Osteoarthritis Initiative and the Johnston County Osteoarthritis Project. Knee OA was defined as definite osteophytes and possible joint space narrowing or total joint replacement in one or both knees. Stage 2 (validation) was performed in the Multicenter Osteoarthritis Study and the Genetics of Osteoarthritis study. We genotyped lead meta-analysis variants (P ≤ 1 × 10^-4 ) from stage 1 and tested the association between these variants and knee OA. We then combined results from all cohorts in a meta-analysis.

Results

Lead variants from stage 1, representing 49 unique loci, were analyzed in stage 2; none met genome-wide significance in the combined analysis of stage 1 and stage 2. We validated 1 locus (rs4867568 near LSP1P3) with nominal significance (P < 0.05), which was also our top finding in the combined meta-analysis (odds ratio [OR] 0.84 [95% confidence interval (95% CI) 0.79-0.91], P = 3.02 × 10^-6 ). We observed nominally significant associations (P < 0.05) with 3 previously reported OA loci: rs143383 in GDF5 (OR 1.12 [95% CI 1.04-1.21], P = 2.13 × 10^-3 ), rs835487 in CHST11 (OR 0.93 [95% CI 0.85-0.99], P = 0.03), and rs8044769 in FTO (OR 1.10 [95% CI 1.03-1.19], P = 6.13 × 10^-3 ).

Conclusion

These findings provide suggestive evidence of a novel knee OA locus and confirm previously reported associations in GDF5, CHST11, and FTO.

Background

Osteoporosis is a major public health problem associated with excess disability and mortality. It is estimated that 50-70% of the variation in osteoporotic fracture risk is attributable to genetic factors. The purpose of this hypothesis-generating study was to identify possible genetic determinants of fracture among African American (AA) women in a GWAS meta-analysis.

Methods

Data on clinical fractures (all fractures except fingers, toes, face, skull or sternum) were analyzed among AA female participants in the Women's Health Initiative (WHI) (N = 8155), Cardiovascular Health Study (CHS) (N = 504), BioVU (N = 704), Health ABC (N = 651), and the Johnston County Osteoarthritis Project (JoCoOA) (N = 291). Affymetrix (WHI) and Illumina (Health ABC, JoCoOA, BioVU, CHS) GWAS panels were used for genotyping, and a 1:1 ratio of YRI:CEU HapMap haplotypes was used as an imputation reference panel. We used Cox proportional hazard models or logistic regression to evaluate the association of ~ 2.5 million SNPs with fracture risk, adjusting for ancestry, age, and geographic region where applicable. We conducted a fixed-effects, inverse variance-weighted meta-analysis. Genome-wide significance was set at P < 5 × 10^- 8.

Results

One SNP, rs12775980 in an intron of SVIL on chromosome 10p11.2, reached genome-wide significance (P = 4.0 × 10^- 8). Although this SNP has a low minor allele frequency (0.03), there was no evidence for heterogeneity of effects across the studies (I² = 0). This locus was not reported in any previous osteoporosis-related GWA studies. We also interrogated previously reported GWA-significant loci associated with fracture or bone mineral density in our data. One locus (SMOC1) generalized, but overall there was not substantial evidence of generalization. Possible reasons for the lack of generalization are discussed.

Conclusion

This GWAS meta-analysis of fractures in African American women identified a potentially novel locus in the supervillin gene, which encodes a platelet-associated factor and was previously associated with platelet thrombus formation in African Americans. If validated in other populations of African descent, these findings suggest potential new mechanisms involved in fracture that may be particularly important among African Americans.

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.

Background

Cut-points to define slow walking speed have largely been derived from expert opinion.

Methods

Study participants (13 589 men and 5043 women aged ≥65years) had walking speed (m/s) measured over 4-6 m (mean ± SD: 1.20 ± 0.27 m/s in men and 0.94 ± 0.24 m/s in women.) Mobility limitation was defined as any self-reported difficulty with walking approximately 1/4 mile (prevalence: 12.6% men, 26.4% women). Sex-stratified classification and regression tree (CART) models with 10-fold cross-validation identified walking speed cut-points that optimally discriminated those who reported mobility limitation from those who did not.

Results

Among 5043 women, CART analysis identified 2 cut-points, classifying 4144 (82.2%) with walking speed ≥0.75 m/s, which we labeled as "fast"; 478 (9.5%) as "intermediate" (walking speed ≥0.62 m/s but <0.75 m/s); and 421 (8.3%) as "slow" (walking speed <0.62 m/s). Among 13 589 men, CART analysis identified 3 cut-points, classifying 10 001 (73.6%) with walking speed ≥1.00 m/s ("very fast"); 2901 (21.3%) as "fast" (walking speed ≥0.74 m/s but <1.00 m/s); 497 (3.7%) as "intermediate" (walking speed ≥0.57 m/s but <0.74 m/s); and 190 (1.4%) as "slow" (walking speed <0.57 m/s). Prevalence of self-reported mobility limitation was lowest in the "fast" or "very fast" (11% for men and 19% for women) and highest in the "slow" (60.5% in men and 71.0% in women). Rounding the 2 slower cut-points to 0.60 m/s and 0.75 m/s reclassified very few participants.

Conclusions

Cut-points in walking speed of approximately 0.60 m/s and 0.75 m/s discriminate those with self-reported mobility limitation from those without.

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.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).