Prevalence, incidence, and re-occurrence risk of musculoskeletal pain in older adults in the United Kingdom: a population-based study

Background

Throughout the literature, pain burden has been assessed by asking different questions, often cross-sectionally, different populations of interest. We know little about pain re-occurrence and how to translate knowledge between pain questions within the population of interest. We aimed to estimate the burden of musculoskeletal pain by estimating prevalence, incidence rates, and re-occurrence risk of back, hand, hip, knee, and foot pain using different questions from UK population-based samples and predict the number of affected individuals in the UK in 2030.

Methods

We used two UK population-representative studies, with two eight-year-apart follow-ups and two pain questions assessing recent pain episodes and often troubled pain when walking. We estimated prevalence, 8-year incidence rates, and 8-year pain re-occurrence risk for women and men aged 50 years and older and the relation between the two pain questions.

Results

Among UK individuals older than 50 years, the prevalence of musculoskeletal pain episode was 20%-50%, and the incidence was 20-40/1,000 person-years, while the prevalence of pain when walking was 10%-25%, and the incidence was 6-12/1,000 person-years. The most prevalent musculoskeletal pain types were back and knee pain; of five women experiencing back or knee pain episodes, three are expected to be often troubled by pain. Hip and foot pain had similar estimates in both questions. Hand pain peaked in women aged 50-65 years. Women had higher prevalence and incidence rates, but men had higher 8-year re-occurrence risk of all types of musculoskeletal pain. Reporting a pain episode was associated with two times higher risk, but often troubled by pain when walking was associated with four to seven times times higher risk of the same pain in 8 years. Women and men with a body mass index (BMI) of ≥27 kg/m2 were twice as likely to experience musculoskeletal pain than those with BMI<27 kg/m2. In 2030, we expect 2-7 million people older than 50 years in the United Kingdom to seek site-specific musculoskeletal pain-focused healthcare.

Conclusions

In individuals older than 50 years, the experience of musculoskeletal pain at least doubles the chance of experiencing it again. Women report musculoskeletal pain more often, but men report more persistent pain. Musculoskeletal pain presents a significant burden to public health.

Effects of body weight and fat mass on back pain - direct mechanical or indirect through inflammatory and metabolic parameters?

BACKGROUND

While reports indicate the association between obesity and back pain, its mechanism is still unclear. Thus, we aimed to investigate the effects of weight and its components on back pain in middle-aged women while considering direct mechanical and indirect effects via inflammatory and metabolic parameters.

METHODS

We used data from the Chingford 1000 Women Study, two follow-ups seven years apart. We assessed effects of weight, body mass index (BMI), total fat mass (TFM), total lean mass (TLM) and total bone mineral density (TBMD), measured by dual-energy X-ray absorptiometry, on back pain episode. We used inflammatory (C-reactive protein, interleukin-6, and tumour necrosis factor-alpha) and metabolic parameters (systolic and diastolic blood pressure, triglyceride, high-density lipoprotein cholesterol, and fasting blood glucose) as mediators of indirect effects. We investigated associations of interest cross-sectionally and longitudinally using binary logistic regression and parallel mediation model.

RESULTS

We included 826 Chingford middle-aged women (mean age=60.7, SD=5.9) from the first used follow-up in cross-sectional and mediation analyses and 645 women that attended the follow-up seven years later, in longitudinal analyses. We found that increased weight was directly associated with increased odds of having back pain episode (OR=1.02; 95% CI 1.01-1.03), similarly as BMI (OR=1.05; 95% CI 1.02-1.08) and TFM (OR=1.03; 95% CI 1.01-1.04) consistently across the cross-sectional and longitudinal models, but not TLM or TBMD. However, we did not find consistent indirect effects of weight or its components through measured inflammatory or metabolic parameters on back pain.

CONCLUSIONS

Our results show that in middle-aged women, weight, BMI and TFM are directly related to back pain, indicating prominence of mechanical loading effect.

Specific body mass index trajectories were related to musculoskeletal pain and mortality: 19-year follow-up cohort

OBJECTIVE

We aimed to study 19-year body mass index (BMI) patterns and their (1) bidirectional relationship with musculoskeletal pain and (2) mortality risk.

STUDY DESIGN AND SETTING

We used data from the Chingford study and group-based trajectory modelling to define 19-year BMI patterns. We investigated whether baseline back, hand, hip, and knee pain (as single- and multi-site) predicted 19-year BMI trajectory, and whether 19-year BMI patterns predicted pain in year 20. We explored BMI trajectories and mortality risk over 25 years (life expectancy).

RESULTS

We included 938 women (mean age: year-1=54, year-20=72) and found seven distinct 19-year BMI trajectories: two normal-weighted (reference), slightly overweight, lower and upper overweight-to-obese, lower and upper obese. BMI patterns capturing the increase overweight-to-obese (BMI 27-34 overtime) were bidirectionally related to knee and multi-site pain. The lower obese pattern (BMI 33-38) was unidirectionally associated with lower limb pain. Women with BMI above 40 had an increased all-cause and cardiovascular mortality risk.

CONCLUSION

For most postmenopausal women, the overweight WHO category was a transition. Two patterns capturing increase overweight-to-obese were mutually related to musculoskeletal pain, i.e., knee and multi-site pain contributed to becoming obese, and trajectories of becoming obese increased the odds of experiencing pain later.

Effect of Heeled Shoes on Joint Symptoms and Knee Osteoarthritis in Older Adults: A 5-Year Follow-Up Study

OBJECTIVE

Our aims were to examine the effects of heeled shoes on incident knee osteoarthritis (OA) and joint pain.

METHODS

We used longitudinal data from the Chingford 1000 Women Study (Chingford Study), a prospective cohort of women aged 50 years or older. Participants with musculoskeletal disorders and/or a history of knee-related injury/surgery were excluded. Participants were followed for up to 5 years for incident outcomes including 1) radiographic knee OA (RKOA) and 2) joint pain (feet, knees, hips, and back). Footwear data, including ever worn heels of 2 inches or more and daytime/evening hours (per week) spent wearing heeled shoes over five decades (ages <20 years, 20-30 years, 30-40 years, and >50 years), were available at Year 10 whereas knee radiographs and joint symptom data were also collected at Year 15. Cumulative time spent wearing heeled shoes was calculated for women reporting ever-use of heeled shoes (≥2 inches). Multiple logistic regression was used to examine the relationship between exposures and outcomes (from Year 10 to Year 15).

RESULTS

A total of 356 women were eligible at Year 10 with a median (interquartile range) age of 60 (56-65) years. Compared with non-use, ever-use of heeled shoes (≥2 inches) was not associated with incident RKOA (1.35; 95% confidence interval: 0.56-3.27). No associations were observed between increasing cumulative time spent wearing heels and incident outcomes.

CONCLUSION

Compared with the non-use of heeled shoes, ever-use of heels (≥2 inches) was not associated with incident RKOA and incident joint symptoms. Further, increasing cumulative time spent wearing heels was not associated with any of our outcomes.

Association of Lumbar Spine Radiographic Changes With Severity of Back Pain-Related Disability Among Middle-aged, Community-Dwelling Women

IMPORTANCE

Previous studies, using mostly cross-sectional data, provide conflicting evidence of an association between lumbar spine radiographic changes and the severity of back pain-related disability. Such conflicting evidence may be associated with widely unnecessary diagnostic imaging of the lumbar spine.

OBJECTIVE

To examine both cross-sectional and longitudinal associations between lumbar spine radiographic changes and the severity of back pain-related disability among middle-aged, community-dwelling women.

DESIGN, SETTING, AND PARTICIPANTS

