Abstract
Humans with cervical spinal cord injury (SCI) often recover voluntary control of elbow flexors and, to a much lesser extent, elbow extensor muscles. The neural mechanisms underlying this asymmetrical recovery remain unknown. Anatomical and physiological evidence in animals and humans indicates that corticospinal and reticulospinal pathways differentially control elbow flexor and extensor motoneurons; therefore, it is possible that reorganization in these pathways contributes to the asymmetrical recovery of elbow muscles after SCI. To test this hypothesis, we examined motor-evoked potentials (MEPs) elicited by transcranial magnetic stimulation over the arm representation of the primary motor cortex, maximal voluntary contractions, the StartReact response (a shortening in reaction time evoked by a startling stimulus), and the effect of an acoustic startle cue on MEPs elicited by cervicomedullary stimulation (CMEPs) on biceps and triceps brachii in males and females with and without chronic cervical incomplete SCI. We found that SCI participants showed similar MEPs and maximal voluntary contractions in biceps but smaller responses in triceps compared with controls, suggesting reduced corticospinal inputs to elbow extensors. The StartReact and CMEP facilitation was larger in biceps but similar to controls in triceps, suggesting enhanced reticulospinal inputs to elbow flexors. These findings support the hypothesis that the recovery of biceps after cervical SCI results, at least in part, from increased reticulospinal inputs and that the lack of these extra inputs combined with the loss of corticospinal drive contribute to the pronounced weakness found in triceps.
SIGNIFICANCE STATEMENT Although a number of individuals with cervical incomplete spinal cord injury show limited functional recovery of elbow extensors compared with elbow flexor muscles, to date, the neural mechanisms underlying this asymmetrical recovery remain unknown. Here, we provide for the first time evidence for increased reticulospinal inputs to biceps but not triceps brachii and loss of corticospinal drive to triceps brachii in humans with tetraplegia. We propose that this reorganization in descending control contributes to the asymmetrical recovery between elbow flexor and extensor muscles after cervical spinal cord injury.
Introduction
Spinal cord injury (SCI) often results in decreased voluntary control of muscles below the level of injury. In individuals with cervical lesions, limited functional recovery has been reported in elbow extensors compared with elbow flexor muscles (Ditunno et al., 1992, 2000; Calancie et al., 2004; McKay et al., 2011). This is interesting considering that the location of corticospinal projections (Penfield and Boldrey, 1937) and motoneurons (Schirmer et al., 2011) that innervate these muscles largely overlap, and that most injuries to the cervical spinal cord affect more than one spinal segment (Schaefer et al., 1989; Benavides et al., 2020). Loss of motoneurons near the injury site (Grumbles and Thomas, 2017), alterations in motor unit recruitment patterns and muscle atrophy (Thomas et al., 1997a,b, 1998; Johanson et al., 2013) might contribute to the poor recovery in triceps brachii. However, it is likely that other factors are involved in the asymmetrical recovery between elbow flexor and extensor muscles after SCI.
Anatomical and electrophysiological studies have shown that elbow flexors and extensors receive differential control from corticospinal and reticulospinal pathways. For example, the terminal distribution of crossed and ipsilateral corticospinal projections to interneurons and motoneurons is larger in biceps compared with triceps brachii (Morecraft et al., 2013). In monkeys (Phillips and Porter, 1964; McKiernan et al., 1998) and humans (Palmer and Ashby, 1992; de Noordhout et al., 1999), elbow flexor motoneurons receive large corticospinal monosynaptic facilitation, whereas triceps motoneurons receive lesser monosynaptic facilitation but strong disynaptic inhibition. Reticulospinal axons are bilaterally distributed with a majority of ipsilateral axons in monkeys and humans (Nathan et al., 1996; Sakai et al., 2009). Stimulation of the reticular formation in monkeys produced biceps facilitation and triceps suppression ipsilaterally (Davidson and Buford, 2006; Davidson et al., 2007; Herbert et al., 2010; Hirschauer and Buford, 2015), consistent with findings in humans showing that ipsilateral cortico-reticulospinal responses are more easily observed in biceps compared with triceps brachii (Ziemann et al., 1999).
Both corticospinal (Oudega and Perez, 2012) and reticulospinal (Asboth et al., 2018) pathways undergo reorganization after SCI, which might then contribute to the asymmetrical recovery reported in elbow extensors compared with elbow flexor muscles. Corticospinal inputs are largely decreased in triceps brachii in humans with incomplete cervical SCI compared with control subjects (Thomas et al., 1997b). Corticospinal lesions in primates led to enhanced monosynaptic and disynaptic reticulospinal facilitation to forearm flexors compared with extensor muscles (Zaaimi et al., 2012), and reticulospinal neurons may have a greater response to regeneration approaches than corticospinal neurons (Vavrek et al., 2007; Zörner et al., 2014). Thus, we hypothesized that reduced corticospinal inputs to triceps brachii and increased reticulospinal inputs to biceps brachii contribute to the asymmetrical recovery observed in people with cervical incomplete SCI.
To test this hypothesis, we examined motor-evoked potentials (MEPs) elicited from the arm representation of the primary motor cortex, maximal voluntary contractions (MVCs), the StartReact response (a shortening in reaction time evoked by a startling stimulus), and the effect of an acoustic startle cue on an MEP elicited by cervicomedullary stimulation (cervicomedullary MEPs [CMEPs]) in biceps and triceps brachii of individuals with and without cervical SCI.
Materials and Methods
Subjects
Eighteen individuals with SCI (mean age 43.3 ± 16.1 years, 2 female; Table 1) and 15 age-matched controls (mean age 40.1 ± 14.2 years, 4 female, p = 0.7) participated in the study. All participants gave informed consent to the experimental procedures, which were approved by the local ethics committee at Northwestern University and performed in accordance with 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 C1-C6 SCI and by the American Spinal Cord Injury Association Impairment Scale (AIS) as AIS A (n = 4), AIS C (n= 10), and AIS D (n= 4). We tested individuals with preservation of motor output in the biceps and triceps brachii to assess the contribution of residual descending motor pathways. All subjects were asked to perform voluntary elbow flexion and extension, and all of them were able to exert voluntary EMG activity in the biceps and triceps muscles. Nine SCI individuals were under antispastic medication (baclofen and/or gabapentin and/or tizanidine; Table 1). However, none of them showed spasticity in their biceps and triceps brachii, as examined with the Modified Ashworth Scale (MAS = 0).
SCI participantsa
EMG recordings
EMG was recorded from the right biceps and triceps brachii in control subjects and from the less affected side in individuals with SCI through bipolar surface electrodes (Ag-AgCl, 10 mm diameter, 1 cm apart) secured on the skin. Signals were amplified, filtered (30-2000 Hz), and sampled at 2 kHz for offline analysis (CED 1401 with Signal software, Cambridge Electronic Design).
Experimental paradigm
During testing, subjects were seated comfortably in an armchair with their arm supported by a custom platform. The tested arm was positioned with the shoulder and the elbow flexed at 90° and strapped to the platform. During MVCs, subjects were instructed to perform an isometric elbow flexion for biceps MVC and an elbow extension for triceps MVC. Subjects performed three brief MVCs for 3-5 s for both elbow flexion and elbow extension, separated by ∼1 min 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. MVCs were tested at the beginning of the session followed by examination of MEP recruitment curves, M-max, StartReact, and CMEPs.
MEP recruitment curves
Transcranial magnetic stimulation (TMS) was delivered over the arm representation of the primary motor cortex from a BiStim2 (Magstim) through a circular coil with a monophasic current waveform. The coil was positioned at the optimal site for eliciting an MEP in the contralateral biceps and triceps muscles during a small tonic voluntary contraction of ∼10% of the MVC (biceps: controls = 7.0 ± 0.5% of the MVC, SCI = 7.1 ± 0.9% of the MVC, U = 124.0, p = 0.7, r = 0.07; triceps: controls = 7.1 ± 0.6% of the MVC, SCI = 7.9 ± 2.4% of the MVC, U = 106.0, p = 0.3, r = 0.2). Participants wore a cap on which the position of the coil was marked to ensure the stability of TMS across measurements. The TMS intensity ranged from below the active motor threshold (AMT) to the intensity needed to produce a maximal MEP (MEP-max), increasing by 5% increments until the maximal stimulator output (MSO) was reached. The AMT was defined as the minimal stimulus intensity needed to elicit 5 of 10 MEPs 100 µV above the EMG background. Ten stimuli (interstimulus interval: 0.2 Hz) were delivered at each intensity to plot the mean peak-to-peak amplitude of the MEP from the nonrectified response against the TMS intensity in each subject (MEP recruitment curve; Fig. 1A). The experimental data were fitted with the following 3-parameter sigmoid function (Devanne et al., 1997; Carroll et al., 2001):
where MEP-max is the maximal amplitude of the MEP, I50 is the TMS intensity producing half MEP-max, and k is the Boltzmann slope parameter (Fig. 1A). Sigmoid fit was individually performed for each subject. The estimated values of each parameter of the sigmoid (MEP-max, I50 and k) were extracted for each participant, and the mean of each parameter was calculated for both control and SCI participants. The mean estimated MEP-max, I50, and k were used to draw the mean estimated sigmoid curve for each group. MEP-max onset latencies were defined as the time when the rectified EMG reached 2×SD calculated over a 100 ms period of the prestimulus activity of the MEP with the largest amplitude in each participant. Percutaneous electrical stimulation of the brachial plexus was delivered (1 ms rectangular electrical stimulus, DS7AH, Digitimer) through a cathode (10-mm-diameter Ag-AgCl electrode) placed over the Erb point and an anode (Ag-AgCl plaque) placed over the shoulder. The M-max was obtained by increasing the stimulus intensity until no increases in the M-wave amplitude were observed. Two pulses were used at each stimulus intensity (interstimulus interval: 0.2 Hz). The M-max was measured as the peak-to-peak amplitude of the nonrectified response, and it was used to normalize MEP values in each participant.
Experimental setup. A, TMS was applied over the arm representation of the primary motor cortex to activate corticospinal neurons projecting directly or indirectly to biceps and triceps brachii motoneurons located between the fifth and the eighth cervical segment (C5-C8) to elicit an MEP. MEP recruitment curve was constructed by plotting the amplitude of the MEP against the TMS intensity and allow to define the MEP-max, I50, AMT, and the slope of the curve. B, During the StartReact response in some trials, an LED was presented with either an auditory acoustic stimulus or a SAS. The StartReact response was measured between the VRT (defined as the time from cue to onset of the EMG burst in the biceps and triceps brachii after the LED presentation), the VART (defined as the time delay between the presentation of the auditory acoustic stimulus and the EMG onset), and the VSRT (defined as the time between the SAS and the EMG onset).
StartReact response
The StartReact response was tested using a previously described paradigm in humans with SCI (Baker and Perez, 2017). Here, participants were asked to observe an LED located ∼1 m in front of the participant's head. When the LED was illuminated (20 ms), individuals were asked to perform an isometric elbow flexion or an isometric elbow extension as fast as possible. We measured the visual reaction time (VRT; Fig. 1B) as the time from cue to onset of the EMG burst in the targeted muscle after the LED presentation. The reaction time was defined as the time point when the rectified EMG signal exceeded 2×SD calculated over a 100 ms period of the prestimulus activity. In some trials, the LED was presented with either a quiet acoustic stimulus (80 dB, 500 Hz, 50 ms) or a startling acoustic stimulus (SAS, 120 dB, 500 Hz, 50 ms) delivered through two audio speakers (T-15, Polk Audio) located right behind the participant's head. The loud intensity evoked a clear startle in control subjects and in some individuals with SCI on initial presentation. The time delay between the presentation of the quiet acoustic stimulus and the onset of the EMG response was defined as the visual-auditory reaction time (VART; Fig. 1B), whereas the time between the SAS and the EMG onset was defined as the visual-startle reaction time (VSRT; Fig. 1B). A familiarization trial consisting of three repetitions of each task responding to the LED was completed at the beginning of each experiment to ensure that participants were able to complete all tasks and to familiarize them with the startling cue (Fisher et al., 2013). In each task, 20 responses (interstimulus interval: 0.2 Hz) were recorded for each condition (VRT, VART, and VSRT) in a randomized order. Data were measured trial by trial. Reaction times exceeding 700 ms were excluded. The VART and VSRT are both mediated via the cochlear nuclei, but only the high-intensity sound of the VSRT activates the reticulospinal pathway (Davis et al., 1982; Brown et al., 1991; Valls-Solé et al., 1999). To estimate changes in the gain of the reticulospinal output, we normalized the data as follows (Baker and Perez, 2017):
where ΔTSR (i.e., shortening of reaction time with a startle stimulus) is the shortening effect of a SAS on the VRT, and ΔTAR (i.e., shortening of reaction time with an auditory stimulus) measures the shortening of reaction time provided by a nonstartling auditory stimulus on the VRT, which presumably does not activate reticulospinal pathways.
CMEPs
The corticospinal tract was stimulated by high-voltage electrical current (200 μs duration, DS7AH, Digitimer) passed between two electrodes (Ag-AgCl, 10 mm diameter) fixed to the skin behind each mastoid process at the cervicomedullary level (Taylor and Gandevia, 2004). The effect of an SAS on CMEP size was tested using previously described methods (Furubayashi et al., 2000). It has been demonstrated that an acoustic startle cue preceding an MEP by 80 ms increases the size of the MEP at a subcortical level (Furubayashi et al., 2000; Kühn et al., 2004; Ilic et al., 2011; Tazoe and Perez, 2017). A test stimulus (CMEP-test) was delivered at an intensity that elicited a CMEP in the biceps or triceps muscles of ∼1 mV during a tonic contraction of the respective muscle at ∼10% of the MVC (biceps: controls = 8.4 ± 1.2% of the MVC, SCI = 8.0 ± 1.1% of the MVC, t(22) = 0.8, p = 0.4, d = 0.3; triceps: controls = 8.4 ± 1.2% of the MVC, SCI = 8.6 ± 1.0% of the MVC, t(22) = 0.4, p = 0.7, d = 0.2). A conditioning stimulus (SAS, 120 dB, 500 Hz, 50 ms; CMEP-startle) was presented 80 ms before a test stimulus through two audio speakers (T-15, Polk Audio) located right behind the participant's head. Three consecutive acoustic stimuli were presented at the beginning of each experiment to ensure that participants were familiarized with the startling cue. Ten responses (interstimulus interval: 0.25 Hz) were recorded in each condition (CMEP-test, CMEP-startle) in a randomized order. Note that similar results were found when we used an interstimulus interval of 0.05 Hz. CMEPs were measured in 12 control and 12 SCI participants randomly selected from each group. The CMEP was measured as the peak-to-peak amplitude of the nonrectified response. The CMEP onset latencies were defined when the rectified EMG reached 2×SD calculated over a 100 ms period of the prestimulus activity.
Data analysis
Normal distribution was tested by the Shapiro-Wilk's test and homogeneity of variances was tested by the Levene's test of equality and Mauchly's test of sphericity. When sphericity could not be assumed, the Greenhouse-Geisser correction statistic was used. Two-way ANOVA was performed to determine the effect of group (control, SCI) and muscle (biceps, triceps) on MVCs, CMEPs, and the reticulospinal gain. Repeated-measures ANOVAs as a mixed model were also performed to determine the effect of group and condition (VRT, VART, VSRT) on reaction times and on the mean EMG activity measured over 100 ms before stimulus and during EMG burst. Condition was used as a repeated-measures factor. Additional repeated-measures ANOVAs were performed on each group separately as needed. The Holm-Sidak adjustment was used to correct all post hoc comparisons. A Student's t test was performed to compare the difference between group or muscle on the voluntary level of contraction exerted during MEP recruitment curves and CMEPs, the MEP-max, slope, AMT, M-max, and between SCI individuals taking and not taking antispastic medications. Pearson correlation coefficient analysis was used as needed. If data were not normally distributed, a log transformation was applied. Normality was retested after the log transformation: if normality was confirmed, a parametric analysis was performed; otherwise, a nonparametric analysis (Mann–Whitney U test) was completed. Statistical analyses were conducted using SigmaPlot (Systat Software), with the significance set at p < 0.05. Group data are presented as mean ± SD.
Results
MVC
Figure 2A illustrates the rectified EMG activity in biceps and triceps during an MVC in a control (blue) and a SCI (orange) participant. The SCI participant showed similar MVC in biceps but smaller MVC in triceps compared with the control subject. Two-way ANOVA showed an effect of group (F(1,62) = 14.6, p < 0.001) and muscle (F(1,62) = 36.4, p < 0.001) and in their interaction (F(1,62) = 8.1, p = 0.006) on MVC (Fig. 2B). Post hoc tests showed that MVCs were reduced in triceps compared with biceps in control (biceps = 1.0 ± 0.5 mV, triceps = 0.4 ± 0.2 mV, t(28) = 2.2, p = 0.03, d = 1.6) and SCI (biceps = 0.8 ± 0.4 mV, triceps = 0.2 ± 0.2 mV, t(34) = 6.6, p < 0.001, d = 1.9) participants. Post hoc tests also showed that biceps MVCs were not different between control and SCI groups (t(31) = 0.7, p = 0.5, d = 0.4), but that triceps MVCs were smaller in SCI compared with controls (t(31) = 4.7, p < 0.001, d = 1.0). Within the SCI group, MVC values remained similar between participants taking and not taking antispastic medications (biceps, t(16) = 1.0, p = 0.3, d = 0.5; triceps t(16) = 0.3, p = 0.8, d = 0.1).
MVCs. A, EMG recorded during biceps and triceps brachii MVCs in a control subject (blue) and a participant with SCI (orange). The SCI participant exhibited similar biceps MVC, but smaller triceps MVC, compared with the control subject. B, Box plot charts represent the group data. Abscissa indicates the groups tested (blue bars represent the control group; orange bars represent the SCI group). Ordinate indicates the MVC (in millivolts). Top and bottom lines of the box indicate the 75th percentile (top quartile) and 25th percentile (bottom quartile), respectively. Line in the box indicates the 50th percentile (median). The two bars extend from the maximum and minimum value. Dot within the box indicates the arithmetical mean. *p < 0.05, ***p < 0.001.
MEP recruitment curves
Figure 3A illustrates the mean biceps brachii MEP recruitment curves in control (blue) and SCI (orange) groups. The x axis of the graph shows the intensity of TMS and the y axis shows the MEP size normalized to the M-max. Control and SCI groups showed similar MEP-max and slope. In biceps brachii, no difference was found for the MEP-max between control (47.0 ± 24.3% of the M-max) and SCI (45.4 ± 25.7% of the M-max, U = 134.0, p = 0.9, r = 0.007) groups. There was also no difference found for the slope of the MEP recruitment curve between control (1.3 ± 0.9) and SCI (1.3 ± 0.9, U = 126.0, p = 0.8, r = 0.06) groups. Since MEPs were normalized to the M-max, we compared these responses across groups. Differences were found for the M-max of the biceps brachii between control (17.7 ± 6.2 mV) and SCI (13.4 ± 5.5 mV, U = 77.0, p = 0.04, r = 0.4) groups. No difference was found in the MEP-max and the slope between control (MEP-max = 8.6 ± 5.5 mV; slope = 0.2 ± 0.2) and SCI (MEP-max = 7.2 ± 5.9 mV, U = 109.0, p = 0.4, r = 0.2; slope = 0.2 ± 0.3, U = 115.0, p = 0.5, r = 0.1) groups, suggesting that normalization procedures did not influence our results. The latency of the MEP-max was delayed in SCI (12.3 ± 1.6 ms) compared with controls (11.3 ± 0.7 ms, t(31) = 2.3, p = 0.03, d = 0.8), consistent with the view that corticospinal drive is impaired in SCI individuals (Thomas et al., 1997b; Calancie et al., 1999; Perez, 2012). No difference was found in the AMT between control (41.5 ± 8.6% of the MSO) and SCI (45.6 ± 5.6% of the MSO, t(31) = 1.6, p = 0.1, d = 0.5) groups. Within the SCI group, no difference was found for the MEP-max (t(16) = 0.9, p = 0.4, d = 0.4), the slope (t(16) = 1.2, p = 0.2, d = 0.6), and the M-max (t(16) = 1.3, p = 0.2, d = 0.6) between participants taking and not taking antispastic medications.
MEP recruitment curves. Mean MEP recruitment curves in biceps (A) and triceps (B) brachii of control (blue) and SCI (orange) groups. The x axis indicates the intensity of TMS and the y axis indicates the MEP size normalized to the M-max. In biceps, the SCI group showed similar MEP-max and slope compare with the control group while in triceps, the SCI group showed reduced MEP-max and slope compare with the control group. Box plot charts represent the group data. Abscissa indicates the groups tested (blue bars represent the control group; orange bars represent the SCI group). Ordinate indicates the MEP-max as a percentage of the M-max (left) and the slope (right). Top and bottom lines of the box indicate the 75th percentile (top quartile) and 25th percentile (bottom quartile), respectively. Line in the box indicates the 50th percentile (median). The two bars extend from the maximum and minimum value. Dot within the box indicates the arithmetical mean. *p < 0.05.
Figure 3B illustrates the mean triceps MEP recruitment curves in control (blue) and SCI (orange) groups. The x axis of the graph shows the intensity of TMS and the y axis shows the MEP size normalized to the M-max. The SCI group showed reduced MEP-max and slope compared with controls. In triceps, we found that the MEP-max was reduced in SCI (14.8 ± 8.7% of the M-max) compared with controls (27.0 ± 20.9% of the M-max, t(31) = 2.3, p = 0.03, d = 0.8). The slope of the recruitment curve was also found reduced in SCI (0.4 ± 0.2) compared with controls (0.8 ± 0.7, t(31) = 2.5, p = 0.02, d = 0.8). Since MEPs were normalized to the M-max, we compared these responses across groups. We found that the M-max of the triceps brachii was reduced in SCI (4.9 ± 2.6 mV) compared with controls (10.6 ± 4.4 mV, U = 32.5, p < 0.001, r = 0.6). We also found that the MEP-max and the slope were reduced in SCI (MEP-max = 0.8 ± 0.6 mV; slope = 0.02 ± 0.02) compared with controls (MEP-max = 3.0 ± 3.1 mV, U = 27.0, p < 0.001, r = 0.7; slope = 0.08 ± 0.07, U = 35.0, p < 0.001, r = 0.6), suggesting that the normalization procedures did not influence our results. We found that the latency of the MEP-max was delayed in SCI (14.4 ± 2.9 ms) compared with controls (11.3 ± 1.1 ms, U = 36.5, p < 0.001, r = 0.6), consistent with the view that corticospinal drive is impaired in SCI individuals. We found that the AMT was increased in SCI (56.0 ± 10.3% of the MSO) compared with controls (47.9 ± 9.9% of the MSO, t(31) = 2.3, p = 0.03, d = 0.8). Within the SCI group, no difference was found for the MEP-max (t(16) = 1.6, p = 0.1, d = 0.7), the slope (t(16) = 1.5, p = 0.2, d = 1.0), and the M-max (t(16) = 1.5, p = 0.2, d = 0.7) between participants taking and not taking antispastic medications.
Because triceps was weaker than biceps in both groups, in an additional control experiment, we asked all control and SCI participants to adjust their contraction to match biceps and triceps EMG levels during the acquisition of the MEP recruitment curves (controls: biceps = 0.03 ± 0.02 mV, triceps = 0.03 ± 0.02 mV, U = 108.0, p = 0.9, r = 0.03; n = 15; SCI: biceps = 0.01 ± 0.01 mV, triceps = 0.01 ± 0.01 mV, U = 150.0, p = 0.7, r = 0.06; n = 18). In the control group, we found no differences for the MEP-max and the slope between biceps (MEP-max = 34.6 ± 23.7% of the M-max; slope = 0.7 ± 0.6) and triceps (MEP-max = 26.8 ± 20.6% of the M-max, U = 87.0, p = 0.3, r = 0.2; slope = 0.8 ± 0.7, U = 111.0, p = 0.9, r = 0.01). In the SCI group, we found a reduction of the MEP-max and the slope in triceps (MEP-max = 14.6 ± 8.6% of the M-max; slope = 0.4 ± 0.2) compared with biceps (MEP-max = 31.6 ± 24.0% of the M-max, U = 94.0, p = 0.03, r = 0.4; slope = 0.8 ± 0.6, U = 95.0, p = 0.03, r = 0.4). This information confirms that triceps was more affected than biceps in the SCI group. Because triceps was weaker in SCI than in control participants, in a second additional experiment, we asked 12 SCI and 12 control participants to match their triceps EMG level during the acquisition of the MEP recruitment curves (controls = 0.03 ± 0.02 mV; SCI 0.02 ± 0.02 mV, U = 41.0, p = 0.1, r = 0.4). We found that the MEP-max and the slope of the recruitment curve were larger in controls (MEP-max = 31.8 ± 20.8% of the M-max; slope = 1.0 ± 0.7) compared with SCI participants (MEP-max = 13.9 ± 9.1% of the M-max, U = 32.0, p = 0.02, r = 0.5; slope = 0.4 ± 0.3, U = 31.0, p = 0.02, r = 0.5). This also confirms our previous results showing that triceps is more affected in SCI compared with controls, supporting the view that the level of voluntary contraction does not affect the MEP-max value (Devanne et al., 1997).
StartReact response
Figure 4A illustrates the mean EMG activity related to the VRT (black), VART (gray), and VSRT (red) in the biceps and triceps of a control subject and a SCI participant. Reaction times in the biceps brachii were similar between control and SCI participants and that reaction times in the triceps were prolonged in all conditions in the SCI compared with the control participant. Notably, in biceps, the shortening in reaction time with a SAS was similar between control and SCI participants but smaller with an acoustic stimulus in SCI participant. In triceps, the shortening of reaction time with an SAS and an acoustic stimulus was similar between control and SCI participants. Repeated-measures ANOVA showed no effect of group (F(1,31) = 1.8, p = 0.2) but an effect of condition (F(2,62) = 273.9, p < 0.001) and in their interaction (F(2,62) = 3.6, p = 0.03) on biceps brachii reaction times. Post hoc tests showed that VRT and VSRT were not different between control (VRT: 185.0 ± 24.1 ms; VSRT: 117.7 ± 24.0 ms) and SCI groups (VRT: 196.0 ± 39.5 ms, t(31) = 0.7, p = 0.5, d = 0.3; VSRT: 126.8 ± 33.8 ms, t(31) = 1.0, p = 0.3, d = 0.3). Post hoc tests also showed that VART was prolonged in SCI (163.8 ± 40.8 ms) compared with controls (140.3 ± 26.9 ms, t(31) = 2.1, p = 0.04, d = 0.7). Repeated-measures ANOVA showed an effect of group (F(1,31) = 7.8, p = 0.009) and condition (F(2,62) = 211.6, p < 0.001) but not in their interaction (F(2,62) = 0.1, p = 0.9) on triceps reaction times. Post hoc tests showed that VRT, VSRT, and VART were prolonged in SCI (VRT: 253.1 ± 61.7 ms; VSRT: 167.8 ± 59.1 ms; VART: 194.1 ± 64.4 ms) compared with controls (VRT: 198.2 ± 28.7 ms, t(31) = 2.7, p = 0.01, d = 1.1; VSRT: 127.0 ± 21.7 ms, t(31) = 2.8, p = 0.008, d = 0.9; VART: 151.1 ± 22.1 ms, t(31) = 2.5, p = 0.02, d = 0.9).
StartReact. A, Mean EMG activity related to the VRT (black), VART (gray), and VSRT (red) in the biceps and triceps brachii of a control subject and a participant with SCI. Reaction time was measured as the time from cue (dashed line) to onset of the EMG burst (arrows). Reaction times in biceps brachii were similar between control and SCI participants, and reaction times in triceps brachii were prolonged in all conditions in the SCI participant compared with control. Notably, in biceps brachii, the shortening of reaction time was similar between control and SCI participants for ΔTSR but smaller in SCI participant for ΔTAR. In triceps brachii, the shortening of reaction time for ΔTSR and ΔTAR was similar between control and SCI participants. B, Box plot charts represent the group data. Abscissa indicates the groups tested (blue bars represent the control group; orange bars represent the SCI group). Ordinate indicates the reticulospinal gain. Top and bottom lines of the box indicate the 75th percentile (top quartile) and 25th percentile (bottom quartile), respectively. Line in the box indicates the 50th percentile (median). The two bars extend from the maximum and minimum value. Dot within the box indicates the arithmetical mean. ***p < 0.001.
Repeated-measures ANOVA showed no effect of group (F(1,31) = 0.1, p = 0.7), condition (F(2,62) = 0.7, p = 0.5), nor in their interaction (F(2,62) = 0.1, p = 0.9) on the mean rectified EMG activity in the biceps muscle measured 100 ms before stimulus presentation. Similarly, repeated-measures ANOVA showed no effect of group (F(1,31) = 0.6, p = 0.5), condition (F(2,62) = 0.04, p = 1.0), nor in their interaction (F(2,62) = 0.4, p = 0.7) on the mean rectified EMG activity in the triceps muscle measured 100 ms before stimulus presentation. Repeated-measures ANOVA showed no effect of group (F(1,31) = 0.01, p = 0.9), but an effect of condition (F(2,62) = 50.1, p < 0.001) and in their interaction (F(2,62) = 4.7, p = 0.01) on the mean rectified EMG activity in the biceps measured over 100 ms after the EMG burst onset. Here in both groups, the mean EMG activity was larger during VSRT (control: 0.30 ± 0.15 mV; SCI: 0.25 ± 0.27 mV) compared with VRT (control: 0.16 ± 0.08 mV, t(14) = 13.0, p < 0.001, dz = 2.9; SCI: 0.16 ± 0.20 mV, t(17) = 7.0, p < 0.001, dz = 0.9). Also, the mean EMG activity was larger during VSRT compared with VART (control: 0.19 ± 0.10 mV, t(14) = 9.8, p < 0.001, dz = 2.9; SCI: 0.21 ± 0.26 mV, t(17) = 4.0, p < 0.001, dz = 1.0). The mean EMG activity was larger during VART compared with VRT in the control (t(14) = 3.2, p = 0.004, dz = 1.4) and SCI groups (t(17) = 3.0, p = 0.005, dz = 1.6). Similarly, repeated-measures ANOVA showed an effect of group (F(1,31) = 5.0, p = 0.03), condition (F(2,62) = 18.4, p < 0.001), and in their interaction (F(2,62) = 7.0, p = 0.002) on the mean rectified EMG activity in the triceps measured over 100 ms after the EMG burst onset. Here in both groups, the mean EMG activity was larger during VSRT (control: 0.17 ± 0.14 mV; SCI: 0.07 ± 0.10 mV) compared with VRT (control: 0.08 ± 0.05 mV, t(14) = 9.2, p < 0.001, dz = 1.5; SCI: 0.04 ± 0.06 mV, t(17) = 11.8, p < 0.001, dz = 2.5). Also, the mean EMG activity was larger during VSRT compared with VART (control: 0.10 ± 0.07 mV, t(14) = 6.8, p < 0.001, dz = 2.9; SCI: 0.05 ± 0.07 mV, t(17) = 6.9, p < 0.001, dz = 3.1). The mean EMG activity was larger during VART compared with VRT in control (t(14) = 2.4, p = 0.02, dz = 2.0) and SCI groups (t(17) = 4.9, p < 0.001, dz = 2.6).
To ensure that any changes in reaction times were related to changes in reticulospinal gain, we compared the ΔTSR/ΔTAR ratios across groups. Two-way ANOVA showed an effect of group (F(1,62) = 6.2, p = 0.02) and muscle (F(1,62) = 7.5, p = 0.008) and in their interaction (F(1,62) = 7.4, p = 0.008) on the ΔTSR/ΔTAR ratios (Fig. 4B). Post hoc tests showed that the reticulospinal gain was not different between biceps (1.6 ± 0.4) and triceps in the control group (1.6 ± 0.3, t(28) = 0.004, p = 1.0, d = 0.0) but was larger in biceps (2.3 ± 0.8) compared with triceps in the SCI group (1.5 ± 0.5, t(34) = 4.0, p < 0.001, d = 1.2). Post hoc tests also showed that the reticulospinal gain in biceps was larger in SCI compared with controls (t(31) = 3.7, p < 0.001, d = 1.1), while no difference was found in triceps between groups (t(31) = 0.2, p = 0.9, d = 0.2). Within the SCI group, the reticulospinal gain remained similar between participants taking and not taking antispastic medication (biceps, t(16) = 0.9, p = 0.4, d = 0.4; triceps, t(16) = 1.7, p = 0.3, d = 0.5).
CMEPs
Figure 5A illustrates the mean CMEP-test in a control subject and a participant with SCI (CMEP-test, black) and the facilitation of the CMEP induced by an SAS (CMEP-startle, gray). The SCI participant exhibited a larger facilitation of the CMEP-startle in biceps but a similar facilitation in triceps compared with the control participant. Two-way ANOVA showed no effect of group (F(1,44) = 0.2, p = 0.6), muscle (F(1,44) = 3.5, p = 0.1), nor in their interaction (F(1,44) = 0.7, p = 0.4) on the amplitude of the CMEP-test (controls: biceps = 1.3 ± 0.6 mV and triceps: 0.9 ± 0.3 mV; SCI: biceps = 1.3 ± 0.7 mV and triceps = 1.1 ± 0.6 mV). However, two-way ANOVA showed an effect of group (F(1,44) = 10.5, p = 0.002), muscle (F(1,44) = 13.1, p < 0.001), and in their interaction (F(1,44) = 11.7, p = 0.001) on the amplitude of the CMEP-startle (Fig. 5B). Post hoc tests showed that the CMEP-startle expressed as a percentage of the CMEP-test was not different between biceps (121.5 ± 8.3% of the CMEP-test) and triceps (120.1 ± 20.0% of the CMEP-test, t(22) = 0.1, p = 0.9, d = 0.09) in the control group, but it was larger in biceps (165.5 ± 38.0% of the CMEP-test) compared with triceps (118.9 ± 13.6% of the CMEP-test, t(22) = 5.0, p < 0.001, d = 1.6) in the SCI group. Post hoc tests also showed that the CMEP-startle was larger in SCI compared with controls in biceps (t(22) = 4.7, p < 0.001, d = 1.6) but not in triceps (t(22) = 0.1, p = 0.9, d = 0.07). We found that the latency of the CMEP test was delayed in biceps and triceps in SCI (biceps: 9.5 ± 1.1 ms; triceps: 9.1 ± 1.2 ms) compared with controls (biceps: 8.2 ± 0.8 ms, t(22) = 3.2, p = 0.004, d = 1.3; triceps: 7.6 ± 1.2 ms, t(22) = 3.1, p = 0.005, d = 1.2), consistent with the view that descending motor drive is impaired in SCI individuals. Within the SCI group, no difference was found on the CMEP-startle facilitation in biceps (t(10) = 0.1, p = 0.9, d = 0.06) and triceps (t(10) = 0.5, p = 0.6, d = 0.3) between participants taking and not taking antispastic medications.
CMEPs. A, Mean CMEP-test (black) and CMEP-startle (gray) elicited in the biceps and triceps brachii of a control subject and a participant with SCI. The SCI participant exhibited a larger facilitation of the CMEP-startle in biceps but to a similar extent in triceps brachii compared with the control subject. B, Graph represents group data. Abscissa indicates the groups tested (blue bars represent the control group; orange bars represent the SCI group). Ordinate indicates the facilitation of CMEP-startle as a percentage of the CMEP-test. Error bars indicate SEM. ***p < 0.001.
Correlations
No correlations were found between any of the physiological outcomes acquired in biceps (MVC, r = 0.1, p = 0.7; MEP-max, r = 0.3, p = 0.3; slope, r = 0.4, p = 0.1; M-max, r = 0.2, p = 0.4; reticulospinal gain, r = −0.5, p = 0.1) and triceps (MVC, r = 0.2, p = 0.4; MEP-max, r = −0.1, p = 0.8; slope, r = 0.01, p = 1.0; M-max, r = 0.2, p = 0.3; reticulospinal gain, r = 0.2, p = 0.5) and the level of injury. Similarly, no correlations were found between physiological outcomes in biceps (MVC, r = −0.2, p = 0.3; MEP-max, r = −0.2, p = 0.4; slope, r = −0.3, p = 0.2; M-max, r = −0.3, p = 0.2; reticulospinal gain, r = 0.2, p = 0.4) and triceps (MVC, r = −0.1, p = 0.8; MEP-max, r = 0.4, p = 0.1; slope, r = 0.4, p = 0.2; M-max, r = −0.3, p = 0.2; reticulospinal gain, r = −0.3, p = 0.3) and the medication intake and the time after injury (biceps: MVC, r = −0.004, p = 1.0; MEP-max, r = 0.1, p = 0.8; slope, r = −0.01, p = 1.0; M-max, r = −0.2, p = 0.5; reticulospinal gain, r = −0.4, p = 0.1 and triceps: MVC, r = 0.02, p = 0.9; MEP-max, r = 0.1, p = 0.7; slope, r = 0.01, p = 1.0; M-max, r = 0.2, p = 0.5; reticulospinal gain, r = −0.1, p = 0.6).
Discussion
In the present study, we show that corticospinal and reticulospinal pathways contribute to a different extent to the control of biceps and triceps in humans with tetraplegia. On the one hand, we found similar MEPs and MVCs but increased StartReact and CMEP facilitation in the biceps in SCI compared with controls. On the other hand, we found smaller MEPs and MVCs and similar StartReact and CMEP facilitation in the triceps in SCI compared with controls. Collectively, these results suggest that the recovery of biceps after cervical SCI results, at least in part, from increased reticulospinal inputs and that the lack of these extra inputs combined with the loss of corticospinal drive might contribute to the pronounced weakness found in triceps. This knowledge might help the design of therapeutic strategies for upper limb function after SCI.
Descending control of elbow flexor and extensor muscles after SCI
Our findings support the view that biceps receive strong corticospinal inputs compared with triceps in people with and without SCI. We found that the amplitude of MEPs (Brouwer and Hopkins-Rosseel, 1997) and maximal contractions (Ditunno et al., 1992, 2000; Murray et al., 2000; Calancie et al., 2004; McKay et al., 2011) were larger in biceps compared with triceps in both groups. We also found that MEPs in triceps were reduced in amplitude and had delayed latencies in SCI compared with controls, as previously shown (Thomas et al., 1997b; Calancie et al., 1999; Perez, 2012). MEPs elicited by TMS over the primary motor cortex provide an index of corticospinal excitability (Lemon et al., 1995; Petersen et al., 2003) and that the corticospinal tract plays a major role in the control of voluntary movement (Lemon, 2008) in monkeys (Fetz and Cheney, 1980) and humans (Day et al., 1987, 1989). An intriguing question is why the amplitude of corticospinal responses in SCI participants was similar to controls in biceps but not in triceps. Motor cortical representations of biceps and triceps are in close locations in monkeys (Kwan et al., 1978) and humans (Penfield and Boldrey, 1937). Indeed, there is a large overlap on the corticospinal sites from which elbow flexor and extensor movements can be evoked (Graziano et al., 2004; Graziano, 2006) and on corticospinal neurons controlling elbow flexor and extensor muscles (Witham et al., 2016). The location of motoneurons innervating biceps and triceps also largely overlaps (Schirmer et al., 2011), and most injuries to the cervical spinal cord affect more than one spinal segment (Schaefer et al., 1989; Benavides et al., 2020). Thus, it is unlikely that, following SCI, corticospinal projections and/or motoneurons innervating triceps were more damaged compared with those innervating biceps. It is also unlikely that our results relate to differences in stimulation parameters. As before (Perez et al., 2014), we found similar MEP thresholds for biceps and triceps in controls. This also agrees with our results in SCI participants showing that, although MEP thresholds were higher in triceps than biceps, smaller MEPs were elicited in triceps.
We favor two possible explanations for our results. The first possibility is that the recovery of responses in biceps after SCI involves reticulospinal compensatory mechanisms. This is supported by our results showing that the StartReact and CMEP facilitation were larger in biceps compared with controls. The StartReact likely reflects an involuntary release of a planned movement when an unexpected startle cue is presented (Brown et al., 1991; Valls-Solé et al., 1995, 1999). Although the shortening in reaction time likely involves subcortical structures, such as the reticular formation (Davis and Gendelman, 1977; Valls-Solé et al., 1999; Carlsen et al., 2003, 2004, 2009), a contribution from cortical pathways cannot be excluded (Marinovic and Tresilian, 2016). Indeed, Dean and Baker (2017) showed that StartReact responses in finger muscles showed a similar degree of muscle fractionation during rapid responses without a startling cue, suggesting that reticulospinal and corticospinal pathways contributed in parallel to the StartReact. If this is the case, why did we not observe the same results in triceps? In triceps, we found that the StartReact response was similar between groups. Stimulation of the pontomedullary reticular formation in monkeys results in a series of bilateral reciprocal patterns with a tendency to facilitate ipsilateral flexors and inhibit ipsilateral extensors and to facilitate contralateral extensors and inhibit contralateral flexors (Davidson and Buford, 2006; Davidson et al., 2007; Herbert et al., 2010; Hirschauer and Buford, 2015). Because reticulospinal axons are bilaterally distributed with a majority of ipsilateral axons in monkeys and humans (Nathan et al., 1996; Sakai et al., 2009), it might not be surprising that ipsilateral effects are predominant. Indeed, corticospinal lesions in primates lead to enhanced monosynaptic and disynaptic reticulospinal facilitation to forearm flexors compared with extensor muscles (Zaaimi et al., 2012). We therefore hypothesized a preferential facilitatory contribution of reticulospinal inputs to biceps compared with triceps motoneurons after SCI. One possibility is that reticulospinal inputs compensate for the loss of corticospinal axons (Pettersson et al., 2000; Zaaimi et al., 2012). This is supported by our results showing that, when a startle cue was given 80 ms before a CMEP, an interval reflecting influences from subcortical pathways (Furubayashi et al., 2000; Tazoe and Perez, 2017), the CMEP amplitude increased in biceps but not in triceps. This is also supported by the fact that these changes were observed in SCI but not in controls, underlying the view that these changes might be related to reorganization after the injury.
The second possibility is that our results relate to differences in the nature of the distribution of corticospinal projections to biceps and triceps. The terminal distribution of crossed and ipsilateral corticospinal projections to interneurons and motoneurons is larger in biceps compared with triceps (Morecraft et al., 2013). In addition, in monkeys (Phillips and Porter, 1964; McKiernan et al., 1998) and humans (Palmer and Ashby, 1992; de Noordhout et al., 1999), biceps motoneurons receive large corticospinal monosynaptic facilitation, whereas triceps motoneurons receive lesser monosynaptic facilitation but strong disynaptic inhibition. Thus, it is possible that responses in biceps in SCI participants reflect the large direct facilitatory corticospinal control. A larger facilitatory effect on biceps can contribute to the increases observed in reticulospinal gain. This agrees with evidence showing that increases in corticospinal excitability can result in shorter reaction time (Benito et al., 2012; Bunday and Perez, 2012; Long et al., 2017). This is also supported by our data showing larger MEP responses in biceps compared with triceps in SCI participants. Overall, our results support the possibility that both processes, either reticulospinal inputs compensate for the loss of corticospinal axon and/or that large direct facilitatory corticospinal control contributed to the results observed in biceps. It is of interest to note that the properties of motoneurons for both muscles differ. Persistent inward currents, which amplify synaptic input and promote self-sustained firing, are higher in the lateral head of the triceps than in the biceps brachii during isometric contractions (Wilson et al., 2015). This is consistent with evidence showing a higher incidence for persistent inward currents in extensor compared with flexor muscles (Hounsgaard et al., 1988; Cotel et al., 2009). Although the loss of descending drive to motoneurons leads to enhanced persistent inward currents (Gorassini et al., 2004), the similar effects of a startling cue in our measurements in the elbow extensor muscle between SCI and control suggest that it is less likely that differences in amplification of synaptic input at the motoneuronal level contributed to our findings.
Functional considerations
Multiple descending motor pathways likely contribute to the recovery of function after SCI. A critical question is how these contributions might take place. For example, we recently showed lesser corticospinal and larger reticulospinal influences to spastic muscles in humans with SCI (Sangari and Perez, 2019), suggesting that the contributions from corticospinal and reticulospinal tracts are imbalanced when spasticity is present (Owen et al., 2017; McPherson et al., 2018; Choudhury et al., 2019). Here, we tested SCI participants who did not have spasticity in biceps or triceps to avoid the effect of this variable on our results. Spasticity is also present to a much lesser extent in upper compared with lower limb muscles (Sköld et al., 1999). Interestingly, we found increased reticulospinal inputs to biceps and decreased corticospinal inputs to triceps. Is this good or bad reorganization? Discussions are still ongoing about the positive or negative nature of reticulospinal reorganization after CNS damage (Brunnstrom, 1970; Karbasforoushan et al., 2019). Our present and previous results (Baker and Perez, 2017) favor the interpretation that reticulospinal reorganization can contribute to functional tasks. This is consistent with evidence in animal models of SCI showing that injured reticulospinal axons form new connections that can contribute to functional recovery (Wakabayashi et al., 2001; Ballermann and Fouad, 2006; Filli et al., 2014; May et al., 2017). Others have also suggested that the recovery of strength in patients with CNS damage involves strengthening of the reticulospinal pathway (Xu et al., 2017; Choudhury et al., 2019). Therapeutically, increased engagement of both corticospinal and reticulospinal inputs might then help the recovery of triceps after SCI; and if spasticity is present, favoring corticospinal inputs might be a preferred approach (Sangari and Perez, 2019). Whether reorganization in reticulospinal pathways is “good” or “bad” might depend on a number of factors, and this might also differ across patient populations, across control of upper and lower limb muscles, and even across muscles of the upper limb as shown in the present study.
Footnotes
M.A.P. was supported by the National Institute of Neurological Disorders and Stroke and the Department of Veterans Affairs. S.S. was supported by the Craig H. Neilsen Foundation.
The authors declare no competing financial interests.
- Correspondence should be addressed to Monica A. Perez at mperez04{at}sralab.org
