Study Participants
This cross-sectional analysis used baseline data from a randomized controlled trial registered at clinicaltrials.gov (NCT02512393) to assess the efficacy of transcranial direct current stimulation (tDCS) in mitigating knee OA pain. Detailed selection criteria and enrollment procedures have been described previously (Ahn et al., 2017b). In brief, at baseline, 40 participants with self-reported knee OA pain (20 NHWs and 20 Asian Americans) were recruited in North Central Florida between September 2015 and August 2016 through local advertisements. Participants were eligible if aged 50–70, if they reported unilateral or bilateral knee OA pain according to American College of Rheumatology criteria, were able to speak and read English, and were willing and able to provide written informed consent before enrollment.
Exclusion criteria ensured that participants did not have concurrent medical conditions that could confound OA-related outcomes or coexisting diseases that could hinder protocol completion. Thus, the following were the exclusion criteria: (1) having undergone prosthetic knee replacement or non-arthroscopic surgery on the affected knee, (2) a serious medical illness, such as uncontrolled hypertension, heart failure, or a recent history of acute myocardial infarction, (3) peripheral neuropathy, (4) systemic rheumatic disorders, such as rheumatoid arthritis, systemic lupus erythematosus, and fibromyalgia, (5) alcohol or substance abuse, (6) cognitive impairment, defined as a Mini-Mental Status Exam score of 23 or lower, (7) a history of brain surgery, tumor, seizure, stroke, or intracranial metal implantation, (8) pregnancy or lactation, and (9) hospitalization for a psychiatric illness within the past year.
Measurement
Knee OA Symptoms: Knee OA pain and physical function
Knee OA pain and physical function were measured using the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) pain and functional subscales, where higher scores indicate greater pain and physical functional disability (Bellamy et al., 1988). The pain subscale includes 5 items on a 5-point Likert scale (0 being none to 4 being extreme) measuring the pain severity during walking, climbing stairs, sleeping, resting, and standing. The participants’ responses to the pain questions were summed up to derive an aggregated score of pain intensity (range 0–20). The functional subscale asks patients to rate the degree of difficulty in accomplishing 17 activities of daily living on a 5-point scale (0 being none to 4 being extreme). The participants’ responses were aggregated to produce a composite score of functional disability (range 0–68). The subscales in WOMAC demonstrate reliability and validity in evaluating knee OA patients (Bellamy & Buchanan, 1986; Roos et al., 1999).
Inflammatory Markers
In this study, we also gathered data regarding the following inflammatory markers: CRP, tumor necrosis factor-alpha (TNF-α), IL-1β, IL-6, and IL-10. Owing to excessively missing data (45.0%–80.0%), we excluded IL-1β and IL-6 from the current analysis. In the original study (blinded for review), blood samples were obtained prior to treatment initiation (i.e., tDCS) on Day 1 and again after the completing the fifth treatment on Day 5. For our analysis, we utilized pre-treatment data acquired on Day 1. Blood was drawn into ethylenediaminetetraacetic acid plasma tubes. Samples were inverted five times and stored on ice until further processing. Within 30 min of being collected, samples were centrifuged at 1,600 × g and 4 ºC for 15 min, aliquoted, and immediately stored in a −80 ºC freezer.
The plasma samples underwent solid-phase extraction using an Oasis™ Hydrophilic–Lipophilic-Balanced (30 mg) 96-well plate along with a vacuum manifold (Waters Corp., Milford, MA, USA), according to the manufacturer’s protocol. Briefly, the plate was conditioned with acetonitrile and equilibrated twice with 0.1% trifluoroacetic acid (TFA) in high-performance liquid chromatography (HPLC)-grade water. Samples were acidified with 1% TFA (1:1) and loaded onto the plate. The plate was washed thrice with 0.1% TFA in HPLC-grade water. The samples were eluted in 60% acetonitrile/40% HPLC-grade water/0.1% TFA and dried in a Savant AES1010 Automatic Environmental SpeedVAC® w/VaporNet Radiant Cover. Thereafter, they were reconstituted using the original sample volume in assay buffer.
Plasma CRP levels were measured in duplicate using enzyme-linked immunosorbent assays, following the manufacturers’ instructions (cat# DCRP00, R&D Systems, Minneapolis, MN; cat# ADI-900-071, Enzo Life Sciences, Inc., Farmingdale, NY, respectively). For CRP, the average intra- and interassay CV values were < 10.0% and < 7.0%, respectively. TNF-α and IL-10 plasma levels were measured in triplicate using a commercial multiplex immunoassay kit (cat# HCYTMAG-60Kl; MilliporeSigma, Burlington, MA) and analyzed using the MILLIPLEX® Analyzer 3.1 xPONENT® System. Data acquisition was accomplished using the same system and data analysis performed via MILLIPLEX® Analyst Software. Intra- and interassay CVs were < 19.0% for all markers.
Statistical Analysis
Descriptive and comparative statistics were employed to determine sample characteristics. As the inflammatory markers were not normally distributed, log transformation was applied to mitigate skewness. When missing data were present (CRP, n = 4), listwise deletion was performed, resulting in a streamlined dataset for analysis. Race-stratified analyses were conducted owing to an observed interaction between race and certain inflammatory markers, such as CRP (data not shown). We examined the relationships between body mass index (BMI) and each inflammatory marker, aiming to circumvent possible collinearity (since adiposity proves to be significantly associated with systemic inflammation), and specifically investigated their associations with knee OA symptoms. Multivariable linear regression models were used in the main analysis of each outcome.
Explanatory variables included age, sex (male vs. female), marital status (partnered vs. unpartnered), BMI (kg/m²), and Kellgren–Lawrence (KL) radiographic grade (0–1 vs. 2–4), all of which potentially affect both inflammation and knee OA symptoms (Hunter et al., 2008; Heidari et al., 2016). Pain catastrophizing and negative affect were also considered because of their potential relation to widespread pain (Dave et al., 2015; Finan et al., 2013). Candidate variables comprised those with p < .200 in bivariate analyses: single-factor analysis of variance in cases of variance equality, the Kruskal–Wallis test for qualitative variables, and the simple linear regression test for quantitative variables. Multivariable analysis based on stepwise selection at an alpha value of 0.05 was conducted to preserve the most relevant variables in the model and distinguish those independently associated with the outcomes. In all multivariable models, systematic adjustment for age was performed. Finally, we conducted diagnostic tests to ensure that the multivariable models satisfied linear regression model assumptions. All statistical analysis was carried out using R Studio version 4.0.2 (www.R-project.org).