This population-based prospective cohort study used data from the Chingford 1000 Women Study. Analyses included data collected from year 6 (1994-1996; physical activity was measured), year 9 (1997-1999; treated as baseline), and year 15 (2003-2005), with a total length of follow-up for longitudinal analyses of 6 years. Data were analyzed from April 17 to November 3, 2020.

EXPOSURES

Primary exposure was lumbar spine radiographic changes, defined using the Kellgren-Lawrence (K-L) grade. Secondary exposures were defined using presence of osteophytes and disc space narrowing. The composite score combined the number of lumbar spine segments with definite changes detected on radiographic images (ie, radiographic changes) (K-L grade ≥2, which means at least definite osteophyte and possible narrowing of disc space are present; osteophyte and disc space narrowing grade ≥1, which means at least mild or definite changes are present).

MAIN OUTCOMES AND MEASURES

Self-reported back pain-related disability measured in years 9 and 15 assessed by the St Thomas disability questionnaire.

RESULTS

Among 650 women (mean [SD] age, 61.3 [5.9] years) in cross-sectional analyses and 443 women (mean [SD] age, 60.6 [6.0] years) in longitudinal analyses, there was no evidence to support an association between higher number of lumbar segments with radiographic changes (K-L grade, osteophytes, and disc space narrowing) and more severe back pain-related disability (eg, cross-sectional analyses using the K-L grade; 1 segment vs 0 segment: adjusted odds ratio, 1.22 [95% CI, 0.76-1.96]). No interactions were found of an association between lumbar spine radiographic changes and the severity of back pain-specific disability with age, body mass index, or smoking status.

CONCLUSIONS AND RELEVANCE

In this cohort of middle-aged, community-dwelling women, there was no evidence to support an association between a higher number of lumbar segments with radiographic changes (K-L grade, osteophytes, and disc space narrowing) and more severe back pain-related disability cross-sectionally or over time. These findings provide further evidence against routinely using diagnostic imaging of the lumbar spine.

Genome-wide association and functional studies identify a role for matrix Gla protein in osteoarthritis of the hand

OBJECTIVE

Osteoarthritis (OA) is the most common form of arthritis and the leading cause of disability in the elderly. Of all the joints, genetic predisposition is strongest for OA of the hand; however, only few genetic risk loci for hand OA have been identified. Our aim was to identify novel genes associated with hand OA and examine the underlying mechanism.

METHODS

We performed a genome-wide association study of a quantitative measure of hand OA in 12 784 individuals (discovery: 8743, replication: 4011). Genome-wide significant signals were followed up by analysing gene and allele-specific expression in a RNA sequencing dataset (n=96) of human articular cartilage.

RESULTS

We found two significantly associated loci in the discovery set: at chr12 (p=3.5 × 10-10) near the matrix Gla protein (MGP) gene and at chr12 (p=6.1×10-9) near the CCDC91 gene. The DNA variant near the MGP gene was validated in three additional studies, which resulted in a highly significant association between the MGP variant and hand OA (rs4764133, Betameta=0.83, Pmeta=1.8*10-15). This variant is high linkage disequilibrium with a coding variant in MGP, a vitamin K-dependent inhibitor of cartilage calcification. Using RNA sequencing data from human primary cartilage tissue (n=96), we observed that the MGP RNA expression of the hand OA risk allele was significantly lowercompared with the MGP RNA expression of the reference allele (40.7%, p<5*10-16).

CONCLUSIONS

Our results indicate that the association between the MGP variant and increased risk for hand OA is caused by a lower expression of MGP, which may increase the burden of hand OA by decreased inhibition of cartilage calcification.

Pain in knee osteoarthritis is associated with variation in the neurokinin 1/substance P receptor (TACR1) gene

BACKGROUND

Substance P (SP) is a pain- and inflammation-related neuropeptide which preferentially binds to the neurokinin receptor 1 (NK₁ ). SP and NK₁ receptors have been implicated in joint pain, inflammation and damage in animal models and human studies of osteoarthritis (OA). The aim of this study was to test if genetic variation at the neurokinin 1 receptor gene (TACR1) is associated with pain in individuals with radiographic knee OA.

METHODS

Participants from the Genetics of OA and Lifestyle study were used for the discovery group (n = 1615). Genotype data for six SNPs selected to cover most variation in the TACR1 gene were used to test for an association with symptomatic OA. Replication analysis was performed using data from the Chingford 1000 Women Study, Hertfordshire Cohort Study, Tasmanian Older Adult Cohort Study and the Clearwater OA Study. In total, n = 1715 symptomatic OA and n = 735 asymptomatic OA individuals were analysed.

RESULTS

Out of six SNPs tested in the TACR1 gene, one (rs11688000) showed a nominally significant association with a decreased risk of symptomatic OA in the discovery cohort. This was then replicated in four additional cohorts. After adjusting for age, gender, body mass index and radiographic severity, the G (minor) allele at rs11688000 was associated with a decreased risk of symptomatic OA compared to asymptomatic OA cases (p = 9.90 × 10^-4 , OR = 0.79 95% 0.68-0.90 after meta-analysis).

CONCLUSIONS

This study supports a contribution from the TACR1 gene in human OA pain, supporting further investigation of this gene's function in OA.

SIGNIFICANCE

This study contributes to the knowledge of the genetics of painful osteoarthritis, a condition which affects millions of individuals worldwide. Specifically, a contribution from the TACR1 gene to modulating pain sensitivity in osteoarthritis is suggested.

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.

Symptoms of Raynaud's phenomenon (RP) in fibromyalgia syndrome are similar to those reported in primary RP despite differences in objective assessment of digital microvascular function and morphology

Symptoms of Raynaud’s phenomenon (RP) are common in fibromyalgia syndrome (FMS). We compared symptom characteristics and objective assessment of digital microvascular function using infrared thermography (and nailfold capillaroscopy where available) in patients with FMS (reporting RP symptoms) and primary RP. We retrospectively reviewed the outcome of microvascular imaging studies and RP symptom characteristics (captured using patient-completed questionnaire at the time of assessment) for patients with FMS (reporting RP symptoms) and patients with primary RP referred for thermographic assessment of RP symptoms over a 2-year period. Of 257 patients referred for thermographic assessment of RP symptoms between 2010 and 2012, we identified 85 patients with primary RP and 43 patients with FMS. There were no differences in RP symptom characteristics between FMS and primary RP (p > 0.05 for all comparisons). In contrast, patients with FMS had higher baseline temperature of the digits (32.1 vs. 29.0 °C, p = 0.004), dorsum (31.9 vs. 30.2 °C, p = 0.005) and thermal gradient (temperature of digits minus temperature of dorsum; +0.0 vs. −0.9 °C, p = 0.03) compared with primary RP. Significant differences between groups persisted following local cold challenge. In primary RP, patient reporting “blue” digits, bi-phasic and tri-phasic RP was associated with lower digital perfusion. In contrast, no associations between skin temperature and RP digital colour changes/phases were identified in FMS. Our findings suggest that symptoms of RP in FMS may have a different aetiology to those seen in primary RP. These findings have potential implications for both the classification of RP symptoms and the management of RP symptoms in the context of FMS. Digital colour changes reported by patients might reflect the degree of digital microvascular compromise in primary RP.

Individuals with high bone mass have an increased prevalence of radiographic knee osteoarthritis

