Closed-loop neuromodulation of spinal sensorimotor circuits controls refined locomotion after complete spinal cord injury

See allHide authors and affiliations

Science Translational Medicine  24 Sep 2014:
Vol. 6, Issue 255, pp. 255ra133
DOI: 10.1126/scitranslmed.3008325


Neuromodulation of spinal sensorimotor circuits improves motor control in animal models and humans with spinal cord injury. With common neuromodulation devices, electrical stimulation parameters are tuned manually and remain constant during movement. We developed a mechanistic framework to optimize neuromodulation in real time to achieve high-fidelity control of leg kinematics during locomotion in rats. We first uncovered relationships between neuromodulation parameters and recruitment of distinct sensorimotor circuits, resulting in predictive adjustments of leg kinematics. Second, we established a technological platform with embedded control policies that integrated robust movement feedback and feed-forward control loops in real time. These developments allowed us to conceive a neuroprosthetic system that controlled a broad range of foot trajectories during continuous locomotion in paralyzed rats. Animals with complete spinal cord injury performed more than 1000 successive steps without failure, and were able to climb staircases of various heights and lengths with precision and fluidity. Beyond therapeutic potential, these findings provide a conceptual and technical framework to personalize neuromodulation treatments for other neurological disorders.


Regulation of dysfunctional neuronal circuits with electrical neuromodulation therapies has broadened the spectrum of treatment options for neurological disorders (1, 2). For example, deep brain stimulation of basal ganglia has become a common medical practice to alleviate motor symptoms of Parkinson’s disease (36). Experimentally, epidural electrical stimulation (EES) and intraspinal stimulation of the lumbosacral and thoracic spinal segments have enabled improved motor control capabilities in rodents and humans with spinal cord injury (SCI) (79) and Parkinson’s disease (10, 11). In these experiments, empirical knowledge and visual observations guided the tuning of electrical neuromodulation parameters. After manual adjustment of pulse width, amplitude, and frequency, stimulation of the brain and spinal cord remained constant. Such an approach is likely suboptimal to mediate maximum therapeutic effects in patients (1, 2). Indeed, preliminary evaluations have suggested that nonregular temporal patterns (5, 6) and closed-loop adjustment (4, 12) of deep brain stimulation augment the therapeutic impact of neuromodulation.

The ability to modulate the parameters and timing of EES of lumbosacral circuits to improve motor control after SCI remains unclear. Various studies have shown that tuning frequency, amplitude, and pulse width of EES can modulate specific aspects of standing, stepping, and isolated movements in both animal models and humans (8, 1315). However, to our knowledge, a comprehensive mapping of the functional relationships between stimulation parameters and limb kinematics that could be leveraged to achieve closed-loop control of EES has not been conducted.

Several limitations have restricted the development of closed-loop control systems for neuromodulation (16). First, efficient control policies require a mechanistic understanding of the relationships between stimulation features, the recruitment of specific sensorimotor circuits, and the resulting functional effects. Second, closed-loop control of neuromodulation relies on the availability of a stable and robust feedback signal to tune stimulation parameters. Third, implementation and optimization of closed-loop control algorithms necessitate advanced technological platforms combining feedback and feed-forward loops that match the natural flow of information in the modulated neural systems. Here, we sought to address these issues simultaneously to modulate lumbosacral circuits through closed-loop EES.

Previous studies using EES to improve standing and walking in severely paralyzed subjects have operated under the assumption that the stimulation primarily increases the general level of excitability within spinal circuits (8) and activates central pattern generators (13). However, EES is also able to recruit specific afferent pathways connected to distinct sensorimotor circuits, which, in turn, generate well-defined responses in leg muscles after each pulse of stimulation (1719). We first capitalized on this understanding of EES mechanisms to establish a mapping between EES parameters, neuromodulation of sensorimotor circuits, and adjustments of leg movements. Second, we developed robust algorithms to extract a useful feedback signal from spatiotemporal characteristics of leg kinematic during locomotion. Third, we established a versatile technological platform with embedded control policies that interfaced feedback signals and feed-forward models operating in real time. These integrated developments allowed us to conceive a closed-loop neuromodulation system that controlled a broad range of leg movements during natural locomotion in paralyzed rats.


Modulation of EES parameters tunes specific aspects of leg movements

We first characterized the impact of EES amplitude, pulse width, and frequency on the modulation of leg kinematic, ground reaction force, and muscle activity in adult rats that received a complete transection of the mid-thoracic spinal cord (fig. S1). After 5 weeks of recovery, all of the tested rats (n = 5) dragged both hindlimbs along the treadmill belt when supported bipedally in a robotic body weight support system (Fig. 1A). The combination of serotonin agonists (5HT1A/2/7) and EES delivered at S1 and L2 spinal segments (Fig. 1B) promoted continuous and alternating stepping movements of the legs. Frequency, amplitude, and pulse width of EES were adjusted to cover the entire range of values that promoted functional movements.

