The genetic contribution to severe post-traumatic osteoarthritis
- Ana M Valdes1,2,
- Sally A Doherty1,
- Kenneth R Muir3,4,
- Margaret Wheeler1,
- Rose A Maciewicz5,
- Weiya Zhang1,
- Michael Doherty1
- 1Academic Rheumatology, Nottingham City Hospital, Nottingham, UK
- 2Department of Twin Research & Genetic Epidemiology, King's College London, St Thomas’ Hospital, London, UK
- 3Health Sciences Research Institute, Warwick Medical School, University of Warwick, Coventry, UK
- 4Institute of Population Health, University of Manchester, Manchester, UK
- 5Respiratory & Inflammation Research Area, AstraZeneca, Charnwood R&D, Loughborough, Leicestershire, UK
- Correspondence to Dr Ana M Valdes, Department of Twin Research & Genetic Epidemiology, Kings College London, St Thomas Hospital campus, London SE1 7EH, UK;
- Received 24 August 2012
- Revised 6 December 2012
- Accepted 29 December 2012
- Published Online First 26 January 2013
Objective To compare the combined role of genetic variants loci associated with risk of knee or hip osteoarthritis (OA) in post-traumatic (PT) and non-traumatic (NT) cases of clinically severe OA leading to total joint replacement.
Methods A total of 1590 controls, 2168 total knee replacement (TKR) cases (33.2% PT) and 1567 total hip replacement (THR) cases (8.7% PT) from 2 UK cohorts were genotyped for 12 variants previously reported to be reproducibly associated with risk of knee or hip OA. A genetic risk score was generated and the association with PT and NT TKR and THR was assessed adjusting for covariates.
Results For THR, each additional genetic risk variant conferred lower risk among PT cases (OR=1.07, 95% CI 0.96 to 1.19; p=0.24) than NT cases (OR 1.11, 95% CI 1.06 to 1.17; p=1.55×10−5). In contrast, for TKR, each risk variant conferred slightly higher risk among PT cases (OR 1.12, 95% CI 1.07 to 1.19; p=1.82×10−5) than among NT cases (OR 1.08, 95% CI 1.03 to 1.1; p=0.00063).
Conclusions Based on the variants reported to date PT TKR cases have at least as high a genetic contribution as NT cases.
OA is divided into ‘primary’ OA, considered to result mainly from genetic and constitutional factors, and ‘secondary’ OA when there is an obvious cause of joint insult.3 A long-recognised risk factor for OA is direct joint injury, such as a subchondral fracture, meniscectomy, anterior cruciate tear or a deforming tibial or femoral fracture that results in poor alignment or leg length shortening.4 ,5 An estimated 12% of symptomatic OA may be attributed to a post-traumatic (PT) cause6 and differences in radiographic patterns have been reported between primary and PT knee OA cases.7 However, two studies that examined the relationship between hand OA and the risk of developing knee OA following meniscectomy found that the presence of hand OA associated with a higher frequency and radiographic severity of post-meniscectomy OA,8 ,9 demonstrating that development of OA following a meniscal tear and subsequent surgery are not necessarily of secondary origin.
Genetic factors are known to influence risk of hip OA, knee OA and generalised OA (see Valdes and Spector10). Recent genetic association studies have identified a number of genetic variants associated with knee or hip OA or with severe large joint OA leading to total joint arthroplasty.11–15 To date, 12 independent variants have been reported to be associated with risk of hip or knee OA or related traits with a statistical significance level of p<1×10−7 with (table 1).
The aim of this study was to examine whether PT OA might have a lower genetic contribution than non-traumatic (NT) idiopathic large joint OA. We investigated the combined role of the published OA-risk genetic variants in two case control studies from the UK focusing on severe OA leading to total knee replacement (TKR) and total hip replacement (THR).
Patients and methods
THR and TKR cases were recruited from two case-control studies in Nottingham, UK: the Nottingham OA Case-Control and the Genetics of Osteoarthritis and Lifestyle (GOAL) study. Controls were obtained from both these studies.
PT OA definition
PT THR was defined as THR in the presence of hip injury and PT TKR as TKR in the presence of knee injury as described in the online supplementary methods section.
Genomic DNA was extracted from peripheral blood leukocytes of affected individuals and controls using standard protocols. Genotyping was carried out by Kbioscience Ltd, Hertfordshire, UK. Single nucleotide polymorphisms (SNPs) were genotyped using the Kompetitive Allele Specific PCR (KASPar) chemistry.
The sum of risk alleles (which can take values 0, 1 or 2 for each variant) as defined in table 1 over all variants was computed for each study participant. Individuals with more than 1 missing genotype for all 12 variants were excluded, otherwise the genotype was imputed using the allele frequency in the population.
Logistic regressions using TKR or THR as outcomes and including age, gender, (BMI) and knee genetic risk score (for TKR) or hip genetic risk score (for THR) were performed. Analyses were carried out including all TKR (or THR) cases, only PT cases (PT-TKR or PT-THR) and only NT cases (NT-TKR or NT-THR). Results are reported as OR and corresponding 95% CI for the GOAL and Nottingham cohorts separately. ORs were also meta-analysed for both cohorts using a fixed effects analysis. Analyses were performed using R V.2.13.1 (http://www.r-project.org/).
The list of variants investigated in this study is shown in table 1 along with the best-reported p value in the literature. The role of injury as a risk factor was assessed in the GOAL study where history of trauma was recorded for cases and controls. After adjustment for age, sex and BMI the risk estimates found were OR=3.27, (95 % CI 2.51 to 4.25; p<1.1×10(−18)) for TKR and OR=5.08 (95 % CI 3.05 to 8.47; p<4.2×10(−10)) for THR. The combined genetic and descriptive statistics risks in each cohort are shown in table 2. These factors were analysed in each cohort to assess their role in defining risk of THR and TKR following trauma or in the absence of injury.
The mean of each of these risk factors stratified by injury status shows, as expected, a significantly younger age and a higher proportion of men in the post-trauma category (see online supplementary table S1).
The knee genetic risk score for TKR and the hip genetic risk score for THR are both significantly related to risk of total joint replacement (table 3). The average OR contributed by each additional variant carried is similar to those reported in the studies that identified these genetic variants (1.11 for hip OA risk variants and 1.09 for knee OA risk). Age and sex contribute differently to risk of PT-TKR and NT-TKR (see online supplementary table S1). We did not observe a significant difference in the risk conferred by knee genetic risk factors, but the OR is slightly higher for PT-TKR cases (OR=1.12) than for non-traumatic cases (OR=1.08). For THR only the gender composition is significantly different between NT and PT cases (see online supplementary table S1). The genetic risk conferred does not achieve statistical significance among PT-THR cases and is lower (OR=1.07) than among NT-THR cases (OR=1.11) where it is highly significant. Results were also computed excluding individuals with missing genotypes (see online supplementary table S2) and are very similar to those using imputed genotypes for missing data.
In this study we report for the first time the role of genetic factors in PT clinically severe OA requiring total joint replacement and compare it to the role that these same genetic factors play among ‘primary’ cases. Our hypothesis was that primary cases would carry a larger number of genetic risk variants than PT OA cases. However, the OR estimates for hip and knee risk scores are not significantly different between PT and NT subjects. Among TKR cases PT cases appear to carry a non-significantly larger proportion of risk variants. Among THR cases a lower contribution of genetic risk factors was seen in PT cases that did not achieve statistical significance.
One possible explanation for the lack of differences between PT and primary cases is that the genome-wide association study (GWAS) analyses that identified these variants have not made a distinction between the two types of patients. Given that PT cases constitute a substantial proportion of all knee OA cases, any genetic risk factors identified are likely to be implicated in both types of patients. This is possible since in the case of hip OA the prevalence of hip injury is much lower and the hip genetic risk load among PT THR cases is lower than in primary cases. However the effect of knee risk variants is as strong or even stronger among PT TKR cases indicating that the genes so far identified to influence risk of knee OA play a role in the pathogenesis of severe PT knee OA. However, the OR for hip genes in the PT group is not significantly different from 1 and with the current sample size and the effect size observed the study is not powered to determine if PT THR has a modest genetic component.
Recent research supports the view that acute joint damage that occurs at the time of an injury initiates a sequence of events that can lead to progressive articular surface damage. Determination of the risk factors affecting joint tissues and their respective roles in disease progression is critical to advances in the treatment of PT OA.4 Our data indicate that among the risk factors to take into account in PT knee OA are genetic influences.
There are some limitations to this study. Prior injuries were self-reported rather than based on medical records, and may be open to recall bias. Also some of the cases classified as PT OA may actually be ‘primary’ cases due to the type of injury experienced. Nevertheless, using the current definitions we found strong significant differences in age and gender between PT OA and NT cases, as reported in other studies.3 ,4 Moreover, the questionnaire applied identified a very high risk due to injury for TKR (OR=3.27) and THR (OR=5.08) in the case control design. Such risks are much higher than all genetic factors combined. Thus, in spite of the caveats, the population classified as PT based on questionnaire data appears enriched for true PT OA cases. Another limitation is that the results apply to UK populations and may not generalise to other ethnic groups, but the genetic variants used are only those reported in Caucasians as having a role in risk of OA.
In conclusion, PT hip OA cases appear to have a lower genetic risk load than NT cases, whereas PT knee OA cases have at least as high a genetic contribution as NT cases. These data support the perspective of OA as common complex condition and further challenge a clear separation between ‘primary’ and ‘secondary’ forms of OA.
Handling editor Tore K Kvien
Contributors All authors contributed to the study design, data interpretation and the final manuscript. In addition, AMV analysed and interpreted the data and prepared manuscript. SAD and MD evaluated study subjects. AMV and MD supervised the study.
Funding Supported by EC-FP7 programme grant 200800 TREAT-OA and 277849 EurHEALTHAgeing. AstraZeneca UK funded the GOAL study sample and data collection.
Competing interests RAM is an employee of Astrazeneca plc.
Ethics approval the research ethics committees of Nottingham City Hospital and North Nottinghamshire.
Provenance and peer review Not commissioned; externally peer reviewed.
Open Access This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/3.0/ and http://creativecommons.org/licenses/by-nc/3.0/legalcode