We previously reported an association between high bone mass (HBM) and a bone-forming phenotype of radiographic hip osteoarthritis (OA). As knee and hip OA have distinct risk factors, in this study we aimed to determine (i) whether HBM is also associated with knee OA, and (ii) whether the HBM knee OA phenotype demonstrates a similar pattern of radiographic features to that observed at the hip. HBM cases (defined by DXA BMD Z-scores) from the UK-based HBM study were compared with unaffected family controls and general population controls from the Chingford and Hertfordshire cohort studies. A single blinded observer graded AP weight-bearing knee radiographs for features of OA (Kellgren-Lawrence score, osteophytes, joint space narrowing (JSN), sclerosis) using an atlas. Analyses used logistic regression, adjusting a priori for age and gender, and additionally for BMI as a potential mediator of the HBM-OA association, using Stata v12. 609 HBM knees in 311 cases (mean age 60.8years, 74% female) and 1937 control knees in 991 controls (63.4years, 81% female) were analysed. The prevalence of radiographic knee OA, defined as Kellgren-Lawrence grade≥2, was increased in cases (31.5% vs. 20.9%), with age and gender adjusted OR [95% CI] 2.38 [1.81, 3.14], p<0.001. The association between HBM and osteophytosis was stronger than that for JSN, both before and after adjustment for BMI which attenuated the ORs for knee OA and osteophytes in cases vs. controls by approximately 50%. Our findings support a positive association between HBM and knee OA. This association was strongest for osteophytes, suggesting HBM confers a general predisposition to a subtype of OA characterised by increased bone formation.

The effect of FTO variation on increased osteoarthritis risk is mediated through body mass index: a Mendelian randomisation study

OBJECTIVE

Variation in the fat mass and obesity-associated (FTO) gene influences susceptibility to obesity. A variant in the FTO gene has been implicated in genetic risk to osteoarthritis (OA). We examined the role of the FTO polymorphism rs8044769 in risk of knee and hip OA in cases and controls incorporating body mass index (BMI) information.

METHODS

5409 knee OA patients, 4355 hip OA patients and up to 5362 healthy controls from 7 independent cohorts from the UK and Australia were genotyped for rs8044769. The association of the FTO variant with OA was investigated in case/control analyses with and without BMI adjustment and in analyses matched for BMI category. A mendelian randomisation approach was employed using the FTO variant as the instrumental variable to evaluate the role of overweight on OA.

RESULTS

In the meta-analysis of all overweight (BMI≥25) samples versus normal-weight controls irrespective of OA status the association of rs8044769 with overweight is highly significant (OR[CIs] for allele G=1.14 [01.08 to 1.19], p=7.5×10(-7)). A significant association with knee OA is present in the analysis without BMI adjustment (OR[CIs]=1.08[1.02 to 1.14], p=0.009) but the signal fully attenuates after BMI adjustment (OR[CIs]=0.99[0.93 to 1.05], p=0.666). We observe no evidence for association in the BMI-matched meta-analyses. Using mendelian randomisation approaches we confirm the causal role of overweight on OA.

CONCLUSIONS

Our data highlight the contribution of genetic risk to overweight in defining risk to OA but the association is exclusively mediated by the effect on BMI. This is consistent with what is known of the biology of the FTO gene and supports the causative role of high BMI in OA.

Osteophytes, enthesophytes, and high bone mass: a bone-forming triad with potential relevance in osteoarthritis

Objective

Previous studies of skeletal remains have suggested that both enthesophytes and osteophytes are manifestations of an underlying bone-forming tendency. A greater prevalence of osteophytes has been observed among individuals with high bone mass (HBM) compared with controls. This study was undertaken to examine the possible interrelationships between bone mass, enthesophytes, and osteophytes in a population of individuals with extreme HBM.

Methods

Cases of HBM (defined according to bone mineral density [BMD] Z scores on dual x-ray absorptiometry) from the UK-based HBM study were compared with a control group comprising unaffected family members and general population controls from the Chingford and Hertfordshire cohort studies. Pelvic radiographs from cases and controls were pooled and evaluated, in a blinded manner, by a single observer, who performed semiquantitative grading of the radiographs for the presence and severity of osteophytes and enthesophytes (score range 0–3 for each). Logistic regression analysis was used to identify significant associations, with a priori adjustment for age, sex, and body mass index.

Results

In this study, 226 radiographs from HBM cases and 437 radiographs from control subjects were included. Enthesophytes (grade ≥1) and moderate enthesophytes (grade ≥2) were more prevalent in HBM cases compared with controls (adjusted odds ratio [OR] 3.00 [95% confidence interval (95% CI) 1.96–4.58], P < 0.001 for any enthesophyte; adjusted OR 4.33 [95% CI 2.67–7.02], P < 0.001 for moderate enthesophytes). In the combined population of cases and controls, the enthesophyte grade was positively associated with BMD at both the total hip and lumbar spine (adjusted P for trend < 0.001). In addition, a positive association between osteophytes and enthesophytes was observed; for each unit increase in enthesophyte grade, the odds of any osteophyte being present were increased >2-fold (P < 0.001).

Conclusion

Strong interrelationships were observed between osteophytes, enthesophytes, and HBM, which may be helpful in defining a distinct subset of patients with osteoarthritis characterized by excess bone formation.

Meta-analysis identifies loci affecting levels of the potential osteoarthritis biomarkers sCOMP and uCTX-II with genome wide significance

BACKGROUND

Research for the use of biomarkers in osteoarthritis (OA) is promising, however, adequate discrimination between patients and controls may be hampered due to innate differences. We set out to identify loci influencing levels of serum cartilage oligomeric protein (sCOMP) and urinary C-telopeptide of type II collagen (uCTX-II).

METHODS

Meta-analysis of genome-wide association studies was applied to standardised residuals of sCOMP (N=3316) and uCTX-II (N=4654) levels available in 6 and 7 studies, respectively, from TreatOA. Effects were estimated using a fixed-effects model. Six promising signals were followed up by de novo genotyping in the Cohort Hip and Cohort Knee study (N = 964). Subsequently, their role in OA susceptibility was investigated in large-scale genome-wide association studies meta-analyses for OA. Differential expression of annotated genes was assessed in cartilage.

RESULTS

Genome-wide significant association with sCOMP levels was found for a SNP within MRC1 (rs691461, p = 1.7 × 10(-12)) and a SNP within CSMD1 associated with variation in uCTX-II levels with borderline genome-wide significance (rs1983474, p = 8.5 × 10(-8)). Indication for association with sCOMP levels was also found for a locus close to the COMP gene itself (rs10038, p = 7.1 × 10(-6)). The latter SNP was subsequently found to be associated with hip OA whereas COMP expression appeared responsive to the OA pathophysiology in cartilage.

CONCLUSIONS

We have identified genetic loci affecting either uCTX-II or sCOMP levels. The genome wide significant association of MRC1 with sCOMP levels was found likely to act independent of OA subtypes. Increased sensitivity of biomarkers with OA may be accomplished by taking genetic variation into account.

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.

Association between cartilage and bone biomarkers and incidence of radiographic knee osteoarthritis (RKOA) in UK females: a prospective study

OBJECTIVE

There is a need to find biochemical markers that would identify people with increased risk of developing radiographic knee osteoarthritis (RKOA). The aim of this study was to evaluate the ability of cartilage and bone biomarkers (cartilage oligomeric matrix protein (COMP), aggrecan, cellular inhibitor of apoptosis protein (cIAP), N-telopeptide-to-helix (NTx)) to predict RKOA incidence in a 10-year follow-up of UK females from the Chingford community study.

METHOD

Joint space narrowing (JSN), osteophytes (OSP) and Kellgren-Lawrence (K/L) grades were scored from radiographs of both knees at study baseline and 10 years later in 1,003 women aged 45-64. Circulating levels of biomarkers and demographic variables were measured at baseline. Statistical association analysis was conducted between the potential predictor factors measured at baseline and documentation of RKOA at 10-year follow-up.

RESULTS