Fig. 1. Modulation of EES frequency tune foot trajectory during locomotion.

(A) Rats with complete SCI were positioned bipedally over a treadmill belt using an adjustable body weight support system. In the inset, markers were positioned over leg joints, including the limb endpoint. (B) EES was applied along the midline of L2 and S1 spinal segments. Serotonin agonists (5HT) were administered systemically before stepping. (C) Representative, color-coded stick diagram decomposition of hindlimb movements during stance (black), drag (cyan), and swing (gray), together with 10 successive trajectories of the limb endpoint. The red dot indicates the maximum step height. The black arrows represent the direction and intensity of the limb endpoint trajectory at swing onset. (D) Principal component (PC) analysis applied on all 147 gait variables (table S1) measured over 10 successive gait cycles during EES frequencies ranging from 20 to 90 Hz, with a 5-Hz increment for a representative rat. Least-square elliptic fits are drawn to emphasize frequency-dependent modulation of gait patterns. The bar graph reports the score for each EES frequency along PC1. Data are means + SEM. a.u., arbitrary units.

Increments of EES frequency led to a graded and consistent modulation of gait patterns in all the tested rats within a range of 20 to 90 Hz (Fig. 1, C and D). To characterize this modulation and identify the most relevant features reflecting tuning of leg movements with EES frequency, we computed 147 kinematic, force, and electromyographic (EMG) variables that provided a comprehensive description of the locomotor state (listed in table S1). We subjected these combined parameters to a PC analysis. Stepping patterns associated with each EES frequency were differentiated along PC1, which explained nearly 50% of the total variance in the data sets (Fig. 1D). Analysis of factor loadings on PC1 revealed that step height was the more robust variable (factor loading, 0.97) to capture the modulation of leg movements during changes in EES frequency (fig. S2). There were linear correlations between EES frequency and step height that were reproducible within a session (Fig. 2A), over weeks (Fig. 2A), and among rats (Fig. 2B), despite idiosyncratic gait patterns. Tuning EES frequency over the useful range of values (shaded areas in Fig. 2A) mediated a threefold increase in foot elevation, which allowed rats to cover step heights ranging from 2.9- to 6.8-cm mean across rats.

Fig. 2. Reproductive modulation of step height with EES frequency within and across days, and between rats.

(A) Relationship between EES frequency and step height for a representative rat tested 30 min apart in the same session 5 weeks post-lesion, and at 7 and 9 weeks post-lesion. The linear regression (back line) was computed over the range of functional EES frequencies (red dots, shaded areas). Each dot represents a gait cycle. (B) Plots reporting the linear regressions from five rats recorded in the same conditions as in (A).

Increase in EES amplitude (Fig. 3, A to C) mediated adjustments of lower limb kinematics similar to those observed with changes in EES frequency, including a significant, albeit limited increase in step height (2.8 to 3.9 cm) (Fig. 3C). Changes in limb kinematics were less systematic when adjusting pulse width (Fig. 3, D to F). Moreover, increase in EES amplitude and pulse width induced coactivation of leg muscles (Fig. 3, A and D), which led to a rapid saturation of stepping movements, and thus limited the relevance of these parameters for neuromodulation strategies.

Fig. 3. EES amplitude and pulse width mediate limited tuning of leg movements during locomotion.

(A) EES amplitude was gradually increased by fixed increments of 60 μA from the onset of EES-induced stepping, until saturation. EES frequency was kept constant at 40 Hz. Stick diagram decomposition of hindlimb movements, successive limb endpoint trajectories, and alternating flexor (Flex) and extensor (Ext) EMG activity are displayed for each experimental condition using the same convention as in Fig. 1. (B) PC analysis was applied on all the gait variables measured in a representative rat during EES amplitude changes shown in (A). The bar graph reports the score for each EES amplitude along PC1 for all tested rats (n = 5). (C) Histogram plots reporting changes in step height in response to modulation of EES amplitude. (D) EES pulse width was gradually increased by fixed increments of 60 μs in the same rat and under the same conditions as (A). (E) PC analysis was applied on all the gait variables measured during EES pulse width changes shown in (D). (F) Histogram plots reporting changes in step height in response to modulation of EES pulse width. *P < 0.05, **P < 0.01, repeated-measures analysis of variance (ANOVA) with Bonferroni post hoc analysis. Data are means ± SEM (n = 5).

EES frequency modulates specific sensorimotor circuits

