You are here

Article Text

Article menu

Download PDFPDF

XML
Osteoarthritis
Risk of venous thromboembolism in knee, hip and hand osteoarthritis: a general population-based cohort study
  1. Chao Zeng1,2,3,4,5,
  2. Kim Bennell6,
  3. Zidan Yang4,
  4. Uyen-Sa D T Nguyen7,
  5. Na Lu8,
  6. Jie Wei2,3,9,
  7. Guanghua Lei1,4,5,
  8. Yuqing Zhang2,3
  1. 1 Department of Orthopaedics, Xiangya Hospital, Central South University, Changsha, China
  2. 2 Division of Rheumatology, Allergy, and Immunology, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, USA
  3. 3 The Mongan Institute, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, USA
  4. 4 Hunan Key Laboratory of Joint Degeneration and Injury, Xiangya Hospital, Central South University, Changsha, China
  5. 5 National Clinical Research Center of Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China
  6. 6 Centre for Health, Exercise and Sports Medicine, The University of Melbourne, Melbourne, Victoria, Australia
  7. 7 Department of Biostatistics and Epidemiology, University of North Texas Health Science Center, School of Public Health, Fort Worth, Texas, USA
  8. 8 Arthritis Research Centre, Richmond, British Columbia, Canada
  9. 9 Health Management Center, Xiangya Hospital, Central South University, Changsha, China
  1. Correspondence to Professor Guanghua Lei, Department of Orthopaedics, Xiangya Hospital, Central South University, Changsha, China; lei_guanghua{at}csu.edu.cn; Professor Yuqing Zhang, Division of Rheumatology, Allergy, and Immunology, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, USA; yzhang108{at}mgh.harvard.edu

Abstract

Objectives Osteoarthritis is a leading cause of immobility and joint replacement, two strong risk factors for venous thromboembolism (VTE). We aimed to examine the relation of knee, hip and hand osteoarthritis to the risk of VTE and investigate joint replacement as a potential mediator.

Methods We conducted three cohort studies using data from The Health Improvement Network. Up to five individuals without osteoarthritis were matched to each case of incident knee (n=20 696), hip (n=10 411) or hand (n=6329) osteoarthritis by age, sex, entry time and body mass index. We examined the relation of osteoarthritis to VTE (pulmonary embolism and deep vein thrombosis) using a multivariable Cox proportional hazard model.

Results VTE developed in 327 individuals with knee osteoarthritis and 951 individuals without osteoarthritis (2.7 vs 2.0 per 1000 person-years), with multivariable-adjusted HR being 1.38 (95% CI 1.23 to 1.56). The indirect effect (HR) of knee osteoarthritis on VTE through knee replacement was 1.07 (95% CI 1.01 to 1.15), explaining 24.8% of its total effect on VTE. Risk of VTE was higher in hip osteoarthritis than non-osteoarthritis (3.3 vs 1.8 per 1000 person-years; multivariable-adjusted HR=1.83, 95% CI 1.56 to 2.13). The indirect effect through hip replacement yielded an HR of 1.14 (95% CI 1.04 to 1.25), explaining 28.1% of the total effect. No statistically significant difference in VTE risk was observed between hand osteoarthritis and non-osteoarthritis (1.5 vs 1.6 per 1000 person-years; multivariable-adjusted HR=0.88, 95% CI 0.67 to 1.16).

Conclusion Our large population-based cohort study provides the first evidence that knee or hip osteoarthritis, but not hand osteoarthritis, was associated with an increased risk of VTE, and such an association was partially mediated through knee or hip replacement.

  • osteoarthritis
  • epidemiology
  • knee
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.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, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