Age and body mass index (BMI), were significant predictors of incidence of RKOA as assessed by K/L and OSP. Considering biomarkers, independent significant association was found between COMP circulating levels and K/L scores (Odd Ratio (OR) = 2.87, 95% Confidence Interval (CI) = 1.19-6.89, P = 0.018). Significant negative association was detected between aggrecan plasma concentrations and JSN, with OR = 0.37 (95% CI 0.15-0.89), P = 0.026.

CONCLUSIONS

Aggrecan and COMP circulating levels contribute to identification of phenotype-specific RKOA incidence. These data suggest potentially protective role of aggrecan in cartilage loss, as measured by JSN. High COMP levels are risk factors for development of RKOA, as assessed by K/L scores.

Neuropathic features of joint pain: a community-based study

OBJECTIVE

Quantitative sensory testing (QST) and questionnaire-based assessments have been used to demonstrate features of neuropathic pain in subjects with musculoskeletal pain. However, their direct relationship has not been investigated in the community. The purpose of this study was to conduct an observational study to describe the characteristics of joint pain and to examine the relationship between QST measures and the PainDETECT Questionnaire (PD-Q).

METHODS

Warm detection, heat pain, and mechanical pain thresholds as well as mechanical pain sensitivity over the sternum were determined and the PD-Q scores were calculated in a cross-sectional study of 462 participants in the Chingford Study. Comparisons were made between subjects with and those without joint pain. Logistic regression modeling was used to describe the association between neuropathic pain features, as determined by the PD-Q score, and each of the QST measures individually, adjusting for age, body mass index, and use of pain-modifying medications.

RESULTS

A total of 66.2% of the subjects reported recent joint pain, with a median average pain severity of 5 of 10. There was increased sensitivity to painful stimuli in the group with pain as compared to the pain-free group, and this persisted after stratification by pain-modifying medication use. While only 6.7% of subjects had possible neuropathic pain features and 1.9% likely neuropathic pain features according to the standard PD-Q thresholds, features of neuropathic pain were common and were present in >50% of those reporting pain of at least moderate severity. Heat pain thresholds and mechanical pain sensitivity were significantly associated with features of neuropathic pain identified using the PD-Q, with an odds ratio (OR) of 0.88 (95% confidence interval [95% CI] 0.79-0.97; P = 0.011) and an OR of 1.24 (95% CI 1.04-1.48; P = 0.018), respectively.

CONCLUSION

QST measures and the PD-Q identified features of neuropathic pain in subjects in this community-based study, with significant overlap between the findings of the two techniques.

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.

The association between hip morphology parameters and nineteen-year risk of end-stage osteoarthritis of the hip: a nested case-control study

Objective

Subtle deformities of the hip joint are implicated in the etiology of osteoarthritis (OA) of the hip. Parameters that quantify these deformities may aid understanding of these associations. We undertook this study to examine relationships between such parameters and the 19-year risk of total hip arthroplasty (THA) for end-stage OA.

Methods

A new software program designed for measuring morphologic parameters around the hip was developed and validated in a reliability study. THA was the outcome measure for end-stage OA. A nested case–control study was used with individuals from a cohort of 1,003 women who were recruited at year 1 in 1989 and followed up to year 20 (the Chingford Study). All hips with THA by year 20 and 243 randomly selected control hips were studied. Pelvis radiographs obtained at year 2 were analyzed for variations in hip morphology. Measurements were compared between the THA case group and the control group.

Results

Patients with THA had a higher prevalence of cam deformity than did their respective controls (median alpha angle 62.4° versus 45.8° [P = 0.001]; mean modified triangular index height 28.5 mm versus 26.9 mm [P = 0.001]) as well as a higher prevalence of acetabular dysplasia (mean lateral center edge angle 29.5° versus 34.3° [P = 0.001]; median extrusion index 0.25 versus 0.185 [P = 0.009]). Logistic regression analyses clustering by subject and adjusting for radiographic hip OA at year 2 showed that these morphologic parameters were still significantly associated with THA by year 20. The alpha angle and lateral center edge angle predicted the risk of THA independently when included in the same model.

Conclusion

This investigation describes measurements that predict the risk of THA for end-stage OA by year 20, independently of the presence of radiographic hip OA at year 2. These measurements can be made on an anteroposterior pelvis radiograph, which is an inexpensive and commonly used clinical method of investigation.

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.

Lower limbs composition and radiographic knee osteoarthritis (RKOA) in Chingford sample--a longitudinal study

Our aim in this longitudinal study was to evaluate to what extent fat and lean tissue mass variations are associated and can predict RKOA in a large sample of British women followed-up over 10 years. Kellgren/Lawrence (K/L), joint space narrowing (JSN) and osteophyte (OSP) grades were scored from radiographs of both knees in 909 middle-aged women from the Chingford registry. Body composition components were assessed using the dual energy X-ray absorptiometry (DXA) method. In cross-sectional analysis, combined effect of age, BMI and leg tissue composition was required for best fitting model explaining variations of K/L scoring and osteophytes at lateral compartment. To explain medial osteophytes, age and BMI were sufficient to generate the best fitting model. In prediction analysis, leg lean mass was the more powerful predictor of K/L, medial osteophytes than BMI. In conclusion, BMI appears to influence the development of knee OA through both fat and/or lean mass, depending on RKOA phenotype.

Effects of age on genetic influence on bone loss over 17 years in women: the Healthy Ageing Twin Study (HATS)

The rate of bone loss varies across the aging period via multiple complex mechanisms. Therefore, the role of genetic factors on bone loss may also change similarly. In this study, we investigated the effect of age on the genetic component of bone loss in a large twin-based longitudinal study. During 17 years of follow-up in TwinsUK and Healthy Ageing Twin Study (HATS), 15,491 hip and lumbar spine dual-energy X-ray absorptiometry (DXA) scans were performed in 7056 twins. Out of these subjects, 2716 female twins aged >35 years with at least two scans separated for >4 years (mean follow-up 9.7 years) were included in this analysis. We used a mixed-effects random-coefficients regression model to predict hip and spine bone mineral density (BMD) values for exact ages of 40, 45, 50, 55, 60, 65, 70, 75, and 80 years, with adjustment for baseline age, weight, height, and duration of hormone replacement therapy. We then estimated heritability of the changes in BMD measures between these age ranges. Heritability estimates for cross-sectional hip and spine BMD were high (ranging between 69% and 88%) at different ages. Heritability of change of BMD was lower and more variable, generally ranging from 0% to 40% for hip and 0% to 70% for spine; between age 40 and 45 years genetic factors explained 39.9% (95% confidence interval [CI], 25%-53%) of variance of BMD loss for total hip, 46.4% (95% CI, 32%-58%) for femoral neck, and 69.5% (95% CI, 59%-77%) for lumbar spine. These estimates decreased with increasing age, and there appeared to be no heritability of BMD changes after the age of 65 years. There was some evidence at the spine for shared genetic effects between cross-sectional and longitudinal BMD. Whereas genetic factors appear to have an important role in bone loss in early postmenopausal women, nongenetic mechanisms become more important determinants of bone loss with advanced age.

Longitudinal study of variation in body mass index in middle-aged UK females

The importance of changing patterns of obesity in society and its implications for public health are well recognized. However, the adult life course of body mass index (BMI) changes in individuals over time is largely unknown and has mostly been extrapolated from cross-sectional studies. The present study examines individual specific variation of BMI during a 15-year follow-up period in a community-based sample of UK females. We attempted to establish whether there is a common, generalized pattern which captures variation in BMI over time. The participants of this study belong to a prospective population cohort of British women studied intensively since 1989: the Chingford Study. The sample originally consisted of 1,003 women aged 45-68 years, who were assessed annually for BMI during follow-up period. Polynomial regression models were used to assess longitudinal BMI variation. We observed a great stability in individual BMI variation during the follow-up period, reflected by high correlations between the baseline BMI and follow-up BMI 10 and 15 years later (r = 0.876, N = 810, and r = 0.824, N = 638, respectively). We also found that three different major age-related patterns in BMI could be clearly identified: no change in 30.6% in 58% it increased and in 11.4% it decreased with age. Thus, our data suggest that individual age-related changes in BMI are very different. Therefore, simply combining all individuals into groups by any other criteria (age, sex, etc.) and overlooking the distinctive patterns of BMI change may lead to biased inferences in epidemiologic and etiologic research of the future.

