Article Text
Abstract
Objectives: To determine whether there is a genetic contribution to knee osteoarthritis (OA) as defined by structural change on plain radiographs.
Design: Sibling study. Comparison of knee OA prevalence in the community with knee OA prevalence among siblings of index cases with knee OA.
Subjects: 490 knee OA index cases listed for total knee replacement and 737 of their siblings aged >40 years, 1729 community subjects aged >40 years.
Main outcome measure: Odds ratio for radiographic knee OA in siblings v community participants.
Results: The age, sex, and knee pain adjusted odds ratios in siblings were 2.9 (95% confidence interval 2.3 to 3.7) for tibiofemoral OA and 1.7 (1.4 to 2.2) for patellofemoral OA. This difference remained after adjustment for important environmental risk factors. The heritability estimate for knee OA was 0.62.
Conclusions: Siblings are at increased risk of knee OA in comparison with the general population. This is likely to be due to genetic factors.
- BMI, body mass index
- CI, confidence interval
- OA, osteoarthritis
- OR, odds ratio
- PFJ, patellofemoral joint
- PF OA, patellofemoral osteoarthritis
- RR, risk ratio
- TFJ, tibiofemoral joint
- TF OA, tibiofemoral osteoarthritis
- TKR, total knee replacement
- heritability
- osteoarthritis
- knees
- sibling study
Statistics from Altmetric.com
- BMI, body mass index
- CI, confidence interval
- OA, osteoarthritis
- OR, odds ratio
- PFJ, patellofemoral joint
- PF OA, patellofemoral osteoarthritis
- RR, risk ratio
- TFJ, tibiofemoral joint
- TF OA, tibiofemoral osteoarthritis
- TKR, total knee replacement
Osteoarthritis (OA) is a common age associated multifactorial condition.1 The knee is the most clinically significant site affected, with an estimated 10% of the population aged >55 experiencing pain and severe disability.2 Few studies have examined in detail the genetic contribution to knee OA. A classic twin study undertaken in the UK has estimated the heritability of knee OA—that is, the proportion of the variance in occurrence of knee OA that may be explained by genetic factors, to be 39%.3 Further analyses have shown that the association between knee OA and high body mass index, which has a heritable component, is unlikely to be mediated by shared genetic factors.4 However, this is the only twin study of knee OA and the study included only women. Another strategy was to examine the risks for total knee replacement (TKR) and total hip replacement in siblings of subjects who had undergone these procedures.5 The relative risks of joint replacement in this UK study were 4.8 at the knee and 1.9 at the hip. However, the confidence intervals (CIs) were wide, these data were not adjusted for important risk factors, and spouses were used as the control population. Two community studies in North America and the Netherlands have investigated familial aggregation and segregation of OA but have found no evidence for a genetic contribution at the knee.6,7 However, the first of these studies was limited by the absence of patellofemoral joint (PFJ) views and the young age of participants, and the second study may have included too few subjects to show an effect. We therefore undertook a sibling study to better estimate the genetic contribution to symptomatic and clinically significant knee OA. Importantly, our design included assessment of all three knee compartments, community subjects as controls, and adjustment for recognised constitutional and environmental risk factors for knee OA.
METHODS
Ethical approval
Approval was obtained from the North Nottinghamshire Research Ethics Committee and the Nottingham City Hospital Research Ethics Committee.
Ascertainment of index cases of knee OA
Index cases were 2857 white patients who had undergone or were awaiting TKR for OA at three Nottinghamshire orthopaedic centres between 1997 and 2000. These people were ascertained from the Trent regional arthroplasty list, an electronic database. All patients undergoing this procedure within the NHS (National Health Service) in the Trent region should be “captured” on this database. Those who had been, or would be, <40 years old at the time of operation were excluded. Those with a self reported diagnosis of inflammatory arthritis, trauma, or joint sepsis were excluded at this stage. Diagnosis was validated by reviewing the preoperative radiographs.
All local index subjects were sent a questionnaire asking for the contact details of all of their siblings, regardless of whether they were thought to have arthritis. Non-responders were re-mailed after 4 weeks.
Ascertainment of siblings
Siblings were identified from the questionnaire responses of the index subjects; half-siblings and those aged <40 were excluded.
To reduce oversampling those with OA, only eligible siblings living locally were asked to attend for metrologist interview and examination and to undergo knee radiographs. It was emphasised that all participants would be welcomed, regardless of whether or not they had arthritis. Those who did not reply were re-mailed after 4 weeks, and the ages of non-responders were obtained by telephone contact. Consent to examine any previous knee radiographs was sought from non-local siblings.
Ascertainment of community participants
Data and radiographs from Nottingham community based knee pain studies performed at our centre between 1995 and 1999 were used. Random samples aged >40 from general practice lists in Nottingham had been surveyed using postal questionnaires.8,9 Subjects were defined as having knee pain if they responded yes to “Have you ever had pain in or around the knee on most days for at least a month?” and “If so, have you experienced any pain during the last year?” In nested case-control studies those with current knee pain were matched for age and sex with those without, and these groups were then asked to undergo knee radiographs. In a nested intervention study only those with knee pain underwent radiographs and clinical examination. Exclusion criteria included TKR (although those awaiting this procedure were included), terminal illness, and incapacity to give informed consent.
Data collection: knee OA risk factors
Postal questionnaires mailed to the community participants and to the siblings were used to collect data on age, sex, presence or absence of knee pain, knee injuries or operations, and smoking history. Metrologists measured body mass index (BMI) in kg/m2 and examined for Heberden’s or Bouchard’s nodes. Reproducibility for the assessment of nodes was assessed by three metrologists re-examining a group of 10 patients, showing a spectrum of hand changes, in random order on the same day.
Radiographic technique: siblings
The majority (556) of sibling radiographs were performed as for the community sample.
An additional currently available 92 sibling radiographs were paired schuss10 and skyline views. The remaining non-local siblings had non-standardised radiographs, as did local siblings who had undergone knee radiographs in the preceding 6 months.
Radiographic technique: community sample
All radiographs were taken under the same standardised conditions, and were weightbearing, fully extended, anteroposterior views (55 kV, 8 mA/s, full scale deflection 100 cm, and skyline 30° flexion views (60 kV, 10 mA/s, full scale deflection 100 cm) after the method of Laurin et al11 using standard Fuji film.
Radiographic assessment: grading
A single observer (RLN) read all the radiographs. Community radiographs and the 556 paired sibling radiographs that had been performed in the same way were mixed and read “blind”. The remaining sibling radiographs could not be “blinded”, and the most recent radiographs for both knees were assessed. If a TKR was evident the film just before surgery was selected.
Radiographs of the tibiofemoral joint (TFJ) and the PFJ (lateral views) were assigned Kellgren and Lawrence grades.12 Individual features of OA at both the TFJ and the PFJ (skyline views) were assessed using an adaptation of a validated line drawing atlas.13 Joint space narrowing in all compartments was graded −1 to 5. Osteophyte was graded 0 to 5 at all sites (medial and lateral femoral trochlear, medial and lateral patella, medial and lateral femoral condyle, and medial and lateral tibia).
Radiographic assessment: reproducibility
Reproducibility was examined half way through the study by re-reading 40 community radiographs 1 month after the first assessment. Radiographs were chosen by an independent researcher to represent a broad range of OA severity.
Case definition
OA was defined on the basis of structural change. The definitions of OA were:
Definite tibiofemoral OA (TF OA)—that is, Kellgren and Lawrence grade ⩾2 at the TFJ
Moderate-severe TF OA—that is, Kellgren and Lawrence grade 3–4 at the TFJ
Definite patellofemoral OA (PF OA)—that is, for skyline radiographs: line atlas osteophyte score ⩾2 plus joint space narrowing score ⩾2; for lateral radiographs: Kellgren and Lawrence grade ⩾2
Moderate-severe PF OA—that is, for skyline radiographs: line atlas osteophyte score ⩾3 plus joint space narrowing score ⩾3; for lateral radiographs: Kellgren and Lawrence grade 3–4.
Data analysis
The TFJ radiograph with the highest Kellgren and Lawrence grade and the right knee PFJ film were analysed for each subject. Reproducibility for OA grading was assessed by the κ statistic14 using STATA for windows.15 Direct standardisation of OA prevalence was by age, sex, and knee pain status using data from a Nottingham survey.8 Odds ratios (ORs) from logistic regression using SPSS were used to examine risk factors for OA.16 Poisson regression analysis was with STATA.15 Heritability was calculated by comparing the liability to OA among siblings with that among the community17 using age standardised prevalence estimates.18 The assumptions were that liability to OA was distributed normally in both groups, and the variances were equal.
RESULTS
Response rate: index cases
Questionnaires were mailed to 2857 index cases. After the second mailing 2016 replies were received, a response rate of 70.6% (fig 1⇓). From these the following were excluded: 133 who did not want to participate, 806 with no living siblings, 32 who did not supply siblings’ contact details, 7 with half-siblings, 30 who had died, and 27 with alternative diagnoses. The final pool of index cases consisted of 981 valid respondents, 490 of whom had valid siblings who participated in the study. Of the later group 460 preoperative radiographs of the index cases confirmed OA, and 30 could not be retrieved but their siblings were included.
Recruitment of index cases and siblings.
Response rate: siblings
Invitations were mailed to 1672 siblings, including 348 who were not local. These 1672 people were derived from the 981 index cases. 1214 replies were received, a response rate of 72.6% (fig 1⇑). From these the following were excluded: 129 declined participation, 12 were not siblings, 19 had died, 7 were too unwell, 5 were aged <40, and 34 had index cases who did not have primary OA. The remaining 1008 valid respondents were invited to undergo knee radiographs. Of these, 787 (78.1%) were subsequently radiographed, 8 were excluded owing to inflammatory arthritis, and radiographs of 737 siblings were retrieved.
Response rate: community sample
The overall response rate to the two postal questionnaires was 70.1% after the second mailing. For participants of the nested case-control studies asked to undergo radiographs, uptake was 92% among participants with knee pain and 77% among participants without knee pain.
Study population characteristics
The 490 index cases comprised 276 women (56.3%, mean age 70.5 years) and 214 men (43.7%, mean age 68.7 years). All the siblings were white, as were 99% of the community group. The 737 siblings were older than the 1729 community participants (table 1⇓) and comprised a higher proportion of men (42.6% v 37.2%, p = 0.012). Knee pain was common among both groups. The mean BMI of male siblings did not significantly differ from that of male community participants. Female siblings were heavier than female community participants, although this difference was small (mean difference 0.8 kg/m2, 95% CI 0.3 to 1.4 kg/m2). Similar proportions of community participants and siblings had never smoked, but a higher proportion of siblings were current smokers. The presence of nodes and history of meniscectomy were more common among siblings.
Participant characteristics
The age-sex distribution of 235 (53%) of the non-responder siblings was known, and this group contained proportionally more men (48.1%). Non-responders were similar in age to participants (women: mean 68.6 years; men: mean 64.9 years).
Reproducibilty of OA assessments
Reproducibility for OA assessments was good. The κ coefficient was 0.84 for Kellgren and Lawrence grades, and ranged from 0.69 to 0.90 for line atlas features.
Reproducibility for assessment of nodes
Intraobserver reproducibility was good, with κ coefficient ⩾0.75 at all sites; interobserver reproducibility was also good, with κ coefficient ⩾0.6 at all sites.
Knee OA prevalence
Knee OA was strongly age associated (figs 2 and 3⇓⇓). A χ2 test for trend for the association between OA and age was highly significant for the siblings and the community sample (p<0.0001 for both samples). TF OA was more common in men, whereas PF OA was more common among women. The prevalence of knee OA was higher in siblings than among the community for all definitions, excepting moderate-severe PF OA in men, for which prevalences were similar.
Tibiofemoral osteoarthritis prevalence in siblings and in the community.
Patellofemoral osteoarthritis prevalence in siblings and in the community.
The age, sex, and knee pain standardised prevalence estimates confirm the excess of knee OA among siblings (table 2⇓).
Age, sex, and knee pain standardised prevalence estimates
Effect of other key risk factors
Adjusting for important knee OA risk factors did not alter the finding of an increased knee OA risk in siblings (table 3⇓). For self reported knee injury and operations, only meniscectomy was associated with OA and hence was used in these calculations.
Odds ratios of OA in siblings
Sensitivity analysis and heritability estimate
Selection of one sibling from each family for analysis did not alter the increased chance of knee OA among siblings, with age, sex, and knee pain adjusted OR = 2.8 (95% CI 2.0 to 3.9).
Exclusion of the 181 siblings with “unblinded” radiographs had no influence on the main finding, with adjusted OR 2.3 (95% CI 1.9 to 3.0). The excess of OA among siblings was maintained when the 23 siblings who had undergone TKR were omitted from analysis, giving adjusted OR 2.4 (95% CI 1.9 to 3.0).
Risk ratios (RRs) were calculated by Poisson regression using age as a surrogate for “genetic” exposure duration. The adjusted RRs were similar to the ORs with RR = 1.6 (95% CI 1.4 to 2.0).
The heritability (h2) of knee OA was estimated at 0.62 overall (95% CI 0.49 to 0.75).
DISCUSSION
This study shows that siblings had double the risk of knee OA compared with community people. This is despite the fact that the community group was not a true random sample, as people with knee pain were overrepresented, and therefore the prevalence of OA among this group is probably an overestimate. The heritability of knee OA, with 62% of disease variance likely to be genetically determined, is greater than previous estimates.3,5 Interestingly, we found that heritability was higher in men (0.78) than in women (0.49). This may explain why our overall estimate is higher than that reported in a previous twin study that was restricted to women (heritability 0.39).3 Of further interest is a study of similar design to ours that reported an increased risk of hip OA in male siblings compared with female siblings.19 The heritability (0.31) reported in the Oxford study was based on the phenotype of either hip or knee replacement rather than the more stringent criterion of radiographic knee OA as in our study.5 It is likely that this difference in phenotype is in part responsible for the difference in heritability estimates. Furthermore, spouses acted as controls in the Oxford study. It is recognised that assortative mating may occur, leading to similarities for disease between spouses, and these may have reduced the heritability estimates in that study.20
The significantly increased risk of OA amongst siblings remained after adjustment for other key risk factors known to influence knee OA: age, sex, BMI, smoking, presence of Heberden’s or Bouchard’s nodes, and meniscectomy. The genetic risk was present in both men and women across a broad age span. It has been suggested that risk factors for OA at the TFJ and PFJ may differ from one another,21 but sibling risk of OA at the TFJ was similar to that at the PFJ. Similarly, there were no real differences in genetic risk at both joints according to radiographic severity of OA. However, the smaller numbers of people with more severe disease may have precluded the detection of any such difference. Sensitivity analyses showed that the increased risk in siblings was robust as it was unaffected by inclusion of only one sibling from each family, exclusion of siblings with TKR, or by statistical methodology.
This observed genetic risk may be generalised to the wider community as efforts were made to ensure that representative samples were studied. Exclusion of young index cases, together with review of preoperative radiographs, ensured that people with premature forms of OA or alternative diagnoses were omitted. Recruitment was restricted to siblings who lived locally. The purpose of this was to reduce the chance of oversampling siblings with knee OA, as it is conceivable that those with knee symptoms and knee OA would be more likely to travel to take part in the study. In addition, emphasis was also placed on recruiting all local siblings whether or not they had a knee problem, with the aim of reducing the chances of people with knee OA being preferentially recruited. The response rates were good, reducing the potential for non-responder bias.
There are several caveats to this study. Firstly, we calculated heritability in order to enable comparison with other studies, but it should be interpreted with care as we were able to adjust only for age. Secondly, although an increased risk of OA was observed in siblings after adjustment for important risk factors, aspects of shared family environment such as sporting and leisure activities or diet, which were not assessed, may contribute to the observed apparent genetic risk. However, associations between OA and these factors are not well established, and would be unlikely to completely explain the increased prevalence of knee OA among siblings. A third caveat is that the twofold increased risk of knee OA relates to siblings of patients with severe disease (who needed TKR) and it is unknown whether the risk may be less marked in families with milder disease. Finally, this study included only white people and there is a need to investigate the role of genetics in OA among other ethnic groups.
This is the first large sibling study of knee OA and adds to the growing body of evidence that there is a genetic contribution to common forms of OA. Identification of the responsible genes will enhance our understanding of OA pathogenesis.
Acknowledgments
We are indebted to all the participants in these studies. We acknowledge the assistance of Alison Carr with statistical advice, and Margaret Ball, Sheila O’Reilly, Kim Thomas, Margaret Wheeler, and Elizabeth Williams with data retrieval.
We are grateful to GlaxoSmithKline, USA for financial support (RLN) and the Arthritis Research Campaign for infrastructure support (ICAC grant D0541).