We next sought to identify the physiological principles underlying the modulation of step height with EES frequency. EES primarily recruits proprioceptive afferent fibers (17, 19, 20). Here, we found that the recruitment of proprioceptive fibers engages monosynaptic and polysynaptic reflex circuits, which elicited well-defined motor responses in leg muscles after each pulse of EES (Fig. 4A). We dissociated these responses functionally based on their respective latencies (medium- and late-latency responses). In turn, phase-dependent modulation of these responses in flexor versus extensor muscles resulted in alternating bursts of EMG activity, enabling continuous locomotion in otherwise paralyzed rats (Fig. 4B and fig. S3, A and B).

Fig. 4. EES frequency modulates specific reflex circuits.

(A) Schematic of neural structures recruited by EES. In, interneuron; Mn, motoneuron; DRG, dorsal root ganglion; Ib, group Ib fibers conveying force feedback; Ia and II, group Ia and group II fibers conveying muscle length–related feedback. A single pulse of EES evokes a medium-latency response (red) and a late-latency response (blue) in leg muscles through activation of proprioceptive segmental circuits. (B) Sequence of EMG activity recorded in the contralateral extensor (medial gastrocnemius) and ipsilateral flexor (tibialis anterior) muscles during continuous stepping at three different EES frequencies. The lower traces show, for each frequency, motor responses evoked in extensor (ext) and flexor (flex) muscles at distinct moments of the gait cycle (i, ii, and iii), but over the same duration. (C) Power spectrum of EMG bursts for flexor and extensor muscles at different stimulation frequencies. The red and gray vertical lines indicate EES frequency and its harmonic multiples, respectively. (D) Color-coded modulation of the amplitude of medium- and late-latency responses in flexor and extensor muscles over the course of a gait cycle for EES frequencies ranging from 20 to 90 Hz, with a 10-Hz increment. The white lines indicate the end of stance for each frequency. (E) Three-dimensional (3D) linear correlation between color-coded EES frequency and normalized EMG amplitude of contralateral extensor and ipsilateral flexor muscles. Each point represents an individual gait cycle. The robustness of the correlation was calculated through PC analysis as the amount of explained variance by PC1. All panels report data recorded in the same representative rat (representative of n = 5).

EES-induced motor responses remained locked to each pulse of stimulation across the entire range of tested EES frequencies (Fig. 4, B and C). Consequently, increment of EES frequency induced a graded increase in the number of medium- and late-latency responses in leg muscles (Fig. 4, C and D, and fig. S3C). At the lower frequencies, both responses independently contributed to the elaboration of muscle activity in both flexor and extensor muscles (fig. S3D). At higher frequencies, the monosynaptic input likely reached motoneurons at the same time as the polysynaptic input from the preceding stimulation. Temporal convergence of monosynaptic and polysynaptic inputs onto motoneurons eventually led to a single, medium-latency response, which accounted for the increase in the amplitude of this response with larger EES frequencies (fig. S3, C and D).

The graded increase in the number and amplitude of motor responses led to a linear enhancement of the overall bilateral EMG activity in extensor and flexor ankle muscles (Fig. 4E and fig. S3D). In turn, the concurrent modulation of ankle extensors in the stance leg and flexor muscles in the swing leg accounted for the progressive tuning of step height when adjusting EES frequency (Fig. 4E). These results indicate that tuning of EES frequency mediates a progressive and reproducible modulation of monosynaptic and polysynaptic reflex circuits, which promotes predictive adjustment of leg movements during locomotion in rats.

Closed-loop EES achieves high-fidelity control of leg movements

Linear tuning of step height with increase in EES frequency opened the intriguing possibility to automate leg movement control in real time through closed-loop modulation of EES frequency. To test this hypothesis, we designed a technological platform that combined (i) online monitoring of kinematics, ground reaction force, and muscle activity in real time; (ii) signal-processing algorithms to reconstruct locomotor states, including limb endpoint trajectory, gait phases, and whole-body position; and (iii) a computational infrastructure to embed control policies in the feedback loop to adjust EES frequency based on current and desired stepping behaviors (Fig. 5A, fig. S4, and movie S1). Optimization of communication protocols and information-processing time allowed the entire loop to operate in less than 20 ms (see Materials and Methods), which was appropriate for real-time control of neuromodulation during gait despite the fast rhythm of stepping in rats (cycle duration, 810 ± 156 ms, SD).

Fig. 5. Real-time control of step height during locomotion.

(A) The architecture of the closed-loop system combined real-time monitoring of kinematic information and automated control policies. After each foot strike, the signal-processing algorithm measured the step height (red dot, hi), which was fed into the control loop. The control structure integrated parallel loops combining feedback correction (PI controller) and a predictive forward model. This predictive component was elaborated from the linear relationship between EES frequency and step height, and was updated at each iteration through an adaptive algorithm (see Materials and Methods). The controller combined the error (ei) correction based on the current reference (ri), and the prediction to generate the adjustment of EES frequency (fi+1) to adjust the current frequency (fi) and thus maintain the maximum step height within the user-defined reference band (ri+1, shaded areas). A saturation function was inserted to constrain EES within the functional range of frequencies. (B) Example of successive step heights and concomitant EES frequencies during a random sequence of reference band changes for one rat. (C) Mean values ± SEM (n = 15 trials) of successive step heights for various combinations of steady and constant-changing reference bands (shaded area).