The natural history of radiographic knee osteoarthritis: a fourteen-year population-based cohort study

OBJECTIVE

To establish the natural history of radiographic knee osteoarthritis (OA) over 14 years in a community-based cohort.

METHODS

We examined women from the Chingford Women's Study, a community-based cohort followed up for more than 14 years. We selected women for whom bilateral radiographs of the knees (with the legs in full extension) were obtained at approximately 5-year intervals. Radiographs were scored for OA in a blinded manner, using Kellgren/Lawrence (K/L) grades. Descriptive statistics and odds ratios (ORs) were used to compare the incidence, worsening, and progression of radiographic knee OA.

RESULTS

A complete radiography series was available for 561 of the original 1,003 subjects enrolled in the study. The median age of these subjects at baseline was 53 years (interquartile range 48-58 years). At baseline, 13.7% of the subjects had radiographic knee OA (K/L grade≥2) in at least one knee, and the prevalence increased to 47.8% by year 15. The annual cumulative incidence of radiographic knee OA was 2.3% between baseline and year 15. The annual rates of disease progression and worsening between baseline and year 15 were 2.8% and 3.0%, respectively. Subjects with a K/L grade of 1 at baseline were more likely to experience worsening by year 15 compared with subjects with a baseline grade of 0 (OR 4.5, 95% confidence interval 2.7-7.4).

CONCLUSION

This is the longest natural history study of radiographic knee OA to date. The results showed relatively low rates for the incidence and progression of radiographic knee OA; more than half of all subjects had no radiographic evidence of knee OA over a 15-year period of time. Subjects with a baseline K/L grade of 1 were more likely than subjects with other baseline K/L grades to experience worsening of knee OA.

A role for PACE4 in osteoarthritis pain: evidence from human genetic association and null mutant phenotype

OBJECTIVES

The aim of this study was to assess if genetic variation in the PACE4 (paired amino acid converting enzyme 4) gene Pcsk6 influences the risk for symptomatic knee osteoarthritis (OA).

METHODS

Ten PCSK6 single nucleotide polymorphisms were tested for association in a discovery cohort of radiographic knee OA (n=156 asymptomatic and 600 symptomatic cases). Meta-analysis of the minor allele at rs900414 was performed in three additional independent cohorts (total n=674 asymptomatic and 2068 symptomatic). Pcsk6 knockout mice and wild-type C57BL/6 mice were compared in a battery of algesiometric assays, including hypersensitivity in response to intraplantar substance P, pain behaviours in response to intrathecal substance P and pain behaviour in the abdominal constriction test.

RESULTS

In the discovery cohort of radiographic knee OA, an intronic single nucleotide polymorphism at rs900414 was significantly associated with symptomatic OA. Replication in three additional cohorts confirmed that the minor allele at rs900414 was consistently increased among asymptomatic compared to symptomatic radiographic knee OA cases in all four cohorts. A fixed-effects meta-analysis yielded an OR=1.35 (95% CI 1.17 to 1.56; p=4.3×10(-5) and no significant between-study heterogeneity). Studies in mice revealed that Pcsk6 knockout mice were significantly protected against pain in a battery of algesiometric assays.

CONCLUSIONS

These results suggest that a variant in PCSK6 is strongly associated with protection against pain in knee OA, offering some insight as to why, in the presence of the same structural damage, some individuals develop chronic pain and others are protected. Studies in Pcsk6 null mutant mice further implicate PACE4 in pain.

Are Twins and Singletons Comparable? A Study of Disease-related and Lifestyle Characteristics in Adult Women

The classic twin study is sometimes described as “the perfect natural experiment” for the investigation of the aetiology of complex disease, but assumptions of the twin design need to be empirically tested if their results are to be considered unbiased and representative of singleton populations. In this study comparisons of disease and prevalence of lifestyle characteristics have been made between twin participants in the St Thomas' Hospital UK adult twin registry, the largest twin volunteer register in the UK for the study of diseases of ageing, and a parallel population-based study of singleton women. The only differences found were for weight, where monozygotic (MZ) twins were lighter and had a smaller variance than dizygotic (DZ) twins and singletons. For the other variables studied, volunteer twins were not found to differ from age-matched singleton women in distribution or prevalence of: bone mineral density, osteoarthritis, blood pressure, hypertensive drug use, height, history of hysterectomy and ovariectomy, menopausal status and current alcohol and overall tobacco consumption. We conclude that the results of twin studies can be generalised to singleton populations for these measures and disease outcomes.

Prevalence of reported knee pain over twelve years in a community-based cohort

OBJECTIVE

To describe the temporal patterns of knee pain in a community-based cohort over 12 years.

METHODS

Data on self-reported knee pain at 4 time points over 12 years were analyzed in participants from the Chingford Women's Study of osteoarthritis (OA) and osteoporosis. Pain status was defined as any pain in the preceding month and pain on most days in the preceding month. This status was used to classify participants according to pain patterns of asymptomatic, persistent, incident, or intermittent pain. Multinomial logistic regression was used to identify baseline predictors for each pain pattern.

RESULTS

Among the 489 women with complete followup data, the median age at baseline was 52 years (interquartile range [IQR] 48-58 years), the median body mass index (BMI) was 24.39 kg/m(2) (IQR 22.46-27.20), and 11.7% of the women had a Kellgren/Lawrence radiographic OA severity grade of ≥2 in at least one knee. Among subjects reporting any pain in the preceding month versus those reporting pain on most days in the preceding month, 9% versus 2% had persistent pain, 24% versus 16% had incident pain, and 29% versus 18% had intermittent pain. A higher BMI was predictive of persistent pain (odds ratio [OR] 1.14, 95% confidence interval [95% CI] 1.04-1.25) and incident pain (OR 1.10, 95% CI 1.02-1.18). The presence of radiographic knee OA was predictive of persistent pain (OR 3.70, 95% CI 1.34-10.28; P = 0.012), and reported knee injury was predictive of both persistent pain (OR 4.13, 95% CI 1.34-12.66; P = 0.013) and intermittent pain (OR 4.25, 95% CI 1.81-9.98; P = 0.001).

CONCLUSION

Significant variability in the temporal fluctuation of self-reported knee pain was seen in this community-based prospective study over a period of 12 years, with few women consistently reporting knee pain at each time point. Distinct baseline predictors for each pain pattern were identified and may explain the observed heterogeneity of self-reported knee pain when pain status is measured at only one time point.

Association between DHEAS and bone loss in postmenopausal women: a 15-year longitudinal population-based study

