Study Participants
This cross-sectional analysis used baseline data from the randomized controlled trial registered at clinicaltrials.gov (NCT02512393) to examine the efficacy of transcranial direct current stimulation on knee OA pain. Detailed selection criteria and enrollment procedures have been documented previously [23]. In summary, at baseline, 40 participants with 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 they were aged 50–70 years, had self-reported unilateral or bilateral knee OA pain as per American College of Rheumatology criteria, could speak and read English, and were willing and able to provide written informed consent before enrollment.
Exclusion criteria ensured participants did not have concurrent medical conditions that could confound OA-related outcomes or coexisting diseases that could impede protocol completion, including (1) prosthetic knee replacement or non-arthroscopic surgery on the affected knee, (2) serious medical illness, such as uncontrolled hypertension, heart failure, or 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 (i.e., Mini-Mental Status Exam score ≤ 23), (7) history of brain surgery, tumor, seizure, stroke, or intracranial metal implantation, (8) pregnancy or lactation, and (9) hospitalization for psychiatric illness within the past year.
Measurement
The collected basic characteristics included age, gender, body mass index (BMI; kg/m²), Kellgren-Lawrence radiographic grade, employment status, marital status, educational attainment, and household income.
Clinical Knee OA Pain
The Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) pain subscale. Average knee pain for the past 48 hours was measured by the pain subscale of the WOMAC, which consisted of 5 items on a 0-4 Likert scale measuring the pain severity during walking, climbing stairs, sleeping, resting, and standing [24]. The participants’ responses to each pain question were summed up to derive an aggregated score for pain intensity (range 0-20). The subscales in WOMAC demonstrate reliability and validity in evaluating knee OA patients [25,26].
Short-form-McGill Pain Questionnaire-2 (SF-MPQ-2). The SF-MPQ-2 has been validated and widely used to assess the multidimensional qualities of pain [27]. It consists of four subscales, including continuous pain (6 items: throbbing pain, cramping pain, gnawing pain, aching pain, heavy pain, tender), intermittent pain (6 items: shooting pain, stabbing pain, splitting pain, electric-shock pain, piercing), neuropathic pain (6 items: hot-burning pain, cold-freezing pain, pain caused by light touch, itching, tingling or pins and needles, numbness), and affective description of pain (4 items: tiring-exhausting, sickening, fearful, punishing-cruel). Each subscale score was computed as the average of answered items, with higher scores indicating greater pain intensity.
Quantitative Sensory Testing
A multimodal quantitative sensory testing (QST) battery was completed for an assessment of pain sensitivity using precisely controlled protocols that elicit pain with thermal and mechanical stimuli. This includes heat pain (i.e., threshold and tolerance), pressure pain threshold (PPT), punctate mechanical pain (PMP), and conditioned pain modulation (CPM). The sequence of heat and mechanical testing was randomized and counterbalanced, while CPM was always administered last to minimize any potential carryover effects. Standardized recorded instructions were provided to all participants.
Thermal testing procedures. All thermal stimuli were delivered using a computer-controlled TSA-II NeuroSensory Analyzer (Medoc Ltd., Ramat Yishai, Israel) to measure heat pain thresholds and heat pain tolerances on both the index knee and the ipsilateral ventral forearm using an ascending method of limits. At each body site, the thermode position was moved between trials to prevent sensitization and/or habituation of cutaneous receptors. Starting from a baseline of 32°C, the thermode temperature increased at a rate of 0.5°C per second until participants responded by pressing a button on a handheld device. Participants were instructed to press the button when heat first becomes painful to assess the heat pain threshold, and when they could no longer tolerate the heat pain to assess heat pain tolerance. Three trials of heat pain threshold were conducted at the first test site, followed by three trials of heat pain tolerance were conducted. Then, three trials each at the second test site were conducted, with a 5-minute rest period between sites. The average of the three trials was calculated for each individual, providing overall heat pain threshold and tolerance temperatures for analysis.
Mechanical testing procedures. Mechanical pain response was measured via two approaches. First, PPT was assessed by applying blunt mechanical pressure to deep tissues (i.e., muscle and joint) via a handheld digital pressure algometer (Wagner, Greenwich, CT, USA). Increasing pressure was applied at a constant rate of 0.3 kgf/cm2 per second to measure the PPT at four sites–– the medial and lateral aspects of the index knee, ipsilateral quadriceps, and trapezius. The order of testing sites was counterbalanced and randomized. For assessing PPT, participants were instructed to inform the experimenter when the sensation “first becomes painful” occurred, and the pressure was recorded. The results of the 3 trials at each body site were averaged for each site, and then these PPTs at four sites were averaged to derive an overall measure of PPT. Second, PMP stimuli evaluated cutaneous mechanical sensitivity on both the index patella and the back of the ipsilateral hand. We used calibrated nylon monofilament that delivered a target force of 300 g to obtain verbal ratings of the pain intensity on a scale of 0 (no pain sensation) to 100 (the most intense pain sensation imaginable) following 10 contacts at the rate of 1 contact per second. An overall score for each site was computed by averaging across two trials.
Conditioned pain modulation (CPM). Ten minutes after assessing the thermal or mechanical pain, the CPM was evaluated. CPM reflects the endogenous pain inhibitory pathway (i.e., descending pain inhibition) also known as the “pain inhibits pain” paradox [28] CPM was assessed by determining the change in PPT on the trapezius, immediately following the immersion of the contralateral hand up to the wrist in the cold water bath (12°C) for 1 minute. The initial pre-immersion PPT measurement was conducted just before placing the hand in the water. Thirty seconds after hand immersion, participants were asked to rate the cold pain intensity (0–100) from the immersed hand followed by the second PPT measurement and were informed to keep their hand in the water bath for as long as tolerable up to 1 minute. After the removal of the hand, the final PPT measurement was taken. This temperature was chosen based on prior experience with middle-aged and older adults with knee OA, where 12°C was found to produce moderate yet tolerable pain for most participants. Water was continually circulated and maintained at a constant temperature by a refrigeration unit (Neslab, Portsmouth, NH, USA). An increase in PPT following cold water immersion demonstrated pain inhibition.
Pain Catastrophizing
The Coping Strategies Questionnaire-Revised (CSQ-R) measures the use of strategies for coping with pain by assessing six domains–– distraction, catastrophizing, ignoring pain sensations, distancing from pain, coping self-statements, and praying. Participants rate how often they use specific strategies on a seven-point Likert scale from 0 (“never”) to 6 (“always”), with higher scores indicating greater usage for each domain. This study used the 6-item catastrophizing subscale, with scores calculated as the mean of the responses. The reliability and validity of the CSQ subscales have previously been shown to be acceptable [29,30].
Negative Affect
The Positive and Negative Affect Schedule (PANAS) includes 20 items that evaluate the frequency of both pleasant and unpleasant emotions individuals experience [31]. The inventory is divided into two subscales, each with 10 items for positive and negative emotions. Negative affect is calculated from the sum of 10 items (afraid, ashamed, distressed, guilty, hostile, irritable, jittery, nervous, scared, and upset), rated on a 5-point scale from 1 (“very slightly or not at all”) to 5 (“extremely”). A lower total negative score (PANAS-) indicates less negative affect (range: 10-50). The PANAS has been validated and demonstrates reliability, with an alpha coefficient range of .84 to .87 for negative affect [31].
Statistical Analyses
Descriptive statistics were used to evaluate participant characteristics. Composite measures for QST were created by calculating z-scores for the heat pain threshold and tolerance at the arm and knee; PPT at the medial and lateral aspects of the index knee, ipsilateral quadriceps, and trapezius; and PMP at the index patella and hand. The z-scores for each pain measure were subsequently averaged across the body sites to yield overall heat pain threshold, heat pain tolerance, PPT, and PMP values for the analyses.
Separate path analytical models were estimated to assess the indirect effects (mediation) of ethnicity (coded 0 for NHWs and 1 for Asian Americans) via pain catastrophizing or negative affect on each clinical and experimental pain measure. The path models facilitated the examination of both direct and indirect effects. Model fit, path coefficient estimates, and 95% highest posterior density (HPD) confidence intervals (CIs; “credibility” in Bayesian terms) for parameter estimates were generated using the Bayesian estimation method in Mplus (version 8.8; Muthén & Muthén, Los Angeles, CA) [32]. Bayesian estimation is advantageous in that it precludes the necessity of the normality assumption in the sampling distribution of estimates and potentially provides more accurate parameters in small-sample cases [33]. Model fit was evaluated using the criteria and methods recommended by Muthén and Asparouhov [34]. Where 95% CIs did not overlap with zero, the effect was considered significant.