We used this real-time platform to implement a closed-loop control policy capable of adjusting step height in real time during locomotion. The continuous flow of kinematic information was discretized per gait cycle and merged into a single variable (step height) that was used to automatically adjust EES frequency during the next gait cycle (Fig. 5A). The controller adjusted stimulation to maintain the maximum elevation of the foot within a reference band, whose upper and lower bounds (±5 mm) were selected to account for the intrinsic variability of step height during continuous locomotion (4.6 ± 0.9 mm, SD). This single-input single-output (SISO) control structure integrated parallel loops combining adaptive feed-forward prediction and feedback correction (Fig. 5A and fig. S4).

The feed-forward model was derived from the linear correlation between step height and EES frequency recorded for each rat at the beginning of each session (Fig. 2). The slope and intercept values of the model were iteratively updated after each gait cycle to account for time-varying changes in stepping behavior. A proportional-integral (PI) feedback loop complemented the feed-forward model predictions to minimize the discrepancy between the desired and actual step height (movie S1). This architecture achieved a high degree of control over a broad range of step heights during both steady and constantly changing locomotor states (Fig. 5, B and C).

To quantify the degree of step height control, we designed a comprehensive suite of testing paradigms, in which the reference band was displaced through a range of fixed increments or decrements after each step, or after a series of steps (Fig. 5C). These tasks tested the performances of the controller during steady-state (fig. S5A) and dynamic (fig. S5, B and C) behaviors over the entire range of step height modulations. Despite these multifaceted challenges, the controller successfully tuned EES frequency to instantly adjust step height to varying reference bands (movie S2). Precision of step height control significantly increased in variability (fig. S5D) when imposing adjustments that reached the physiological limits of the studied biomechanical system.

Closed-loop EES prevents rapid fatigue during continuous locomotion

We next sought to evaluate the ability of closed-loop EES to improve key aspects of gait performance. Human patients with SCI exhibit a rapid exhaustion of locomotor output during continuous stepping on a treadmill (21), which is also observed in injured rats (22). During noncontrolled EES (40 Hz), rats with complete SCI showed a progressive decrease in the amplitude of muscle activity during stepping (Fig. 6, A and B), which led to a graded decrease in step height until the animals collapsed (Fig. 6C).

Fig. 6. Closed-loop EES compensates for exhaustion.

(A) Stick diagram decomposition of hindlimb movements and extensor muscle activity for both noncontrolled and controlled conditions at three different time points: beginning, middle, and end of the noncontrolled session. Conventions are the same as in Fig. 1C. The gray- and cyan-shaded areas indicate the duration of the stance and drag phases, respectively. The reference (gray band) was set on the basis of the average step height measured across rats. (B) Histogram plots reporting the durations of stepping until collapse and within the reference band, and the amplitude of extensor muscle activity (medial gastrocnemius, MG) for the three time points shown in (A) during noncontrolled and controlled conditions for all the rats (n = 3 rats). (C) Representative examples of successive step heights and EES frequency modulation during a continuous sequence of stepping under noncontrolled and controlled conditions. The vertical dotted line corresponds to the saturation time when EES frequency reached the maximum allowed value for the controller, that is, 95 Hz (Fig. 5A). Each data point corresponds to a gait cycle. Data are means ± SEM. *P < 0.05, **P < 0.01, Kruskal-Wallis test.

Instead, closed-loop control of step height doubled the duration of stepping before collapse compared to noncontrolled EES (Fig. 6C). The controller continuously increased EES frequency over time, which limited the decrease in the amplitude of muscle activity (Fig. 6B), and enabled rats to maintain step height within the constant reference band for extended durations despite exhaustion (Fig. 6C and movie S3). Under controlled conditions, 100% of the tested rats (n = 3) performed at least 1000 successive steps without failure over the duration of the testing session.

Closed-loop EES enables locomotion across combinations of staircases

We lastly aimed to illustrate the practical impact of closed-loop EES to restore walking during more natural locomotor conditions. For this, we positioned the rats overground in a robotic interface that moved the rats forward along a runway while providing constant-force support in the vertical (25% of body weight) and mediolateral (rigid, 100%) directions (Fig. 7A). We interposed staircases of various heights (1.3 to 3.5 cm), lengths (12 to 25 cm), and numbers (one to four stairs) along the runway to test the ability of closed-loop EES to modulate step height to accommodate leg movements during stair climbing. The reference band was set at a constant height during locomotion along the horizontal surface until crossing a predefined distance from the staircase (1.5 cm) where the reference band was raised to an appropriate height (3 cm above staircase height) to pass the staircase. During noncontrolled EES, rats tumbled against and failed to pass the lower staircase (Fig. 7A). In contrast, closed-loop EES enabled the rats (n = 3) to climb all of the tested staircase combinations successfully in 99 of 100 successive attempts (P < 0.001 versus noncontrolled EES, Kruskal-Wallis test) (Fig. 7, A and B, and movie S4).