Our aim was to examine the association between serum dehydroepiandrosterone sulfate (DHEAS) at baseline and BMD change at the femoral neck (FN) and lumbar spine (LS) in postmenopausal women during a 15-year follow-up. All participants were from the Chingford Study. BMD at the FN and LS were measured eight times during the 15-year follow-up by dual-energy X-ray absorptiometry. DHEAS at baseline was measured using radioimmunoassay. Data on height, weight, and hormone-replacement therapy (HRT) status were obtained at each visit. Multilevel linear regression modeling was used to examine the association between longitudinal BMD change at the FN and LS and DHEAS at baseline. Postmenopausal women (n = 1,003) aged 45-68 years (mean 54.7) at baseline were included in the study. After adjustment for baseline age, estradiol, HRT, and BMI, BMD at the FN decreased on average 0.49% (95% CI 0.31-0.71%) per year; and the decline was slowed down by 0.028% per squared year. Increase of DHEAS (each micromole per liter) was associated with 0.49% less bone loss at the FN (95% CI 0.21-0.71%, P = 0.001). However, this strong association became slightly weaker over time. Similar but weaker results were obtained for LS BMD. Our data suggest that high serum DHEAS at baseline is associated with less bone loss at both FN and LS and this association diminishes over time. The nature of the association is unclear, but such an association implies that, in managing BMD loss, women might benefit from maintaining a high level of DHEAS.

Does obesity predict knee pain over fourteen years in women, independently of radiographic changes?

OBJECTIVE

To examine longitudinal patterns in body mass index (BMI) over 14 years and its association with knee pain in the Chingford Study.

METHODS

We studied a total of 594 women with BMI data from clinic visits at years (Y) 1, 5, 10, and 15. Knee pain at Y15 was assessed by questionnaire. Associations between BMI over 14 years and knee pain at Y15 were examined using logistic regression.

RESULTS

BMI significantly increased from Y1 to Y15 (P < 0.0005) with medians (interquartile ranges) of 24.5 kg/m(2) (22.5-27.2 kg/m(2) ) and 26.5 kg/m(2) (23.9-30.1 kg/m(2) ), respectively. At Y15, 45.1% of subjects had knee pain. A greater BMI at Y1 (odds ratio [OR] 1.34, 95% confidence interval [95% CI] 1.05-1.69), at Y15 (OR 1.34, 95% CI 1.10-1.61), and change in BMI over 15 years (OR 1.40, 95% CI 1.00-1.93) were significant predictors of knee pain at Y15 (P < 0.05). BMI change was associated with bilateral (OR 1.61, 95% CI 1.05-1.76, P = 0.024) but not unilateral knee pain (OR 1.22, 95% CI 0.73-1.76, P = 0.298). The association between BMI change and knee pain was independent of radiographic knee osteoarthritis (OA). The strength of association between BMI and knee pain at Y15 was similar during followup measurements.

CONCLUSION

Over 14 years, a higher BMI predicts knee pain at Y15 in women, independently of radiographic knee OA. When adjusted, the association was significant in bilateral, not unilateral, knee pain, suggesting alternative pathologic mechanisms may exist. The longitudinal effect of BMI on knee pain at Y15 is equally important at any time point, which may assist reducing the population burden of knee pain.

Large scale replication study of the association between HLA class II/BTNL2 variants and osteoarthritis of the knee in European-descent populations

Osteoarthritis (OA) is the most common form of arthritis and a major cause of disability. This study evaluates the association in Caucasian populations of two single nucleotide polymorphisms (SNPs) mapping to the Human Leukocyte Antigen (HLA) region and deriving from a genome wide association scan (GWAS) of knee OA in Japanese populations. The frequencies for rs10947262 were compared in 36,408 controls and 5,749 knee OA cases from European-descent populations. rs7775228 was tested in 32,823 controls and 1,837 knee OA cases of European descent. The risk (major) allele at rs10947262 in Caucasian samples was not significantly associated with an odds ratio (OR)  = 1.07 (95%CI 0.94 -1.21; p = 0.28). For rs7775228 the meta-analysis resulted in OR = 0.94 (95%CI 0.81-1.09; p = 0.42) for the allele associated with risk in the Japanese GWAS. In Japanese individuals these two SNPs are in strong linkage disequilibrium (LD) (r(2) = 0.86) with the HLA class II haplotype DRB1*1502 DQA1*0103 DQB1*0601 (frequency 8%). In Caucasian and Chinese samples, using imputed data, these SNPs appear not to be in LD with that haplotype (r(2)<0.07). The rs10947262 and rs7775228 variants are not associated with risk of knee OA in European descent populations and they do not appear tag the same HLA class II haplotype as they do in Japanese individuals.

The Ile585Val TRPV1 variant is involved in risk of painful knee osteoarthritis

Objective

To assess if a coding variant in the gene encoding transient receptor potential cation channel, subfamily V, member 1 (TRPV1) is associated with genetic risk of painful knee osteoarthritis (OA).

Methods

The Ile585Val TRPV1 variant encoded by rs8065080 was genotyped in 3270 cases of symptomatic knee OA, 1098 cases of asymptomatic knee OA and 3852 controls from seven cohorts from the UK, the USA and Australia. The genetic association between the low-pain genotype Ile–Ile and risk of symptomatic and asymptomatic knee OA was assessed.

Results

The TRPV1 585 Ile–Ile genotype, reported to be associated with lower thermal pain sensitivity, was associated with a lower risk of symptomatic knee OA in a comparison of symptomatic cases with healthy controls, with an odds ratio (OR) of 0.75 (95% CI 0.64 to 0.88; p=0.00039 by meta-analysis) after adjustment for age, sex and body mass index. No difference was seen between asymptomatic OA cases and controls (OR=1.02, 95% CI 0.82 to 1.27 p=0.86) but the Ile–Ile genotype was associated with lower risk of symptomatic versus asymptomatic knee OA adjusting for covariates and radiographic severity (OR=0.73, 95% CI 0.57 to 0.94 p=0.0136). TRPV1 expression in articular cartilage was increased by inflammatory cytokines (tumour necrosis factor α and interleukin 1). However, there were no differences in TRPV1 expression in healthy and arthritic synovial tissue.

Conclusions

A genotype involved in lower peripheral pain sensitivity is significantly associated with a decreased risk of painful knee OA. This indicates a role for the pro-nociceptive gene TRPV1 in genetic susceptibility to symptomatic knee OA, which may also be influenced by a role for this molecule in cartilage function.

Recommendations for standardization and phenotype definitions in genetic studies of osteoarthritis: the TREAT-OA consortium

OBJECTIVE

To address the need for standardization of osteoarthritis (OA) phenotypes by examining the effect of heterogeneity among symptomatic (SOA) and radiographic osteoarthritis (ROA) phenotypes.

METHODS

Descriptions of OA phenotypes of the 28 studies involved in the TREAT-OA consortium were collected. We investigated whether different OA definitions result in different association results by creating various hip OA definitions in one large population based cohort (the Rotterdam Study I (RSI)) and testing those for association with gender, age and body mass index using one-way ANOVA. For ROA, we standardized the hip-, knee- and hand ROA definitions and calculated prevalence's of ROA before and after standardization in nine cohort studies. This procedure could only be performed in cohort studies and standardization of SOA definitions was not feasible at this moment.

RESULTS

In this consortium, all studies with SOA phenotypes (knee, hip and hand) used a different definition and/or assessment of OA status. For knee-, hip- and hand ROA five, four and seven different definitions were used, respectively. Different hip ROA definitions do lead to different association results. For example, we showed in the RSI that hip OA defined as "at least definite joint space narrowing (JSN) and one definite osteophyte" was not associated with gender (P =0.22), but defined as "at least one definite osteophyte" was significantly associated with gender (P=3×10(-9)). Therefore, a standardization process was undertaken for ROA definitions. Before standardization a wide range of ROA prevalence's was observed in the nine cohorts studied. After standardization the range in prevalence of knee- and hip ROA was small.

CONCLUSION

Phenotype definitions influence the prevalence of OA and association with clinical variables. ROA phenotypes within the TREAT-OA consortium were standardized to reduce heterogeneity and improve power in future genetics studies.

Large-scale meta-analysis of interleukin-1 beta and interleukin-1 receptor antagonist polymorphisms on risk of radiographic hip and knee osteoarthritis and severity of knee osteoarthritis

