Advertisement
Research Article| Volume 24, ISSUE 4, P597-604, April 2016

Download started.

Ok

Long term use of analgesics and risk of osteoarthritis progressions and knee replacement: propensity score matched cohort analysis of data from the Osteoarthritis Initiative

Open ArchivePublished:November 09, 2015DOI:https://doi.org/10.1016/j.joca.2015.11.003

      Summary

      Objectives

      To determine the association between the long-term use of analgesics and progression of osteoarthritis (OA) as evidenced by up to 3-years follow-up worsening of radiographic Kellgren–Lawrence (KL) grade and incidence of knee replacement (KR).

      Design

      Using nearest neighbor matching of the propensity scores with caliper in the Osteoarthritis Initiative (OAI) cohort, 173 index (Analgesic +) and 173 referent (Analgesic −) subjects were included. Analgesic + and − subjects had analgesics in all and none of their visits, respectively. Analgesic + and − subjects were balanced in their demographics, baseline, first, second and third year body mass index (BMI), Western Ontario and McMaster (WOMAC) total score, Physical and Mental health summary scales (SF-12), Physical Activity Scale for the Elderly (PASE) and Charleston Comorbidity Scale. Analgesic + and − subjects were also matched for baseline radiographic KL grade. Interval increase in the KL grade and incidence of KR were defined as the outcome.

      Results

      Included subjects had average 6.5 years of follow-up. By the third year, 44 subjects had an interval increase in the KL grade; 29 in Analgesic + and 15 among Analgesic − subjects (P = 0.024). By the eighth-year, 41 subjects had their first KR; 29 in Analgesic + and 12 among Analgesic − subjects (P = 0.005). Hazard Ratio (HR) of OA progression and KR for Analgesic + subjects was 1.91 (1.02–3.57) and 2.57 (1.31–5.04), respectively.

      Conclusions

      Long-term use of analgesics may be associated with radiographic progression of knee OA and increased risk of future KR.

      Keywords

      Introduction

      Background and rationale

      Knee Osteoarthritis (OA) is one of the most prevalent chronic joint disorders in the United States
      • Lawrence R.C.
      • Felson D.T.
      • Helmick C.G.
      • Arnold L.M.
      • Choi H.
      • Deyo R.A.
      • et al.
      Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II.
      . Knee OA has a progressive course, which is associated with gradually worsening symptoms such as chronic pain, stiffness and restricted range of motion
      • Felson D.T.
      • Zhang Y.
      • Hannan M.T.
      • Naimark A.
      • Weissman B.N.
      • Aliabadi P.
      • et al.
      The incidence and natural history of knee osteoarthritis in the elderly. The Framingham Osteoarthritis Study.
      . Knee replacement (KR) is recognized as the ultimate outcome in the progression of OA
      • Raynauld J.P.
      • Martel-Pelletier J.
      • Haraoui B.
      • Choquette D.
      • Dorais M.
      • Wildi L.M.
      • et al.
      Risk factors predictive of joint replacement in a 2-year multicentre clinical trial in knee osteoarthritis using MRI: results from over 6 years of observation.
      . Various medical treatments can be offered prior to KR in order to improve symptoms and function
      • Bruyere O.
      • Cooper C.
      • Pelletier J.P.
      • Branco J.
      • Luisa Brandi M.
      • Guillemin F.
      • et al.
      An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
      . Several guidelines have been published to standardize the recommendations regarding the medical therapy and use of analgesics in OA patients
      • Bruyere O.
      • Cooper C.
      • Pelletier J.P.
      • Branco J.
      • Luisa Brandi M.
      • Guillemin F.
      • et al.
      An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
      • McAlindon T.E.
      • Bannuru R.R.
      • Sullivan M.C.
      • Arden N.K.
      • Berenbaum F.
      • Bierma-Zeinstra S.M.
      • et al.
      OARSI guidelines for the non-surgical management of knee osteoarthritis.
      • Jevsevar D.S.
      • Brown G.A.
      • Jones D.L.
      • Matzkin E.G.
      • Manner P.A.
      • Mooar P.
      • et al.
      The American Academy of Orthopaedic Surgeons evidence-based guideline on: treatment of osteoarthritis of the knee.
      • Hochberg M.C.
      • Altman R.D.
      • April K.T.
      • Benkhalti M.
      • Guyatt G.
      • McGowan J.
      • et al.
      American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee.
      . Analgesics are the largest group of medications that are prescribed for OA patients. A variety of analgesics including non-steroidal anti-inflammatory drugs (NSAIDS), Cyclooxygenase II (COX II) inhibitors, acetaminophen, salicylates and narcotic analgesics are used for symptom relief
      • Ding C.
      Do NSAIDs affect the progression of osteoarthritis?.
      • Reijman M.
      • Bierma-Zeinstra S.M.
      • Pols H.A.
      • Koes B.W.
      • Stricker B.H.
      • Hazes J.M.
      Is there an association between the use of different types of nonsteroidal antiinflammatory drugs and radiologic progression of osteoarthritis? The Rotterdam Study.
      • Martel-Pelletier J.
      • Roubille C.
      • Abram F.
      • Hochberg M.C.
      • Dorais M.
      • Delorme P.
      • et al.
      First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
      • Wolfe F.
      • Michaud K.
      • Burke T.A.
      • Zhao S.Z.
      Longer use of COX-2-specific inhibitors compared to nonspecific nonsteroidal antiinflammatory drugs: a longitudinal study of 3639 patients in community practice.
      • Wilcox C.M.
      • Cryer B.
      • Triadafilopoulos G.
      Patterns of use and public perception of over-the-counter pain relievers: focus on nonsteroidal antiinflammatory drugs.
      • Adatia A.
      • Rainsford K.D.
      • Kean W.F.
      Osteoarthritis of the knee and hip. Part II: therapy with ibuprofen and a review of clinical trials.
      • Towheed T.E.
      • Maxwell L.
      • Judd M.G.
      • Catton M.
      • Hochberg M.C.
      • Wells G.
      Acetaminophen for osteoarthritis.
      • Fendrick A.M.
      • Greenberg B.P.
      A review of the benefits and risks of nonsteroidal anti-inflammatory drugs in the management of mild-to-moderate osteoarthritis.
      • Pavelka K.
      A comparison of the therapeutic efficacy of diclofenac in osteoarthritis: a systematic review of randomised controlled trials.
      • Bjordal J.M.
      • Klovning A.
      • Ljunggren A.E.
      • Slordal L.
      Short-term efficacy of pharmacotherapeutic interventions in osteoarthritic knee pain: a meta-analysis of randomised placebo-controlled trials.
      • Fuzier R.
      • Serres I.
      • Bourrel R.
      • Palmaro A.
      • Montastruc J.L.
      • Lapeyre-Mestre M.
      Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
      . Selective and non-selective NSAIDS, acetaminophen and narcotic analgesics are the most commonly used analgesics
      • Wilcox C.M.
      • Cryer B.
      • Triadafilopoulos G.
      Patterns of use and public perception of over-the-counter pain relievers: focus on nonsteroidal antiinflammatory drugs.
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      . Despite short-term pain relief, the long-term effect of various analgesics on OA progression and outcome is not clear
      • Bjordal J.M.
      • Klovning A.
      • Ljunggren A.E.
      • Slordal L.
      Short-term efficacy of pharmacotherapeutic interventions in osteoarthritic knee pain: a meta-analysis of randomised placebo-controlled trials.
      • Bijlsma J.W.
      • Berenbaum F.
      • Lafeber F.P.
      Osteoarthritis: an update with relevance for clinical practice.
      . Previous in vitro and animal studies have reported that use of analgesics and NSAIDs has deleterious effects on the structural progression of OA
      • Mastbergen S.C.
      • Jansen N.W.
      • Bijlsma J.W.
      • Lafeber F.P.
      Differential direct effects of cyclo-oxygenase-1/2 inhibition on proteoglycan turnover of human osteoarthritic cartilage: an in vitro study.
      • Gencosmanoglu B.E.
      • Eryavuz M.
      • Dervisoglu S.
      Effects of some nonsteroidal anti-inflammatory drugs on articular cartilage of rats in an experimental model of osteoarthritis.
      . Results of recent observational studies were more controversial. While some reports have confirmed animal studies on negative effects, others have claimed a protective role, especially for the selective COX-II inhibitors
      • Wolfe F.
      • Michaud K.
      • Burke T.A.
      • Zhao S.Z.
      Longer use of COX-2-specific inhibitors compared to nonspecific nonsteroidal antiinflammatory drugs: a longitudinal study of 3639 patients in community practice.
      • Mastbergen S.C.
      • Marijnissen A.C.
      • Vianen M.E.
      • Zoer B.
      • van Roermund P.M.
      • Bijlsma J.W.
      • et al.
      Inhibition of COX-2 by celecoxib in the canine groove model of osteoarthritis.
      • El Hajjaji H.
      • Marcelis A.
      • Devogelaer J.P.
      • Manicourt D.H.
      Celecoxib has a positive effect on the overall metabolism of hyaluronan and proteoglycans in human osteoarthritic cartilage.
      • Deeks J.J.
      • Smith L.A.
      • Bradley M.D.
      Efficacy, tolerability, and upper gastrointestinal safety of celecoxib for treatment of osteoarthritis and rheumatoid arthritis: systematic review of randomised controlled trials.
      .

      Recent studies and pre-specified hypothesis

      Recently, Lapane et al. evaluated the use of NSAIDs on symptoms and OA progression among subjects with knee OA
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      . Interestingly, long-term use of NSAIDs was associated with decreased joint space narrowing (JSN). In contrast, another first line analysis of the Osteoarthritis Initiative (OAI) cohort by Pelletier et al. reported a lower joint space width among analgesic + subjects
      • Martel-Pelletier J.
      • Roubille C.
      • Abram F.
      • Hochberg M.C.
      • Dorais M.
      • Delorme P.
      • et al.
      First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
      . In addition, patients with OA have great variations in their use of analgesics
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      • Conaghan P.G.
      • Peloso P.M.
      • Everett S.V.
      • Rajagopalan S.
      • Black C.M.
      • Mavros P.
      • et al.
      Inadequate pain relief and large functional loss among patients with knee osteoarthritis: evidence from a prospective multinational longitudinal study of osteoarthritis real-world therapies.
      . Many are using more than one analgesic at a time and switches are common among common types of analgesics
      • Fuzier R.
      • Serres I.
      • Bourrel R.
      • Palmaro A.
      • Montastruc J.L.
      • Lapeyre-Mestre M.
      Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      .

      Objectives

      In this study, our objective was to investigate the long-term use of analgesics in association with radiographic progression of OA and subsequent KR as the ultimate structural outcome of knee OA.

      Methods

      Study design, setting and participants

      For this analysis, we used the data from the OAI. The OAI is a cohort of 4796 subjects with or at risk of OA, which are regularly followed on an annual basis. Datasets are available for public access through the OAI website at the http://www.oai.ucsf.edu. Details of the study participants, enrollments, evaluations and follow-ups are also available for public access from the OAI website. The primary hypothesis of the study was to find the maximum number of the index (Analgesic +) and referent (Analgesics −) subjects which have the most similar follow-up assessments except their exposure to analgesics.

      Exposure: long-term use of analgesics

      Our definition of Analgesic + was documented use of analgesics in all available follow-ups. Physician-prescribed medications were adjudicated and were present in the OAI datasets. Following a recently published study, we used prescription analgesics as it's more reflective of its regular use especially in clinically relevant doses (compared to non-prescription analgesics)
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      . For comparison, we aimed at selecting our referent subjects from those with no report of prescription analgesics in neither of their follow-up visits. Thus, we checked the Medication Inventory Form (MIF) of each subject. Out of the total 4796 participants, 2654 subjects had analgesics in their MIF in “some” of the visits. The remaining 2142 subjects had analgesics in their MIF in “all”, or “none” of the follow-up visits.

      Outcomes: radiographic Kellgren and Lawrence (KL) grade and KR

      We had two distinct outcomes of interest. First, increase in the follow-up KL grade compared to the matched baseline KL grade (only between grade changes) as a measure of radiographic OA progression
      • Gossec L.
      • Jordan J.M.
      • Lam M.A.
      • Fang F.
      • Renner J.B.
      • Davis A.
      • et al.
      Comparative evaluation of three semi-quantitative radiographic grading techniques for hip osteoarthritis in terms of validity and reproducibility in 1404 radiographs: report of the OARSI-OMERACT Task Force.
      . Second, KR was chosen as the ultimate outcome of OA
      • Raynauld J.P.
      • Martel-Pelletier J.
      • Haraoui B.
      • Choquette D.
      • Dorais M.
      • Wildi L.M.
      • et al.
      Risk factors predictive of joint replacement in a 2-year multicentre clinical trial in knee osteoarthritis using MRI: results from over 6 years of observation.
      . Following the previous study
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      , the increase in KL grade was evaluated up to the third year of follow-up. The higher KL grade of the two knees was selected as the variable of interest for each participant
      • Barr A.J.
      • Dube B.
      • Hensor E.M.
      • Kingsbury S.R.
      • Peat G.
      • Bowes M.A.
      • et al.
      The relationship between clinical characteristics, radiographic osteoarthritis and 3D bone area: data from the Osteoarthritis Initiative.
      . KR incidence was assessed throughout the 8 years of follow-up using the latest updated datasets that were present online (Datasets were accessed at November 2014). For the purpose of consistency and to avoid issues that may arise from inconsistent or unstandardized readings, we used OAI's central readings of X-rays for KL grade.

      Matching variables

      In order to define matched referents, for every index, we found one matched referent using nearest neighbor matching of the propensity scores with calipers. For every index, we selected the one best referent who matched with the index in not only the baseline, but the follow-up variables as well (except radiographic KL grade, which was matched only for the initial baseline values). Therefore, we attempted to find subjects with a similar course of OA while having different exposures (long term use of analgesic vs no analgesics: Analgesic +/−). Variables that could possibly affect the probability of using analgesics were considered in finding the matching referents
      • Brookhart M.A.
      • Schneeweiss S.
      • Rothman K.J.
      • Glynn R.J.
      • Avorn J.
      • Sturmer T.
      Variable selection for propensity score models.
      . Choice of the possible confounding variables was based on previous OA risk prediction models which were validated on the OAI cohort's population
      • Raynauld J.P.
      • Martel-Pelletier J.
      • Haraoui B.
      • Choquette D.
      • Dorais M.
      • Wildi L.M.
      • et al.
      Risk factors predictive of joint replacement in a 2-year multicentre clinical trial in knee osteoarthritis using MRI: results from over 6 years of observation.
      • Zhang W.
      • McWilliams D.F.
      • Ingham S.L.
      • Doherty S.A.
      • Muthuri S.
      • Muir K.R.
      • et al.
      Nottingham knee osteoarthritis risk prediction models.
      • Riddle D.L.
      • Kong X.
      • Jiranek W.A.
      Factors associated with rapid progression to knee arthroplasty: complete analysis of three-year data from the osteoarthritis initiative.
      • Kerkhof H.J.
      • Bierma-Zeinstra S.M.
      • Arden N.K.
      • Metrustry S.
      • Castano-Betancourt M.
      • Hart D.J.
      • et al.
      Prediction model for knee osteoarthritis incidence, including clinical, genetic and biochemical risk factors.
      . Demographic variables, body mass index (BMI), Western Ontario and McMaster Questionnaire (WOMAC) total score, Physical Summary Scale (PSS) for the Medical Outcome Study (MOS) 12 item short form (SF-12) health survey, mental summary scale (MSS) for the MOS 12 item short form (SF-12) health survey, Physical Activity Scale for the Elderly (PASE) score, and Charlson comorbidity score included the matching variables. Index and referent subjects (Analgesic + and Analgesic −) were selected to have initial matched KL grades, to compare follow-up changes in radiographic KL grades. The greater WOMAC score of the two knees was selected as the variable of interest for each participant
      • Baum T.
      • Joseph G.B.
      • Arulanandan A.
      • Nardo L.
      • Virayavanich W.
      • Carballido-Gamio J.
      • et al.
      Association of magnetic resonance imaging-based knee cartilage T2 measurements and focal knee lesions with knee pain: data from the Osteoarthritis Initiative.
      . After excluding 384 subjects with missing values in either the demographics, KL grades or the Charleston comorbidity scales, 1758 subjects remained. No participant had missing values in neither the exposure nor the outcome. Following exclusion of the above incomplete index cases, 125 of the remaining 61,405 variable-values, which were accounting for 0.002 of all values, were missing. These 125 missing values were imputed using expectation maximization (EM) technique, which provided us with one completely imputed dataset (no remaining missing values)
      • Dong Y.
      • Peng C.Y.
      Principled missing data methods for researchers.
      • Greenland S.
      • Finkle W.D.
      A critical look at methods for handling missing covariates in epidemiologic regression analyses.
      .

      Propensity score matching

      The imputed dataset was used for the purpose of matching. In order to determine referent subjects (Analgesic −) for each index (Analgesics+), we searched the dataset for the one nearest neighbor, non-replaced match
      • Austin P.C.
      A comparison of 12 algorithms for matching on the propensity score.
      . The index (Analgesic +) and referent (Analgesic −) subjects had analgesics in all and none of their visits, respectively. To select the best available referents, propensity score matching was utilized. A 1:1 greedy matching algorithm was employed similar to the method described by Rosenbaum et al.
      • Rosenbaum P.R.
      • Rubin D.B.
      The bias due to incomplete matching.
      . For every participant in the index (Analgesic +), one best-matched referent (Analgesic −) participant was selected. Using greedy matching algorithm, we found the first best match pair and then moved to find the next best match in a hierarchical manner. Best match was defined as subjects with the highest digit match on their propensity scores, which were calculated using logistic regression analyzes to predict the probability of exposure (long-term use of analgesics). A caliper of 0.1 was employed in finding the highest digit match of the propensity scores. As mentioned earlier, out of the total 4796 participants, 2654 subjects had analgesics in their MIF in “some” of the visits. The remaining 2142 subjects had analgesics in their MIF in “all”, or “none” of the follow-up visits; 384 subjects were excluded due to having missing values in either the demographics, KL grades or the Charleston comorbidity scales, and 1758 subjects remained. By exclusion of subjects outside common support of both groups and by considering the caliper, of 187 index subjects (Analgesic +) 14 index cases were excluded and best matches were selected for 173 Analgesic + index subjects, resulting in 173 Analgesic − referents. This method of propensity matching and use of calipers has been extensively reviewed by Austin et al.
      • Austin P.C.
      A comparison of 12 algorithms for matching on the propensity score.
      and is demonstrated to induce minimal, if any, bias for the comparisons. Fig. 1 shows the inclusion of subjects from the 4796 OAI participants. Table I demonstrates the characteristics of the index and referent subjects (Analgesic + and Analgesic −). There was no significant difference in any of the included variables in the matching process between the two groups neither in the average, standardized average or average ranks.
      Figure thumbnail gr1
      Fig. 1Flowchart outlining inclusion of the study subjects.
      Table ICharacteristics of the Analgesic +/− cohorts as means and standard differences
      VariableAnalgesic +Analgesic −Standard mean differenceP value
      Demographics
       Age61.1 ± 8.7461.0 ± 8.62.0030.98
       Gender (Female)0.68 ± 0.470.74 ± 0.44−.1380.19
       History of knee injury0.40 ± 0.490.44 ± 0.50−.0700.51
       Family history of KR0.15 ± 0.360.12 ± 0.32.0980.34
       Yearly income category3.73 ± 1.093.76 ± 1.01−.0270.80
      BMI
       Screening BMI28.9 ± 4.8128.6 ± 4.59.0670.50
       Twelfth month BMI29.0 ± 4.7428.8 ± 4.59.0480.64
       Twenty-fourth month BMI29.1 ± 5.0728.7 ± 4.59.0680.50
       Thirty-sixth month BMI29.1 ± 5.2328.6 ± 4.52.0890.37
      WOMAC total score (Reference range: 0–96)
       Enrollment WOMAC17.6 ± 16.517.8 ± 16.2−.0120.91
       Twelfth month WOMAC17.9 ± 16.416.9 ± 17.3.0570.58
       Twenty-fourth month WOMAC17.6 ± 17.016.8 ± 16.4.0410.68
       Thirty-sixth month WOMAC18.2 ± 16.217.3 ± 16.3.0550.58
      PSS for the MOS 12 item short form (SF-12) health survey (Reference range: 0–100)
       Enrollment PSS46.6 ± 9.3647.5 ± 8.87−.0890.38
       Twelfth month PSS45.3 ± 9.2145.9 ± 9.82−.0640.54
       Twenty-fourth month PSS45.2 ± 10.346.1 ± 9.52−.0790.43
       Thirty-sixth month PSS44.4 ± 9.5545.1 ± 9.90−.0650.53
      MSS for the MOS 12 item short form (SF-12) health survey (Reference range: 0–100)
       Enrollment MSS53.3 ± 9.7153.3 ± 8.19−.0070.94
       Twelfth month MSS53.2 ± 8.7953.9 ± 8.51−.0850.39
       Twenty-fourth month MSS53.0 ± 9.4353.4 ± 8.83−.0370.72
       Thirty-sixth month MSS53.3 ± 9.7353.5 ± 8.57−.0190.84
      PASE score (Reference range: 0–361)
       Enrollment PASE158 ± 78.3164 ± 83.4−.0740.50
       Twelfth month PASE159 ± 81.2168 ± 87.7−.1000.34
       Twenty-fourth month PASE145 ± 78.4149 ± 78.1−.0480.66
       Thirty-sixth month PASE143 ± 75.5146 ± 75.0−.0390.71
      Charlson Comorbidity Score (Weighted scores range: 1–6)
       Enrollment Comorbidity score0.40 ± 0.700.35 ± 0.84.0720.53
       Twenty-fourth month Comorbidity score0.50 ± 0.830.50 ± 1.10.0001.00
      Radiographic KL Score (Reference range: 0–4)
       Enrollment KL grade1.82 ± 1.251.92 ± 1.21−.0850.43
       Overall Propensity Score0.20 ± 0.160.20 ± 0.16.0110.90
      Yearly income categories: 1. Less than $10K; 2. $10K to <$25K; 3. $25K to <$50K; 4. $50K to <$100K; 5. $100K or greater.

      Comparison of index and referent subjects (Analgesic + and Analgesic −)

      Cumulative numbers of subjects with an increase in KL grade and/or KRs were compared between the two groups using Chi-square (X2). Fisher's exact test was used whenever the expected number of events was less than five in a two by two table. For a clear description of the number of participants with OA progression throughout the follow-up, year by year analyzes were performed. Number of subjects with OA progression at each time point were calculated from the total of 173 baseline subjects. Cox regression analysis was employed to extract the Hazard Ratio (HR) of increase in KL grade (OA progression) and KR associated with long-term use of analgesics. Additional covariate adjustment was performed using the previously obtained propensity scores
      • Austin P.C.
      • Grootendorst P.
      • Normand S.L.
      • Anderson G.M.
      Conditioning on the propensity score can result in biased estimation of common measures of treatment effect: a Monte Carlo study.
      • Austin P.C.
      Double propensity-score adjustment: a solution to design bias or bias due to incomplete matching.
      . This approach has shown to derivate the results toward null hypothesis in selected reports
      • Austin P.C.
      • Grootendorst P.
      • Normand S.L.
      • Anderson G.M.
      Conditioning on the propensity score can result in biased estimation of common measures of treatment effect: a Monte Carlo study.
      • Hade E.M.
      • Lu B.
      Bias associated with using the estimated propensity score as a regression covariate.
      . However, we took a conservative stance (avoiding type I error) and calculated the adjusted HRs as well. A two-tailed P value of less than 5% was considered significant. Analyzes were performed using SPSS (v.18, Chicago, IL) and the R platform (v. 2.8.1). We note that the aim of this investigation was to examine an association (rather than depicting a cause and effect graph) while dealing with the imbalance in a series of known determinants; using a relatively matched sample of subjects, in the OAI cohort.

      Results

      Participants and descriptive data

      The mean age was 61.1 ± 8.74 and 61.0 ± 8.62 among the Analgesic + and Analgesic − subjects, respectively. Sixty-eight percent of Analgesic + and seventy-four percent of Analgesic − subjects were females. Table I demonstrates the standardized mean difference of the matching variables between Analgesic + and Analgesic − subjects. The standardized mean difference is used to represent the size of each variable's difference, relative to its variability, between Analgesic + and Analgesic − groups. The standardized mean difference is calculated by dividing the mean difference in each variable's value by the standard deviation of that variable in the studied group. There was no statistical difference in the demographics, BMI, WOMAC score, PSS, MSS and comorbidity scales of the two groups; neither at the baseline, first year, second year or third year follow-ups (Table I). The mean KL grades at the baseline radiographs were 1.82 ± 1.25 and 1.92 ± 1.21 for Analgesic + and Analgesic − subjects, respectively (P value = 0.43) (Table I). There was no significant difference in the distribution of the subjects among different KL grades (P values > 0.05). For both the Analgesic + and Analgesic − subjects, median and inter-quartile ranges of KL grades were 2 (1–3). The mean KL grade at the third year of follow-up was 1.94 ± 1.29 and 2.00 ± 1.24 for Analgesic + and Analgesic − subjects, respectively (P-value: 0.61). Table II demonstrates the details about the type of analgesics as reported in the follow-up reports. As it is evident in Table II, many Analgesic + subjects receive more than one type of analgesics. Moreover, patients frequently change the type of their analgesics in long follow-ups.
      Table IIDetails of the prescription analgesics during the years of follow-up
      BaselineTwelfth monthTwenty-fourth monthThirty-sixth monthForty-eighth monthSeventy-second month
      Type of analgesics
       NSAID839897969388
       COX-II inhibitor594337353635
       Acetaminophen262529273028
       Salicylate121212969
       Narcotics373446465253
       Other analgesics182221204846
      Number of analgesics
       1128129124127123124
       2293031332023
       3151116121816
       41321129
       5000001

      Outcome data and main results

      At the end of third year, 16.76% (29/173) of Analgesic + subjects, and 8.67% (15/173) of the Analgesic − subjects had OA progression in terms of increased KL grades (P value = 0.024). By the eighth year of follow-up, 16.76% (29/173) and 6.94% (12/173) subjects had KR among the Analgesic + and Analgesic − subjects, respectively (P value = 0.005). Long-term use of analgesics was associated with 1.91 (Adjusted HR 95% CI: 1.02–3.57) and 2.57 (Adjusted HR 95% CI: 1.31–5.04) times greater risk of OA progression and KR (Table III, Table IV). Long-term use of analgesics was associated with significantly higher risk of KR even when controlled for the effect of increase in the KL grade in the first 3 years (HR: 2.30 (1.16–4.54); P-value: 0.017). In other words, use of analgesics was associated with KR among subjects with and without KL grade increase in the first 3 years.
      Table IIICumulative number of subjects with increase in the radiographic KL scores and KR according to the study group (Analgesic + vs Analgesic −)
      Analgesic + (n = 173)Analgesic − (n = 173)P value
      Increase in KL grade
       Baseline
      Analgesic + and Analgesic − subjects were matched for their baseline KL scores.
      N/A
      NA: Not applicable.
      N/A
       First year1390.378
       Second year18110.174
       Third year29150.024
      KR
       Baseline00
       First year410.216
       Second year830.219
       Third year1340.044
       Fourth year1860.011
       Fifth year2090.033
       Sixth year27100.003
       Seventh year29110.002
       Eighth year29120.005
      Analgesic + and Analgesic − subjects were matched for their baseline KL scores.
      NA: Not applicable.
      Table IVRisk of increase in the KL score and KR according to the study group (Analgesic + vs Analgesic −)
      Crude HRP valueAdjusted HRP value
      Risk of Increase in KL grade (3 years)
      Analgesic + vs Analgesic −1.93

      1.04–3.61
      0.0391.93

      1.03–3.60
      0.038
      Risk of KR (8 years)
      Analgesic + vs Analgesic −2.54

      1.29–4.97
      0.0072.57

      1.31–5.04
      0.006
      Adjustment is performed for the derived propensity score.
      We performed an additional analysis of the recently released eighth year follow-up readings to further test our hypothesis. In the recently updated OAI project, reports of the KL grading on 264 (out of 346; 76.3% of included subjects) subjects were not available for the eighth year of follow-up. Thus, we continued our analysis knowing that our study power may be diminished due to small sample size. However, results of the available eighth year KL grades showed that 16 (out of 43) subjects in the Analgesics + and seven (out of 39) subjects in the Analgesic − group had OA progression (increase in KL grade compared to the baseline radiography). The comparison was significant with a P value (using Chi–square test) of 0.05. Other determinants had no significant difference among the Analgesic + and − group: Age (58.35 ± 8.44 vs 57.23 ± 7.63; P-value: 0.53), Gender (72% vs 85%; P value: 0.17), BMI (27.06 ± 4.44 vs 26.39 ± 3.86; P-value: 0.47), PASE (164.91 ± 73.74 vs 166.63 ± 77.91; P-value: 0.92), WOMAC (9.32 ± 9.71 vs 12.16 ± 15.10; P-value: 0.32) and baseline KL grade (0.26 ± 0.44 vs 0.33 ± 0.48; P-value: 0.45).
      We also investigated the OAI cohort for the availability of data on the use of steroid injections in the symptomatic knees. Such information is collected on a yearly basis in the annual follow-ups (from baseline to the seventh year follow-up). Subjects are asked about steroid injections within the last 6 months. In the whole OAI cohort (n = 4796), only 1.94% (baseline) to 4.45% (seventh year – the last available data on steroid injections) of subjects had reports of using steroid injections. In this study (n = 346), number of subjects who had steroid injections was 3.2% (baseline) and 6.4% (seventh year – the last available data on steroid injection), respectively; which was not statistically significant among Analgesic + and Analgesic − subjects (P > 0.05). Moreover, WOMAC scores for pain, stiffness, and disability were analyzed separately to investigate whether Analgesic + and Analgesic − subjects have significantly different sub-scores, or not. Our analyzes revealed that none of the sub-scores had neither statistically (all P values > 0.10) nor clinically significant difference among Analgesic + and Analgesic − subjects (Supplementary Table 1).

      Discussion

      Key results

      In this study, our findings suggest that the increased risk of KR in the index (Analgesic +) subjects may be at least partly explained by OA progression. Subjects with long-term use of analgesics were more likely to undergo KR in comparison to subjects with similar OA severity (OA severity was similar between the Analgesic + and Analgesic − subjects at the baseline). The most important discussion point is that the indication for KR is made clinically and is not only based on radiographs. If the patient has knee OA but is clinically asymptomatic, he/she does not receive joint replacement while a patient with persistent pain and analgesics consumption will receive joint replacement much earlier. Interestingly, other studies show that prescription analgesic consumption may even increase after the KR surgery
      • Fuzier R.
      • Serres I.
      • Bourrel R.
      • Palmaro A.
      • Montastruc J.L.
      • Lapeyre-Mestre M.
      Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
      . These findings suggest new determinants of analgesic consumption, regardless of OA severity and progression. In addition to their analgesic effects, different analgesics might be associated with different pathophysiologic effects. Local anti-inflammatory effects (in case of ibuprofen)
      • Adatia A.
      • Rainsford K.D.
      • Kean W.F.
      Osteoarthritis of the knee and hip. Part II: therapy with ibuprofen and a review of clinical trials.
      , favorable and unfavorable morphologic changes in cartilage tissues (in case of tiaprofenic and diclofenac)
      • Gencosmanoglu B.E.
      • Eryavuz M.
      • Dervisoglu S.
      Effects of some nonsteroidal anti-inflammatory drugs on articular cartilage of rats in an experimental model of osteoarthritis.
      and effects on hyaluronan, proteoglycans and their metabolism (in case of celecoxib)
      • El Hajjaji H.
      • Marcelis A.
      • Devogelaer J.P.
      • Manicourt D.H.
      Celecoxib has a positive effect on the overall metabolism of hyaluronan and proteoglycans in human osteoarthritic cartilage.
      , are of the most investigated mechanisms. Besides pathophysiologic explanations, exploring the differences in the level of satisfaction and expectations of patients in decision making for KR may potentially explain the differences between Analgesic + and Analgesic − subjects
      • Becker R.
      • Doring C.
      • Denecke A.
      • Brosz M.
      Expectation, satisfaction and clinical outcome of patients after total knee arthroplasty.
      . Furthermore, especially with regard to the incidence of KR, economic discrepancies and impact of healthcare access, information and education are of the deterministic factors in the management of OA; as expressed by the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO)
      • Bruyere O.
      • Cooper C.
      • Pelletier J.P.
      • Branco J.
      • Luisa Brandi M.
      • Guillemin F.
      • et al.
      An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
      .
      Our results are similar to previous reports that highlighted the deleterious effects of analgesics, particularly non-selective NSAIDs, in the course of OA
      • Ding C.
      Do NSAIDs affect the progression of osteoarthritis?.
      • Reijman M.
      • Bierma-Zeinstra S.M.
      • Pols H.A.
      • Koes B.W.
      • Stricker B.H.
      • Hazes J.M.
      Is there an association between the use of different types of nonsteroidal antiinflammatory drugs and radiologic progression of osteoarthritis? The Rotterdam Study.
      • Fendrick A.M.
      • Greenberg B.P.
      A review of the benefits and risks of nonsteroidal anti-inflammatory drugs in the management of mild-to-moderate osteoarthritis.
      • Pavelka K.
      A comparison of the therapeutic efficacy of diclofenac in osteoarthritis: a systematic review of randomised controlled trials.
      • Gencosmanoglu B.E.
      • Eryavuz M.
      • Dervisoglu S.
      Effects of some nonsteroidal anti-inflammatory drugs on articular cartilage of rats in an experimental model of osteoarthritis.
      . In this study, Analgesic + subjects with similar WOMAC scores (including pain, stiffness and disability) and initial radiographic KL grade had a higher risk of OA progression and KR. Recently, Lapane et al. evaluated the use of NSAIDs on symptoms and OA progression among subjects with knee OA
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      . They analyzed the OAI progression cohort using marginal structural modeling to control for the time variant variables including pain, body weight, physical and mental health status. Interestingly, long-term use of NSAIDs was associated with decreased JSN, albeit not statistically significant. In contrast, another first line analysis of the OAI cohort by Pelletier et al. reported a lower joint space width among analgesic + subjects
      • Martel-Pelletier J.
      • Roubille C.
      • Abram F.
      • Hochberg M.C.
      • Dorais M.
      • Delorme P.
      • et al.
      First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
      . In addition, patients with OA have great variations in their use of analgesics
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      • Lapane K.L.
      • Yang S.
      • Driban J.B.
      • Liu S.H.
      • Dube C.E.
      • McAlindon T.E.
      • et al.
      Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
      • Conaghan P.G.
      • Peloso P.M.
      • Everett S.V.
      • Rajagopalan S.
      • Black C.M.
      • Mavros P.
      • et al.
      Inadequate pain relief and large functional loss among patients with knee osteoarthritis: evidence from a prospective multinational longitudinal study of osteoarthritis real-world therapies.
      and may use or switch to more than one analgesic at a time
      • Fuzier R.
      • Serres I.
      • Bourrel R.
      • Palmaro A.
      • Montastruc J.L.
      • Lapeyre-Mestre M.
      Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      . Compared to previous studies, we investigated the effect of the long-term analgesic use on changes in follow-up radiographic KL grade measurements as an objective metric for OA progression. Subjects included in our study had NSAIDs among their prescriptions for the majority of the follow-up period as well. Certain differences in the choice of exposure (Analgesics in general vs class specific drug effects; long-term exposure vs short term exposure) study design (only including subjects with balanced symptoms vs including subjects with improved symptoms) and outcome (KR as a clinical decision based on patients' preference; radiographic KL grade that includes joint space, sclerosis, osteophytes and joint deformity vs JSN) contribute to the different findings among different studies.
      We observed a meaningful difference, between Analgesic + and − subjects, for OA progression and KR. Likewise, a previous report of the Rotterdam Study indicated an increased risk of OA progression in subjects who were using diclofenac for more than 180 days
      • Reijman M.
      • Bierma-Zeinstra S.M.
      • Pols H.A.
      • Koes B.W.
      • Stricker B.H.
      • Hazes J.M.
      Is there an association between the use of different types of nonsteroidal antiinflammatory drugs and radiologic progression of osteoarthritis? The Rotterdam Study.
      . Our results may further expand this observation. Even with similar pain, stiffness and activity profiles, those who achieve these profiles with the aid of analgesics, accompanied an increased risk of OA progression in the OAI cohort; whereas, structural joint damages were also shown to predict the risk of future KR
      • Roemer F.W.
      • Kwoh C.K.
      • Hannon M.J.
      • Hunter D.J.
      • Eckstein F.
      • Wang Z.
      • et al.
      Can structural joint damage measured with MR imaging be used to predict knee replacement in the following year?.
      . In other words, having similar WOMAC score while using analgesics may further contribute to the difference between the Analgesic + and Analgesic − subjects. Aside from the theoretical independent contribution of the analgesic drugs to progression of OA, our findings may also speculate the potential role of pain perception in protection from hazardous positions, activities and weight bearings
      • Stehling C.
      • Lane N.E.
      • Nevitt M.C.
      • Lynch J.
      • McCulloch C.E.
      • Link T.M.
      Subjects with higher physical activity levels have more severe focal knee lesions diagnosed with 3T MRI: analysis of a non-symptomatic cohort of the osteoarthritis initiative.
      • Mansournia M.A.
      • Danaei G.
      • Forouzanfar M.H.
      • Mahmoodi M.
      • Jamali M.
      • Mansournia N.
      • et al.
      Effect of physical activity on functional performance and knee pain in patients with osteoarthritis : analysis with marginal structural models.
      that can accelerate OA progression. Despite having balanced Analgesic + and − populations, this study does not explain why some of the subjects were prescribed NSAIDs over the others. We could show that the Analgesic + and − groups were balanced in their symptoms (pain, stiffness and disability). Whether index subjects were less tolerant to pain than referents and whether that was the reason why they chose KR more often, is not elucidated in this study. Even beyond OA progression, Analgesic + subjects were more frequent candidates for KR; further suggesting their analgesic seeking behavior and their need to accept new alternatives interventions in addition to analgesics for the treatment of painful OA. In this regard, Collins et al. argued that subjects with symptomatic knee OA (KL grade ≥ 2) had little changes in their pain; during their 6-year follow-up study of the OAI cohort. Depression, medical comorbidities, and lower education (among others) were suggested to be associated with pain trajectories
      • Collins J.E.
      • Katz J.N.
      • Dervan E.E.
      • Losina E.
      Trajectories and risk profiles of pain in persons with radiographic, symptomatic knee osteoarthritis: data from the osteoarthritis initiative.
      . On the other hand, continued use of conventional analgesics (such as NSAIDs) is not necessarily associated with complete pain relief either
      • Conaghan P.G.
      • Peloso P.M.
      • Everett S.V.
      • Rajagopalan S.
      • Black C.M.
      • Mavros P.
      • et al.
      Inadequate pain relief and large functional loss among patients with knee osteoarthritis: evidence from a prospective multinational longitudinal study of osteoarthritis real-world therapies.
      . While a great portion of OA patients are frequent users of analgesics
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      still many patients suffer from inadequate pain relief (IPR), which is also a common reason for impaired quality of life (QoL)
      • Conaghan P.G.
      • Peloso P.M.
      • Everett S.V.
      • Rajagopalan S.
      • Black C.M.
      • Mavros P.
      • et al.
      Inadequate pain relief and large functional loss among patients with knee osteoarthritis: evidence from a prospective multinational longitudinal study of osteoarthritis real-world therapies.
      .

      Limitations

      This study had several limitations. Because of the number of available cases we could not draw meaningful comparisons for different types of analgesics. Table II demonstrates the details of analgesics used in each follow-up. The same is true for the number of the analgesics that were used by each participant at each time point (each follow-up) and during the entire follow-up period. Our study design made it impossible to investigate the number of analgesics in relation to OA progression and KR. Future investigations, using a larger number of Analgesic + subjects, should be designed to address this issue. Considering the WHO three ladder categorization of analgesics, we used all three categories in our definition of analgesic use
      • Bruyere O.
      • Cooper C.
      • Pelletier J.P.
      • Branco J.
      • Luisa Brandi M.
      • Guillemin F.
      • et al.
      An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
      • Hochberg M.C.
      • Altman R.D.
      • April K.T.
      • Benkhalti M.
      • Guyatt G.
      • McGowan J.
      • et al.
      American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee.
      . While different analgesics may have heterogeneous anti-inflammatory and disease-modifying effects, they all share the common property of reducing pain perception
      • Bruyere O.
      • Cooper C.
      • Pelletier J.P.
      • Branco J.
      • Luisa Brandi M.
      • Guillemin F.
      • et al.
      An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
      • Hochberg M.C.
      • Altman R.D.
      • April K.T.
      • Benkhalti M.
      • Guyatt G.
      • McGowan J.
      • et al.
      American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee.
      • Martel-Pelletier J.
      • Roubille C.
      • Abram F.
      • Hochberg M.C.
      • Dorais M.
      • Delorme P.
      • et al.
      First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
      • Wilcox C.M.
      • Cryer B.
      • Triadafilopoulos G.
      Patterns of use and public perception of over-the-counter pain relievers: focus on nonsteroidal antiinflammatory drugs.
      • Bjordal J.M.
      • Klovning A.
      • Ljunggren A.E.
      • Slordal L.
      Short-term efficacy of pharmacotherapeutic interventions in osteoarthritic knee pain: a meta-analysis of randomised placebo-controlled trials.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      • Trijau S.
      • Avouac J.
      • Escalas C.
      • Gossec L.
      • Dougados M.
      Influence of flare design on symptomatic efficacy of non-steroidal anti-inflammatory drugs in osteoarthritis: a meta-analysis of randomized placebo-controlled trials.
      . In addition to the common use of analgesic medications in OA patients, there is a large variation in physician's recommendations and patients' consistency in their choice of medications
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      . Moreover, many subjects use multiple analgesics simultaneously and commonly switch between the types of analgesics during the course of symptomatic OA
      • Fuzier R.
      • Serres I.
      • Bourrel R.
      • Palmaro A.
      • Montastruc J.L.
      • Lapeyre-Mestre M.
      Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
      • Maserejian N.N.
      • Fischer M.A.
      • Trachtenberg F.L.
      • Yu J.
      • Marceau L.D.
      • McKinlay J.B.
      • et al.
      Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
      • Kingsbury S.R.
      • Hensor E.M.
      • Walsh C.A.
      • Hochberg M.C.
      • Conaghan P.G.
      How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
      • Gore M.
      • Sadosky A.B.
      • Leslie D.L.
      • Tai K.S.
      • Emery P.
      Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
      . For all the aforementioned reasons, it was not practical to find adequate samples of subjects with the long-term use of any specific analgesic who also share a similar course of OA. Moreover, this study compares the extremes ends (no use vs persistent and long term use) of a spectrum of subjects using analgesic medications. These findings may not be achieved in subjects with minimal or short-term use of analgesics. We believe that future studies should reveal whether there is a significant difference among different analgesics for OA progression owing to their differential actions. Another limitation is the lack of analgesic seeking behavioral metrics in patients who are on analgesics. Finally, we note that especial caution should be made when applying our conclusions to other, possibly community dwelling populations. The analgesic + population in our study had a history of long-term analgesic use. Meanwhile, many other subjects in the OAI cohort may have used analgesics only in some of the follow-up years [Fig. 1], and thus were excluded from this analysis. In fact, in a trade off, we achieved finely tuned comparisons by assuming a condition that subjects had used analgesics in all of their follow-up visits.
      Considering recent studies elucidating the effect of different analgesics on different aspects of OA progression warrants further investigations using advanced imaging modalities (including Magnetic Resonance Imaging – MRI) for more sensitive and specific detection of the early signs of knee OA. Larger studies may unveil the specific KL grading component that may be affected in Analgesic + subjects (sclerosis vs JSN vs osteophytes). Our next planned study is to investigate the MRI findings of Analgesic + and Analgesic − subjects, looking for signs of subtle cartilage damage, bone marrow lesions, synovitis and other structural knee damages that may not be apparent in knee radiographs. The ultimate question is whether individuals with more severe degenerative changes on MRI take more analgesics or (vice versa) that analgesics lead to progression of OA (the direction we proposed in this investigation).

      Conclusions

      In conclusion, our results indicate that, after controlling the effect of potential time-varying confounders, long-term use of analgesics can be associated with radiographic progression of knee OA and increased risk of future KR.

      Authors' contributions

      All authors have made a substantial contribution in the conception, design or acquisition of data, interpretation of data, and critically reviewing the draft of the article and approving the final version. In addition, NHN had a role in the design, analysis and drafting the manuscript. AG and FWR had an additional role in the interpretation of data, critical review of the manuscript and provision of study materials. JE had an additional role for statistical revision, intellectual revision of the manuscript and interpretation of the results. BZ and SD had an additional role in the conception of the study, provision of the analyzes and final preparation of the manuscript.

      Conflicts of interest

      Nima Hafezi-Nejad has no conflicts of interest. Ali Guermazi is president and shareholder of Boston Imaging Core Lab, LLC. He is a consultant to Genzyme, MerckSerono, TissueGene and OrthoTrophix. Frank W Roemer is CMO and shareholder of Boston Imaging Core Lab, LLC. John Eng and Bashir Zikria declare no conflicts of interest. Shadpour Demehri has grants from GERRAF 2014–2016; Carestream Health Inc. 2013–2015 for Cone – Beam CT clinical trial. He is a consultant to Toshiba Medical Systems.

      Acknowledgments

      The OAI is a public-private partnership comprised of five contracts (N01-AR-2-2258; N01-AR-2-2259; N01-AR-2-2260; N01-AR-2-2261; N01-AR-2-2262) funded by the National Institutes of Health, a branch of the Department of Health and Human Services, and conducted by the OAI Study Investigators. Private funding partners include Merck Research Laboratories; Novartis Pharmaceuticals Corporation, GlaxoSmithKline; and Pfizer, Inc. Private sector funding for the OAI is managed by the Foundation for the National Institutes of Health. This manuscript was prepared using an OAI public use data set and does not necessarily reflect the opinions or views of the OAI investigators, the NIH, or the private funding partners.

      Appendix A. Supplementary data

      The following is the supplementary data related to this article:

      References

        • Lawrence R.C.
        • Felson D.T.
        • Helmick C.G.
        • Arnold L.M.
        • Choi H.
        • Deyo R.A.
        • et al.
        Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II.
        Arthritis Rheum. 2008; 58: 26-35
        • Felson D.T.
        • Zhang Y.
        • Hannan M.T.
        • Naimark A.
        • Weissman B.N.
        • Aliabadi P.
        • et al.
        The incidence and natural history of knee osteoarthritis in the elderly. The Framingham Osteoarthritis Study.
        Arthritis Rheum. 1995; 38: 1500-1505
        • Raynauld J.P.
        • Martel-Pelletier J.
        • Haraoui B.
        • Choquette D.
        • Dorais M.
        • Wildi L.M.
        • et al.
        Risk factors predictive of joint replacement in a 2-year multicentre clinical trial in knee osteoarthritis using MRI: results from over 6 years of observation.
        Ann Rheum Dis. 2011; 70: 1382-1388
        • Bruyere O.
        • Cooper C.
        • Pelletier J.P.
        • Branco J.
        • Luisa Brandi M.
        • Guillemin F.
        • et al.
        An algorithm recommendation for the management of knee osteoarthritis in Europe and internationally: a report from a task force of the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis (ESCEO).
        Semin Arthritis Rheum. 2014; 44: 253-263
        • McAlindon T.E.
        • Bannuru R.R.
        • Sullivan M.C.
        • Arden N.K.
        • Berenbaum F.
        • Bierma-Zeinstra S.M.
        • et al.
        OARSI guidelines for the non-surgical management of knee osteoarthritis.
        Osteoarthr Cartil. 2014; 22: 363-388
        • Jevsevar D.S.
        • Brown G.A.
        • Jones D.L.
        • Matzkin E.G.
        • Manner P.A.
        • Mooar P.
        • et al.
        The American Academy of Orthopaedic Surgeons evidence-based guideline on: treatment of osteoarthritis of the knee.
        J Bone Jt Surg Am. 2013; 95 (2nd edn.): 1885-1886
        • Hochberg M.C.
        • Altman R.D.
        • April K.T.
        • Benkhalti M.
        • Guyatt G.
        • McGowan J.
        • et al.
        American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee.
        Arthritis Care Res (Hoboken). 2012; 64: 465-474
        • Ding C.
        Do NSAIDs affect the progression of osteoarthritis?.
        Inflammation. 2002; 26: 139-142
        • Reijman M.
        • Bierma-Zeinstra S.M.
        • Pols H.A.
        • Koes B.W.
        • Stricker B.H.
        • Hazes J.M.
        Is there an association between the use of different types of nonsteroidal antiinflammatory drugs and radiologic progression of osteoarthritis? The Rotterdam Study.
        Arthritis Rheum. 2005; 52: 3137-3142
        • Martel-Pelletier J.
        • Roubille C.
        • Abram F.
        • Hochberg M.C.
        • Dorais M.
        • Delorme P.
        • et al.
        First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
        Ann Rheum Dis. 2015 Mar; 74 ([Epub 2013 Dec 13]): 547-556https://doi.org/10.1136/annrheumdis-2013-203906
        • Wolfe F.
        • Michaud K.
        • Burke T.A.
        • Zhao S.Z.
        Longer use of COX-2-specific inhibitors compared to nonspecific nonsteroidal antiinflammatory drugs: a longitudinal study of 3639 patients in community practice.
        J Rheumatol. 2004; 31: 355-358
        • Wilcox C.M.
        • Cryer B.
        • Triadafilopoulos G.
        Patterns of use and public perception of over-the-counter pain relievers: focus on nonsteroidal antiinflammatory drugs.
        J Rheumatol. 2005; 32: 2218-2224
        • Adatia A.
        • Rainsford K.D.
        • Kean W.F.
        Osteoarthritis of the knee and hip. Part II: therapy with ibuprofen and a review of clinical trials.
        J Pharm Pharmacol. 2012; 64: 626-636
        • Towheed T.E.
        • Maxwell L.
        • Judd M.G.
        • Catton M.
        • Hochberg M.C.
        • Wells G.
        Acetaminophen for osteoarthritis.
        Cochrane Database Syst Rev. 2006; : CD004257
        • Fendrick A.M.
        • Greenberg B.P.
        A review of the benefits and risks of nonsteroidal anti-inflammatory drugs in the management of mild-to-moderate osteoarthritis.
        Osteopath Med Prim Care. 2009; 3: 1
        • Pavelka K.
        A comparison of the therapeutic efficacy of diclofenac in osteoarthritis: a systematic review of randomised controlled trials.
        Curr Med Res Opin. 2012; 28: 163-178
        • Bjordal J.M.
        • Klovning A.
        • Ljunggren A.E.
        • Slordal L.
        Short-term efficacy of pharmacotherapeutic interventions in osteoarthritic knee pain: a meta-analysis of randomised placebo-controlled trials.
        Eur J Pain. 2007; 11: 125-138
        • Fuzier R.
        • Serres I.
        • Bourrel R.
        • Palmaro A.
        • Montastruc J.L.
        • Lapeyre-Mestre M.
        Analgesic drug consumption increases after knee arthroplasty: a pharmacoepidemiological study investigating postoperative pain.
        Pain. 2014; 155: 1339-1345
        • Maserejian N.N.
        • Fischer M.A.
        • Trachtenberg F.L.
        • Yu J.
        • Marceau L.D.
        • McKinlay J.B.
        • et al.
        Variations among primary care physicians in exercise advice, imaging, and analgesics for musculoskeletal pain: results from a factorial experiment.
        Arthritis Care Res (Hoboken). 2014; 66: 147-156
        • Kingsbury S.R.
        • Hensor E.M.
        • Walsh C.A.
        • Hochberg M.C.
        • Conaghan P.G.
        How do people with knee osteoarthritis use osteoarthritis pain medications and does this change over time? Data from the Osteoarthritis Initiative.
        Arthritis Res Ther. 2013; 15: R106
        • Gore M.
        • Sadosky A.B.
        • Leslie D.L.
        • Tai K.S.
        • Emery P.
        Therapy switching, augmentation, and discontinuation in patients with osteoarthritis and chronic low back pain.
        Pain Pract. 2012; 12: 457-468
        • Bijlsma J.W.
        • Berenbaum F.
        • Lafeber F.P.
        Osteoarthritis: an update with relevance for clinical practice.
        Lancet. 2011; 377: 2115-2126
        • Mastbergen S.C.
        • Jansen N.W.
        • Bijlsma J.W.
        • Lafeber F.P.
        Differential direct effects of cyclo-oxygenase-1/2 inhibition on proteoglycan turnover of human osteoarthritic cartilage: an in vitro study.
        Arthritis Res Ther. 2006; 8: R2
        • Gencosmanoglu B.E.
        • Eryavuz M.
        • Dervisoglu S.
        Effects of some nonsteroidal anti-inflammatory drugs on articular cartilage of rats in an experimental model of osteoarthritis.
        Res Exp Med (Berl). 2001; 200: 215-226
        • Mastbergen S.C.
        • Marijnissen A.C.
        • Vianen M.E.
        • Zoer B.
        • van Roermund P.M.
        • Bijlsma J.W.
        • et al.
        Inhibition of COX-2 by celecoxib in the canine groove model of osteoarthritis.
        Rheumatology (Oxford). 2006; 45: 405-413
        • El Hajjaji H.
        • Marcelis A.
        • Devogelaer J.P.
        • Manicourt D.H.
        Celecoxib has a positive effect on the overall metabolism of hyaluronan and proteoglycans in human osteoarthritic cartilage.
        J Rheumatol. 2003; 30: 2444-2451
        • Deeks J.J.
        • Smith L.A.
        • Bradley M.D.
        Efficacy, tolerability, and upper gastrointestinal safety of celecoxib for treatment of osteoarthritis and rheumatoid arthritis: systematic review of randomised controlled trials.
        BMJ. 2002; 325: 619
        • Lapane K.L.
        • Yang S.
        • Driban J.B.
        • Liu S.H.
        • Dube C.E.
        • McAlindon T.E.
        • et al.
        Effects of prescription non-steroidal anti-inflammatory agents on symptoms and disease progression among patients with knee osteoarthritis.
        Arthritis Rheumatol. 2014;
        • Martel-Pelletier J.
        • Roubille C.
        • Abram F.
        • Hochberg M.C.
        • Dorais M.
        • Delorme P.
        • et al.
        First-line analysis of the effects of treatment on progression of structural changes in knee osteoarthritis over 24 months: data from the osteoarthritis initiative progression cohort.
        Ann Rheum Dis. 2015; 74: 547-556
        • Conaghan P.G.
        • Peloso P.M.
        • Everett S.V.
        • Rajagopalan S.
        • Black C.M.
        • Mavros P.
        • et al.
        Inadequate pain relief and large functional loss among patients with knee osteoarthritis: evidence from a prospective multinational longitudinal study of osteoarthritis real-world therapies.
        Rheumatology (Oxford). 2014;
        • Gossec L.
        • Jordan J.M.
        • Lam M.A.
        • Fang F.
        • Renner J.B.
        • Davis A.
        • et al.
        Comparative evaluation of three semi-quantitative radiographic grading techniques for hip osteoarthritis in terms of validity and reproducibility in 1404 radiographs: report of the OARSI-OMERACT Task Force.
        Osteoarthr Cartil. 2009; 17: 182-187
        • Barr A.J.
        • Dube B.
        • Hensor E.M.
        • Kingsbury S.R.
        • Peat G.
        • Bowes M.A.
        • et al.
        The relationship between clinical characteristics, radiographic osteoarthritis and 3D bone area: data from the Osteoarthritis Initiative.
        Osteoarthr Cartil. 2014; 22: 1703-1709
        • Brookhart M.A.
        • Schneeweiss S.
        • Rothman K.J.
        • Glynn R.J.
        • Avorn J.
        • Sturmer T.
        Variable selection for propensity score models.
        Am J Epidemiol. 2006; 163: 1149-1156
        • Zhang W.
        • McWilliams D.F.
        • Ingham S.L.
        • Doherty S.A.
        • Muthuri S.
        • Muir K.R.
        • et al.
        Nottingham knee osteoarthritis risk prediction models.
        Ann Rheum Dis. 2011; 70: 1599-1604
        • Riddle D.L.
        • Kong X.
        • Jiranek W.A.
        Factors associated with rapid progression to knee arthroplasty: complete analysis of three-year data from the osteoarthritis initiative.
        Jt Bone Spine. 2012; 79: 298-303
        • Kerkhof H.J.
        • Bierma-Zeinstra S.M.
        • Arden N.K.
        • Metrustry S.
        • Castano-Betancourt M.
        • Hart D.J.
        • et al.
        Prediction model for knee osteoarthritis incidence, including clinical, genetic and biochemical risk factors.
        Ann Rheum Dis. 2014; 73: 2116-2121
        • Baum T.
        • Joseph G.B.
        • Arulanandan A.
        • Nardo L.
        • Virayavanich W.
        • Carballido-Gamio J.
        • et al.
        Association of magnetic resonance imaging-based knee cartilage T2 measurements and focal knee lesions with knee pain: data from the Osteoarthritis Initiative.
        Arthritis Care Res (Hoboken). 2012; 64: 248-255
        • Dong Y.
        • Peng C.Y.
        Principled missing data methods for researchers.
        Springerplus. 2013; 2: 222
        • Greenland S.
        • Finkle W.D.
        A critical look at methods for handling missing covariates in epidemiologic regression analyses.
        Am J Epidemiol. 1995; 142: 1255-1264
        • Austin P.C.
        A comparison of 12 algorithms for matching on the propensity score.
        Stat Med. 2014; 33: 1057-1069
        • Rosenbaum P.R.
        • Rubin D.B.
        The bias due to incomplete matching.
        Biometrics. 1985; 41: 103-116
        • Austin P.C.
        • Grootendorst P.
        • Normand S.L.
        • Anderson G.M.
        Conditioning on the propensity score can result in biased estimation of common measures of treatment effect: a Monte Carlo study.
        Stat Med. 2007; 26: 754-768
        • Austin P.C.
        Double propensity-score adjustment: a solution to design bias or bias due to incomplete matching.
        Stat Methods Med Res. 2014;
        • Hade E.M.
        • Lu B.
        Bias associated with using the estimated propensity score as a regression covariate.
        Stat Med. 2014; 33: 74-87
        • Becker R.
        • Doring C.
        • Denecke A.
        • Brosz M.
        Expectation, satisfaction and clinical outcome of patients after total knee arthroplasty.
        Knee Surg Sports Traumatol Arthrosc. 2011; 19: 1433-1441
        • Roemer F.W.
        • Kwoh C.K.
        • Hannon M.J.
        • Hunter D.J.
        • Eckstein F.
        • Wang Z.
        • et al.
        Can structural joint damage measured with MR imaging be used to predict knee replacement in the following year?.
        Radiology. 2014; : 140991
        • Stehling C.
        • Lane N.E.
        • Nevitt M.C.
        • Lynch J.
        • McCulloch C.E.
        • Link T.M.
        Subjects with higher physical activity levels have more severe focal knee lesions diagnosed with 3T MRI: analysis of a non-symptomatic cohort of the osteoarthritis initiative.
        Osteoarthr Cartil. 2010; 18: 776-786
        • Mansournia M.A.
        • Danaei G.
        • Forouzanfar M.H.
        • Mahmoodi M.
        • Jamali M.
        • Mansournia N.
        • et al.
        Effect of physical activity on functional performance and knee pain in patients with osteoarthritis : analysis with marginal structural models.
        Epidemiology. 2012; 23: 631-640
        • Collins J.E.
        • Katz J.N.
        • Dervan E.E.
        • Losina E.
        Trajectories and risk profiles of pain in persons with radiographic, symptomatic knee osteoarthritis: data from the osteoarthritis initiative.
        Osteoarthr Cartil. 2014; 22: 622-630
        • Trijau S.
        • Avouac J.
        • Escalas C.
        • Gossec L.
        • Dougados M.
        Influence of flare design on symptomatic efficacy of non-steroidal anti-inflammatory drugs in osteoarthritis: a meta-analysis of randomized placebo-controlled trials.
        Osteoarthr Cartil. 2010; 18: 1012-1018