https://doi.org/10.1136/annrheumdis-2020-217782

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Key messages

    What is already known about this subject?

    • Osteoarthritis (OA) is a leading cause of immobility and joint replacement, two strong risk factors for venous thromboembolism (VTE); however, no study has examined the association between OA and risk of VTE.

    • To date, patients with OA have hardly been considered a high-risk population for VTE.

    What does this study add?

    • Our novel findings showed that the risk of incident VTE increased approximately 40% among individuals with knee OA and 80% among individuals with hip OA compared with those without OA.

    • An increased risk of VTE in knee or hip OA was partly mediated through knee or hip replacement.

    How might this impact on clinical practice or future developments?

    • If our findings are confirmed by other studies, it will help identify a large proportion of high-risk population for VTE given that knee and hip OA are common among older individuals and have a growing prevalence due to the global ageing population and obesity epidemic.

    Introduction

    Osteoarthritis (OA) is a common joint disorder, often affecting the knee, hip and hand, and is a leading cause of pain and disability.1 The WHO estimated that approximately 80% of patients with knee or hip OA have movement limitation and 25% are unable to perform their major daily-living activities.2

    Venous thromboembolism (VTE) is a common thromboembolic condition and potentially fatal disorder. Approximately 30% of patients with pulmonary embolism (PE) die within the first year after PE diagnosis.3 4 Strong risk factors for VTE include immobility and surgery.3 Thus, one would expect that OA, especially knee or hip OA, may be associated with a high risk of VTE through immobility and knee or hip joint replacement, a common surgery for advanced OA disease. Conversely, compared with knee and hip OA, hand OA is less likely to lead to immobility or lower-limb joint replacement; thus, individuals with hand OA may not experience a higher risk of VTE as those with knee or hip OA when compared with the general population. Although OA is a common disease and VTE is a potentially fatal condition, to our knowledge, no study has specifically examined the association between OA and risk of VTE.

    We conducted three matched cohort studies to compare the risk of VTE, as well as PE and deep vein thrombosis (DVT) separately, among patients with knee, hip or hand OA with that among individuals without OA in a general population. In addition, we assessed the extent to which the association of knee or hip OA with VTE was mediated by knee or hip replacement.

    Methods

    Data source

    We used data from The Health Improvement Network (THIN), which is a prospectively gathered, anonymised medical database derived from the records of general practitioners (GPs) in the UK. It draws approximately 17 million participants from 790 general practices and is representative of the UK population in terms of demographics and medical conditions. The Read classification system is used to code specific diagnoses,5 whereas a dictionary based on the Multilex classification system is used to code drugs.6 The validity of THIN for use in clinical and epidemiological research studies has been demonstrated in a previous study.7

    Study design and cohort definition

    Eligible participants included individuals aged 50 years or older8 with ≥1 year of continuous enrolment with the general practice from January 2000 through December 2017. We conducted three matched cohort studies to compare the risk of VTE among individuals with clinically diagnosed incident knee, hip or hand OA only, respectively, compared with up to five individuals without OA (ie, without any site OA), matched by same age, sex, calendar year of entry (ie, the index date or date of incident OA diagnoses) and body mass index (BMI) (±0.5 kg/m2). Individuals with OA that occurred at multiple sites (eg, knee, hip and hand), or without indicating specific site of OA, or individuals without BMI data were not eligible. We used Read codes to define knee, hip and hand OA (online supplementary table S1) according to previous studies using THIN database.9–12 This approach has been preferred in validation studies, as opposed to other approaches, such as medical visits, referrals or prescription records.13–15 Individuals were excluded from the analysis if they had knee or hip replacement prior to the index date (date of incident OA diagnosis), had neither GP visit nor specialist referral within 1 year prior to the index date, or had a history of VTE prior to the index date. Individuals who had been prescribed non-steroidal anti-inflammatory drugs, paracetamol or opioids within 1 year prior to the index date were also excluded as the potential prevalent OA cases may introduce selection bias.16 We conducted two additional matched cohort studies to compare the risk of VTE in the incident knee or hip OA cohort with that in the incident hand OA cohort.

    Assessment of outcomes

    The primary outcome was VTE (the combined endpoint of PE and DVT), with secondary outcomes of PE and DVT, separately. An individual was considered to have developed the outcome of interest if he or she had a recorded Read code of PE or DVT (online supplementary table S1) and received anticoagulant treatment.17 In addition, since VTE is potentially fatal and some individuals might have died before receiving anticoagulation, we also included individuals with a recorded code of PE or DVT but without a prescription for anticoagulant or autopsy results if there was a fatal outcome within 1 month of the VTE diagnosis.18 19 This VTE definition has a confirmation rate of 94% in the UK Clinical Practice Research Datalink,17 a similar database to THIN.20

    Assessment of covariates

    Sociodemographic and anthropometric characteristics, lifestyle factors, comorbidities, as well as medication use were obtained from the database before the index date (see table 1).

    View this table:
    Table 1

    Baseline characteristics of individuals in the incident OA cohorts and their matched non-OA cohorts

    Assessment of potential mediator

    Because major surgery is a strong risk factor for VTE and a previous study has found that knee or hip replacement was associated with a 10-fold to 11-fold increased risk of VTE during 6 months of follow-up, but not thereafter,21 we considered knee or hip replacement a potential mediator in the causal pathways between knee or hip OA and VTE (online supplementary figure S1). We also considered other major surgeries (ie, surgeries other than knee and hip replacement) a potential mediator. The putative mediators (knee or hip replacement and other major surgeries) were assessed within 6 months prior to the development of VTE, death, disenrolment from a GP practice participating in THIN or the end of follow-up.21

    Statistical analysis

    We estimated that approximately 22 000, 12 000 and 7000 individuals with incident knee, hip and hand OA, with 110 000, 60 000 and 35 000 matched individuals without OA, respectively, will be identified from 2000 to 2017. Assuming the risk of incident VTE among matched individuals without OA is 1.0%, with type 2 error of 0.05, our study will have 80% power to detect a relative risk of 1.21, 1.29 and 1.39, respectively.

    We compared baseline characteristics between each incident OA cohort (knee, hip and hand) and its comparative cohort.22 Person-years of follow-up for each individual were computed as the amount of time from the index date to the date of the following events: VTE, death, disenrolment from a GP practice participating in THIN or the end of follow-up (31 December 2017), whichever occurred first. We plotted the cumulative incidence curves for VTE, PE and DVT in each cohort accounting for the competing risk of death,23 and fitted a Cox proportional hazards model to estimate HRs accounting for matched study design. In the multivariable Cox proportional hazard models, we adjusted for potential confounders (table 1). We tested the proportional hazards assumption using the Kolmogorov supremum test.24 Since some individuals without OA could develop OA during the follow-up, we conducted a sensitivity analysis by censoring their follow-up at the time when they developed OA.

    We performed mediation analyses to examine the extent to which any association between knee or hip OA and VTE (total effect) was mediated through knee or hip replacement. Using marginal structural models, we estimated the natural indirect effect (the effect of knee or hip OA on the risk of VTE mediated through knee or hip replacement) and the natural direct effect (the effect of knee or hip OA on the risk of VTE not mediated through knee or hip replacement), while controlling for potential confounders.25 Additionally, we estimated the percentage of the total effect of OA on the risk of VTE that was mediated via knee or hip replacement.26 We took the same approach to estimate the mediation effect via other major surgeries.

    For all analyses, missing values for smoking status and alcohol consumption were imputed using IVEware for SAS V.9.4, which is based on chained equations.27 To minimise random error, we imputed 20 data sets and then combined estimates from these data sets.28 29 Detailed information is described in online supplementary text, and the results from the corresponding complete case analyses are presented in online supplementary table S2.

    Results

    The flow charts depicting the selection process of individuals are shown in online supplementary figures S2-4. Overall, the characteristics of the knee, hip and hand OA cohorts were similar to their matched non-OA cohorts, respectively, with standardised differences ≤0.1 (table 1).30

    The risk of VTE was higher in the knee OA cohort than in its matched non-OA cohort (figure 1A). As shown in table 2, 327 new cases of VTE occurred in the knee OA cohort (2.7 per 1000 person-years) and 951 cases in the non-OA cohort (2.0 per 1000 person-years) during the follow-up period. The rate difference (RD) of VTE between the knee OA and non-OA cohorts was 0.7 (95% CI 0.4 to 1.0) per 1000 person-years and the multivariable-adjusted HR was 1.38 (95% CI 1.23 to 1.56) during the entire follow-up period. Since the proportional hazard assumption was violated (p<0.05), follow-up time was divided into ≤2, 4, 6, 8 and 10 years. The HR over 2 years was 1.60 (95% CI 1.29 to 1.98), 4 years was 1.54 (95% CI 1.32 to 1.81), 6 years was 1.49 (95% CI 1.30 to 1.71), 8 years was 1.41 (95% CI 1.24 to 1.61) and 10 years was 1.41 (95% CI 1.25 to 1.59). The risk of knee replacement was 10-fold greater among the knee OA cohort than its matched non-OA cohort (13.3% vs 1.3%), and the risk of VTE was 13-fold greater among those with knee replacement than those without knee replacement within 6 months prior to the occurrence of VTE, death, disenrolment or the end of follow-up (15.8% vs 1.2%). The indirect effect (HR) of knee OA on VTE through knee replacement was 1.07 (95% CI 1.01 to 1.15), explaining 24.8% of its total effect on VTE. The indirect effect (HR) of knee OA on VTE through other major surgeries was 1.01 (95% CI 0.94 to 1.07). Individuals in the knee OA cohort also experienced a higher risk of PE (figure 1B) or DVT (figure 1C) than those in the non-OA cohort. Sensitivity and complete case analyses did not change the results materially (table 3 and online supplementary table S2).

    Figure 1
    Figure 1

    Cumulative incidence of venous thromboembolism (A), pulmonary embolism (B) and deep vein thrombosis (C) in 20 696 individuals with incident knee OA as compared with 81 137 individuals without OA, matched by age, sex, entry time and BMI. BMI, body mass index; OA, osteoarthritis.

    View this table:
    Table 2

    Association between incident OA and risk of venous thromboembolism (primary analysis)

    The hip OA cohort also had a higher risk of VTE (3.3 per 1000 person-years) than its matched non-OA cohort (1.8 per 1000 person-years) (figure 2A). The RD of VTE between the hip OA and non-OA cohorts was 1.5 (95% CI 1.0 to 2.0) per 1000 person-years. The matched and multivariable-adjusted HRs for VTE were 1.83 (95% CI 1.57 to 2.13) and 1.83 (95% CI 1.56 to 2.13), respectively. The risk of hip replacement was 20-fold greater among the hip OA cohort than its matched non-OA cohort (36.0% vs 1.8%), and the risk of VTE was 12-fold greater among those with hip replacement than those without hip replacement within 6 months prior to the occurrence of VTE, death, disenrolment or the end of follow-up (14.9% vs 1.2%). The indirect effect through hip replacement yielded an HR of 1.14 (95% CI 1.04 to 1.25), explaining 28.1% of its total effect on VTE. The indirect effect (HR) of hip OA on VTE through other major surgeries was 1.01 (95% CI 0.92 to 1.10). Compared with the non-OA cohort, the hip OA cohort was also at a higher risk of either PE (figure 2B) or DVT (figure 2C). The results from the sensitivity and complete case analyses did not change materially (table 3 and online supplementary table S2).

    Figure 2
    Figure 2

    Cumulative incidence of venous thromboembolism (A), pulmonary embolism (B) and deep vein thrombosis (C) in 10 411 individuals with incident hip OA as compared with 41 594 individuals without OA, matched by age, sex, entry time and BMI. BMI, body mass index; OA, osteoarthritis.

    No statistically significant difference in the risk of VTE was observed between the hand OA cohort (1.5 per 1000 person-years) and its matched non-OA cohort (1.6 per 1000 person-years) (figure 3A). Compared with the non-OA cohort, the multivariable-adjusted HR was 0.88 (95% CI 0.67 to 1.16) (table 2). The hand OA cohort did not experience a higher risk of either PE (figure 3B) or DVT (figure 3C). The results from the sensitivity and complete case analyses did not change materially (table 3 and online supplementary table S2). In addition, the risk of VTE was higher in the incident knee OA cohort (multivariable-adjusted HR=1.46, 95% CI 1.05 to 2.02) and the hip OA cohort (multivariable-adjusted HR=1.87, 95% CI 1.07 to 3.26) than that in the hand OA cohort.

    Figure 3
    Figure 3

    Cumulative incidence of venous thromboembolism (A), pulmonary embolism (B) and deep vein thrombosis (C) in 6329 individuals with incident hand OA as compared with 25 206 individuals without OA, matched by age, sex, entry time and BMI. BMI, body mass index; OA, osteoarthritis.

    View this table:
    Table 3

    Association between incident OA and risk of venous thromboembolism (sensitivity analysis)

    Discussion

    Our novel findings showed that the risk of incident VTE increased approximately 40% among individuals with knee OA and 80% among individuals with hip OA compared with those without OA. An increased risk of VTE in knee or hip OA was partly mediated through knee or hip replacement.

    Comparison with previous studies

    To date, many studies have examined the association between OA and risk of cardiovascular events; results, however, have been inconclusive.31–51 Some studies reported that OA was associated with an increased risk of cardiovascular events,31 32 39 44 46 while others failed to confirm it.34 37 38 42 51 In addition, several studies found that knee or hip OA, but not hand OA, was associated with an increased risk of cardiovascular events, and immobility may be a potential explanation for such site-specific difference.40 41 47 Indeed, results from the Wuchuan Osteoarthritis Study showed that the relation of symptomatic knee OA to all-cause mortality was almost entirely (96.8%) mediated through its effect on walking disability.52 Results from the Osteoarthritis Initiative showed that 22.4% of the effect of symptomatic knee OA on all-cause mortality was mediated through disability in physical function and 26.5% was mediated through impaired physical component summary scores of quality of life.53 A similar finding was also reported in the Study of Osteoporotic Fractures, and the percentage of the total effect of hip OA on all-cause mortality mediated by physical function was 42.9%.32 A few studies also showed that hand OA may increase the risk of cardiovascular events, and OA-related systemic inflammation may be the potential explanation.35 36 49

    Possible explanations

    Several explanations may account for our findings. First, immobilisation or prolonged physical inactivity could substantially increase the risk of VTE through inducing muscular and diaphragm dysfunction, venous stasis, and a hypercoagulability state.3 54 Thus, OA-related immobility may play a key role in the development of VTE. Second, knee or hip replacement is associated with a substantially increased risk of VTE.21 In the current study, 24.8% of the effect of knee OA on the risk of VTE was mediated through knee replacement and 28.1% of the hip OA effect on VTE was mediated through hip replacement, suggesting knee and hip replacement surgery indeed played an important role in the association between knee or hip OA and risk of VTE.

    Strengths and limitations

    First, we matched individuals with OA with individuals without OA on age, sex, BMI and entry time, and adjusted for major potential confounders for VTE in the analyses. The matched and multivariable-adjusted effect estimates were consistent, suggesting that residual confounders should not greatly alter the observed associations. Second, we compared the risk of VTE among the incident instead of prevalent OA cases with those of individuals without OA, to minimise potential selection bias.16 Third, no apparent mediation effect was observed through other major surgeries except knee and hip replacement, indicating that the knee or hip OA cohort was comparable with its matched non-OA cohort in potential confounders. Finally, immobility and surgery, including knee or hip replacement, are the two major risk factors for VTE,3 21 lending credibility to our observations that the increased risk of VTE was only among patients with knee or hip OA, but not hand OA.

    Limitations of our study should be acknowledged. First, we used Read codes to define knee, hip and hand OA, and it is possible that some individuals with mild symptomatic OA may not seek care from their GPs and could be included in the comparison cohort. Such misclassification may bias our effect estimates in either way (ie, either towards or away from the null). However, when we censored the follow-up time for individuals who were diagnosed with OA during the follow-up, the results did not change materially. Second, the residual confounding still cannot be ruled out, although we used several approaches to control for potential confounding bias. For example, in THIN database, severity for most diseases was not collected. If severity of other comorbidities in the OA cohort differs from that in the non-OA cohort, this may lead to potential bias. Third, excluding patients with OA who had no BMI information or who cannot find matched individuals without OA may limit generalisability of our findings. Fourth, we cannot directly examine whether the association between OA and risk of VTE was mediated through immobility owing to lack of detailed information of immobility and physical activity level in THIN.

    Clinical implications

    VTE is a condition that can result in substantial morbidity and mortality; however, to date, the cause is still unknown in a third to a half of VTE cases.55 Our study showed that the incidence rates of VTE among patients with knee (2.7 per 1000 person-years) or hip (3.3 per 1000 person-years) OA were similar to that among patients with rheumatoid arthritis (3.3 per 1000 person-years) from the THIN database.18 Although the RD of VTE between the knee or hip OA and non-OA cohorts was relatively small (0.7 or 1.5 per 1000 person-years), the present study, if confirmed by other studies, will help identify a large proportion of high-risk population for VTE given that knee and hip OA are much more common than rheumatoid arthritis and have a growing prevalence due to the global ageing population and obesity epidemic.1

    Clinicians managing patients with knee or hip OA, especially those with movement limitations, could then have a high index of suspicion for VTE, particularly in the presence of non-specific symptoms and signs. Our findings also help shed light on potential mechanisms linking OA to VTE and may guide development of evidence-based preventive strategies for VTE. For example, we postulate that immobility due to knee or hip OA is likely to play an important role linking knee and hip OA with an elevated risk of VTE; thus, strategies aimed at encouraging patients with knee or hip OA to actively participate in physical activities seem justified.56–61 In addition, as joint replacement surgery, a clinically relevant and cost-effective treatment for end-stage OA, continuously grows worldwide, more efforts (eg, early mobilisation following surgery) are still needed to further reduce the risk of VTE after surgery.62 63

    Conclusion

    Our study provides the first evidence that knee or hip OA, but not hand OA, was associated with an increased risk of VTE, and such an association was partially mediated through knee or hip replacement. Since OA is a prevalent disease and VTE is a serious condition, further studies are needed to confirm our findings and to understand biological mechanisms, so that appropriate prevention, monitoring and treatment strategies can be developed to reduce the disease burden in patients with knee or hip OA.

    References

    1. ↵
      2. Hunter DJ ,
      3. Bierma-Zeinstra S
      . Osteoarthritis. Lancet 2019;393:1745–59.doi:10.1016/S0140-6736(19)30417-9 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31034380
      OpenUrlPubMed
    2. ↵
      2. Brooks PM
      . Impact of osteoarthritis on individuals and society: how much disability? social consequences and health economic implications. Curr Opin Rheumatol 2002;14:573–7.doi:10.1097/00002281-200209000-00017 pmid:http://www.ncbi.nlm.nih.gov/pubmed/12192258
      OpenUrlCrossRefPubMedWeb of Science
    3. ↵
      2. Di Nisio M ,
      3. van Es N ,
      4. Büller HR
      . Deep vein thrombosis and pulmonary embolism. Lancet 2016;388:3060–73.doi:10.1016/S0140-6736(16)30514-1 pmid:http://www.ncbi.nlm.nih.gov/pubmed/27375038
      OpenUrlPubMed
    4. ↵
      2. Tritschler T ,
      3. Kraaijpoel N ,
      4. Le Gal G , et al
      . Venous thromboembolism: advances in diagnosis and treatment. JAMA 2018;320:1583–94.doi:10.1001/jama.2018.14346 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30326130
      OpenUrlCrossRefPubMed
    5. ↵
      2. Chisholm J
      . The read clinical classification. BMJ 1990;300:1092. doi:10.1136/bmj.300.6732.1092 pmid:http://www.ncbi.nlm.nih.gov/pubmed/2344534
      OpenUrlFREE Full Text
    6. ↵
      1. First Databank
      . Multilex. Available: http://www.folbhealth.co.uk/solutions/multilex/ [Accessed 8 Nov 2017].
    7. ↵
      2. Lewis JD ,
      3. Schinnar R ,
      4. Bilker WB , et al
      . Validation studies of the health improvement network (THIN) database for pharmacoepidemiology research. Pharmacoepidemiol Drug Saf 2007;16:393–401.doi:10.1002/pds.1335 pmid:17066486
      OpenUrlCrossRefPubMedWeb of Science
    8. ↵
      2. Heit JA
      . Epidemiology of venous thromboembolism. Nat Rev Cardiol 2015;12:464–74.doi:10.1038/nrcardio.2015.83 pmid:http://www.ncbi.nlm.nih.gov/pubmed/26076949
      OpenUrlCrossRefPubMed
    9. ↵
      2. Zeng C ,
      3. Dubreuil M ,
      4. LaRochelle MR , et al
      . Association of tramadol with all-cause mortality among patients with osteoarthritis. JAMA 2019;321:969–82.doi:10.1001/jama.2019.1347 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30860559
      OpenUrlPubMed
    10. ↵
      2. Dubreuil M ,
      3. Louie-Gao Q ,
      4. Peloquin CE , et al
      . Risk of myocardial infarction with use of selected non-steroidal anti-inflammatory drugs in patients with spondyloarthritis and osteoarthritis. Ann Rheum Dis 2018;77:1137–42.doi:10.1136/annrheumdis-2018-213089 pmid:http://www.ncbi.nlm.nih.gov/pubmed/29674321
      OpenUrlAbstract/FREE Full Text
    11. ↵
      2. Neogi T ,
      3. Li S ,
      4. Peloquin C , et al
      . Effect of bisphosphonates on knee replacement surgery. Ann Rheum Dis 2018;77:92–7.doi:10.1136/annrheumdis-2017-211811 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28988204
      OpenUrlAbstract/FREE Full Text
    12. ↵
      2. Misra D ,
      3. Lu N ,
      4. Felson D , et al
      . Does knee replacement surgery for osteoarthritis improve survival? the jury is still out. Ann Rheum Dis 2017;76:140–6.doi:10.1136/annrheumdis-2016-209167 pmid:http://www.ncbi.nlm.nih.gov/pubmed/27190096
      OpenUrlAbstract/FREE Full Text
    13. ↵
      2. Lix L ,
      3. Yogendran M ,
      4. Burchill C
      . Defining and validating chronic diseases: an administrative data approach, 2006. Available: http://www.umanitoba.ca/centres/mchp/reports.htm [Accessed 29 Sept 2019].
    14. ↵
      2. Lix L ,
      3. Yogendran M ,
      4. Mann J
      . Defining and validating chronic diseases: an administrative data approach. An update with ICD-10-CA, 2008. Available: https://www.researchgate.net/publication/242284353 [Accessed 29 Sept 2019].
    15. ↵
      2. Prieto-Alhambra D ,
      3. Judge A ,
      4. Javaid MK , et al
      . Incidence and risk factors for clinically diagnosed knee, hip and hand osteoarthritis: influences of age, gender and osteoarthritis affecting other joints. Ann Rheum Dis 2014;73:1659–64.doi:10.1136/annrheumdis-2013-203355 pmid:http://www.ncbi.nlm.nih.gov/pubmed/23744977
      OpenUrlAbstract/FREE Full Text
    16. ↵
      2. Choi HK ,
      3. Nguyen U-S ,
      4. Niu J , et al
      . Selection bias in rheumatic disease research. Nat Rev Rheumatol 2014;10:403–12.doi:10.1038/nrrheum.2014.36 pmid:http://www.ncbi.nlm.nih.gov/pubmed/24686510
      OpenUrlCrossRefPubMed
    17. ↵
      2. Huerta C ,
      3. Johansson S ,
      4. Wallander M-A , et al
      . Risk factors and short-term mortality of venous thromboembolism diagnosed in the primary care setting in the United Kingdom. Arch Intern Med 2007;167:935–43.doi:10.1001/archinte.167.9.935 pmid:http://www.ncbi.nlm.nih.gov/pubmed/17502535
      OpenUrlCrossRefPubMedWeb of Science
    18. ↵
      2. Choi HK ,
      3. Rho Y-H ,
      4. Zhu Y , et al
      . The risk of pulmonary embolism and deep vein thrombosis in rheumatoid arthritis: a UK population-based outpatient cohort study. Ann Rheum Dis 2013;72:1182–7.doi:10.1136/annrheumdis-2012-201669 pmid:http://www.ncbi.nlm.nih.gov/pubmed/22930596
      OpenUrlAbstract/FREE Full Text
    19. ↵
      2. Ogdie A ,
      3. Kay McGill N ,
      4. Shin DB , et al
      . Risk of venous thromboembolism in patients with psoriatic arthritis, psoriasis and rheumatoid arthritis: a general population-based cohort study. Eur Heart J 2018;39:3608–14.doi:10.1093/eurheartj/ehx145 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28444172
      OpenUrlCrossRefPubMed
    20. ↵
      2. Carbonari DM ,
      3. Saine ME ,
      4. Newcomb CW , et al
      . Use of demographic and pharmacy data to identify patients included within both the Clinical Practice Research Datalink (CPRD) and The Health Improvement Network (THIN). Pharmacoepidemiol Drug Saf 2015;24:999–1003.doi:10.1002/pds.3844 pmid:26213344
      OpenUrlPubMed
    21. ↵
      2. Lu N ,
      3. Misra D ,
      4. Neogi T , et al
      . Total joint arthroplasty and the risk of myocardial infarction: a general population, propensity score-matched cohort study. Arthritis Rheumatol 2015;67:2771–9.doi:10.1002/art.39246 pmid:http://www.ncbi.nlm.nih.gov/pubmed/26331443
      OpenUrlPubMed
    22. ↵
      2. Mamdani M ,
      3. Sykora K ,
      4. Li P , et al
      . Reader's guide to critical appraisal of cohort studies: 2. assessing potential for confounding. BMJ 2005;330:960–2.doi:10.1136/bmj.330.7497.960 pmid:http://www.ncbi.nlm.nih.gov/pubmed/15845982
      OpenUrlFREE Full Text
    23. ↵
      2. Gooley TA ,
      3. Leisenring W ,
      4. Crowley J , et al
      . Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. Stat Med 1999;18:695–706.doi:10.1002/(SICI)1097-0258(19990330)18:6<695::AID-SIM60>3.0.CO;2-O pmid:http://www.ncbi.nlm.nih.gov/pubmed/10204198
      OpenUrlCrossRefPubMedWeb of Science
    24. ↵
      2. Lin DY ,
      3. Wei LJ ,
      4. Ying Z
      . Checking the COX model with cumulative sums of martingale-based residuals. Biometrika 1993;80:557–72.doi:10.1093/biomet/80.3.557
      OpenUrlCrossRefWeb of Science
    25. ↵
      2. Lange T ,
      3. Vansteelandt S ,
      4. Bekaert M
      . A simple unified approach for estimating natural direct and indirect effects. Am J Epidemiol 2012;176:190–5.doi:10.1093/aje/kwr525 pmid:http://www.ncbi.nlm.nih.gov/pubmed/22781427
      OpenUrlCrossRefPubMedWeb of Science
    26. ↵
      2. VanderWeele TJ
      . Mediation analysis: a practitioner's guide. Annu Rev Public Health 2016;37:17–32.doi:10.1146/annurev-publhealth-032315-021402 pmid:http://www.ncbi.nlm.nih.gov/pubmed/26653405
      OpenUrlCrossRefPubMed
    27. ↵
      2. van Buuren S
      . Flexible imputation of missing data. 2nd Edition. Boca Raton: Chapman & Hall/CRC, 2018.
    28. ↵
      2. Rubin DB
      . Multiple imputation for nonresponse in surveys. New York: J Wiley, 1987.
    29. ↵
      2. Sterne JAC ,
      3. White IR ,
      4. Carlin JB , et al
      . Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ 2009;338:b2393. doi:10.1136/bmj.b2393 pmid:http://www.ncbi.nlm.nih.gov/pubmed/19564179
      OpenUrlFREE Full Text
    30. ↵
      2. Franklin JM ,
      3. Rassen JA ,
      4. Ackermann D , et al
      . Metrics for covariate balance in cohort studies of causal effects. Stat Med 2014;33:1685–99.doi:10.1002/sim.6058 pmid:http://www.ncbi.nlm.nih.gov/pubmed/24323618
      OpenUrlCrossRefPubMed
    31. ↵
      2. Atiquzzaman M ,
      3. Karim ME ,
      4. Kopec J , et al
      . Role of nonsteroidal antiinflammatory drugs in the association between osteoarthritis and cardiovascular diseases: a longitudinal study. Arthritis Rheumatol 2019;71:1835–43.doi:10.1002/art.41027 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31389178
      OpenUrlPubMed
    32. ↵
      2. Barbour KE ,
      3. Lui L-Y ,
      4. Nevitt MC , et al
      . Hip osteoarthritis and the risk of all-cause and disease-specific mortality in older women: a population-based cohort study. Arthritis Rheumatol 2015;67:1798–805.doi:10.1002/art.39113 pmid:http://www.ncbi.nlm.nih.gov/pubmed/25778744
      OpenUrlPubMed
    33. ↵
      2. Cleveland RJ ,
      3. Alvarez C ,
      4. Schwartz TA , et al
      . The impact of painful knee osteoarthritis on mortality: a community-based cohort study with over 24 years of follow-up. Osteoarthritis Cartilage 2019;27:593–602.doi:10.1016/j.joca.2018.12.008 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30583096
      OpenUrlPubMed
    34. ↵
      2. Corsi M ,
      3. Alvarez C ,
      4. Callahan LF , et al
      . Contributions of symptomatic osteoarthritis and physical function to incident cardiovascular disease. BMC Musculoskelet Disord 2018;19:393. doi:10.1186/s12891-018-2311-4 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30414614
      OpenUrlPubMed
    35. ↵
      2. Haara MM ,
      3. Manninen P ,
      4. Kröger H , et al
      . Osteoarthritis of finger joints in Finns aged 30 or over: prevalence, determinants, and association with mortality. Ann Rheum Dis 2003;62:151–8.doi:10.1136/ard.62.2.151 pmid:http://www.ncbi.nlm.nih.gov/pubmed/12525385
      OpenUrlAbstract/FREE Full Text
    36. ↵
      2. Haugen IK ,
      3. Ramachandran VS ,
      4. Misra D , et al
      . Hand osteoarthritis in relation to mortality and incidence of cardiovascular disease: data from the Framingham heart study. Ann Rheum Dis 2015;74:74–81.doi:10.1136/annrheumdis-2013-203789 pmid:http://www.ncbi.nlm.nih.gov/pubmed/24047870
      OpenUrlAbstract/FREE Full Text
    37. ↵
      2. Hoeven TA ,
      3. Leening MJG ,
      4. Bindels PJ , et al
      . Disability and not osteoarthritis predicts cardiovascular disease: a prospective population-based cohort study. Ann Rheum Dis 2015;74:752–6.doi:10.1136/annrheumdis-2013-204388 pmid:http://www.ncbi.nlm.nih.gov/pubmed/24385204
      OpenUrlAbstract/FREE Full Text
    38. ↵
      2. Holbrook TL ,
      3. Wingard DL ,
      4. Barrett-Connor E
      . Self-Reported arthritis among men and women in an adult community. J Community Health 1990;15:195–208.doi:10.1007/BF01350257 pmid:http://www.ncbi.nlm.nih.gov/pubmed/2365841
      OpenUrlCrossRefPubMed
    39. ↵
      2. Hsu P-S ,
      3. Lin H-H ,
      4. Li C-R , et al
      . Increased risk of stroke in patients with osteoarthritis: a population-based cohort study. Osteoarthritis Cartilage 2017;25:1026–31.doi:10.1016/j.joca.2016.10.027 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28300652
      OpenUrlPubMed
    40. ↵
      2. Kendzerska T ,
      3. Jüni P ,
      4. King LK , et al
      . The longitudinal relationship between hand, hip and knee osteoarthritis and cardiovascular events: a population-based cohort study. Osteoarthritis Cartilage 2017;25:1771–80.doi:10.1016/j.joca.2017.07.024 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28801210
      OpenUrlPubMed
    41. ↵
      2. Kluzek S ,
      3. Sanchez-Santos MT ,
      4. Leyland KM , et al
      . Painful knee but not hand osteoarthritis is an independent predictor of mortality over 23 years follow-up of a population-based cohort of middle-aged women. Ann Rheum Dis 2016;75:1749–56.doi:10.1136/annrheumdis-2015-208056 pmid:http://www.ncbi.nlm.nih.gov/pubmed/26543059
      OpenUrlAbstract/FREE Full Text
    42. ↵
      2. Liu R ,
      3. Kwok WY ,
      4. Vliet Vlieland TPM , et al
      . Mortality in osteoarthritis patients. Scand J Rheumatol 2015;44:70–3.doi:10.3109/03009742.2014.922213 pmid:http://www.ncbi.nlm.nih.gov/pubmed/25179456
      OpenUrlCrossRefPubMed
    43. ↵
      2. Mendy A ,
      3. Park J ,
      4. Vieira ER
      . Osteoarthritis and risk of mortality in the USA: a population-based cohort study. Int J Epidemiol 2018;47:1821–9.doi:10.1093/ije/dyy187 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30169829
      OpenUrlPubMed
    44. ↵
      2. Nüesch E ,
      3. Dieppe P ,
      4. Reichenbach S , et al
      . All cause and disease specific mortality in patients with knee or hip osteoarthritis: population based cohort study. BMJ 2011;342:d1165. doi:10.1136/bmj.d1165 pmid:http://www.ncbi.nlm.nih.gov/pubmed/21385807
      OpenUrlAbstract/FREE Full Text
    45. ↵
      2. Rahman MM ,
      3. Kopec JA ,
      4. Anis AH , et al
      . Risk of cardiovascular disease in patients with osteoarthritis: a prospective longitudinal study. Arthritis Care Res 2013;65:1951–8.doi:10.1002/acr.22092 pmid:http://www.ncbi.nlm.nih.gov/pubmed/23925995
      OpenUrlPubMed
    46. ↵
      2. Tsuboi M ,
      3. Hasegawa Y ,
      4. Matsuyama Y , et al
      . Do musculoskeletal degenerative diseases affect mortality and cause of death after 10 years in Japan? J Bone Miner Metab 2011;29:217–23.doi:10.1007/s00774-010-0214-z pmid:http://www.ncbi.nlm.nih.gov/pubmed/20711854
      OpenUrlCrossRefPubMed
    47. ↵
      2. Turkiewicz A ,
      3. Kiadaliri AA ,
      4. Englund M
      . Cause-Specific mortality in osteoarthritis of peripheral joints. Osteoarthritis Cartilage 2019;27:848–54.doi:10.1016/j.joca.2019.02.793 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30797945
      OpenUrlPubMed
    48. ↵
      2. Veronese N ,
      3. Cereda E ,
      4. Maggi S , et al
      . Osteoarthritis and mortality: a prospective cohort study and systematic review with meta-analysis. Semin Arthritis Rheum 2016;46:160–7.doi:10.1016/j.semarthrit.2016.04.002 pmid:http://www.ncbi.nlm.nih.gov/pubmed/27179749
      OpenUrlPubMed
    49. ↵
      2. Veronese N ,
      3. Stubbs B ,
      4. Solmi M , et al
      . Osteoarthristis increases the risk of cardiovascular disease: data from the osteoarthritis initiative. J Nutr Health Aging 2018;22:371–6.doi:10.1007/s12603-017-0941-0 pmid:http://www.ncbi.nlm.nih.gov/pubmed/29484350
      OpenUrlPubMed
    50. ↵
      2. Veronese N ,
      3. Trevisan C ,
      4. De Rui M , et al
      . Association of osteoarthritis with increased risk of cardiovascular diseases in the elderly: findings from the Progetto Veneto Anziano study cohort. Arthritis Rheumatol 2016;68:1136–44.doi:10.1002/art.39564 pmid:http://www.ncbi.nlm.nih.gov/pubmed/26714002
      OpenUrlPubMed
    51. ↵
      2. Watson DJ ,
      3. Rhodes T ,
      4. Guess HA
      . All-Cause mortality and vascular events among patients with rheumatoid arthritis, osteoarthritis, or no arthritis in the UK general practice research database. J Rheumatol 2003;30:1196–202.pmid:http://www.ncbi.nlm.nih.gov/pubmed/12784389
      OpenUrlAbstract/FREE Full Text
    52. ↵
      2. Liu Q ,
      3. Niu J ,
      4. Li H , et al
      . Knee symptomatic osteoarthritis, walking disability, NSAIDs use and all-cause mortality: population-based Wuchuan osteoarthritis study. Sci Rep 2017;7:3309. doi:10.1038/s41598-017-03110-3 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28607349
      OpenUrlPubMed
    53. ↵
      2. Wang Y ,
      3. Nguyen U-SDT ,
      4. Lane NE , et al
      . Knee osteoarthritis, potential mediators, and risk of all-cause mortality: data from the osteoarthritis initiative. Arthritis Care Res 2020. doi:doi:10.1002/acr.24151. [Epub ahead of print: 21 Jan 2020].pmid:http://www.ncbi.nlm.nih.gov/pubmed/31961495
    54. ↵
      2. Pottier P ,
      3. Hardouin JB ,
      4. Lejeune S , et al
      . Immobilization and the risk of venous thromboembolism. A meta-analysis on epidemiological studies. Thromb Res 2009;124:468–76.doi:10.1016/j.thromres.2009.05.006 pmid:http://www.ncbi.nlm.nih.gov/pubmed/19481781
      OpenUrlCrossRefPubMed
    55. ↵
      2. Kearon C ,
      3. Ageno W ,
      4. Cannegieter SC , et al
      . Categorization of patients as having provoked or unprovoked venous thromboembolism: guidance from the SSC of ISTH. J Thromb Haemost 2016;14:1480–3.doi:10.1111/jth.13336 pmid:http://www.ncbi.nlm.nih.gov/pubmed/27428935
      OpenUrlCrossRefPubMed
    56. ↵
      2. Lear SA ,
      3. Hu W ,
      4. Rangarajan S , et al
      . The effect of physical activity on mortality and cardiovascular disease in 130 000 people from 17 high-income, middle-income, and low-income countries: the PURE study. Lancet 2017;390:2643–54.doi:10.1016/S0140-6736(17)31634-3 pmid:http://www.ncbi.nlm.nih.gov/pubmed/28943267
      OpenUrlCrossRefPubMed
    57. ↵
      2. Mok A ,
      3. Khaw K-T ,
      4. Luben R , et al
      . Physical activity trajectories and mortality: population based cohort study. BMJ 2019;365:l2323.doi:10.1136/bmj.l2323 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31243014
      OpenUrlAbstract/FREE Full Text
    58. ↵
      2. Ekelund U ,
      3. Tarp J ,
      4. Steene-Johannessen J , et al
      . Dose-Response associations between accelerometry measured physical activity and sedentary time and all cause mortality: systematic review and harmonised meta-analysis. BMJ 2019;366:l4570.doi:10.1136/bmj.l4570 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31434697
      OpenUrlAbstract/FREE Full Text
    59. ↵
      1. National Clinical Guideline Centre (UK)
      . Osteoarthritis: care and management in adults, 2014. Available: https://www.ncbi.nlm.nih.gov/books/NBK248069/pdf/Bookshelf_NBK248069.pdf [Accessed 8 Nov 2019].
    60. ↵
      2. Bannuru RR ,
      3. Osani MC ,
      4. Vaysbrot EE , et al
      . OARSI guidelines for the non-surgical management of knee, hip, and polyarticular osteoarthritis. Osteoarthritis Cartilage 2019;27:1578–89.doi:10.1016/j.joca.2019.06.011 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31278997
      OpenUrlCrossRefPubMed
    61. ↵
      2. Kolasinski SL ,
      3. Neogi T ,
      4. Hochberg MC , et al
      . 2019 American College of Rheumatology/Arthritis Foundation guideline for the management of osteoarthritis of the hand, hip, and knee. Arthritis Rheumatol 2020;72:220–33.doi:10.1002/art.41142 pmid:http://www.ncbi.nlm.nih.gov/pubmed/31908163
      OpenUrlPubMed
    62. ↵
      2. Price AJ ,
      3. Alvand A ,
      4. Troelsen A , et al
      . Knee replacement. Lancet 2018;392:1672–82.doi:10.1016/S0140-6736(18)32344-4 pmid:http://www.ncbi.nlm.nih.gov/pubmed/30496082
      OpenUrlCrossRefPubMed
    63. ↵
      2. Ferguson RJ ,
      3. Palmer AJ ,
      4. Taylor A , et al
      . Hip replacement. Lancet 2018;392:1662–71.doi:10.1016/S0140-6736(18)31777-X pmid:http://www.ncbi.nlm.nih.gov/pubmed/30496081
      OpenUrlPubMed
    View Abstract

    Footnotes

    • Handling editor Josef S Smolen

    • Contributors YZ and GL had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. GL and YZ are joint corresponding authors. All authors have read, provided critical feedback on intellectual content and approved the final manuscript. Concept and design: GL, YZ, CZ. Acquisition, analysis or interpretation of data: all authors. Drafting of the manuscript: CZ, GL, YZ. Critical revision of the manuscript for important intellectual content: all authors. Statistical analysis: ZY, NL, JW, YZ. Obtained funding: CZ, KB, JW, GL. Administrative, technical or material support: CZ, GL, YZ. Supervision: GL, YZ.

    • Funding This work was supported by the National Natural Science Foundation of China (81772413, 81702207, 81702206, 81930071), the National Health and Medical Research Council Principal Research Fellowship, the National Key Research and Development Project (2018YFB1105705), and the Key Research and Development Program of Hunan Province (2018SK2070, 2018SK2071).

    • Competing interests None declared.

    • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

    • Patient consent for publication Not required.

    • Ethics approval This study received approval from the medical ethical committee at Xiangya Hospital, with waiver of informed consent. This study was approved by the THIN Scientific Review Committee (20SRC003). THIN is a registered trademark of Cegedim SA in the UK and other countries. Reference made to the THIN database is intended to be descriptive of the data asset licensed by IQVIA. This work uses de-identified data provided by patients as part of their routine primary care.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data availability statement No data are available. While we are unable to share the data sets, code lists are available on request.