OBJECTIVE

To clarify the role of common genetic variation in the Interleukin-1β (IL1B) and Interleukin-1R antagonist (IL1RN) genes on risk of knee and hip osteoarthritis (OA) and severity of knee OA by means of large-scale meta-analyses.

METHODS

We searched PubMed for articles assessing the role of IL1B and IL1RN polymorphisms/haplotypes on the risk of hip and/or knee OA. Novel data were included from eight unpublished studies. Meta-analyses were performed using fixed- and random-effects models with a total of 3595 hip OA and 5013 knee OA cases, and 6559 and 9132 controls respectively. The role of ILRN haplotypes on radiographic severity of knee OA was tested in 1918 cases with Kellgren-Lawrence (K/L) 1 or 2 compared to 199 cases with K/L 3 or 4.

RESULTS

The meta-analysis of six published studies retrieved from the literature search and eight unpublished studies showed no evidence of association between common genetic variation in the IL1B or IL1RN genes and risk of hip OA or knee OA (P>0.05 for rs16944, rs1143634, rs419598 and haplotype C-G-C (rs1143634, rs16944 and rs419598) previously implicated in risk of hip OA). The C-T-A haplotype formed by rs419598, rs315952 and rs9005, previously implicated in radiographic severity of knee OA, was associated with reduced severity of knee OA (odds ratio (OR)=0.71 95%CI 0.56-0.91; P=0.006, I(2)=74%), and achieved borderline statistical significance in a random-effects model (OR=0.61 95%CI 0.35-1.06 P=0.08).

CONCLUSION

Common genetic variation in the Interleukin-1 region is not associated with prevalence of hip or knee OA but our data suggest that IL1RN might have a role in severity of knee OA.

Selenium concentration and speciation in biofortified flour and bread: Retention of selenium during grain biofortification, processing and production of Se-enriched food

The retention and speciation of selenium in flour and bread was determined following experimental applications of selenium fertilisers to a high-yielding UK wheat crop. Flour and bread were produced using standard commercial practices. Total selenium was measured using inductively coupled plasma-mass spectrometry (ICP-MS) and the profile of selenium species in the flour and bread were determined using high performance liquid chromatography (HPLC) ICP-MS. The selenium concentration of flour ranged from 30ng/g in white flour and 35ng/g in wholemeal flour from untreated plots up to >1800ng/g in white and >2200ng/g in wholemeal flour processed from grain treated with selenium (as selenate) at the highest application rate of 100g/ha. The relationship between the amount of selenium applied to the crop and the amount of selenium in flour and bread was approximately linear, indicating minimal loss of Se during grain processing and bread production. On average, application of selenium at 10g/ha increased total selenium in white and wholemeal bread by 155 and 185ng/g, respectively, equivalent to 6.4 and 7.1μg selenium per average slice of white and wholemeal bread, respectively. Selenomethionine accounted for 65-87% of total extractable selenium species in Se-enriched flour and bread; selenocysteine, Se-methylselenocysteine selenite and selenate were also detected. Controlled agronomic biofortification of wheat crops for flour and bread production could provide an appropriate strategy to increase the intake of bioavailable selenium.

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⁻⁹), 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.

Genetic variation in the SMAD3 gene is associated with hip and knee osteoarthritis

OBJECTIVE

Smad3 (or, MADH3) is a key intracellular messenger in the transforming growth factor beta signaling pathway. In mice, Smad3 deficiency accelerates growth plate chondrocyte maturation and leads to an osteoarthritis (OA)-like disease. We undertook this study to investigate the role of genetic variation in SMAD3 in the risk of large-joint OA in humans.

METHODS

Ten tag single-nucleotide polymorphisms (SNPs) in the SMAD3 gene region were tested in a discovery set: 313 patients who had undergone total knee replacement, 214 patients who had undergone total hip replacement, and 520 controls from the UK. The SNP associated with both hip and knee OA was subsequently genotyped in 1,221 controls and 1,074 cases from 2 cohorts of patients with hip OA and 2,537 controls and 1,575 cases from 4 cohorts of patients with knee OA.

RESULTS

A SNP (rs12901499) mapping to intron 1 of SMAD3 was associated with both knee and hip OA (P < 0.0022 and P < 0.021, respectively) in the discovery set. In all study cohorts, the major allele (G) was increased among OA patients relative to controls. A meta-analysis for knee OA yielded an odds ratio (OR) of 1.22 (95% confidence interval [95% CI] 1.12-1.34), P < 7.5 x 10(-6). For hip OA, the OR was 1.22 (95% CI 1.09-1.36), P < 4.0 x 10(-4). No evidence for heterogeneity was found (I(2) = 0%).

CONCLUSION

Our data indicate that genetic variation in the SMAD3 gene is involved in the risk of both hip OA and knee OA in European populations, confirming the results from animal models on the potential importance of this molecule in the pathogenesis of OA.

Free-radical carbon-carbon bond formation in organic synthesis

Organic chemists have begun to use intra- and intermolecular free-radical addition reactions to develop useful synthetic transformations. Carboncentered radicals can form bonds with electron-rich or electron-deficient alkenes, allenes, and acetylenes. Radical addition reactions can also be used to construct hindered carbon-carbon bonds. These characteristics, as well as the large number of functional groups that tolerate free-radical conditions contribute to the importance of such reactions in organic synthesis.

Serum branched-chain amino acid to histidine ratio: a novel metabolomic biomarker of knee osteoarthritis

OBJECTIVE

There is a pressing need to develop reliable molecular biomarkers in osteoarthritis. The aim of the study was to identify novel serum biomarkers for osteoarthritis using a metabolomics approach.

METHODS

A two-stage study design was utilised. 123 knee osteoarthritis cases and 299 controls were selected from the TwinsUK cohort as a discovery sample. 76 knee osteoarthritis cases and 100 controls from the Chingford Study were used as replication. Knee osteoarthritis was defined as either radiographic, medically diagnosed or total knee joint replacement due to primary osteoarthritis. All the subjects were unrelated white women. Their serum samples were assessed for targeted metabolite profiling by electrospray ionisation tandem mass spectrometry using the AbsoluteIDQ kit. 163 serum metabolites were assessed and their concentrations obtained. The ratios of metabolite concentrations as proxies for enzymatic reaction rates were calculated and tested for the association with knee osteoarthritis. Significance was assessed after adjustment for multiple testing (Bonferroni method) and potential confounders.

RESULTS

In the discovery stage, the authors identified 14 ratios significantly associated with knee osteoarthritis with p < or = 1.9 x 10(-6). Two of these 14 ratios were successfully confirmed in the replication stage-the ratios of valine to histidine and xleucine to histidine, with p=0.002. The significance remained after adjustment for age and body mass index.

CONCLUSION

This is the first serum-based metabolomic study of osteoarthritis in humans. The branched-chain amino acids to histidine ratio has potential clinical use as an osteoarthritis biomarker and shows the clinical potential of metabolomics.

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.

Interleukin-6 is a significant predictor of radiographic knee osteoarthritis: The Chingford Study

OBJECTIVE

There is a great need for identification of biomarkers that could improve the prediction of early osteoarthritis (OA). We undertook this study to determine whether circulating levels of interleukin-6 (IL-6), tumor necrosis factor alpha (TNFalpha), and C-reactive protein (CRP) can serve as useful markers of radiographic knee OA (RKOA) in a normal human population.

METHODS