Fig. 7. Real-time control of locomotion across staircases.

(A) Rats with complete SCI were attached to a robotic postural interface that provided adjustable, constant-force support in the vertical and mediolateral directions while moving the animal forward at a constant speed (9 cm/s). The reference band was displaced when crossing a virtual line located 1.5 cm away from each staircase. Stick diagram decomposition of hindlimb movements is shown together with EES frequency modulation during a representative execution along a staircase. Changes in vertical ground reaction forces (GRFs) when stepping onto the force plate are displayed at the bottom. Shaded areas indicate the occurrence of stance. The pie plots report the percentage of trials during which the rats (n = 3 rats) tumbled, touched, or passed the staircase across 100 successive repetitions under noncontrolled and controlled conditions. (B) Polynomial fit of step height values during locomotion along staircases of various heights and numbers. The concurrent vertical ground reaction forces measured during the step preceding and coinciding with the staircase are displayed at the bottom. Automated tuning of EES frequency significantly augmented vertical ground reaction forces proportional to the displacement of the reference band (**P < 0.01, ***P < 0.005, Kruskal-Wallis test for each staircase). Data are averages (n = 3 rats, 20 trials per rat) and SEM (shaded areas). The pie plots report the percentage of success for the first and second staircases over 20 successive trials.

Automated tuning of EES frequency during the stance phase preceding the staircase significantly augmented vertical ground reaction forces proportional to the displacement of the reference band (Fig. 7B), which generated the adequate foot trajectory to overcome the staircase. The same precision of step height control was obtained during sequences of staircases, which required a continuous adjustment of EES frequency to pass the successive stairs while stabilizing upward and downward movements of the body after each step (fig. S6). These combined results show the ability of the developed neuromodulation algorithms to achieve real-time adjustment of leg kinematic and ground reaction forces with adequate precision to restore complex walking behaviors in paralyzed rats.


Open-loop regulation of spinal sensorimotor circuits with neuromodulation therapy improves motor control after SCI (7, 8, 14, 23). Here, we demonstrate that closed-loop neuromodulation yielded superior control of leg movement compared to continuous stimulation, enabling rats with complete SCI to perform complex locomotor tasks with precision and fluidity. Although the current demonstration is in rodents, research platforms for human spinal cord stimulation (16) and a newly developed robotic interface (24) provide the technological infrastructure to achieve immediate translation of our control algorithms into neuromodulation therapies to facilitate robot-assisted training in humans with SCI.

EES of lumbosacral circuits facilitates locomotion through the recruitment of proprioceptive afferent fibers, which engages two complementary mechanisms (17). First, the propagation of the tonic neural drive elicited along myelinated afferent fibers activates central pattern-generating networks (25, 26). Second, each stimulation pulse elicits a medium- and late-latency response in leg muscles, which corresponds to the activation of monosynaptic and polysynaptic reflex circuits, respectively (13, 17, 19). Central pattern-generating networks modulate the gain of these reflex circuits according to the state of the locomotor system (19, 27). In turn, the successive motor responses compose the bursts of leg muscle activity during locomotion (17, 19). Here, we found that spinally evoked motor responses remain locked to the stimulation rate over a large range of frequencies, enabling fine-tuning of muscle activity through linear adjustment of stimulation. Likewise, frequency-dependent modulation of intraspinal microstimulation delivered in cervical segments near-linearly adjusts the amount of grasp force in anesthetized nonhuman primates (28). This principle underlies natural muscle force production, which primarily relies on the modulation of motoneuron firing rate (29).

We leveraged this mechanistic framework to achieve real-time control of lower limb kinematics and vertical ground reaction forces during continuous locomotion. This control policy allowed paralyzed rats to produce a broad range of foot trajectories during continuous locomotion, and to climb combinations of staircases. Frequency-dependent tuning functions were consistent across days and across rats, which supported the design of algorithms with robust performances. This treatment paradigm engages central pattern-generating networks (30). In addition, each pulse elicits monosynaptic and polysynaptic motor responses, which are modulated according to the phase of the movement. Because the human spinal cord exhibits similar responses to EES (23, 31), our approach—frequency-dependent tuning of muscle activity—could translate into a neuromodulation therapy to improve motor control and rehabilitation in patients with SCI.

