Transcutaneous Auricular Vagus Nerve Stimulation as a Treatment for Neuropathic Pain Following Spinal Cord Injury

Recruitment feasibility will be assessed by documenting the number and proportion of eligible individuals who consent to participate and the total time required to recruit the target sample of 32 participants. Recruitment rate will be calculated as the percentage of eligible candidates who agree to participate. These data will be summarized descriptively for the overall study population.

Retention will be assessed as the proportion of enrolled participants who complete the full 30-day intervention and post-intervention assessment. Attrition will be recorded as the number and proportion of participants who withdraw prior to study completion. Reasons for withdrawal or loss to follow-up will be documented where available. Retention and attrition rates will be summarized and compared descriptively between the active taVNS and sham taVNS groups.

Adherence will be assessed using automated usage data recorded by the taVNS device and associated smartphone application. Daily adherence will be calculated as the percentage of the prescribed 4-hour daily stimulation period completed (actual stimulation time divided by prescribed stimulation time × 100). Overall adherence will be summarized as the mean daily adherence across the 30-day intervention period. Adherence will be compared descriptively between the active taVNS and sham taVNS groups.

Participant acceptability will be assessed at the end of the intervention using the Acceptability of Intervention Measure (AIM). Acceptability scores will be summarized descriptively and compared between the active taVNS and sham taVNS conditions.

All adverse events which occur during the intervention will be recorded during the twice weekly phone calls and/or study visits and compared between the active taVNS and sham taVNS conditions. A description of the adverse event, number of events, and total number of participants affected will be recorded.

Blinding success will be assessed after completion of the intervention by asking participants and investigators to indicate which treatment group they believe the participant was assigned to (active taVNS or sham taVNS). Blinding effectiveness will be quantified using the James Blinding Index, calculated separately for participants and investigators. Blinding outcomes will be summarized descriptively.

Neuropathic pain severity and interference will be assessed using the Neuropathic Pain Symptoms Inventory (NPSI), a self-reported questionnaire with total scores ranging from 0 to 100, where higher scores indicate greater neuropathic pain severity and pain-related interference. Change scores will be calculated as post-intervention values minus baseline values. Mean changes and variability will be summarized and compared descriptively between the active taVNS and sham taVNS groups.

Neuropathic pain severity and pain-related interference will be assessed using the International Spinal Cord Injury Pain Basic Data Set (ISCIPBDS), version 3, which includes numeric rating scales ranging from 0 to 10, where higher scores indicate greater pain intensity and greater pain-related interference. Change scores will be calculated as post-intervention values minus baseline values. Mean changes and variability will be summarized and compared descriptively between the active taVNS and sham taVNS groups.

Plasma concentrations of C-reactive protein (CRP) will be measured from fasting blood samples collected at baseline and at the end of the intervention. Concentrations will be reported in milligrams per liter (mg/L), with higher values indicating greater systemic inflammation. Change values will be calculated as post-intervention values minus baseline values and summarized descriptively between the active taVNS and sham taVNS groups.

Plasma concentrations of interleukin-1 beta (IL-1β) will be measured from fasting blood samples collected at baseline and at the end of the intervention. Concentrations will be reported in picograms per milliliter (pg/mL), with higher values indicating greater pro-inflammatory activity. Change values will be calculated as post-intervention values minus baseline values and summarized descriptively between the active taVNS and sham taVNS groups.

Plasma concentrations of interleukin-6 (IL-6) will be measured from fasting blood samples collected at baseline and at the end of the intervention. Concentrations will be reported in picograms per milliliter (pg/mL), with higher values indicating greater pro-inflammatory activity. Change values will be calculated as post-intervention values minus baseline values and summarized descriptively between the active taVNS and sham taVNS groups.

Plasma concentrations of interferon gamma (IFN-γ) will be measured from fasting blood samples collected at baseline and at the end of the intervention. Concentrations will be reported in picograms per milliliter (pg/mL), with higher values indicating greater pro-inflammatory activity. Change values will be calculated as post-intervention values minus baseline values and summarized descriptively between the active taVNS and sham taVNS groups.

Plasma concentrations of tumor necrosis factor alpha (TNF-α) will be measured from fasting blood samples collected at baseline and at the end of the intervention. Concentrations will be reported in picograms per milliliter (pg/mL), with higher values indicating greater pro-inflammatory activity. Change values will be calculated as post-intervention values minus baseline values and summarized descriptively between the active taVNS and sham taVNS groups.

Vagal tone will be assessed using heart rate variability (HRV) derived from a 5-minute resting electrocardiogram recording. Root mean square of successive differences (RMSSD) and high-frequency power will be calculated using standard HRV analysis procedures. Change in HRV measures from baseline to post-intervention will be summarized descriptively and compared between the active taVNS and sham taVNS groups.

Quality of life will be assessed using selected short forms from the Spinal Cord Injury-Quality of Life (SCI-QOL) measurement system, including Pain Interference, Depression, and Positive Affect and Well-Being domains. Scores will be calculated as standardized T-scores. Change in quality-of-life scores from baseline to post-intervention will be summarized and compared descriptively between treatment groups.