RKOA data were obtained from the cohort of the Chingford Study, a prospective population-based study of healthy, middle-aged British women. The RKOA-affected status of the subjects was assessed using the Kellgren/Lawrence (K/L) grade as determined on radiographs obtained at baseline (n = 908) and at 10 years and 15 years thereafter. Serum levels of CRP, IL-6, and TNFalpha were assayed at 5, 8, and 15 years, using high-sensitivity commercial assays. A K/L grade of >or=2 in either knee was used as the outcome measure. Statistical analyses included analysis of variance for repeated measurements and logistic regression models, together with longitudinal modeling of dichotomous responses.

RESULTS

During 15 years of followup, the prevalence of RKOA (K/L grade >or=2) increased from 14.7% to 48.7% (P < 0.00001 versus baseline). The body mass index (BMI) and circulating levels of CRP and IL-6 were consistently and significantly higher in subjects diagnosed as having RKOA. When multiple logistic regression was applied to the data, the variables of older age (P = 3.93 x 10(-5)), higher BMI at baseline (P = 0.0003), and increased levels of IL-6 at year 5 (P = 0.0129) were determined to be independent predictors of the appearance of RKOA at year 10. The results were fully confirmed using longitudinal modeling of repeated measurements of the data obtained at 3 visits. The odds ratio for RKOA in subjects whose IL-6 levels were in the fourth quartile of increasing levels (versus the first quartile) was 2.74 (95% confidence interval 1.94-3.87).

CONCLUSION

This followup study showed that individuals were more likely to be diagnosed as having RKOA if they had a higher BMI and increased circulating levels of IL-6. These results should stimulate more work on IL-6 as a potential 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.

Further genetic evidence suggesting a role for the RhoGTPase-RhoGEF pathway in osteoporosis

Osteoporosis is a highly heritable trait that appears to be influenced by multiple genes. Genome-wide linkage studies have highlighted the chromosomal region 3p14-p21 as a quantitative trait locus for BMD. We have previously published evidence suggesting that the ARHGEF3 gene from this region is associated with BMD in women. The product of this gene activates the RHOA GTPase, the gene for which is also located within this region. The aim of this study was to evaluate the influence of genetic polymorphism in RHOA on bone density in women. Sequence variation within the RHOA gene region was determined using 9 single nucleotide polymorphisms (SNPs) in a discovery cohort of 769 female sibs. Of the 9 SNPs, one was found to be monomorphic with the others representing 3 distinct linkage disequilibrium (LD) blocks. Using FBAT software, significant associations were found between two of these LD blocks and BMD Z-score of the spine and hip (P=0.001-0.036). The LD block tagged by the SNP rs17595772 showed maximal association, with the more common G allele at rs17595772 associated with decreased BMD Z-score. Genotyping for rs17595772 in a replication cohort of 780 postmenopausal women confirmed an association with BMD Z-score (P=0.002-0.036). Again, the G allele was found to be associated with a reduced hip and spine BMD Z-score. These results support the implication of the RhoGTPase-RhoGEF pathway in osteoporosis, and suggest that one or more genes in this pathway may be responsible for the linkage observed between 3p14-p21 and BMD.

Meta-analysis of genome-wide scans for human adult stature identifies novel Loci and associations with measures of skeletal frame size

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value = 6.1×10⁻⁸ and rs910316 in TMED10, P-value = 1.4×10⁻⁷) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value = 3×10⁻⁷ and rs849141 in JAZF1, P-value = 3.2×10⁻¹¹). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value = 4×10⁻⁵ and rs6817306 in LCORL, P-value = 4×10⁻⁴), hip axis length (including rs6830062 at LCORL, P-value = 4.8×10⁻⁴ and rs4911494 at UQCC, P-value = 1.9×10⁻⁴), and femur length (including rs710841 at PRKG2, P-value = 2.4×10⁻⁵ and rs10946808 at HIST1H1D, P-value = 6.4×10⁻⁶). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.

The first genetic association studies of adult height have confirmed a role of many common variants in influencing human height, but to date, the genetic basis of differences between different skeletal components of height have not been addressed. Here, we take advantage of recent technical and methodological advances to examine the role of common genetic variants on both height and skeletal components of height. By examining nearly 20,000 individuals from the UK and the Netherlands, we provide statistically significant evidence that 17 genomic regions are associated with height, including four novel regions. We also examine, for the first time, the association of these 17 regions with skeletal size measurements of spine, femur, and hip axis length, a measurement of hip geometry known to influence the risk of osteoporotic fractures. We find that some height loci are also associated with these skeletal components, although further statistical tests will be required to verify if these genetic variants act differentially on the individual skeletal measurements. The knowledge generated by this and other studies will not only inform the genetics of human quantitative variation, but will also lead to the potential discovery of many medically important polymorphisms.

Effective measurement of knee alignment using AP knee radiographs

The gold standard for measuring knee alignment is mechanical axis determined using full-limb radiographs (FLR). Measurement of joint alignment using antero-posterior (AP) knee radiographs is more accessible, economical and involves less radiation exposure to the patient compared with using full-limb radiographs. The aim of this study was to compare and assess the reproducibility of knee joint axial alignment on full-limb radiographs and conventional AP knee radiographs. Knee alignment was measured in 40 subjects (80 knees) from the TwinsUK registry. Measurement of mechanical knee alignment was from FLR, and anatomic knee alignment from weight-bearing AP knee radiographs. Reproducibility was assessed by intra-class correlation coefficients and kappa statistics. Reproducibility of knee alignment for both methods was good, with intra-observer ICC's of 0.99 for both FLR and AP radiographs. The mean alignment angle on FLR was 178.9 degrees (SD 2.1, range 173-183 degrees ), and 179.0 degrees (SD 2.1, range 173-185 degrees ) on AP films. 58.8% of knees on FLR and 66.3% on AP films were of varus alignment. Good correlations were seen between results for FLR and AP radiographs, with ICC ranging from 0.87-0.92 for left and right knees, and kappa statistics of 0.65-0.74. Standard AP knee radiographs can be used to measure knee alignment with good reproducibility, and provide comparable results to those obtained from FLR. This will facilitate measurement of knee alignment in existing cohort studies to assess malalignment as a risk factor of incident OA, and in clinical practice.

Association of the DVWA and GDF5 polymorphisms with osteoarthritis in UK populations

BACKGROUND

Variants in the growth differentiation factor 5 (GDF5) and in the double von Willebrand factor A (DVWA) have recently been reported to be associated with osteoarthritis (OA) in Asian populations.

OBJECTIVE

To assess the role of such variants in OA susceptibility in two independent UK samples of Caucasian origin.

METHODS

Polymorphisms rs11718863 and rs7639618 (DVWA) and rs143383 (GDF5) were genotyped in 999 patients with knee OA, 843 patients with hip OA and 1166 controls from two UK studies from Nottingham and Chingford.

RESULTS

In agreement with previous reports, the major allele at rs143383 (GDF5) was associated with a higher risk of knee OA in our samples (OR(MH) = 1.29, 95% CI 1.14 to 1.47 p = 8 x 10(-5)). Conversely, the major allele at the DVWA SNP rs7639618, which increased risk in Asians, was not associated with a risk of knee OA, (OR(MH) = 0.88, 95% CI 0.74 to 1.04; p = 0.12). A meta-analysis of Asian and UK knee OA data indicated highly significant heterogeneity (I(2) = 92%, Q = 48.5, p = 7 x 10(-10)) and no significant association with knee OA using a random effects meta-analysis (OR(DL) = 1.18, 95% CI 0.86 to 1.63; p = 0.309).

CONCLUSIONS

These data confirm that the GDF5 variant is consistently associated with the risk of knee OA. Considerable ethnic variation in allele frequencies at the DVWA gene was found and no significant association was found in UK samples or by combining UK and Asian samples. The results suggest that the effect of DVWA amino acid changes on tubulin binding is unlikely to influence the risk of OA in Caucasians.