Previous closed-loop neuromodulation therapies for restoring locomotion using electrical stimulation of muscles (32, 33) or spinal structures (9) have sought to elaborate the motor command extrinsically through finite-state control of discrete stimulation sequences. This type of control policies enabled stepping-like movements, but was not able to modulate muscle activity to produce movement over a broad range of natural gait trajectories. Moreover, rapid muscle fatigue prevented locomotor execution over extended durations (32, 33), which limited the use of this therapeutic strategy for rehabilitation in humans with SCI (34). Instead, our neuromodulation therapy exploits the intrinsic capacities of sensorimotor circuits embedded in the spinal cord to generate coordinated movements. Thus, we superimposed a simple control policy onto the ongoing motor program to achieve continuous, high-fidelity control of leg movements over extended periods of time.

Our control policy combined feedback signals and feed-forward models that relied on a global parameter capturing the targeted locomotor behavior. For this purpose, we extracted spatiotemporal features from lower limb endpoint trajectories. There is evidence that central nervous system circuits may elaborate an explicit representation of limb endpoint kinematics (35). Moreover, neurobehavioral experiments have suggested that the limb endpoint is the primary variable used to coordinate locomotion in animal models (36) and humans (37). We thus monitored and controlled step height, which was the most robust parameter to account for frequency-dependent modulation of gait. In addition, lower limb endpoint is a neurally relevant variable to encode and organize leg movements. This type of control policy establishes a framework to design a corticospinal neuroprosthesis through which neuronal modulations from the motor cortex will directly modulate electrical spinal cord stimulation to adjust lower limb kinematics and improve locomotor recovery (38). Current wearable and implantable motion sensors can generate the necessary feedback to incorporate automated closed-loop control policies into neuromodulation therapies for use outside laboratory environments (39). For this purpose, industrial development of clinically viable neuromodulation platforms has established the appropriate neurotechnology to personalize stimulation algorithms based on feedback signals (16).

Current knowledge of human spinal cord properties in response to electrical neuromodulation (23, 31) suggests that the control policies developed in our study could translate into clinical applications for facilitating neurorehabilitation (8, 23). Existing technological platforms for closed-loop neuromodulation of the human nervous system (16) can readily be interfaced with real-time kinematic feedback to tune EES in clinical settings. Moreover, our results in rats (40) compelled the development of a robotic body weight support system that can provide patients with adjustable, multidirectional support during overground locomotion within a large workspace (24). Thus, our developed algorithms can be implemented to facilitate rehabilitation during robot-assisted locomotion on a treadmill and overground in patients with SCI. Closed-loop neuromodulation of spinal circuits has the potential to improve motor control and prevent rapid fatigue during robot-assisted training in individuals with incomplete SCI, such as ASIA-C (American Spinal Injury Association grade C) patients, which are essential factors to promote activity-dependent plasticity and recovery with rehabilitation (8, 41).

In rats with complete SCI, continuous weight-bearing stepping only emerged in the presence of agonists to specific serotonergic receptor subtypes. Monoaminergic replacement therapies are necessary to mediate locomotor permissive states of spinal circuits deprived of any source of modulation (42). Consequently, closed-loop neuromodulation therapies in humans with severe injuries will be contingent on the development of safe and efficacious monoaminergic replacement strategies, which are under exploration (43), but do not yet exist for general use in people. However, open-loop electrical neuromodulation of spinal circuits instantly restored fine control over joint-specific movements in patients with chronic motor paralysis (23). These results suggest that the therapeutic effects of neuromodulation therapies that we demonstrated in rodent models (7) also applied in humans with severe SCI.

Neuromodulation therapies of the brain and spinal cord target distinct neural structures, but share common principles, including the growing need for closed-loop regulation (1, 2, 44). Our results provide a conceptual and technical framework for the broad development of closed-loop neuromodulation therapies to improve function after neurological disorders.


Study design

The study was segregated into two experimental groups. Characterization of relationships between EES frequency and gait pattern modulation was conducted on a group of five rats that participated in all of the testing sessions over a period of 3 weeks, starting 5 weeks post-lesion. After comprehensive analysis of EES characterization and neuromodulation algorithm development, we generated an additional group of three rats that were tested daily over a period of 4 weeks, starting 5 weeks post-lesion. All experiments involving real-time control of EES were carried out with those three rats. This sample size was selected on the basis of the robustness of EES modulation effects and to obtain enough statistical power to demonstrate significant differences between closed- and open-loop neuromodulation. All the tested rats are presented in the study. No statistical outliners were excluded. One animal from the second group had to be excluded from PC analysis owing to a mechanical failure in one of the four EMG wires; however, this rat could be included for analysis in all other figures where unilateral EMG signals were sufficient for data analysis. All measurements were obtained using objective readouts with high-precision equipment. Blinding during data acquisition and analysis was not possible because of the obvious effects of EES modulation.

