Abstract
An increasing number of studies supports the view that transcutaneous electrical stimulation of the spinal cord (TESS) promotes functional recovery in humans with spinal cord injury (SCI). However, the neural mechanisms contributing to these effects remain poorly understood. Here we examined motor-evoked potentials in arm muscles elicited by cortical and subcortical stimulation of corticospinal axons before and after 20 min of TESS (30 Hz pulses with a 5 kHz carrier frequency) and sham-TESS applied between C5 and C6 spinous processes in males and females with and without chronic incomplete cervical SCI. The amplitude of subcortical, but not cortical, motor-evoked potentials increased in proximal and distal arm muscles for 75 min after TESS, but not sham-TESS, in control subjects and SCI participants, suggesting a subcortical origin for these effects. Intracortical inhibition, elicited by paired stimuli, increased after TESS in both groups. When TESS was applied without the 5 kHz carrier frequency both subcortical and cortical motor-evoked potentials were facilitated without changing intracortical inhibition, suggesting that the 5 kHz carrier frequency contributed to the cortical inhibitory effects. Hand and arm function improved largely when TESS was used with, compared with without, the 5 kHz carrier frequency. These novel observations demonstrate that TESS influences cortical and spinal networks, having an excitatory effect at the spinal level and an inhibitory effect at the cortical level. We hypothesized that these parallel effects contribute to further the recovery of limb function following SCI.
SIGNIFICANCE STATEMENT Accumulating evidence supports the view that transcutaneous electrical stimulation of the spinal cord (TESS) promotes recovery of function in humans with spinal cord injury (SCI). Here, we show that a single session of TESS over the cervical spinal cord in individuals with incomplete chronic cervical SCI influenced in parallel the excitability cortical and spinal networks, having an excitatory effect at the spinal level and an inhibitory effect at the cortical level. Importantly, these parallel physiological effects had an impact on the magnitude of improvements in voluntary motor output.
- corticospinal
- intracortical inhibition
- neurophysiology
- neuroplasticity
- spinal cord injury
- spinal networks
Introduction
An increasing number of studies support the view that transcutaneous electrical stimulation of the spinal cord (TESS) promotes functional recovery in humans with spinal cord injury (SCI). TESS has shown improvements in sensory and motor function when it is applied alone or in combination with other therapies (Hofstoetter et al., 2013; Gerasimenko et al., 2015; Gad et al., 2017, 2018; Inanici et al., 2018). Improved cardiovascular (Phillips et al., 2018) and bladder (Horst et al., 2013; Radziszewski, 2013) functions and reduced spasticity (Hofstoetter et al., 2014, 2020) have also been observed. While it is clear that single TESS pulses generate spinally mediated responses activating spinal networks that project transsynaptically to motoneurons (Milosevic et al., 2019), the effects of repeated use of TESS on CNS pathways remain unknown. A more complete understanding of the neural mechanisms of TESS-induced recovery is likely to enhance rehabilitation effectiveness.
Electrophysiological (Minassian et al., 2007a; Hofstoetter et al., 2018) and modeling (Rattay et al., 2000; Capogrosso et al., 2013) studies suggest that a single TESS pulse recruits large-to-medium-diameter proprioceptive and cutaneous afferents within posterior rootlets/roots of different spinal segments depending on the stimulus intensity. Most studies in humans with SCI aiming to achieve therapeutic effects use TESS currents between 5 and 50 Hz with intensities fluctuating between 10 and 200 mA using a high carrier frequency between 5 and 10 kHz for several minutes to hours (see references above). In humans, a prolonged period of repeated stimulation of afferent fibers in a peripheral nerve with frequencies between 5 and 50 Hz usually increases the net excitability of descending cortical–motor neuron pathways, at least in part by reducing the activity of GABAergic cortical mechanisms (Kaelin-Lang et al., 2002; Golaszewski et al., 2012). Note that TESS uses a high carrier frequency or a “Russian current”. Russian currents are alternating currents in a frequency range from ∼100 Hz to 10 kHz that are delivered in bursts, with a burst frequency within a “physiological” range (up to ∼100 Hz; Ward et al., 2009). It is thought that this high carrier frequency is beneficial for improving muscle strength (Selkowitz, 1985, 1989) and for suppressing the sensitivity of pain receptors (Ward and Robertson, 1998a,b). High-frequency electrical stimulation of the spinal cord releases serotonin in the dorsal horn (Linderoth et al., 1993), which may depress ascending nociceptive transmission (Fürst, 1999; Millan, 2002). Indeed, stimulation of muscle afferents at frequencies ∼100 Hz has an inhibitory effect on motor cortex (Mima et al., 2004). We hypothesized that TESS has a predominant excitatory effect on spinal networks and predominantly inhibitory effects at the cortical level, both of which could serve to fine-tune the net output that could contribute to further improving functional recovery.
To test our hypothesis, we evoked motor potentials in arm muscles by cortical and subcortical stimulation of corticospinal axons before and after 20 min of TESS (30 Hz pulses with a 5 kHz carrier frequency) applied between C5 and C6 spinous processes in humans with and without chronic incomplete cervical SCI. Structural magnetic resonance imaging (MRI) was used to localize the extent of the spinal cord lesion and to ensure that TESS was applied over injured segments.
Materials and Methods
Participants.
Seventeen individuals with SCI participants (mean age, 43.1 ± 14.0 years; 4 females) and 15 age-matched control subjects (mean age, 36.5 ± 17.5 years; 6 females; p = 0.3) participated in the study (Table 1). All participants gave informed consent to experimental procedures, which were approved by the local ethics committee at the University of Miami in accordance with guidelines established in the Declaration of Helsinki. Participants with SCI had a chronic injury (≥1 year) and were classified using the International Standards for Neurological Classification of Spinal Cord Injury examination as having a C4–C6 SCI and by the American Spinal Cord Injury Association Impairment Scale (AIS) as AIS A (n = 4), AIS B (n = 3), AIS C (n = 2), or AIS D (n = 8). Five SCI individuals were taking antispastic medication (baclofen and/or gabapentin and/or tizanidine; Table 1).
Spinal cord injury participants
Study design.
All subjects participated in the main experiment, in which a 2 d crossover design was used to compare the effect of 20 min of TESS or sham-TESS, applied between C5 and C6 spinous processes, on motor-evoked potentials (MEPs) elicited by transcranial magnetic stimulation (TMS) and by cervicomedullary electrical stimulation (CMEPs) in the biceps brachii. Participants, but not experimenters, were blinded to the interventions. The duration of stimulation was comparable to that used in other plasticity protocols targeting the spinal cord in humans (Taylor and Martin, 2009; Bunday and Perez, 2012). The following tests were also completed in additional sessions in a subset of subjects: (1) we examined the effect of 20 min of TESS on MEPs and CMEPs elicited in the triceps brachii (control subjects = 10; SCI participants = 11) and first dorsal interosseous (control subjects = 8; SCI participants = 8) muscles. This was tested separately because different coil locations and stimulus intensities are needed to elicit MEPs and CMEPs in the triceps brachii and the first dorsal interosseous than in the biceps brachii; (2) because MEP amplitude increased after TESS, but not sham-TESS, we examined the effect of 20 min of TESS on short-interval intracortical inhibition (SICI; control subjects = 14; SCI participants = 17) in a separate session to assess the contribution of intracortical pathways to the observed changes in MEP amplitude; and (3) to assess the effect of TESS without the 5 kHz carrier frequency (TESS without 5 kHz), we examined the effects of 20 min of TESS without 5 kHz on MEPs and CMEPs and on SICI in separate sessions (control subjects = 8; SCI participants = 8). Subjects tested on these sessions were randomly selected.
Electromyographic (EMG) recordings.
Electromyograms were recorded from the biceps brachii muscle of the right side in control subjects and from the less affected hand in individuals with SCI through surface electrodes secured to the skin over the belly of each muscle (Ag-AgCl, 10 mm diameter). To examine whether the effects of TESS spread to other muscles, we also recorded from the triceps brachii and the first dorsal interosseous. The signals were amplified, filtered (30–2000 Hz), and sampled at 10 kHz for offline analysis (CED 1401 with Signal software, Cambridge Electronic Design). Subjects performed three brief elbow flexion maximal voluntary isometric contractions (MVCs) for 3–5 s, separated by ∼30 s of rest. The maximal mean EMG activity of the rectified response generated for 1 s during each MVC was analyzed, and the highest value of the three trials was used. The MVC with the biceps brachii during elbow flexion was larger in control subjects (1.4 ± 0.5 mV) compared with SCI participants (0.8 ± 0.4 mV; p = 0.007).
Experimental setup.
Subjects were seated comfortably in an armchair with both arms placed on a custom platform with the elbow and shoulder flexed at 90° (Fig. 1A). Individuals randomly received TESS and sham-TESS on different sessions separated by 2–3 d. During the TESS session, a custom-made five-channel constant current stimulator (BioStim-5, BioMedical Life Systems; Grishin et al., 2017) was used to noninvasively stimulate cervical spinal networks. Stimulation was delivered between the spinous process of C5–C6 using a 3.2-cm-diameter hydrogel adhesive electrode as cathode and two 7.5 × 12 cm rectangular electrodes (ValuTrode Cloth, Axelgaard Manufacturing) located bilaterally over the iliac crests as anodes. We used blocks of five biphasic pulses (each pulse of 200 μs duration) at a frequency of 30 Hz (every 32.3 ms) for 20 min (Fig. 1B). Thus, each block was filled with a 5 kHz carrier frequency. The stimulus intensity from a single biphasic pulse was defined as the minimal intensity required to induce a root-evoked potential in the biceps brachii ≥50 μV peak-to-peak amplitude above the background EMG in 5 of 10 consecutive trials in the relaxed muscle [control subjects = 77.5 ± 10.8 mA (range, 55–90 mA); SCI participants = 79.2 ± 12.2 mA (range, 55–95 mA); p = 0.6]. During the sham-TESS session, the electrodes were located in the same position, and the intensity of stimulation was set as the same intensity as used in TESS sessions and gradually decreased down to 0 in ∼1 min. Similar sham stimulation procedures have been used successfully in previous neuromodulation studies using TESS to provide the initial sensation of stimulation without the subsequent effects (Murray et al., 2018; Awosika et al., 2019). During TESS without 5 kHz carrier frequency (tested on subjects listed as 1–8 on Table 1), single biphasic pulses (200 μs duration) were delivered at a 30 Hz frequency (every 33.1 ms) for 20 min using the same intensity as described above. This paradigm was tolerated by all subjects, allowing us to match the same stimulation intensity as used in the original TESS protocol. By comparing the effects of TESS with and without the 5 kHz carrier frequency, we can examine whether the additional pulses provided with the carrier frequency are causing any effects of physiological outcomes. Since changes in pulse duration can affect the relative recruitment of sensory axons, we maintained the waveform pulse duration consistent across paradigms (Paillard, 1955; Panizza et al., 1989). The following measurements were taken at rest before, immediately after, and up to 75 min after the end of each stimulation period: MEPs, CMEPs, SICI, and functional outcomes.
Experimental setup. A, Participants were comfortably seated in a customized chair during TESS or sham-TESS for 20 min. TESS was delivered using a surface electrode on the back of the neck between C5 and C6 spinous processes segments (cathode) and a surface electrode in each anterior crest of the hip bone (anode) using a custom-made five-channel stimulator (BioStim5). Electrophysiological and behavioral outcomes were tested Pre, immediately after, and 15, 30, 45, 60, and 75 min after the end of the stimulation or sham stimulation period. B, Schematic representation of the type of current used during TESS. We used five biphasic pulses at 5 kHz with each biphasic pulse lasting for 200 μs. The middle scheme shows the blocks of five biphasic pulses passed at a 30 Hz frequency. The bottom part of the schematic shows the number of pulses delivered in 1 s.
MEPs.
Transcranial magnetic stimuli were delivered from a Magstim 200 Stimulator (Magstim Company) through a figure-of-eight coil (custom-made loop diameter, 7 cm; type number SP15560) with a monophasic current waveform. TMS was delivered to the optimal scalp position for activation of the biceps brachii, triceps brachii, and first dorsal interosseous muscles as needed. The coil was held tangentially to the scalp to induce currents in the brain flowing in a posterior–anterior direction and held to the head of the subject with coil holder (Manfrotto) with the head firmly secured to a headrest to limit head movements. TMS stimuli were delivered at an intensity to generate MEPs of ∼3–5% of the maximal motor response (M-max) in the biceps brachii (control subjects = 4.4 ± 1.5% of the M-max; SCI participants = 5.0 ± 1.1% of the M-max; p = 0.9), triceps brachii (control subjects = 2.9 ± 1.9% of the M-max; SCI participants = 3.1 ± 2.8% of the M-max; p = 0.4), and the first dorsal interosseous (control subjects = 4.5% ± 0.2 of the M-max; SCI participants = 3.6 ± 1.3% of the M-max; p = 0.1) muscles. Peripheral nerve stimulation was delivered to the brachial plexus at the supraclavicular fossa to test the M-max in the biceps brachii (200 μs pulse duration; model DS7AH, Digitimer). The anode was placed at the Erb's point, and the cathode was placed over the acromion for biceps and triceps brachii M-max, and at the ulnar nerve at the wrist for the first dorsal interosseous (biceps brachii: control subjects = 21.7 ± 4.7 mV; SCI participants = 18.2 ± 4.8 mV; p = 0.02; triceps brachii: control subjects = 14.7 ± 3.5 mV; SCI participants = 11.8 ± 3.6 mV; p = 0.03; first dorsal interosseous: control subjects = 22.7 ± 1.1 mV; SCI participants = 13.7 ± 4.3 mV; p < 0.001). TMS pulses were delivered at 4 s intervals (0.25 Hz). Twenty MEPs were averaged at each time point before and after each protocol, and peak-to-peak MEP amplitude was measured.
CMEPs.
CMEPs were elicited by simulation after elicited at the cervicomedullary junction by using high-voltage electrical current (100 μs duration; DS7AH, Digitimer) passed between adhesive Ag-AgCl electrodes fixed to the skin behind the mastoid process. The stimulation intensity was set at an intensity to generate CMEPs of ∼3–5% of the M-max in the biceps brachii (control subjects = 4.6 ± 1.3% of the M-max; SCI participants = 5.2 ± 1.0% of the M-max; p = 0.1), triceps brachii (control subjects = 3.5 ± 2.6% of the M-max; SCI participants = 4.4 ± 2.7% of the M-max; p = 0.2), and the first dorsal interosseous (control subjects = 4.5 ± 0.2% of the M-max; SCI participants = 3.6 ± 1.3% of the M-max; p = 0.1) muscles. We monitored the stimulation to ensure that it was below the intensity required to activate the peripheral nerve directly by increasing the intensity and observing a decrease in latency. The latency of CMEPs was shorter than the MEPs elicited by TMS in both groups in all muscles tested (biceps brachii, control subjects: MEPs = 13.0 ± 0.8 ms; CMEPs = 8.4 ± 0.5 ms; p < 0.001; SCI participants: MEPs = 13.7 ± 0.8 ms; CMEPs = 8.5 ± 1.1 ms; p < 0.001; triceps brachii, control subjects: MEPs = 13.4 ± 0.9 ms; CMEPs = 8.6 ± 0.4 ms; p < 0.001; SCI participants: MEPs = 14.4 ± 1.1 ms; CMEPs = 8.7 ± 1.1 ms; p < 0.001; first dorsal interosseous, control subjects: MEPs = 22.2 ± 0.3 ms; CMEPs = 17.7 ± 0.9 ms; p < 0.001; SCI participants: MEPs = 25.5 ± 2.7 ms; CMEPs = 20.0 ± 1.8 ms; p < 0.001). Five to 10 CMEPs were averaged at each time point before and after each protocol for each muscle on separate sessions, and peak-to-peak CMEP amplitude was measured.
SICI.
We tested SICI in biceps brachii using a previously described method (Kujirai et al., 1993). A conditioning stimulus (CS) was set at an intensity needed to elicit ∼50% of SICI, which corresponded to ∼80% of resting motor threshold. The same stimulation intensity was used before and after the TESS protocol. The test stimuli (TS) was set at an intensity needed to elicit an MEP that equals ∼3–5% of the M-max (control subjects = 4.6 ± 1.1% of the M-max; SCI participants = 4.8 ± 1.4% of the M-max; p = 0.3). The delay between CS and TS was 2.5 ms. SICI was calculated by expressing the amplitude of conditioned MEP as a percentage of the test MEP amplitude [% = (conditioned MEP × 100)/(test MEP)]. Two sets of 15 test MEPs and 15 conditioned MEPs were tested at each time point.
Functional outcomes.
In the SCI group, we tested gross and fine motor functions (jar opening and water bottle test, and key test and coin test, respectively) using subcomponents of the Graded and Redefined Assessment of Strength, Sensibility and Prehension (GRASSP) test. During jar opening, subjects were asked to open a jar lid with a tested hand while holding the jar (7 cm in diameter and 9 cm in height) with the nontested hand in the shortest time possible. For the water bottle test, subjects were asked to lift a 6-cm-diameter bottle, filled with ∼200 ml of water, from the table and pour the water into a cup. During the key test, subjects were asked to pick up a key from the table, insert it in a lock, and give a full turn in the key slot. The coin test consisted in subjects picking up four coins, one by one, and positioning them into a slot as fast as possible. The order of three subcomponents of test was randomized in each subject and across subjects. Measurements were taken 30 min after the end of each session. The distance and position between each subject's hand and the apparatus was recorded and maintained constant for preassessments and postassessments.
MRI.
Structural MRI was used to verify in SCI participants damage to C5–C6 spinal segments (n = 16). The images were acquired on a TIM Trio 3 T system using a 20-channel birdcage head/neck coil (Siemens). A T2-weighted turbo spin-echo sequence (TR = 3 s echo time, TE = 96 ms), with sagittal in-plane resolutions of 0.6 and 3 mm slice thickness, was used to evaluate spinal cord sections. Sagittal sections were used to determine the extension of the lesions in reference to the vertebral bodies, and each segment was limited by the position of the intervertebral discs and confirmed by two observers. The caudal extension of the injury was marked along the extension of each vertebra taking as reference superior and inferior borders. The position of C5–C6 was marked, and the distance from the center of the cord at this location to the surface of the skin as well as the distance to the subcutaneous (SC) fat was measured (Fig. 2A). Note that in SCI participants TESS was applied over injured spinal segments (Fig. 2B). One participant was unable to obtain the spinal cord MRI due to claustrophobia. Note that on this participant with a C4 neurological level of injury (participant 17 in Table 1) TESS electrodes were located between C5 and C6 spinous processes.
Data analysis.
Normal distribution was tested by the Shapiro–Wilk test and homogeneity of variances by the Mauchly's test of sphericity. When sphericity could not be assumed, the Greenhouse–Geisser correction statistic was used. Repeated-measures ANOVA was performed to determine the effect of the group (control subjects and SCI participants), time [before stimulation (Pre), and at 15, 30, 45, 60, and 75 min after stimulation], and intervention (TESS and sham-TESS) on the amplitude of MEPs, CMEPs, and mean rectified background EMG and the effect of group and time on SICI. Repeated-measures ANOVA was also performed to determine the effect of the group, time, and intervention frequency (TESS, sham-TESS, and TESS without 5 kHz) on MEPs, CMEPs, SICI, and mean rectified background EMG; and the effect of intervention frequency and task (coin, key, and bottle) on the time to task completion. Note that statistical analysis was completed using the raw data, and when the data were normalized to the baseline (using 100% values before the intervention) on each variable. Because the same results were found in both analyses, in the text we present the results using the raw data and in the figures we showed the data normalized to the baseline for clarity. Bonferroni post hoc tests were used to test for significant comparisons. Independent t tests were used to compare age, M-max, MVC, and stimulation intensities (TMS and electrical stimulation). t Test-corrected values were used when equal variances were not assumed. Pearson correlation analysis was used as needed. Significance was set at p < 0.05. Group data are presented as the mean ± SD in the text.
Results
MRI
Figure 2A shows sagittal slices obtained from the cervical level of three representative participants with SCI. Measurements were taken to measure the depth from the surface of the skin to the spinal canal. Blue arrows show the distance from the canal toward the beginning of the SC space. This part is composed mainly of bone forming the laminae and the spinous process (isointense T2-weighted signal) as well as paravertebral muscles (hypointense T2-weighted signal). Red arrows show the distance from the spinal canal toward the surface of the skin (D) and included the cutaneous layer and fat tissue underneath (hyperintense T2-weighted signal). The yellow dotted lines show that the top and bottom borders of the injury and the arrow points to the C5–C6 intervertebral space. The distance from the spinal cord to the surface of the skin was variable from participant to participant, the thickness of the fat layer was calculated by subtracting D − SC (SC = 39.5 ± 6.4 mm; D = 57.6 ± 13.2 mm; fat = 18.3 ± 12.4 mm). We found a moderate but not significant correlation between SC and the amount of current needed to observe a threshold root-evoked potential in the biceps brachii muscle (r = 0.42, p = 0.1). No correlations were found between D (r = 0.01 p = 0.1) or fat (r = 0.06 p = 0.8) and the amount of current needed to observe a threshold root-evoked potential in the biceps brachii muscle. Note that 13 participants showed injuries that extended more than one cervical segment, and all injuries were localized between C4 and C7 spinal segments (Fig. 2B).
Structural MRI. A, Three sagittal slices from representative participants with cervical SCI. The blue arrows show the distance from the canal toward the beginning of the SC space. The red arrows show the distance from the spinal canal toward the surface of the skin, including the cutaneous layer and fat tissue underneath. The yellow dotted lines show the top and bottom borders of the injury, and arrows point to the C5–C6 intervertebral space. Individual distances and segmental distributions of the injuries are shown in B.
CMEPs
Repeated-measures ANOVA showed an effect of intervention (F(1,29) = 91.5, p < 0.001) and time (F(5,145) = 20.1, p < 0.001), but not of group (F(1,29) = 2.4, p = 0.13) or in their interaction (F(5,145) = 0.9, p = 0.5) on biceps brachii CMEP amplitude. Additionally, there was an effect of intervention × time (F(5,145) = 21.9, p < 0.001) without any other interactions. Our analysis revealed that CMEP amplitude increased after TESS in both groups at 15 min (control subjects = 141.2 ± 54.0%, p < 0.05; SCI participants = 150.0 ± 48.4%, p < 0.01), 30 min (control subjects = 156.9 ± 33.0%, p < 0.001; SCI participants = 166.9 ± 53.2%, p < 0.001), 45 min (control subjects = 188.9 ± 49.4%, p < 0.001; SCI participants = 189.6 ± 52.9%; p < 0.001), 60 min (control subjects = 181.4 ± 41.4%, p < 0.001; SCI participants = 186.7 ± 59.2%, p < 0.001), and 75 min (control subjects = 153.2 ± 53.6%, p < 0.001; SCI participants = 147.1 ± 41.2%, p < 0.05) compared with baseline. In contrast, CMEP amplitude remained unchanged after sham TESS in both groups at 15 min (98.7 ± 13.5%; p = 0.4), 30 min (93.3 ± 11.3%; p = 0.1), 45 min (96.5 ± 9.3%; p = 0.1), 60 min (97.5 ± 9.4%; p = 0.2), and 75 min (100.0 ± 13.3%; p = 0.9) compared with baseline. Note that CMEPs were facilitated after TESS but not sham-TESS in all control and SCI participants, and representative examples and individual data are shown on Figure 3, A and B (top panels and lower graphs).
CMEPs. A, B, Averaged CMEP traces in the biceps brachii muscle in a representative control subject (A) and a participant with SCI (B) before and after TESS and sham-TESS. Each waveform represents the average of 10 CMEPs. Graphs show group (control subjects, n = 15; SCI, n = 17) and individual data. The ordinate shows the CMEP amplitude as a percentage of the CMEP at baseline (percentage of baseline) in control subjects and SCI participants before and after TESS (control subjects, green circles; SCI participants, orange circles) and sham-TESS (control subjects and SCI participants, black circles). The abscissa shows the time of measurements (Pre, and 15, 30, 45, 60, and 75 min after each protocol). Note that in graphs showing data from individual participant postmeasurement data showed the average from 15 to 75 min. Error bars indicate SDs. *p < 0.05.
We observed similar increases in CMEP amplitude after TESS in triceps brachii [effect of time (F(1.6,31.7) = 4.9, p < 0.01) not Group (F(1,20) = 1.9, p = 0.2) or in their interaction (F(1.6,31.7) = 1.1, p = 0.3)] and first dorsal interosseous [effect of time (F(5,70) = 6.1, p < 0.001) not group (F(1,14) = 4.1, p = 0.1) or in their interaction (F(5,70) = 1.3, p = 0.3)] muscles. The amplitude of CMEPs increased on average in the triceps brachii in control subjects (166.6 ± 46.5%; range, 100.6–243.6%) and SCI participants (150.9 ± 43.3%; range, 104.8–255.2%) compared with baseline for up to 75 min. Similarly, CMEPs in the first dorsal interosseous muscle were facilitated on average in control subjects (162.1 ± 47.8%; range, 113.3–243.9%) and SCI participants (130.6 ± 22.7%; range, 101.3–160.2%) compared with baseline measurements for up to 75 min. Note that CMEP amplitude in all muscles was similar at baseline in both groups (biceps brachii: control subjects = 4.6 ± 1.3% of the M-max; SCI participants = 5.2 ± 1.0% of the M-max; p = 0.1; triceps brachii: control subjects = 3.5 ± 2.6% of the M-max; SCI participants = 4.4 ± 2.8% of the M-max; p = 0.2; first dorsal interosseous: control subjects = 5.4 ± 0.8% of the M-max; SCI = 5.9 ± 0.7% of M-max; p = 0.1).
MEPs
Repeated-measures ANOVA revealed no effect of intervention (F(1,30) = 0.1, p = 0.8), time (F(5,150) = 1.4, p = 0.2), group (F(1,30) = 0.7, p = 0.4), or in their interaction (F(5,150) = 0.3, p = 0.9) on biceps brachii MEP amplitude. We observed no changes in MEP amplitude after TESS or sham-TESS up to 75 min in both groups tested (Fig. 4A,B). Similarly, we observed no effect of time (F(1.5,30.6) = 1.0, p = 0.4), group (F(1,20) = 0.1, p = 0.7), or in their interaction (F(1.5,30.6) = 1.9, p = 0.2) on triceps brachii MEP amplitude. In addition, there was an effect of group (F(1,12) = 11.3, p < 0.01) but not of time (F(5,60) = 1.2, p = 0.3) or in their interaction (F(5,60) = 0.4, p = 0.8) on first dorsal interosseous MEP amplitude.
MEPs. A, B, Averaged MEP traces elicited by TMS in the biceps brachii muscle in a representative participant before and after each protocol. Each waveform represents the average of 20 MEPs. Graphs show group (control subjects, n = 15; SCI participants, n = 17) and individual data. The ordinate shows the MEP amplitude as a percentage of the MEP at baseline (percentage of baseline) in control subjects and SCI participants before and after TESS (control subjects, green circles; SCI participants, orange circles) and sham-TESS (control subjects and SCI participants, black circles). The abscissa shows the time of measurements (Pre, and 15, 30, 45, 60, and 75 min after each protocol). Error bars indicate SDs. *p < 0.05.
SICI
Repeated-measures ANOVA showed an effect of time (F(5,145) = 15.5, p < 0.001) and group (F(1,29) = 10.8, p < 0.01) but not in their interaction (F(5,145) = 1.1, p = 0.4) on SICI. We found that SICI increased after 15 min (control subjects = 55.9 ± 30.9%, p < 0.001; SCI participants = 54.6 ± 39.1%, p = 0.02), 30 min (control subjects = 44.9 ± 37.8%, p < 0.001; SCI participants = 54.0 ± 40.0%, p < 0.01), 45 min (control subjects = 44.3 ± 34.1%, p < 0.001; SCI participants = 53.7 ± 24.0%, p < 0.05), 60 min (control subjects: 44.0 ± 25.1%, p < 0.001; SCI participants = 61.7 ± 39.0%, p < 0.05), and 75 min (control subjects = 55.5 ± 25.3%, p < 0.001; SCI participants = 58.6 ± 33.4%, p < 0.05) after TESS compared with baseline in control and SCI participants. Note that SICI was increased in all of subjects after the TESS protocol in the control group (14 of 14) and in most cases in SCI participants (13 of 15) when comparing baseline measurements with the average of measurements between 15 and 75 min after TESS (Fig. 5A,B, graphs).
SICI. A, B, Averaged traces showing SICI measurements in the biceps brachii muscle in a representative control subject (A) and SCI participant (B) before and after TESS. Each waveform represents the average of 15 test MEPs (black traces) and 15 conditioned MEPs (gray traces). Arrows at the beginning of each trace indicate the TS and CS used during testing. Graphs show group data (control subjects, n = 15; SCI participants, n = 17) and individual data. The ordinate shows the conditioned MEP as a percentage of the test MEP at baseline (percentage of test response) in control subjects and SCI participants before and after TESS (control subjects, green circles; SCI participants, orange circles). The abscissa shows the time of measurements (Pre, and 15, 30, 45, 60, and 75 min after each protocol). In graphs showing data from individual participants, postmeasurements show the average from 15 to 75 min. Error bars indicate SDs. *p < 0.05.
Effects without the carrier frequency
Repeated-measures ANOVA showed an effect of intervention frequency (F(2,28) = 40.1, p < 0.001) and time (F(5,70) = 17.1, p < 0.001), but not of group (F(1,14) = 0.1, p = 0.9) on biceps brachii CMEP amplitude. There was also an effect of intervention frequency × time (F(4.4,61.1) = 1.3, p < 0.001). Note that CMEP amplitude increased after the TESS and TESS without 5 kHz at 15 min (p < 0.02), 30 min (p < 0.001), 45 min (p < 0.001), 60 min (p < 0.001), and 75 min (p < 0.01) compared with baseline. In contrast, no changes in CMEP amplitude were observed after the sham-TESS session (p > 0.2; Fig. 6A,B, left column). In addition, repeated-measures ANOVA showed an effect of intervention frequency (F(1.2,16.7) = 48.8, p < 0.001) and time (F(5,70) = 10.7, p < 0.001) but not group (F(1,14) = 0.2, p = 0.6) on biceps brachii MEP amplitude. There was an effect of intervention frequency × time (F(3.5,49.6) = 10.9, p < 0.001). Notably, MEP amplitude increased after the TESS without 5 kHz at 15 min (p < 0.001), 30 min (p < 0.001), 45 min (p < 0.001), 60 min (p < 0.001), and 75 min (p < 0.001) compared with baseline but not after the TESS session (p > 0.2 at all time points) and sham-TESS session (p > 0.2 at all time points; Fig. 6A,B, middle column). These results together suggest that the 5 kHz carrier frequency had a suppressive effect on corticospinal excitability.
Effect of different currents. A–B, Graphs show group data (A. control subjects, n = 8; B. SCI participants, n = 8) for CMEPs MEPs and SICI before and after the TESS, sham-TESS, and TESS applied without the 5 kHz carrier frequency (TESS w/o 5 kHz) in control subjects and SCI participants, respectively. The ordinate shows CMEP and MEP amplitude as a percentage of CMEPs and MEPs at baseline (percentage of baseline) and SICI as a percentage of SICI at baseline (percentage of baseline) in control subjects and SCI participants before and after TESS (light gray circles), sham-TESS (black circles), and TESS without 5 kHz (dark gray circles). The abscissa shows the time of measurements (Pre, and 15, 30, 45, 60, and 75 min after each protocol). Error bars indicate SDs. *p < 0.05.
Therefore, we measured the effect of all these conditions on SICI. Repeated-measures ANOVA showed an effect of intervention frequency (F(1.4,21.1) = 26.1, p < 0.001) but not time (F(5,75) = 2.0, p = 0.1) or group (F(1,15) = 3.0, p = 0.1) on SICI. There was an effect of intervention frequency × time (F(4.9,74.3) = 4.4, p < 0.001). We found that SICI decreased after the TESS at 15 min (p < 0.001), 30 min (p < 0.001), 45 min (p < 0.001), 60 min (p < 0.001), and 75 min (p < 0.001) compared with baseline but not after the TESS without 5 kHz session (p > 0.2 at all time points) and sham-TESS session (p > 0.6 at all time points; Fig. 6A,B, right column), supporting the view the 5 kHz carrier frequency had a suppressive effect on cortical activity.
Functional outcomes
Repeated-measures ANOVA revealed an effect of intervention frequency (F(2,18) = 5.5, p < 0.05), task (F(2,18) = 8.8, p < 0.01), time (F(1,9) = 12.9, p < 0.01), and in their interaction (F(4,36) = 7.7, p < 0.001) on the time to complete the task (Fig. 7A,B). Additionally, there was an effect of intervention frequency × time (F(2,18) = 12.2, p < 0.001) and task × time (F(2,18) = 7.7, p < 0.01). Our results indicate that the mean time to complete all tasks decreased after TESS (coin = 38.1 ± 5.6%, p < 0.001; key = 41.9 ± 8.2%, p < 0.001; and bottle = 47.6 ± 6.3%; p < 0.001) and TESS without 5 kHz (coin = 22.3 ± 10.3%, p < 0.01; key = 23.6 ± 10.3%, p < 0.01; and bottle = 27.0 ± 7.4%, p < 0.001) compared with sham-TESS. Notably, the decrease in time was larger after TESS compared with TESS without 5 kHz in the coin (p < 0.01), key (p < 0.01), and bottle (p < 0.001) tasks. In contrast, no differences were observed in the time to complete tasks before and after sham-TESS (p = 0.1).
Functional outcomes. A, B, Images showing the three subcomponents of the GRASSP (coint, key, and water bottle tests, respectively) tested in SCI participants (n = 10) before and ∼75 min after sham-TESS, TESS w/o 5 kHz, and TESS. Graphs show individual (left side, A) and group (right side, B). The ordinate shows the time to complete each task as a percentage of the time needed at baseline (percentage of baseline) before and after TESS (light gray circles), sham-TESS (black circles), and TESS without 5 kHz (dark gray circles). The abscissa shows the subcomponents of the GRASSP tested (ccoint, key, and water bottle tests). Error bars indicate SDs. *p < 0.05.
Discussion
Our novel findings indicate that TESS can change the excitability of cortical and spinal neuronal networks in uninjured and SCI humans. The size of CMEPs but not MEPs increased in proximal and distal arm muscles for 75 min after TESS, but not sham-TESS, in control subjects and SCI participants. Paired stimuli-evoked intracortical inhibition increased after TESS in both groups. These results together suggest that TESS can affect cortical and spinal networks in parallel. When TESS was applied without the 5 kHz carrier frequency, both subcortical and cortical motor-evoked responses were facilitated without changing intracortical inhibition, suggesting that the 5 kHz carrier frequency contributed to the cortical inhibitory effects. Notably, improvements in hand and arm function after TESS was greatest when TESS was used with the 5 kHz carrier frequency. We hypothesized that parallel cortical and subcortical TESS effects can facilitate further recovery of limb function following SCI and suggest that knowledge of the neural basis of TESS-induced recovery may help to maximize its use in rehabilitation.
Cortical and subcortical TESS aftereffects
Accumulating evidence supports the view that TESS contributes to improve sensory and motor function in humans with SCI (Hofstoetter et al., 2013; Gerasimenko et al., 2015; Gad et al., 2017, 2018; Inanici et al., 2018). Here, for the first time, we examined the effects of TESS applied for several minutes, using parameters currently used in clinical settings, on the excitability of CNS pathways. We found that 20 min of TESS between C5 and C6 spinous processes increased the amplitude of CMEPs but not MEPs elicited by TMS over the primary motor cortex in proximal and distal arm muscles in control subjects and SCI participants. CMEPs likely result from motoneuron activation by descending volleys elicited by excitation of corticospinal axons (Rothwell et al., 1994), reflecting changes in the efficacy of cortico–motoneuronal synapses or motoneuron excitability (Ugawa et al., 1991; Taylor and Gandevia, 2004). CMEP size increases markedly during a voluntary contraction, without changing its latency, which could reflect the excitability of the motoneuron pool linked to monosynaptic connections from corticospinal axons (Ugawa et al., 1995; Taylor et al., 2002). Thus, it is possible that neuromodulation of spinal networks by TESS boosts the excitability of motoneurons and/or cortico–motoneuronal synapses. This is consistent with evidence from electrophysiological (Minassian et al., 2007b; Hofstoetter et al., 2018) and modeling (Rattay et al., 2000; Capogrosso et al., 2013) studies, indicating that a single TESS pulse recruits large-to-medium-diameter proprioceptive and cutaneous afferents within posterior rootlets/roots of different spinal segments. This is also consistent with results from other paradigms showing that 20 min of stimulation targeting the spinal cord is sufficient time to change the excitability of spinal networks (Taylor and Martin, 2009; Bunday and Perez, 2012). The lack of changes in the amplitude of MEPs, while the amplitude of CMEPs increased, was surprising because collision studies using electrical (Ugawa et al., 1991) and magnetic stimulation over the primary motor cortex (Berardelli et al., 1991; Gandevia et al., 1999) indicate that the two stimuli activate similar corticospinal axons. A possibility is that TESS activated inhibitory cortical circuits projecting onto corticospinal neurons. Animal studies showed that afferent input can influence motor cortical activity via dense intracortical projections between the primary motor cortex and the somatosensory cortex (Goldring et al., 1970). Indeed, stimulation of afferent fibers from the somatosensory cortex can produce both excitation and inhibition in motor cortical cells (Porter et al., 1990). To examine the contribution of cortical networks, we assessed SICI by using a randomized paired-pulse TMS paradigm (Kujirai et al., 1993). Studies using pharmacological agents and epidural recordings indicated that SICI reflects the activation of GABA inhibitory circuits in the primary motor cortex (Di Lazzaro et al., 2000). We found that SICI increased after TESS in both groups, supporting the view that increased intracortical inhibition could represent a mechanism contributing to the lack of changes in MEP amplitude. In agreement, studies in humans using TMS over the primary motor cortex showed that electrical stimulation of afferent fibers innervating upper limb muscles suppresses corticospinal excitability at interstimulus intervals (ISIs) of ∼20–40 ms (Lei and Perez, 2017) and facilitates corticospinal excitability at ISIs of ∼50–100 ms (Devanne et al., 2009). Afferent inputs reach the primary motor cortex in humans no more than 4 ms later than the primary somatosensory cortex (Goldring et al., 1970). Thus, even if after SCI conduction velocity in sensory axons is slower (Tan et al., 2007) and the latency of somatosensory-evoked potentials increased (Ozdemir and Perez, 2018), the timing of arrival of afferent information to the primary somatosensory cortex could be sufficient to exert effects on intracortical interneurons projecting to corticospinal neurons controlling upper limb muscles. In addition, all our participants were able to perceive sensory stimuli in dermatomes of the arm and hand tested, indicating the presence of an anatomical substrate to influence cortical inhibitory networks.
In most studies aiming to achieve therapeutic effects in SCI subjects with TESS, a frequency of 5–50 Hz with a high carrier frequency between 5 and 10 kHz has been applied for several minutes to hours. To further understand the origin of the cortical inhibitory effects induced by TESS, we applied TESS without the 5 kHz carrier frequency. Here, we found that both MEPs and CMEPs were facilitated without changing intracortical inhibition, suggesting that the 5 kHz carrier frequency contributed to the cortical inhibitory effects. This agrees with evidence showing that repeated stimulation of afferent fibers in a peripheral nerve with frequencies between 5 and 50 Hz usually increases corticospinal excitability (Kaelin-Lang et al., 2002; Golaszewski et al., 2012). Why did the 5 kHz carrier frequency contribute to the cortical inhibitory effects? High carrier frequencies or Russian currents are alternating currents normally within the kilohertz frequency, which are delivered in bursts, with a burst frequency within a physiological range (usually up to ∼100 Hz; Ward, 2009). It is thought that this high carrier frequency is beneficial for improving muscle strength (Selkowitz, 1985, 1989) and for suppressing the sensitivity of pain receptors (Ward and Robertson, 1998a,b). High-frequency electrical stimulation of the spinal cord releases serotonin into the dorsal horn (Linderoth et al., 1993), which may depress ascending nociceptive transmission (Fürst, 1999; Millan, 2002). This is also consistent with evidence showing that electrical stimulation of skin afferents at frequencies of ∼100 Hz has inhibitory effects on corticospinal excitability (Mima et al., 2004).
It is interesting that TESS increased inhibition within the primary motor cortex while at the same time improved motor performance in SCI participants. Improved performance after repeated practice and motor learning has been associated with increases in intracortical inhibition in a few cases (Perez et al., 2007). For example, in this study, motor practice involved training participants to improve balanced activation of antagonistic muscles during a cocontraction task. Muscle coactivation ratios between antagonistic muscles of the arm are impaired in humans with SCI compared with control participants (Stahl et al., 2015; Calabro and Perez, 2016). Our observation that TESS applied over C5–C6 spinous processes had similar physiological effects on biceps and triceps brachii muscles suggests that TESS reached neuronal networks controlling antagonistic muscles. This also agrees with evidence showing that electrical stimulation of afferents innervating a flexor muscle could affect intracortical networks contributing to control of the antagonistic extensor muscle (Bertolasi et al., 1998). Thus, changes in muscle coactivation might be a possible factor contributing to our results.
Functional significance
Our results indicate that a single 20 min session of TESS at rest is sufficient to change the excitability in neuronal pathways and motor performance as other neuromodulatory paradigms targeting the spinal synapses (Bunday and Perez, 2012) or the primary motor cortex (Roy et al., 2010) in humans with SCI. Previous findings have shown that the effect of TESS on the excitability of neuronal pathways and motor performance are generally similar when measurements are taking before and after the intervention with or without the cervical stimulation after as well as during cervical stimulation in humans (Gad et al., 2018) and rats (Alam et al., 2017). Benefits in motor performance have been reported in humans with SCI when other neuromodulation paradigms, using electrical and/or magnetic stimulation, are used before (Lu et al., 2016) or concurrently (Hoffman and Field-Fote, 2007; Potter-Baker et al., 2018) with motor training. Some aspects of the timing of the application of TESS relative to performing the motor task, however, might be important to consider. The responsiveness of the motoneuron pool during voluntary activity also decreases following SCI (Vastano and Perez, 2020), which might benefit from extra inputs recruited by TESS. Our results showing effects on muscles innervated by segments distal to the applied stimulation site suggest that TESS effects extend below targeted segments. This is consistent with previous finding showing that TESS effects spread to adjacent proximal and distal motor pools to the stimulated site after a single pulse (Sayenko et al., 2015). It is likely that the effects of TESS on the interneurons and motor pools projecting to distal hand muscles are the consequence of stimulating afferents from more distal segments, as well as via the high degree of propriospinal intersegmental connectivity and/or to leakage of the current to lower segments due to electrode settings. Thus, a clearer understanding on the neural connectivity underlying the TESS-induced recovery that has been observed may help to maximize rehabilitative strategies.
Footnotes
M.A.P. received funding from National Institute of Neurological Disorders and Stroke, the Veterans Administration, and the Walkabout Foundation. H.L. received funding from the European Research Council under the European Union's Horizon 2020 Research and Innovation Programme.
The authors declare no competing financial interests.
- Correspondence should be addressed to Monica A. Perez at mperez04{at}sralab.org