Study population
We performed analyses using publicly available data from the OAI, a multicenter prospective cohort of KOA registered under ClinicalTrials.gov with the identifier NCT00080171. The OAI recruited a total of 4796 participants from four distinct centers: University of Pittsburgh, University of California, University of Maryland and Johns Hopkins University Joint Center, the Ohio State University, and Memorial Hospital/Brown University, all of whom were diagnosed with KOA or exhibited high-risk KOA features. The specific details of this cohort's protocol have been previously well-documented16.
In this study, we used the knee as the primary unit of analysis and included knees that had a complete assessment of quadriceps strength, effusion-synovitis, and knee pain at baseline, as well as at 1- or 2-year follow-ups from the OAI. It is possible that two included knees may originate from the same participant. The specific selection flowchart is presented in Fig. 1.
The characteristics of the study population were assessed across various dimensions, including age, sex, race, body mass index (BMI), radiographic osteoarthritis (ROA), knee alignment, physical activity, depressive symptoms, quadriceps strength, effusion-synovitis assessments, and knee pain at baseline. Upon enrollment, age, sex, race, and BMI were recorded. ROA was evaluated using the Kellgren-Lawrence grading (KLG) system, which is based on X-ray radiographs. To account for the potential influence of knee alignment on quadriceps strength within the knee joint, the degree of alignment for each knee was measured using the femorotibial angle (FTA)17,18. Physical activity and depressive symptoms were assessed using the Physical Activity Scale for the Elderly (PASE) and the Center for Epidemiologic Studies Depression Scale (CESD), respectively. Isometric quadriceps strength was measured at baseline, while effusion-synovitis and knee pain were evaluated at baseline, 1-year, and 2-year follow-ups. The Magnetic Resonance Imaging Osteoarthritis Knee Score (MOAKS) and the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) were employed to quantify the levels of effusion-synovitis and knee pain, respectively. All of these baseline characteristics of included participants were considered as covariates in our analyses to account for their potential influence on final results.
Quadriceps strength
The measurement of quadriceps strength was collected at baseline. According to the protocol established by the OAI, the isometric measurement of quadriceps strength for each knee was conducted using the Good Strength Chair (Metitur Oy, Jyvaskyla, Finland)19. Prior to the testing, all research technicians underwent training followed a standardized protocol, and the test-retest reliability yielded range of 0.88 to 0.92. Participants were seated in the Good Strength Chair, with the tested knee fixated at a 60° angle. Straps were used to secure the pelvis and thigh, stabilizing the participants' trunk and lower limb. A load cell, situated on a lever arm, was attached approximately 2 cm proximal to the calcaneus20. Before the formal test, participants completed two warm-up trials at 50% of their maximum effort. During the test, three repetitions were performed with the participant's maximum effort, and the peak force (measured in Newtons) was recorded as the maximum quadriceps strength. To account for variations in participant weight and height, the maximum quadriceps strength was normalized using body mass and leg length, resulting in units of Nm/kg.
Effusion-synovitis
The assessment of effusion-synovitis was collected at baseline, 12-month and 24-month follow-ups, utilizing magnetic resonance imaging (MRI). The structure of each knee was evaluated using identical 3.0 Tesla Trio MRI systems (Siemens) at each OAI clinical center. The specific MRI pulse sequence parameters outlined in the OAI protocol were previously published elsewhere21.
Two fellowship-trained musculoskeletal radiologists, A.G. and F.W.R., who had over 10 years of experience in evaluating knee structure, performed the reading of the MRIs using the MOAKS22. They were blinded to the participants' clinical information. For this study, our primary focus was on the changes in effusion-synovitis within the initial two-year period. Therefore, semiquantitative assessments were conducted at baseline, 12-month, and 24-month follow-ups. In accordance with the MOAKS system, the distention of the synovial cavity within the intercondylar region, referred to as effusion-synovitis, was categorized into four levels: level 0 (none, representing a physiological amount), level 1 (small, less than 33% maximum distention), level 2 (medium, between 33% and 66% maximum distention), and level 3 (large, more than 66% maximum distention). The effusion-synovitis assessment demonstrated a high level of inter-rater reliability (0.95, 95% CI 0.61-1.00). Changes in effusion-synovitis were defined as the difference between the follow-up score and the baseline score.
Knee pain
The severity of knee pain at each visit (baseline, 12- and 24-month follow-ups) was evaluated using the WOMAC pain subscale. This subscale assesses knee pain during five daily activities over a past seven-day period, including lying, sitting, standing, stair climbing, and walking on a flat surface. The cumulative score of each activity is then combined to derive the total points of the WOMAC scale, which spans from 0 to 2023. Higher scores indicate more severe pain. The WOMAC pain scale has been rigorously validated for its reliability and validity24. Changes in knee pain were determined by subtracting the WOMAC pain score at baseline from the WOMAC pain score at follow-ups.
Statistical analysis
The baseline characteristics of the study population were summarized using mean ± SD or percentages. Linear regression models were employed to analyze the association between quadriceps strength or effusion-synovitis with WOMAC pain. Adjustments were made in the models for sex, age, BMI, race, KLG, alignment, PASE score, CESD score, and WOMAC pain at baseline. Longitudinal regression analyses were conducted using generalized estimating equations (GEE) to account for within-individual correlations between the two knees. The results were reported as beta coefficients along with their respective 95% confidence intervals (95% CIs).
A mediation analysis was conducted using the Mediation Package (version 4.5.0) in R (version 4.2.2) to examine the mediating role of effusion-synovitis in the relationship between quadriceps weakness and worsening knee pain. The robustness was ensured by employed 2000 bootstraps for analysis. The total effect of quadriceps weakness on knee pain comprised both a natural direct effect and a natural indirect effect that was mediated by effusion-synovitis changes. To estimate the extent of mediation in this relationship, the mediation proportion, along with a 95% confidence interval, was calculated. The mediation model was adjusted for sex, age, race, BMI, KLG, alignment, CESD score, PASE score, and WOMAC knee pain at baseline. Furthermore, a sensitivity analysis was conducted to account for potential confounding factors, specifically the diagnosis of radiographic osteoarthritis at baseline. In the sensitivity analysis, knees were categorized as with ROA (KLG ≥ 2) or without ROA at baseline, considering the potential influence of ROA on changes in knee pain.
Regression analyses were performed using IBM SPSS (version 27.0), and mediation analyses were conducted using R software (version 4.2.2). A statistical significance level of P < 0.05 was used.