Animals and animal care

All procedures and surgeries were approved by the Veterinarian Office Vaud, Switzerland. The experiments were conducted on adult female Lewis rats (~200 g of body weight, Centre d’Elevage R. Janvier). Rats were housed individually on a 12-hour light-dark cycle, with access to food and water ad libitum.

Surgical procedure and post-surgical care

All surgical procedures have been described in detail previously (7, 40) and are elaborated in the Supplementary Methods

Locomotor training

Rats were attached to an automated, servo-controlled body weight support system (Robomedica), and positioned over a motorized treadmill belt for bipedal locomotion. The animals were trained every other day, 20 min per session for 4 weeks, starting 8 days post-SCI. Locomotor training was enabled by the combination of EES at S1 and L2, and systemic administration of agonists to 5HT1A/7 (8-OH-DPAT, 0.05 to 0.1 mg/kg body weight) and 5HT2A/C (quipazine, 0.2 to 0.3 mg/kg body weight) receptor subtypes (42). Treadmill belt speed was set at 9 cm/s.

Kinematic, ground reaction force, reflex analysis, and EMG recording procedures as well as data and PC analyses are described in the Supplementary Materials.

Real-time platform for closed-loop adjustment of EES

The platform was implemented within a multithreaded C++ code (Visual Studio 2010, Microsoft) running on a quad core Microsoft Windows 7 computer. Stimulation patterns were applied via an RZ5 processing unit (Tucker-Davis Technologies) connected to an MS16 Stimulus Isolator (Tucker-Davis Technologies). The integrated Vicon recording system generated raw 3D positions of the markers, which were imported into the C++ environment in (soft) real time via Ethernet using the DataStream SDK software. We developed a custom algorithm for online interpolation of missing markers through triangulation and relabeling of each marker to the appropriate joint landmark. All signals were filtered online using least mean squares adaptive filters. Specific gait events such as foot strike, toe-off, or predetermined timings during each gait cycle were automatically detected using online kinematic classification of limb endpoint trajectories. These timings triggered controller calculations and model updates (fig. S3). To evaluate the delay of the entire platform, we sent concomitant pulses at opposite ends of the loop, and calculated the discrepancy between both signals after each iteration, for more than 100 successive cycles.

Control structure for closed-loop adjustment of EES

We modeled stepping patterns as a discrete-time, SISO system (stimulation frequency/step height). The continuous flow of kinematic information was discretized by gait cycles using foot strike events, which separated consecutive gait cycles and defined the set of sampling times K = {k1, k2kn}. These events triggered the update of the stimulation frequency, which was applied using zero-order hold control until the next event. The maximum step height during each gait cycle (defined as hi = max(ht), t ∈ [ki−1, ki]) served as feedback signal for closed-loop corrections. Stimulation values were set identically for both electrodes and effectively acted as a unique value. The controller used an adaptive input-output forward model built from the linear correlations between stimulation frequency and step height. Regression estimates (slope Embedded Image and intercept Embedded Image) were initially set for each subject based on offline analyses. These values were iteratively updated at every gait cycle using online adaptive filters (least mean squares):Embedded Imagewhere μ = 0.1 is the forgetting factor, and ei is the error between recorded and desired outputs at time i. This feed-forward control loop enabled to predict and anticipate the necessary changes in stimulation output. It was combined with a PI feedback loop to correct for unmodeled discrepancies. Values for the proportional (kp) and integral (ki) corrections were set manually for each rat on the first experimental session (kp = 0.2 to 0.35, ki = 0.1 to 0.15), and kept constant during the duration of the recordings. The final stimulation frequency was computed on the basis of the combined feedback and feed-forward components to track a user-defined reference value ri+1 for the maximum foot height during the next gait cycle:Embedded Imagewhere Ei is the cumulative error at time i, and is iteratively calculated as Embedded Image. The reference included a dead band (db = 5 mm), which was included above and below the reference value to account for the inherent variability in stepping when calculating tracking errors: Embedded Image. Final stimulation values were applied within the range of values efficient to promote functional locomotor movements, spanning from 20 Hz onset to 95 Hz saturation (movie S1). The amplitude of the current delivered through each electrode was calibrated at the beginning of each experimental session to obtain optimal facilitation of locomotion at baseline frequency based on visual observation. The real-time platform and control structure are represented and further explained in fig. S4.

Treadmill paradigm

Rats were tested in the same conditions as training. Serotonin agonists were administered 10 min before recording onset. Various experimental conditions were tested during each session that lasted 60 min. The first series of experiments evaluated the functional effects of varying amplitude, pulse width, and frequency of EES. For each session, optimal features for the three parameters of EES to promote robust stepping patterns were searched through visual observations. After this preparatory phase, two features were maintained constant, whereas the tested feature was randomly modulated with fixed increments. A minimum of 10 gait successive cycles were recorded for each experimental condition.

The second series of experiments aimed to evaluate the robustness of the developed control policies. Two testing paradigms were designed to challenge real-time control of foot height. During the first paradigm, termed step increment, the reference band was displaced by a fixed increment (0.5 to 2.5 cm) or decrement (0.5 to 1.5 cm), and maintained at this new reference for six consecutive steps. During the second paradigm, termed triangular waveform, the reference band was gradually displaced by a fixed increment (0.5 to 1.5 cm) or decrement (0.5 to 2.5 cm) at each step until the maximum (7.5 cm) or minimum (3.5 cm) limit of the locomotor system was met. When reaching the limit, the direction of the reference band displacement was inversed (fig. S5).

Locomotor exhaustion

Kinematic and EMG activities were recorded during continuous stepping on a treadmill until exhaustion of locomotor movements. For each rat, recordings stopped when the animal was not able to move the limb forward; both hindlimbs were dragging along the treadmill belt. On the first day, stepping was enabled by serotonin agonists and continuous EES applied overlying L2 and S1 spinal segments at a fixed frequency of 40 Hz. On the second day, the controller was set to modulate EES frequency to maintain step height within a reference band fixed at 4.5 cm. These experiments took place at 7 weeks post-SCI, when substantial exhaustion of locomotor activity occurs clearly (22). The duration of stepping within the reference band, and until exhaustion of locomotion, was measured for each rat and condition.

Staircase paradigm

Rats were positioned in a robotic support system (40) that provided a constant force against the direction of gravity (25% of body weight) while moving the rat forward at a fixed velocity (9 cm/s) along a horizontal runway (Fig. 7). Staircases of various heights (1.3 to 3.5 cm), lengths (12 to 25 cm), and numbers (one or four stairs) were interposed along the runway. For each rat, 100 (single) or 20 (successive) trials were performed along each type of staircase paradigms under noncontrolled or controlled conditions. The percentage of trials during which the rat tumbled, touched, or passed the entire staircase was evaluated offline using video recordings.


All data are reported as means ± SEM, unless otherwise specified. Repeated-measures ANOVAs were used to evaluate differences between normally distributed data (Kolmogorov-Smirnov test) from the various experimental conditions. The two-tailed Kruskal-Wallis test was used for nonparametric evaluations. Post hoc differences were assessed using the Bonferroni test. The statistical α level P < 0.05 was considered significant.



Fig. S1. Completeness of the SCI.

Fig. S2. Modulation of EES frequency tunes multiple aspects of gait patterns.

Fig. S3. Modulation of EES frequency tunes mono- and polysynaptic responses in flexor and extensor muscle during locomotion.

Fig. S4. Real-time monitoring and control platform.

Fig. S5. High-fidelity control of complex foot trajectory.

Fig. S6. Execution along staircases of various heights and lengths.

Table S1. List of computed kinematic, kinetic, and EMG parameters.

Movie S1. Real-time monitoring and control platform.

Movie S2. Closed-loop neuromodulation achieves high-fidelity control of leg movements.

Movie S3. Closed-loop neuromodulation prevents rapid fatigue during continuous locomotion.

Movie S4. Closed-loop neuromodulation enables locomotion across staircases.


  1. Acknowledgments: We thank D. Borton for helpful advice on the manuscript, S. Duis and K. Bartholdi for help with animal care, A. Vuaridel for data processing, and J. Gandar and A. Kuck for their contribution to real-time platform development. Funding: A Starting Grant from the European Research Council (ERC 261247, Walk Again), the European Community’s Seventh Framework Program (CP-IP 258654, NeuWALK), the International Paraplegic Foundation, the Bertarelli Foundation, and funding from the National Center of Competence in Research (NCCR) in Robotics and program of the Swiss National Science Foundation (SpineRepair). Author contributions: N.W., E.M.M., S.R., and M.B. performed the in vivo experiments. N.W., E.M.M., and G.C. analyzed the data. N.W. and E.M.M. conducted the statistical analysis. P.M. implanted the epidural electrodes and placed the SCI, whereas N.W. and G.C. inserted the EMG electrodes. N.W., E.M.M., J.D., M.M., S.M., and G.C. conceived the experiments. N.W., E.M.M., and G.C. prepared the figures with the help of the other authors. N.W., E.M.M., and G.C. wrote the manuscript, and all the authors contributed to its editing. G.C. supervised all aspects of the work. Competing interests: Ongoing European patent application (no. 13191003.6) on the stimulation algorithms. The authors declare no other competing interests. Data and materials availability: Raw data and custom-developed C++ and TDT (Tucker-Davis Technologies) codes will be made available through material transfer agreement upon request to G.C.
View Abstract

Navigate This Article