Abstract
A spinal cord injury (SCI) causes immediate and sustained hemodynamic instability that threatens neurological recovery and impacts quality of life. Here, we establish the clinical burden of chronic hypotensive complications due to SCI and expose the ineffective treatment of these complications with conservative measures. To address this clinical burden, we developed a purpose-built implantable system based on biomimetic epidural electrical stimulation (EES) of the spinal cord that immediately triggered robust pressor responses. The system durably reduced the severity of hypotensive complications in people with SCI, removed the necessity for conservative treatments, improved quality of life, and enabled superior engagement in activities of daily living. Central to the development of this therapy was the head-to-head demonstration in the same participants that EES must target the last three thoracic segments, and not the lumbosacral segments, to achieve the safe and effective regulation of blood pressure in people with SCI. These findings establish a path forward to treat the underappreciated, treatment-resistant hypotensive complications due to SCI.
Introduction
A spinal cord injury (SCI) disrupts the communication between supraspinal vasomotor centers and the regions of the spinal cord that regulate blood pressure1,2. The consequence is pernicious exposure to disabling and potentially life-threatening blood pressure instability that threatens neurological recovery3–8. A prominent feature of this hemodynamic instability is a spectrum of hypotensive complications that include orthostatic hypotension, resting hypotension, postprandial hypotension, and exertional hypotension2. These complications are clinically managed with conservative measures such as uncomfortable abdominal binders, pressure stockings, high sodium diets, and slow-acting pharmacological pressor agents2. Despite these measures, the majority of people who present with hemodynamic instability due to SCI remain chronically hypotensive, which impacts their quality of life and reduces engagement in social and professional activities2,9–12. This symptomatology9,10,13, coupled with the unexplained prevalence of stroke and heart disease in this population, compelled health institutions to encourage the development of new approaches to manage hemodynamic instability in people with SCI14,15.
Recent case studies documented transient increases in blood pressure in response to epidural electrical stimulation (EES) applied over the lumbosacral spinal cord of people with SCI16–21. Since the lumbosacral spinal cord contains a paucity of sympathetic preganglionic neurons that are required for the regulation of blood pressure1,22, these serendipitous observations appeared to be incongruent with the anatomy of the autonomic nervous system22,23. Consequently, the mechanisms underlying these pressor responses remain enigmatic, casting doubt that observations of blood pressure elevations triggered by EES applied over the lum-bosacral spinal cord could translate into a safe and efficacious treatment to manage hemodynamic instability in people with SCI.
To address the unsatisfying management of hemodynamic instability in people with SCI, we previously developed a preclinical model of hemodynamic instability due to SCI with the aim of studying the mechanisms through which EES applied over the spinal cord can lead to increases in blood pressure1,24. We found that EES depolarizes large-diameter afferents where they bend to enter the spinal cord through the dorsal root entry zones1,25. The recruitment of these afferents leads to the activation of splanchnic neurons via the activation of Vsx2-expressing neurons that directly connect to sympathetic preganglionic neurons1,25. Importantly, the amplitude of pressor responses followed a Gaussian distribution wherein the largest responses were elicited when EES was applied over the lower thoracic spinal cord. This distribution correlated with the density of sympathetic preganglionic neurons, in line with the known anatomy of the autonomic nervous system22. We called this region the hemo-dynamic hotspot. We translated this understanding into a neuromodulation strategy operating in closed-loop that adjusts EES amplitudes delivered over the hemodynamic hotspot to regulate blood pressure in real time. This strategy maintained blood pressure stability during transient, varying, and sustained hemodynamic challenges in rat and nonhuman primate models of SCI, and in one human participant1.
Here, we aimed to establish the burden of hypotension and assess the effectiveness of current treatments in large populations of people living with SCI in order to justify new treatments involving neurosurgical interventions. We then resolved the ambiguity on the optimal location of the spinal cord to apply EES to treat hemodynamic instability and exploited this understanding to develop a purpose-built implanted system that harnesses the mechanisms through which EES regulates blood pressure. We finally validate the long-term efficacy of this new system and resulting improvements in quality of life, paving the way towards pivotal clinical trials.
Results
Underappreciated prevalence of medically-refractory hypotensive complications
EES is a promising therapy to address the lack of satisfying management options for hemodynamic instability due to SCI. However, EES requires a neurosurgical intervention that must be weighed against the risks and benefits of the procedure. Establishing this balance requires an understanding of the prevalence, symptomatology, and effectiveness of current management strategies. While hypotension is a recognized medically refractory complication of SCI2, these epidemiological factors have not been quantified conclusively.
To address this knowledge gap, we analyzed the Spinal Cord Injury Community Survey (SCICS)26,27, which includes self-reported information on symptoms of orthostatic hypotension and demographic information in 1,479 individuals living with chronic SCI (Supplementary Fig. 1a). This analysis revealed that 78% of individuals with tetraplegia had been told by a medical practitioner that they have orthostatic hypotension. Of this subset of individuals who were diagnosed, 28% of them were being treated for orthostatic hypotension, yet 91% still experienced symptoms (Fig. 1a). Orthostatic hypotension was still present in half the people with paraplegia (Supplementary Fig. 1b), and analysis of only individuals with complete (ASIA A) injuries led to comparable results (Supplementary Fig. 1c). Together these results suggested that conservative treatments for orthostatic hypotension are insufficient to manage the ongoing symptomatology.
This analysis compelled us to characterize the pattern of hypotension-related symptoms due to SCI. For this purpose, we surveyed a diverse sample of 254 individuals with SCI from 26 different countries, which included responses to the Autonomic Dysfunction Following Spinal Cord Injury (ADFSCI) rating scale. When considering individuals with tetraplegia, we found that every individual who responded to the survey experienced symptoms of hypotension through-out the day. These symptoms were primarily characterized by lightheadedness, dizziness, fatigue, blurred vision, and weakness throughout the day (Fig. 1b). As expected, the frequency and severity of hypotension-related symptoms were higher in people with tetraplegia compared to people with paraplegia, and in those with complete SCI compared to those with incomplete SCI (Supplementary Fig. 1e).
To link these outcomes to objective physiological quantifications, we conducted a formal 10-minute tilt table test in a subset of these individuals (Fig. 1c, Supplementary Fig. 2a-d and Supplementary Table 2). We found that individuals who met the criteria for orthostatic hypotension during the tilt-table test reported significantly more hypotensive symptoms (Fig. 1d-f and Supplementary Fig. 2e-i).
These results established that i) the vast majority of individuals living with tetraplegia experience severe and persistent symptoms of hypotension; ii) hypotension is accompanied by quantifiable symptomatology that increases with the level and completeness of injury, iii) the ADFSCI is a valid self-reported measure of hypotension, and iv) the management of hypotension is insufficient to reduce the symptomatology.
We concluded that these outcomes, combined with the long-term consequences of medically refractory hypotension, justified the consideration of therapeutic strategies involving surgical interventions to manage hypotension due to SCI.
EES is more effective when applied over lower thoracic than lumbosacral segments
Extensive preclinical evidence in mice, rats, and nonhuman primates as well as clinical case studies demonstrated that EES reduces the severity of orthostatic hypotension1,16–19,28. However, the optimal location to deliver EES has remained controversial. We previously showed that EES activates large-diameter afferents where they bend to enter the spinal cord through the dorsal root entry zones and that this recruitment engages sympathetic preganglionic neurons indirectly via Vsx2-expressing excitatory neurons nested in the intermediate lamina of the stimulated region of the spinal cord1,25 (Fig. 2a). We also confirmed an enrichment of these sympathetic preganglionic neurons in the lower thoracic spinal cord (Fig. 2a), and showed that this enrichment coincided with the maximal pressor response to EES. This location was thus termed the hemodynamic hotspot1. Despite this evidence, many studies reported unexpected pressor responses during the application of EES over the lumbosacral spinal cord16–19. The incongruence of these observations with the anatomy of the autonomic nervous system22,23 compelled us to perform a head-to-head comparison in the same participants to establish the optimal location to deliver EES for the regulation of blood pressure in humans with SCI.
Consequently, we initiated two clinical studies in Canada and Switzerland to compare the effect of EES applied over the lower thoracic versus lumbosacral spinal cord on hemodynamic stability (ClinicalTrials.gov ID: NCT04994886, NCT05044923). We enrolled 6 individuals with medically-refractory orthostatic hypotension secondary to their SCI (Supplementary Fig. 3a and Table 1). The six participants included in these clinical study met the established clinical criteria for orthostatic hypotension since none of these participants were able to sustain the orthostatic challenge imposed by verticalization on a tilt table test, showing an average drop of 50 mmHg in systolic blood pressure within 3 minutes (Supplementary Fig. 3b-f). This drop far exceeds the international criteria for orthostatic hypotension, set at a 20 mmHg drop of systolic blood pressure within 3 minutes29.
Since only repurposed technologies originally designed to alleviate neuropathic pain were available to target both locations (Fig. 2b and Supplementary Fig. 3g), we sought to optimize the surgical implantation of two separate paddle leads over the lower thoracic spinal cord and the lumbosacral spinal cord. To facilitate this placement, we generated computational models that we personalized with computed to-mography and magnetic resonance imaging sequences that we previously optimized to enable the identification of the dorsal root entry zones targeted by EES30 (Fig. 2b and Supplementary Fig. 4). We then injected these coordinates into our preoperative planning software that we consolidated with intraoperative imaging of the spinal column. Under general anesthesia, we performed laminotomies to enable the insertion of both leads over the planned anatomical locations. We then fine-tuned the final position of both leads based on the muscle responses to EES that are expected from the known rostrocaudal distribution of motor neuron pools innervating trunk and leg muscles30 (Supplementary Fig. 3h-j).
Once the location of both leads was secured, we applied EES (120 Hz, 1-10 mA, 300 μs) over each pair of electrodes on each row of both paddles to assess pressor responses that we monitored by measuring beat-by-beat blood pressure through an arterial line. We mapped pressor responses to post-operative reconstructions of the spinal cord regions activated by each pair of electrodes. As previously established in rat and nonhuman primate models of SCI1, we found that pressor responses followed a Gaussian distribution, wherein the peak pressor responses were centered over the penultimate segment of the thoracic spinal cord. In contrast, modulation of blood pressure was poorly effective when EES was applied over the lumbosacral segments (Fig. 2c-e and Supplementary Fig. 3k-l).
Postoperative computed tomography acquisitions confirmed that both leads were positioned over the planned locations (Fig. 2f). Intraoperative assessments confirmed that the optimal location to elicit pressor responses with EES was centered over the lower thoracic spinal cord—the hemodynamic hotspot. Therefore, we next verified whether EES applied over the hemodynamic hotspot was superior to EES applied over the lumbosacral spinal cord to stabilize blood pressure during formal tilt-table testing post-operatively, since previous case studies16–19 reported a mitigation of orthostatic hypotension when EES was applied over the lumbosacral segments (Fig. 2g-i and Supplementary Fig. 5a-g).
We conducted comprehensive mapping sessions to identify optimal configurations of electrodes and stimulation parameters that could trigger pressor responses when EES was delivered over the lumbosacral segments (Supplementary Fig. 5a-d). Despite this mapping, we failed to identify parameters of EES that could trigger robust and reproducible pressor responses without eliciting spasms in leg muscles. EES applied over the lumbosacral segments only marginally reduced the severity of orthostatic hypotension (Fig. 2g-i and Supplementary Fig. 5e-g).
This failure contrasted with the robust pressor responses elicited when EES was applied over the lower thoracic spinal cord, which enabled all the participants to improve their tolerance to orthostatic challenges induced by the verticalization during tilt-table tests in the absence of undesirable sensations (Fig. 2g-i and Supplementary Fig. 5a-g).
These results demonstrate that EES applied over the lumbosacral spinal cord failed to consistently trigger meaningful pressor responses, as expected based on the incongruence between this location and the anatomy and physiology of the autonomic nervous system. Instead, EES applied over the hemodynamic hotspot, wherein an enrichment of sympathetic preganglionic neurons involved in the modulation of blood pressure reside, consistently triggered robust pressor responses that enabled all six patients to tolerate severe orthostatic challenges.
These results provide a definitive answer to the debated optimal location of EES for effective treatment of hemodynamic instability in people with SCI. In addition, recent preclinical experiments demonstrated that long-term application of EES over the lumbosacral spinal cord augmented the density of maladaptive synaptic connections from Vsx2-expressing neurons of the lumbosacral to Vsx2-expressing neurons of the lower thoracic spinal cord, which reinforces the severity of autonomic dysreflexia25. Given the discrepancy in safety and efficacy of EES applied over the hemodynamic hotspot versus lumbosacral segments, these results establish the path forward on how to deliver EES to treat hemodynamic instability due to SCI.
Requirements for a system to restore hemodynamic stability after SCI
Having identified the optimal location over which EES must be delivered to regulate blood pressure, we next sought to develop an implanted system that leverages this understanding with the goal of establishing a therapy to improve the management of hemodynamic instability in people with SCI. We reasoned that such a system must combine the following features.
First, the implantable neurostimulation platform should include the capacity to adjust EES in closed-loop. Indeed, while continuous EES is sufficient to stabilize blood pressure in static conditions, stabilization of blood pressure during dynamic changes in body orientation requires constant adjustment of EES amplitudes1. Moreover, constant adjustment of EES would be ideal in acute care settings wherein blood pressure lability impacts the neurological recovery of patients5–8. Second, the paddle lead must include an optimal configuration of electrodes to target the ensemble of dorsal root entry zones innervating the hemodynamic hotspot, while avoiding the undesired recruitment of sensory axons ascending in the dorsal columns. Since functional and anatomical experiments in mice, rats, nonhuman primates, and six humans revealed that the correct location to deliver EES is centered over the last three thoracic segments of the spinal cord, we concluded that the paddle lead must cover the entire extent of this hemodynamic hotspot.
Third, the hardware and software components of the system must incorporate practical interfaces to support the rapid configuration of EES by physicians, and enable patients to operate the therapy safely and independently despite accessibility limitations.
We aimed to develop a system that fulfilled these requirements and to establish the safety and efficacy of each component embedded in this new purpose-built implantable system through a series of stepwise clinical validations.
Purpose-built neurostimulation platform to restore hemodynamic stability after SCI
We aimed to develop a clinical-grade implantable neurostimulation platform that addresses the limitations of existing systems and meets the identified requirements for the safe and effective management of hemodynamic stability in people with SCI.
Since the possibility to adjust EES parameters in closed-loop was a desired feature of the platform, we designed an implantable neurostimulation platform that continuously communicates with a wearable hub through near field magnetic induction (NFMI) (Fig. 3a). This hub can acquire external signals, including continuous measurements of blood pressure, and embeds computing power to support closed-loop control algorithms. This configuration enables closed-loop control of EES parameters with latencies as low as 25 ms. We programmed the firmware to support any configuration of anodes and cathodes, and to resolve the limited range of programmable EES parameters in current neurostimulation platforms.
The hub communicates wirelessly with a touch-screen tablet and a smartwatch, from which therapeutic programs can be selected and parameterized (Fig. 3a). We optimized these user interfaces to comply with usability assessments and user needs. We also endowed the user interface with practical features that reduce the impact of impaired arm and hand functions for operating the therapy independently and safely. In addition, customizable algorithms enable real time adjustment of set-point amplitudes based on hemodynamic requirements that are well known to differ throughout the day, such as after meals and during sleep and exercise2,31–33.
After meeting the requirements for all applicable standards for active medical device development, we tested the preliminary safety and efficacy of this new platform under the protocol of a new clinical study conducted in Switzerland (ClinicalTrials.gov ID: NCT05111093) (Supplementary Fig. 6a). We enrolled 4 new participants who presented with debilitating medically-refractory hypotension (Supplementary Fig. 6b-f and Table 1) and implanted the purpose-built neurostimulation platform, but maintained every other aspect of our preoperative and operative procedures, including the implantation of the same repurposed paddle lead over the hemodynamic hotspot as in our previous experiments (Supplementary Fig. 3g-h).
In all the participants, the implanted neurostimulation platform met the requirements to deliver EES waveforms that elicit robust pressor responses (Supplementary Fig. 6i-l) that enabled all the participants to improve their tolerance to orthostatic challenges (Fig. 3b-d and Supplementary Fig. 6m-s). As we previously reported in rats and nonhuman primates1, we confirmed in all four participants that the magnitude of pressor responses scaled linearly with increasing amplitudes of EES (Fig. 3f and Supplementary Fig. 7a). This linear relationship suggested that blood pressure could be regulated in real time based on a proportional-integral-differential controller regulating EES amplitudes to maintain blood pressure within user-defined levels.
To test this possibility, we leveraged the closed-loop capabilities and versatile software environment of the new neurostimulation platform to implement this proportional-integral-differential controller (Fig. 3e). We then evaluated the performance of this controller during formal and dynamic tilt-tests. As expected based on previous experiments (Fig. 2g-i), we found that the delivery of continuous EES decreased the severity of orthostatic hypotension during formal tilt-tests wherein the amplitude of EES remained constant (Fig. 3b-d). However, dynamic tilt-testing sessions exposed the limitations of continuous EES, since pseudorandom changes in body orientation with respect to the direction of gravity led to dynamic perturbations of blood pressure that could not be managed satisfactorily with continuous EES. In contrast, closed-loop control of EES led to rapid and precise updates to EES amplitudes that translated into finely-tuned levels of blood pressure during dynamic tilt-tests (Fig. 3g and Supplementary Fig. 7b). Closed-loop control of EES applied over the hemodynamic hotspot maintained blood pressure levels within the desired target range for extended periods of time despite transient and dynamic changes in postural orientation (Figure 3g, Supplementary Fig. 7b and Supplementary Fig. 11n). Quantifications also revealed the superiority of EES controlled in closed-loop over both continuous EES and alternative closed-loop approaches relying on inertial measurement units or manual adjustments of EES (Fig. 3h, Supplementary Fig. 7c and Supplementary Fig.11o).
These results confirmed that the newly developed neurostimulation platform fulfilled the technical requirements to modulate blood pressure with EES, and established the possibility of achieving precise control of blood pressure in closed-loop despite transient, varying, and sustained hemodynamic challenges in people with orthostatic hypotension due to SCI.
A new paddle lead targets the hemodynamic hotspot
We next aimed to develop a clinical-grade paddle lead that combines the necessary features to target the large-diameter afferents projecting to the hemodynamic hotspot while avoiding the undesired side effects of non-selective stimulation.
Our mapping experiments showed that the hemodynamic hotspot is centered over the last three segments of the lower thoracic spinal cord. Moreover, the same experiments revealed that EES applied over multiple segments leads to cumulative increases in the amplitude of pressor responses compared to when the same amplitude of EES is applied over a single spinal segment (Fig. 4a).
We thus reasoned that the optimal paddle lead to regulate blood pressure with EES must maximize the coverage of the last three thoracic segments of the spinal cord. However, our analysis of the 10 individuals who were implanted with the repurposed paddle revealed that the coverage of these segments only reached 56% over the targeted region (Fig. 4b and Supplementary Fig. 8a). Moreover, 6 of the 16 electrodes were located over the midline of the spinal cord. These electrodes elicited undesired muscle spasms in the legs through antidromic volleys elicited along sensory axons ascending the dorsal columns (Supplementary Fig. 8b). This analysis emphasized the necessity to design a new paddle lead with a configuration of electrodes that meets the requirements to modulate blood pressure with EES.
To guide the informed design of this lead, we leveraged our digital library of human thoracic spinal cords that were reconstructed from magnetic resonance imaging acquisitions and combined these reconstructions with anatomical values from the literature30,34–39. These anatomical quantifications uncovered the requirements to arrange two columns of equally distributed electrodes that could cover 100% of the last three segments of the human thoracic spinal cord in more than 95% of the population (Fig. 4b). In addition, we conducted computer simulations to identify the optimal mediolateral location of electrodes that would target the dorsal root entry zones innervating the lower thoracic spinal cord while avoiding the undesired recruitment of the sensory axons in the dorsal columns. These simulations combined with quantifications from our anatomical models of the spinal cord revealed that the optimal tradeoff for the placement of these electrodes within the anatomical constraints of the spinal column was 3.7 mm lateral to the dorsal root entry zones (Fig. 4c and Supplementary Fig. 8c-d). We embedded the identified configuration of electrodes into the design of a new paddle lead that we fabricated with conventional medical-grade technologies and validated for biocompatibility (Fig. 4d).
The implanted system that regulates hemodynamics after SCI
Having developed the complete implantable system with the required features to regulate blood pressure with EES, we amended our existing clinical study (ClinicalTrials.gov ID: NCT05111093) to interface the new paddle lead with the purpose-built neurostimulation platform (Supplementary Fig. 9a).
We enrolled 3 new participants who presented with debilitating hypotension (Supplementary Fig. 9b-f and Table 1) and implanted them with the complete new system following the same procedure as described above. The extended length of the lead and optimal coverage of the dorsal root entry zones innervating the hemodynamic hotspot simplified the preoperative planning and intraoperative placement of the lead (Supplementary Fig. 9g-h). Indeed, the extended length avoided compromises on whether the rostral or caudal regions of the hemodynamic hotspot had to be targeted with shorter leads, while the optimized mediolateral electrode locations ensured the selective recruitment of the left versus right dorsal root entry zones (Fig. 4c and Supplementary Fig. 8c-d). These improved features of the paddle lead minimized the need for intraoperative fine-tuning of the paddle lead position. Updates to the personalized computational models of the spinal cord with postoperative imaging data confirmed that the lead covered the entire extent of the hemodynamic hotspot and predicted selective activation of the dorsal root entry zones in all participants (Fig. 5b).
The optimized location of the electrodes also eliminated the need for extensive and cumbersome mapping sessions to identify the optimal configurations and parameters to elicit pressor responses40. Instead, minimal adjustments of prede-termined electrode configurations (three pairs of electrodes covering the hemodynamic hotspot) and parameters (120 Hz frequency, 0.3 ms pulse duration, up to 20 mA amplitude) were sufficient to trigger robust, safe, and well-controlled increases in blood pressure during the first postoperative stimulation session (Supplementary Fig. 9i-l).
As observed in the previous 10 participants, EES applied over the hemodynamic hotspot led to pressor responses that scaled linearly with increasing amplitudes of EES (Supplementary Fig. 9u) and reduced the severity of orthostatic hypotension during a formal 10-minute tilt-table test (Fig 5c and Supplementary Fig. 9m-t). We next assessed the possibility of controlling blood pressure in closed-loop during formal and dynamic tilt-table testing sessions (Fig. 5d). We found that the complete purpose-built system maintained the desired blood pressure levels with high-fidelity despite varying orthostatic challenges (Fig. 5e).
Independent validation of the implanted system
We next asked whether our therapeutic strategy could be deployed independently by a center that had never been exposed to the implanted system and therapy. To assess this scalability, we conducted a new clinical study in the Netherlands (ClinicalTrials.gov ID: NCT05941819) (Supplementary Fig. 10a).
A neurosurgeon who had not been previously involved in the studies followed our procedures to implant the complete system in a participant who presented with severe orthostatic hypotension (Supplementary Fig. 10b-f and Table 1). As observed in all the previous participants, intraoperative delivery of EES over the hemodynamic hotspot triggered robust pressor responses. Postoperative imaging confirmed the location of the lead over the hemodynamic hotspot (Supplementary Fig. 10g-h).
During the first postoperative session, we asked a physician without prior experience with EES to program the therapy based on anatomical principles to determine effective electrode configurations and parameters to modulate blood pressure with the system. The non-expert was able to program the therapy in only a few minutes. As early as the first session postoperatively, EES delivered over the hemodynamic hotspot in a seated position elicited robust pressor responses that increased the systolic blood pressure from 86 mmHg to 105 mmHg (Supplementary Fig. 10i). With EES, the participant was able to complete the 10-min tilt-table test, which was not possible without EES (Supplementary Fig. 10j-n). Closed-loop control of EES further improved the maintenance of blood pressure despite orthostatic challenges of varying severities (Supplementary Fig. 10o-q).
Independent use of the system for blood pressure management
We concluded that the safety and efficacy profile of the new purpose-built implantable system and therapy were sufficient to enable the participants (n = 8 participants) to use the therapy independently at home, and thus to evaluate the long-term consequences of this therapy for the management of blood pressure after SCI. To support this evaluation, the participants were sent home with stimulation programs that were tailored to their individual pattern of hypotensive symptoms by a trained physician.
Long-term efficacy of EES applied over the hemodynamic hotspot
Although the various participants included in our clinical studies were implanted with successively improved technologies, this diverse cohort of female and male participants spanning 18 to 74 years of age (Supplementary Table 1) across 3 sites provided the opportunity to conduct a preliminary evaluation of the safety and effectiveness of EES applied over the hemodynamic hotspot to manage hemodynamic instability in people living with chronic SCI.
When evaluated up to two years after implantation, participants still exhibited severe orthostatic hypotension during a formal tilt table test when EES was turned off. They also continued exhibiting robust hemodynamic stability when EES was turned on, which allowed them to complete the tilttable test (Fig 6a-b and Supplementary Fig. 11f-k). More-over, formal evaluations confirmed the persistence of blood pressure stabilization evaluations during constant-varying orthostatic challenges when EES was controlled in closed-loop (Supplementary Fig. 11m-o). Pressor responses during tilt test with EES were associated with increased concentrations of circulating norepinephrine (Supplementary Fig. 11l), augmented splanchnic volume (Supplementary Fig. 12i), and increased sympathetic nerve activity (Supplementary Fig. 12j), which confirmed the activation of the sympathetic nervous system when EES was applied over the hemodynamic hotspot.
Every tested participant continued demonstrating controllable pressor responses when EES was delivered over the hemodynamic hotspot while seated in their wheelchair, corresponding to the conditions of their daily life (Fig. 6c and Supplementary Fig. 11a-e). Pressor responses in seated conditions coincided with an increase in cerebral blood flow (Fig. 6d).
These objective physiological outcomes also concurred with a reduction in self-reported measures of orthostatic hy-potension severity (Fig. 6e). These results showed that EES applied over the hemodynamic hotspot led to an immediate and durable reduction in the severity of hypotension in this cohort of individuals with SCI.
Effective management of blood pressure improves quality of life
Our clinical outcomes thus far highlighted the durability of EES applied over the hemodynamic hotspot to manage blood pressure instability in people with SCI, but whether the resulting long-term maintenance of blood pressure stability translated into an improvement in quality of life and an increased possibility to participate in activities of daily living remained unclear.
To address this question, we asked whether the participants relied on the therapy at home to manage their blood pressure throughout their daily life. We quantified the use of the therapy from the stimulation logs, and found that the participants turned the stimulation on throughout their waking hours, with increased usage to facilitate verticalization in the morning routine and to reduce postprandial hypotension after meals (Supplementary Fig. 12a). We confirmed the improvement in postprandial hypotension during in-clinic evaluations with controlled food intake (Supplementary Fig. 12b). The daily use of the therapy was accompanied by the medically-supervised cessation of all other conservative treatment options that participants had been relying on to manage hypotension, including abdominal binders and midodrine (Supplementary Fig. 12d). However, these improvements only manifested in the presence of the stimulation, since none of the participants demonstrated a reduction in the severity of orthostatic hypotension during formal tilt-tests when EES was turned off (Fig. 6a-b and Supplementary Fig. 11g-k). Evaluations of system usability41 confirmed that the participants could operate the therapy safely and effectively despite impairements in arm and hand movements (Supplementary Fig. 12c).
These multifaceted ameliorations translated into improvements in quality of life, as measured by formal outcome measures (Fig. 6e). Quantified increases in quality of life coincided with other subjective improvements. For example, one participant reported that she had regained the energy to re-initiate graduate school, and two other participants felt they were once again capable of returning to work (Supplementary Fig. 12k).
The beneficial impact of the therapy also extended to other aspects of life that could be potentially influenced by improved management of blood pressure and autonomic functions. For example, longitudinal outcome measures revealed improvements in sleep quality, bowel management, spasticity, and engagement in social activities (Fig. 6f-i and Supplementary Fig. 12e).
Since the participants regained the ability to sustain verticalization for extended periods of time when the system was turned on, they could engage in rehabilitation exercises that involved standing or walking in a motorized exoskeleton42 (Fig. 6i and Supplementary Fig. 12f). These tasks were not accessible with the system turned off since they experienced hypotensive symptoms in less than 2 minutes during verticalization.
The last three thoracic segments contain motor neurons innervating the trunk musculature. Consequently, EES applied over these segments activated trunk muscles. We reasoned that this added benefit of applying EES over thoracic segments could be leveraged to configure stimulation programs that concomitantly restore hemodynamic stability and improve neurological functions dependent on trunk and abdominal musculature. We mapped the optimal locations of EES to activate specific ensembles of trunk and abdominal muscles, which allowed us to customize stimulation programs to reinforce the activation of the trunk and abdominal musculature, aiming at improving trunk stability and coughing (Supplementary Fig. 12g-h). Quantifications of these neurological functions in well-controlled laboratory settings confirmed that these stimulation programs led to significant improvements (Fig. 6i and Supplementary Fig. 12g-h).
Participants leveraged these new stimulation programs to support activities of daily living. For example, one participant turned on trunk stability stimulation programs to facilitate the upright position of his trunk during meals, which freed his hands to manipulate his fork and knife without experiencing instability (Fig. 6f). Another participant used these programs to engage in recreational skiing activities since both his blood pressure and trunk stability were sufficient to enjoy these rides safely (Fig. 6i).
Safety of the newly-developed system and therapy
Despite decades of experience in the management of neuropathic pain that have established the safety profile of EES, we systematically monitored for device- and procedure-related adverse events throughout the clinical studies.
All the participants have been operating the system to regulate their hemodynamic stability daily for up to two years. They were asked to report any safety event, including during home-use. All the serious and non-serious adverse device effects (S/ADE) reported across the 5 clinical studies for each participant are quantified and described in Table 2 and Supplementary Tables 3 - 5.
Across the 8 participants implanted with the newly-developed system or components of this system, no serious adverse events related to the device were reported. Only one serious adverse event involving an uncontrolled increase in blood pressure due to the inadvertent removal of a urinary probe occurred while the system was turned off (Supplementary Table 5). Uncontrolled increases in blood pressure during EES occurred only once across hundreds of testing sessions conducted over a period that spanned several months and involved the delivery of EES across a broad range of parameters. This non-serious device-related adverse event occurred while targeting the trunk musculature, but ceased as soon as EES was switched off.
Discussion
We conducted two epidemiological analyses and five complementary clinical studies involving implantable medical devices that i) established the burden of hypotension complications in people with SCI, ii) resolved the ambiguity on the optimal location of the spinal cord to apply EES for restoring hemodynamic stability, iii) enabled the stepwise development of a purpose-built system harnessing the understanding of the mechanisms through which EES regulates hemo-dynamics in people with SCI, and iv) validated the preliminary safety and efficacy of this system to reduce hypotensive complications throughout the day and thus improve quality of life.
Before considering the development of a therapy based on EES to restore hemodynamic stability after SCI, we thought that it was essential to resolve the controversy on the optimal location to restore hemodynamics with EES16–21. We unequivocally demonstrate the superiority of EES applied to the hemodynamic hotspot compared to EES applied to the lumbosacral spinal cord. This result was not surprising since attempts to control hemodynamics by stimulating the lumbar region of the spinal cord are incongruent with the basic and well-conserved anatomy of the ancestral autonomic nervous system. To date, the mechanisms underlying the fortuitous observations of pressor responses with EES applied over the lumbosacral spinal cord remain unknown. We surmised that two mechanisms may explain these observations.
The first mechanism that we previously exposed in pre-clinical models is the off-target activation of the dorsal roots innervating the last thoracic segment as they exit the spinal column through the T12 intervertebral foramen1, which roughly coincides with upper lumbar spinal cord segments. Our mapping experiments demonstrate that the optimal coverage of the last three segments of the thoracic spinal cord, especially the penultimate thoracic segment, is essential to maximize pressor responses to EES. Consequently, the inefficient recruitment of only the last thoracic root is unlikely to be an effective and reproducible methodology to manage hemodynamic instability across the population of people with SCI.
The second mechanism, which is the most concerning for the safety of people with SCI, is the risk that EES applied over the lower aspect of the lumbosacral spinal cord directly activates the pathways that cause life-threatening autonomic dysreflexia2,43. Sensory information originating from the overstimulated afferents of the bladder and bowel are the most common triggers of autonomic dysreflexia in people with SCI. Afferents from these visceral organs project to the spinal cord through the lower lumbar and upper sacral dorsal roots, which are directly recruited when applying EES over the lumbosacral spinal cord, especially with configurations targeting the lower lumbar dorsal roots18,21,40. Since no other neuronal populations that contribute to the control of blood pressure exist in this region off the spinal cord1,22, the only explanation for pressor responses elicited by the application of EES over the lower lumbar and sacral spinal cord is the activation of the pathways that can cause autonomic dysreflexia. Indeed, preclinical evidence demonstrated that daily application of EES over the lower lumbosacral segments in rodents with complete upper thoracic SCI augments the density of axonal projections from Vsx2-expressing neurons in the lumbosacral spinal cord onto Vsx2-expressing neurons in the lower thoracic spinal cord25. These neurons are embedded in the neuronal architecture that causes autonomic dysreflexia. Consequently, this aberrant exacerbation of the neuronal architecture responsible for autonomic dysreflexia doubled the severity of pressor responses that characterize this life-threatening condition25. This mechanism may explain the reported improvements of resting hypotension and orthostatic hypotension in individuals with SCI who underwent extensive stimulation of the lower lumbosacral segments.
Together, robust results obtained in a total of 14 participants in 5 clinical studies across 3 countries conducted at 3 independent clinical centers with independent clinical teams demonstrate the preliminary safety and efficacy of EES applied over the hemodynamic hotspot to treat hemodynamic instability in people living with chronic SCI. These findings establish a path forward to treat the underappreciated, treatment-resistant hypotensive complications due to SCI.
Methods
Epidemological analysis
Spinal Cord Injury Community Survey (SCICS). Ethical approval was obtained from an independent ethics board (Veritas Independent Review Board) and the Research Ethics Board of Université Laval (principal investigator’s institution). Ethical approval from local research ethics boards was also obtained to recruit from spinal cord injury centers across Canada. Individuals with SCI (n = 1,479) across Canada were recruited using a national consumer awareness campaign and provided written informed consent to participate26,27. The survey consisted of a series of variables identified by healthcare and service providers, researchers, as well as individuals with Study participants. Fourteen individuals who had suffered a SCI, including demographics, secondary health complications, comorbidities, SCI-related needs, healthcare utilization, community participation, quality of life, as well as overall health ratings26,27. Participants were asked how often they experienced symptoms such as light-headedness and/or dizziness consistent with symptoms of orthostatic hypotension in the past 12 months. Responses were ranked on a 6-point ordinal severity scale ranging from zero (i.e., “Never”) to 5 (i.e., “Every day”). Participants were also asked if they received or sought out treatment concerning these symptoms on a two-point scale (i.e., “Yes” or “No”), along with the degree to which it limited activities from zero (i.e., “Never”) to 5 (i.e., “Every day”). Participants were also asked if they had experienced specific problems, such as heart disease, in the past 12 months. Participants’ American Spinal Injury Association Impairment Scale (ASIA) was estimated using responses to questions about lesion level and sensorimotor/mobility capabilities. A binary approach was used for the evaluation of outcome variables including the level of injury (i.e., cervical SCI vs. non-cervical SCI), the severity of injury (i.e., functionally complete vs. incomplete), presence of orthostatic hypotension (i.e., yes vs. no), and orthostatic hypotensive symptoms (i.e., yes vs. no). For variables ranked on a 6-point ordinal scale, lower scores (i.e., 0-3) were categorized as “No” and higher scores (i.e., 4-5) were categorized as “Yes.”
Autonomic dysfunction in a sample of individuals living with SCI. 254 participants were surveyed with the ADFSCI questionnaire (see Autonomic Dysfunction Following Spinal Cord Injury) to measure the severity and frequency of orthostatic hypotension in a diverse sample of individuals with SCI. From this population, 17 individuals were pseudo-randomly selected to participate in a single formal, 10-minute tilt-table test (see Orthostatic challenge with tilt-table test). Hemodynamic measurements were taken via finger plethysmography (see Post-operative hemodynamic monitoring), and data was analogously acquired with the LabChart Software (ADInstruments, Dunedin, New Zealand) and upsampled to 1000 Hz.
Clinical trials
Study designs and objectives. All experiments were carried out as part of 5 clinical safety (primary objective) and preliminary efficacy (secondary objectives) trials: STIMO-HEMO (NCT04994886, CHUV, Lausanne, Switzerland), HEMO (NCT05044923, University of Calgary, C algary, Canada), HemON (NCT05111093, CHUV, Lausanne, Switzerland), Amendement of HemON (NCT05111093, CHUV, Lausanne, Switzerland) and HemON-NL (NCT05941819, Sint Maartenskliniek, Ubbergen, The Netherlands). All 5 trials received approval by the local ethical committees and national competent authorities. A summary of information related to each study is provided in Supplementary Table 1. All the participants signed a written informed consent at enrollment in the trial. All participants had the option to indicate consent for the publication of identifiable images or videos. All surgical and experimental procedures were performed at the investigational hospital sites, including the neurosurgery department of Lausanne University Hospital (CHUV) in Lausanne, Switzerland, the neurosurgery department of Foothills Medical Center in Calgary, Canada, and the neurosurgery department of the Radboud University Medical Center in Nijmegen, The Netherlands. The study involved eligibility and baseline assessments conducted before the surgery, the neuro-surgical implantation of an investigational device, a postoperative period during which EES protocols were configured, and long-term follow-up periods, including independent use of EES at home. To date, 13 individuals completed the intensive part of the study.
Study participants. Fourteen individuals who had suffered a traumatic SCI participated in one of the 5 studies. Demographic data and neurological status were evaluated according to the International Standards for Neurological Classification of Spinal Cord Injury44, and are reported in Table 1 and Supplementary Data 8. For the STIMO-HEMO and HEMO clinical trials, participants were eligible if they were aged between 18 and 70 years old, had a radiologically confirmed SCI located between C3 and T6 (inclusive), were neurologically classified as A or B on the American Spinal Cord Injury Association (ASIA) Impairment Scale (AIS), had stable medical, physical and psychological condition under the discretion of the investigators, had the injury for at least 1 year, did not have any required spinal instrumentation within the past 6 months, were able and willing to attend all scheduled appointments, and had confirmed orthostatic hypotension and autonomic dysreflexia. Exclusion criteria included individuals in any emergency situations, participants with any diseases and conditions that would increase the morbidity and mortality of the spinal cord injury surgery, participants were unable to withhold antiplatelet, or anticoagulation agents perioperatively, had history of myocardial infarction or cere-brovascular events within the past 6 months, had a current, or anticipated need for opiod pain medications, had pain that would prevent full participation in the rehabilitation program in the judgment of the investigators, had a clinically significant mental illness in the judgment of the investigators, received botulinum toxin injections within 6 months of enrollment, presence of significant pressure ulcers, recurrent urinary tract infection refractory to antibiotics, un-healed spinal fractures, known or suspected drug or alcohol abuse, presence of indwelling baclofen or insulin pump, and any other conditions that would make the subject unable to participate in testing in the judgment of the investigators. Women pregnant during screening, or breastfeeding were also excluded.
For the HemON clinical trial, eligible participants were 18 years of age or older, had a SCI located between C3 and T6 (inclusive), had their injury for at least 1 month, had stable medical, physical, and psychological conditions as considered by the investigators, spoke the French or English language to interact with the study team, available to participate in good faith and attend scheduled appointments per study conditions, and had confirmed orthostatic hypotension. Exclusion criteria include SCI related to neurodegenerative disease, diseases and conditions that would increase the morbidity and mortality of the surgery, the inability to withhold antiplatelet, or anti-coagulation agents perioperatively, history of myocardial infarction or cerebrovascular events within the past 6 months, other conditions that would make the subject unable to participate in testing in the judgment of the investigators, clinically significant mental illness in the judgment of the investigators, having received botulinum toxin vesical and non-vesical injections within 3 months of enrollment, presence of significant pressure ulcers, recurrent urinary tract infections refractory to antibiotics, presence of indwelling baclofen or insulin pumps, other clinically significant concomitant disease states including but not limited to renal failure, hepatic dysfunction, cardiovascular disease, etc., inability to follow study procedures due to language problems, psychological disorders, or dementia, participation in another study with investigational drugs within the 30 days preceding and during the present study, and any relation to the investigator, his/her family members, employees, and other dependent persons. Women who are pregnant or breastfeeding, lack safe contraception if with childbearing capacity, or had the intention to become pregnant during the course of the study were excluded.
For the HemON-NL trial, eligible participants were 18 years of age or older had a SCI located between C3 and T6 (inclusive), had a traumatic SCI, had their injury for at least 1 month, were neurologically classified as AIS A, B, C or D, had stable medical, physical and psychological conditions as considered by the investigators, spoke the Dutch or English language to interact and understand the study team, available to participate in good faith and attend scheduled appointments per study conditions, had if needed continuous support from personal caregiver in daily life during the study visit (including independent transport) and had confirmed orthostatic hypotension. Exclusion criteria include SCI related to diseases and conditions that would increase the morbidity and mortality of the surgery, diseases and conditions that would require regular magnetic resonance imaging (MRI), the inability to perform an MRI due to metal, magnetic or electrical device in the body (e.g. oral implant with magnet, metal splinter, neurostimulator, artificial heart valve, clips, stents etc.) as assessed by the MRI form of Sint Maartenskliniek, the inability to withhold antiplatelet, or anticoagulation agents perioperatively, history of myocardial infarction or cerebrovascular events within the past 6 months, other conditions that would make the subject unable to participate in testing in the judgment of the investigators, clinically significant mental illness in the judgment of the investigators, having received botulinum toxin vesical and non-vesical injections within 3 months of enrollment, presence of significant pressure ulcers, recurrent urinary tract infections refractory to antibiotics, presence of indwelling baclofen or insulin pumps, other clinically significant concomitant disease states including but not limited to renal failure, hepatic dysfunction, cardiovascular disease, etc., inability to follow study procedures due to language problems, psychological disorders, or dementia, participation in another study with investigational drugs within the 30 days preceding and during the present study, and any relation to the investigator, his/her family members, employees, and other dependent persons. Women who are pregnant or breastfeeding, lack safe contraception if with childbearing capacity, or had the intention to become pregnant during the course of the study were excluded.
Safety outcomes. All the feasibility clinical trials outlined in Supplementary Table 1 evaluated the safety and preliminary efficacy of the investigational system and/or therapy to restore hemodynamic stability. The primary outcome measures assess the occurrences of serious adverse events and adverse events that are deemed related or possibly related to the investigational system or study procedure during the course of studies. All reportable safety events have been reported per local regulation to the respective authorities (Health Research Ethics Board at the University of Calgary, Ethics Committee CER-VD, and competent authorities Swissmedic in Switzerland and METC Oost-Nederland in Nijmegen). Data safety monitoring boards have been appointed in all studies and consulted on a yearly basis or adHoc in case of reportable events. All adverse device effects for the 5 clinical trials are summarized in Supplementary Tables 3-5. Two serious adverse events related to the device leading to explantation occurred during the course of the HEMO study due to postoperative infections. Device explantations were reported to local authorities and to the manufacturer. Both were resolved without sequelae following explantation surgery. In studies evaluating the safety of the new implanted device, one reportable serious adverse event related to study procedure occurred (subpubic catheter came out during a rehabilitation session while the stimulation was turned off).
Imaging Data. All the participants underwent structural 3T magnetic resonance imaging (MRI) and computed tomography (CT) imaging of the thoracolumbar spine to enable the elaboration of a personalized anatomical model of the spine and guide the preoperative planning. MRI acquisitions were conducted in a supine position with both arms on the sides. Prior to image acquisition, shim boxes were applied to correct for magnetic field inhomogeneities. All sequences were performed without gadolinium-based contrast agent administration. In the STIMO-HEMO and HemON clinical trials, MRI was performed on a Magnetom PrismaFit (Siemens Healthineers, Erlangen, Germany) with 18-channel body and 32-channel spine array coils. The standard MRI protocol comprised the following four pulse sequences: a) 2D sagittal T2-weighted turbo spin-echo (repetition time (TR), 3080 msec; echo time (TE), 98 msec; voxel size, 0.6 × 0.6 × 3 mm3); b) T2-weighted SPACE (Sampling Perfection with Application-optimized Contrasts using different flip angle Evolution) sequence (TR, 1500 msec; TE, 135 msec; interpolated voxel size, 0.4 x 0.4 x 0.8 mm3); c) 3D axial T2-weighted SPACE with ZOOMit (dynamic excitation pulses to achieve selective/zoomed field-of-view) software (TR, 2500 msec; TE, 106 msec; interpolated voxel size, 0.3 × 0.3 × 0.5 mm3); and d) 3D coronal T2-weighted TrueFISP (True Fast Imaging with Steady state Precession) (TR, 6.04 msec; TE, 3.02 msec; interpolated voxel size, 0.3 × 0.3 × 0.6 mm3). The total scan time was less than 35 minutes over-all. The HemON trial MRI protocol consisted of a subset of these sequences (no 3D axial T2-weighted SPACE with ZOOMit as this trial does not include an implant at the lumbosacral level). In the HEMO trial, a similar protocol was performed on a Discovery MR750 (GE Healthcare, Waukesha, WI, USA) with an 8-channel cervical thoracic lumbar coil. The following two sequences were performed:a) 3D coronal T2-weighted FIESTA-C (Fast Imaging Employing Steady-state Acquisition) (TR, 7.65 msec; TE, 3.48 msec; voxel size, 0.4×0.4×0.4 mm3); and b) 3D sagittal T2-weighted CUBE sequence (TR, 2000 msec; TE, 90 msec; interpolated voxel size, 0.8×1.25×0.8 mm3). The total scan time was 84 min overall. Finally, in the HemON-NL trial, the standard protocol was translated to an Ingenia Omega (Philips Healthcare, Best, The Netherlands) with a 32-channel coil. The protocol comprised the following two sequences: a) sagittal T2-weighted FFE (Fast Field Echo) (TR, 3.34 msec; TE, 8.1 msec; interpolated voxel size, 0.75×0.75×3 mm3); and b) coronal Balanced FFE (TR, 7.0 msec; TE, 3.5 msec; inter-polated voxel size, 0.75×0.75×3 mm3). The total duration of the acquisition was approximately 20 minutes overall.
CT scans were acquired in a supine position with both arms on the sides. No contrast agents were administered. In the STIMO-HEMO and HemON clinical trials, CT of the trunk including the pelvis was performed on a 256-detector row CT system (Revolution Apex, GE Healthcare) using the following acquisition parameters: tube potential, 120 kVp; tube current, 115-200 mA, by enabling automatic tube current modulation; gantry revolution time, 0.35 sec; beam collimation, 128 × 0.625 mm; pitch, 0.508. CT images were reconstructed with a field of view of 150 × 150 and 400 × 400 mm2 for the spine and pelvis, respectively; at a section thickness/interval of 1.25/0.625 and 1.25/1.25 mm, respectively, with both smooth (Standard) and sharp (Bone plus) convolution kernels and a deep learning image reconstruction algorithm (TrueFidelity, GE Health-care; medium strength level) for image noise reduction. The nominal voxel size was 0.29 × 0.29 × 0.625 and 0.78 × 0.78 × 1.25 mm3 for the spine and pelvis, respectively. In the HEMO trial, a Revolution Discovery scanner (GE, Boston, MA, USA) was used to collect Standard and Bone Plus reconstructions (pitch 0.531:1, slice thickness 1.25 mm, interval: 1.25, 120 kV, 115-200 mA, rotation time: 0.4 s) of the thoracolumbar spine. In the HemON trial, an Ingenuity Core 128 scanner (Philips) was used to collect standard iDose and IMR reconstructions.
Personalized anatomical models of the spine We trained two different nnU-net networks45 on a labeled data set of 15 healthy participants from a previous study30. The first network was trained on high-resolution MRI to segment the cerebrospinal fluid, the white matter, and spinal root tissues. The second network was trained on low-resolution MRI data to segment vertebral bodies, intervertebral discs, and the spinal canal. We trained the networks on 12 participants and tested on 3 subjects. We achieved high dice scores for each tissue (cerebrospinal fluid: 95.50; white matter: 90.55; roots: 71.88; vertebrae: 89.90; intervertebral discs: 91.22). When applying the networks to a new data set, we corrected the segmentation manually if needed. The segmentation of the roots was used as guidance to trace each spinal root individually. The rostral location of each root was constrained within the white matter. For each spinal quadrant (dorsal/ventral left/right), the end-of-root within the white matter points was cut 1 mm within the white matter and a spline was made with them. These end-of-root lines were then guide-lines to force rootlets to attach to them, forcing them to end within the white matter. The procedure to generate the rootlets was described previously30. By forcing each rootlet to end in the end-of-roots guidelines, we made sure the rootlets followed the shape of the spine, thereby increasing robustness to scoliosis. We cut the white matter with respect to the rootlet divisions and obtained henceforth the spinal levels.
We used an anatomically informed framework trained on the VerSE datasets46 to segment the different vertebrae on CT imaging47. The authors of the framework showed a dice score of 91.04 on a large test set. Having both segmentations from MRI and CT spaces, we automatically placed labels on each vertebrae segmentation on MRI data. We then co-registered the CT data to the MRI space automatically using the Spinal Cord Toolbox48.
A computer-aided design (CAD) model of the 16 electrodes composing the paddle lead was positioned over the targeted dorsal root entry zones. Several visualizations of the 3D personalized anatomical model were then generated and shared with the clinical and neurosurgical team to guide preoperative surgical planning. Finally, a virtual paddle lead was placed over the selected location of the spinal cord and, and then projected onto a large screen in the operation room to guide the neurosurgeon in the insertion of the lead.
Neurosurgical intervention
The participants were put under protocolized propofol-based total intravenous general anesthesia and were placed in a prone position. Preoperative surgical planning informed the neurosurgeon about the vertebral entry-level and predicted optimal paddle lead position. Based on this knowledge, lateral and anteroposterior fluoroscopy x-rays were performed intra-operatively to identify the location of the laminotomies. A midline skin incision of approximately 5 cm on the back was performed, the fascia opened and the muscles retracted bilaterally. Excision of the midline ligamentous structures and a laminotomy at the desired entry level enabled the insertion of the paddle lead below the vertebra. For participants of the STIMO-HEMO and HEMO trials, a second skin incision or extended opening caudally was made and a second laminotomy was performed in the lumbar region based on the preoperative planning to allow for the insertion of the paddle lead over the lumbosacral spinal cord. The paddle lead(s) (Specify 5-6-5, Medtronic Inc, Minneapolis, MN, USA or ARCIM Thoracic Lead, ONWARD Medical N.V, Eindhoven, Netherlands) are inserted and placed over the midline of the dura-mater and advanced rostrally to the target position under repeated fluoroscopies. Electrophysiological recordings were acquired using standard neuro monitoring systems (Cascade Elite, Cadwell Industries, Kennewick, USA or ISIS Xpress, Inomed Medizintechnik, Emmendingen, Germany). Single-pulses of EES were delivered at 0.5 Hz and increasing amplitude to elicit muscle responses that were recorded from subdermal (Neuro-line Twisted Pair Subdermal, 12 x 0.4 mm, Ambu A/S, Ballerup, Denmark) or intramuscular needle electrodes (Inomed SDN electrodes, 40 x 0.45 mm, Inomed Medizintechnik, Emmendingen, Germany). The symmetry of the responses and the expected rostrocaudal distribution of muscle responses informed the corrections for lateral and rostrocaudal positioning. When the paddles deviated from a straight midline position, small additional laminotomies were made to remove bony protrusions and guide the paddle to a midline placement. Once the final position was determined, the paddle lead(s) was (were) anchored to the muscular fascia. In the STIMO-HEMO, HemON, and HemON-NL trials, the back opening was temporarily closed and the participants were put in lateral decubitus. Subsequently, the back incision was reopened and an abdominal incision of about 5 cm was made to implant each pulse generator in a subcutaneous pocket. In the HEMO trial, incisions of about 5 cm were made bilaterally in the upper buttocks region and subcutaneous pockets were created. The cables of the paddle lead were then tunneled between the back opening and subcutaneous pockets to be connected to the implantable pulse generator (Intellis, Medtronic Inc, Minneapolis, MN, USA or ARCIM IPG, ONWARD Medical N.V., Eindhoven, Netherlands). The implantable pulse generator(s) was (were) implanted in the subcutaneous pocket(s) and all incisions were finally closed. In the STIMO-HEMO trial, both paddle lead cables were tunneled from the same flank side and tunneling was performed between the 2 abdominal subcutaneous pockets.
Technological framework
Stimulation optimization
The configuration of the electrodes for EESXX programs was guided by the understanding of the mechanisms through which EES regulates blood pressure (Fig. 2a). This configuration involved three steps. First, an intraoperative mapping was conducted to identify the relevant rows of electrodes to target the hemodynamic hotspot1 and to elicit the largest pressor response. Second, postoperative imaging was used to update the anatomical models and thus estimate the location of the electrodes that maximized recruitment of the dorsal root entry zones projecting to the hemodynamic hotspot. Third, a single postoperative session was conducted to quantify pressor responses elicited by each row of electrodes. Both monopolar and multipolar configurations of electrodes were tested to ensure the maximal activation of the hemodynamic hotspot. Stimulation frequency was set at 120 Hz, since we established that this frequency maximized pressor responses1. Multiple electrode configurations that each targeted a specific region of the hemodynamic hotspot. Each electrode configuration was delivered with 2 ms delay between each pulse, forming a traveling wave over the spinal cord. In the STIMO-HEMO and HEMO trials, pulse-width was mapped between 100 μs - 500 μs. Five out of six participants had a pulse width of 300 μs, while one participant (P3) had a pulse width of 500 μs. For the HemON and HemON-NL trials, the pulse width was set at 300 μs. The amplitude was set based on the increase in the current delivered through each electrode configuration until the systolic pressure increased by 20 mmHg, the diastolic pressure increased by 10 mmHg, or the patient reported any discomfort such as muscle contractions, or unpleasant sensation such as strong tingling. To validate the combined effect of multiple electrode configurations, we compared the amplitude of the pressor responses when EES was delivered with a configuration targeting a single level of the hemodynamic hotspot, versus all the levels embedded in the hemodynamic hotspot. All other stimulation parameters were maintained unchanged.
Intraoperative hemodynamic monitoring
Intraoperative beat-to-beat hemodynamic data was obtained from a pressure transducer connected to an arterial line cannula placed in the distal radial artery. The arterial line system was zeroed to atmospheric pressure. Arterial cannulation was performed either by palpation or with ultra-sound guidance under sterile conditions. Adjustments in the anesthetic depth as well as administration of vasoactive medications or significant volumes of intravenous fluids were avoided during hotspot testing. All data was recorded using VitalRecorder software (VitalDB, Department of Anesthesiology and Pain Medicine, Seoul National University College of Medicine, Seoul, Korea) at a sampling frequency of 125 Hz. Alternatively, the arterial line was recorded with LabChart24 (ADInstruments, Dunedin, New Zealand).
Postoperative hemodynamic monitoring
Beat-to-beat blood pressure and heart rate were acquired using finger plethysmography (Finometer, Finapres Medical Systems; Amsterdam, Netherlands). Beat-by-beat blood pressure was calibrated to brachial artery blood pressure collected using an arm cuff embedded and synchronized with the Finometer49–53. Brachial arterial pressure was sampled at 200 Hz, while the systolic, diastolic, and mean arterial pressure were extracted from the calibrated arterial pressure at 1 Hz. The heart rate was also sampled at 1 Hz. Raw data and automatically extracted hemodynamic parameters were saved by and exported from the Finometer.
Off-label investigational system
The investigational system used in the STIMO-HEMO and HEMO clinical trials consisted of a set of CE-marked, FDA-approved medical devices used off-label. Two implantable pulse generators (lntellis™ with AdaptiveStim™, Medtronic Inc, Minneapolis, MN, USA) were connected to a paddle lead (Specify™ 5-6-5 SureScan™ MRI, Medtronic Inc, Minneapolis, MN, USA). These implants are indicated for chronic pain management. A tablet application with a wireless communicator device (Intellis™ clinician programmer, Medtronic Inc, Minneapolis, MN, USA) was used by the clinical team to set up the system and optimize the stimulation parameters. A remote control device and transcutaneous charger device (Patient Programmer and Recharger, Medtronic Inc, Minneapolis, MN, USA) were used by the participants to charge the implanted pulse generators and to turn EES programs on and off during their daily life, as well as adapt parameters of EES programs defined by the clinical team.
Purpose-built investigational system to treat hemodynamic instability
The investigational system used in the HemON and HemON-NL clinical trials were designed and built for treating hemodynamic instability based on the understanding of the mechanisms through which EES regulates blood pressure. A step-wise approach was followed. All participants were implanted with a purpose-built implantable pulse generator ARCIM (ONWARD Medical N.V., Eindhoven, Netherlands) that communicates with a purpose-built ecosystem of control devices. The first 4 participants of the HemON clinical trial were implanted with the off-label paddle lead (Specify™ 5-6-5 SureScan™ MRI, Medtronic Inc, Minneapolis, MN, USA) as the participants of STIMO-HEMO and HEMO clinical trials, whereas all other participants were implanted with a newly-designed novel, purpose-built paddle lead (ARCIM Thoracic Lead, ONWARD Medical N.V., Eindhoven, Netherlands).
Purpose-built implantable pulse generator and communication ecosystem to treat hemodynamic instability
The purpose-built ARCIM implantable pulse generator is a new 16-channel neurostimulation platform that was specifically developed to regulate neurological functions with EES. The ARCIM implantable pulse generator controls and delivers current-controlled stimulation pulses that can be predefined by the user, or controlled in real time. The system consists of a hermetically sealed, biocompatible can that surrounds the electrical components and a rechargeable battery that enables the delivery of electrical current. The external envelope is composed of two main components: the header containing the connector block that enables connection with 2 8-contact lead connectors and 2 coils for charging and communication, and the can with a rechargeable battery and electronics circuits. The ARCIM implantable pulse generator was developed according to all applicable standards for medical device development. Conventional biomedical technologies were used to fabricate the system and extensive bench and invivo testing was performed to verify its performance.
The implantable pulse generator was implanted subcutaneously in a pocket created in the abdomen. The ARCIM communication Hub ensured exchange of information with the implanted pulse generator using Near Field Magnetic Induction (NFMI), and also serves as a wireless charging station for the battery embedded in the implanted pulse generator. The communication hub is worn on a belt over or in proximity to the location of the implanted pulse generator. This Hub hosts a Bluetooth Low Energy (BLE) chip that enables fast, reliable wireless communication with external sensors or programmers. This versatility enables the acquisition of varying control signals and rapid exchange of information between the communication hub and the implanted pulse generator, which together establish the conditions for closed-loop control of stimulation parameters with latencies as low as 25ms between the generation of the stimulation command and the execution of this command. The Hub communicates with external programmers such as the ARCIM Clinician Programmer onto which an Android App designed for clinicians allows the configuration of EES programs and evaluation of the system. When EES programs are deemed safe for personal use, the Clinician Programmer provides the possibility to make this stimulation program available to the patient.
The patients, or their caregivers, can control the system through the ARCIM Personal Programmer. This Android Watch application allows users to select, start, and stop EES programs, as well as to tune stimulation amplitudes within predefined safety limits determined by the clinicians.
Device errors, paddle lead impedances, and daily stimulation utilization were extracted from usage logs across all devices.
Furthermore, the Clinician Programmer includes an application programming interface, termed the ARCIM API, that enables other programming softwares to control the stimulation, e.g. for closed-loop control of the stimulation.
All devices and softwares are adherent to the applicable standards and their performance was extensively tested. The entire system received the equivalent of an investigational device exemption from the competent Swiss authorities.
Purpose-built paddle lead
The ARCIM Thoracic Lead developed by ONWARD Medical is a new 16-electrode paddle lead that incorporates a configuration of electrodes optimized to achieve the selective recruitment of the dorsal root entry zones projecting to the last 3 segments of the lower thoracic spinal cord. The configuration of the electrodes was determined based on a digital library of human thoracic spinal cords that we elaborated from CT and MRI acquisitions (Supplementary Fig. 4). The average length and width of spinal segments were calculated from our library and existing literature30,34–39. The combination of these 2 datasets generated global averages and standard deviations for each level of the lower thoracic spinal cord, which determined the total length of the hemo-dynamic hotspot. We created a normal distribution with the calculated global means and standard deviations and found that a length of 64.2 mm accurately represents the combined length of the 3 targeted spinal segments for over 95% of the population. This new paddle lead was designed following all applicable standards for medical implant development and fabricated using conventional biomedical technologies. Extensive bench and in-vivo testing was conducted to validate the mechanical, electrical, and biocompatibility robustness of the paddled lead. The equivalent of an investigational device exemption was granted by the competent Swiss competent authorities.
System for closed-loop control of blood pressure
The purpose-built investigational system was embedded in a framework for closed-loop control of systolic blood pressure based on real time monitoring of hemodynamic parameters. For this purpose, changes in systolic blood pressure were measured with the Finometer (see Postoperative hemodynamic monitoring). This signal was acquired using an analog-to-digital converter (Delsys Trigno Avanti Analog, Delsys Inc., Natick, Massachusetts) that wirelessly streamed the data to a sensor base station (Delsys Trigno Research+, Delsys Inc., Natick, Massachusetts). The 1 Hz systolic blood pressure measured with the Finometer was up-sampled to 1,257 Hz by the Trigno system, which was connected through USB to a Windows computer onto which a medical-grade software developed by EPFL was running. This software30 enables closed-loop control of the ARCIM system through the ARCIM API. The specific closed-loop algorithm used for the control of systolic blood pressure was implemented in a multi-threaded Python application that interfaces with the GDrive+ software via TCP/UDP communication enabling the reception of systolic blood pressure signals and delivery of stimulation updates.
Algorithm for closed-loop control of blood pressure
The systolic blood pressure signal was smoothed with an exponential weighted moving average over a 3-second window. This elaborated signal was sent as an input to a Proportional-Integral-Derivative (PID) controller. To minimize steady-state errors and minimize system saturation, the integral term was penalized over time when the stimulation remained unchanged at the minimum or maximum amplitude of EES. Additionally, a forgetting factor on the integral term reduces error accumulation. A minimum (0 mA) and maximum amplitude (participant dependent) amplitude range was set to ensure the comfort and safety of the participant. The controller target was defined based on average blood pressure measurements from 10-minute tilt-table tests with EES (see Postoperative hemodynamic monitoring). A target range of +/-3 mmHg was set around the user-defined target to account for intrinsic blood pressure variability. If the change in EES amplitude was greater than 0.5 mA, the controller incrementally ramped in steps of 500 ms per 0.5 mA for additional comfort. To calibrate the closed-loop control of EES, we found a linear operating range of stimulation amplitudes and systolic blood pressure. The minimum was the amplitude at which the first increase in systolic blood pressure occurred, and the maximum was the last amplitude at which the systolic blood pressure stopped increasing, or when the patient reported abdominal muscle contractions. Each participant remained in a seated position as the amplitude of EES increased by 1 mA every 1 minute to find the minimum (0 mA) and maximum parameters of stimulation (participant dependent). A regression between the average systolic blood pressure per stimulation amplitude resulted in a linear relationship between the systolic blood pressure and the amplitude of stimulation. The closed-loop was tuned manually using the heuristic Zieger-Nichols method54. The initial proportional constant was the slope of the regression derived from the calibration phase. The proportional, integral, and derivative constants were tuned over one 70-degree tilt. The same parameters were then used across days to test dynamic changes in the severity of orthostatic challenge. The dynamic changes in tilts consisted of a pseudo-random sequence of 8 changes in tilt angle spanning mild (20 degrees) to severe orthostatic challenge (80 degrees). Each tilt position was held for 60 to 120 seconds. The same sequence was run with continuous EES to enable comparison. An inertial measurement unit-based control algorithm was also developed using the same framework. The angle of the tilt table was derived from changes in 3D acceleration of an inertial measurement unit placed onto the tilt table (Delsys Trigno Avanti, Delsys Inc, Natick, Massachusetts). Data from the inertial measurement unit was recorded at 148 Hz. Accelerometer values in three dimensions were used to compute angular changes. Ranges of tilt positions were mapped to supine, seated, and standing positions, and corresponding amplitudes of EES.
Computer simulations
Numerical simulations were conducted using the computational life sciences platform Sim4Life (ZMT Zurich MedTech AG, Switzerland). First, a canonical spinal cord model was developed based on average anatomical measurements of the participant spinal cord anatomies. It was then used to simulate EES exposure throughout the spinal cord anatomy, as well as induced electrophysiological responses30. The dosimetric simulations involved finite element discretization using Sim4Life, which resulted in approximately 5 million volume elements (voxels). Each voxel was attributed dielectric properties corresponding to the assigned tissue types55 (i.e., spinal cord, spinal roots, cerebrospinal fluid, epidural fat, and vertebral bone). Anisotropic conductivity maps were assigned to the white matter and rootlets, leveraging existing Sim4Life functionalities (either diffusion-simulation- or DTI-imaging-based). By precomputing a suitable basis of electric potential distributions (one per epidural stimulation electrode contact) for the canonical bioelectrical spinal model, the electrical fields resulting from arbitrary EES pulses can rapidly be obtained30. For that purpose, the ohmic-current-dominated electro-quasistatic approximation of Maxwell’s equations was solved56. In addition to providing an average anatomical geometry, the spine model also encompassed trajectories of axons entering the spinal cord through the dorsal roots and ascending in the dorsal columns. Nerve fiber trajectories were generated inside the rootlets using Sim4Life’s IM-Safe tool, while the fibers in the dorsal column were uniformly distributed. Sim4Life includes functionality to estimate fiber recruitment using the generalized, Green’s function-based activating function, or detailed electrophysiological modeling. For the present work, the classic activating function (proportional to the discrete second derivative of the extracellular potential56) was used to assess relative stimulability for a given pulse shape, while electrophysiological simulations with the NEURON simulator incorporated in Sim4Life57 permitted conversion to absolute thresholds. For the latter, the predefined diameter-dependent spatially extended non-linear node (SENN) fiber model underlying low-frequency exposure safety standards was used. By assigning fiber diameters according to histology-based probability distributions, the percentage of recruited fibers per root could be estimated and the recruitment level across the different spinal segments and within the dorsal column could be predicted. Activation levels were normalized relative to maximal activation (i.e., recruitment of all dorsal root fibers). Various methodologies have been developed to optimize selective stimulation, including neural network backpropagation and hybrid surrogate modeling (SUMO) - multi-objective genetic algorithm (MOGA) optimization.
Electromyographic recordings
Electromyographic (EMG) activity of selected muscles was acquired at a 1,259kHz sample rate using the 16-channel wireless Delsys Trigno sensors (Delsys Inc, Natick, MA, USA) with bipolar surface electrodes. Skin hairs were shaved and abrasive gel was used (Nuprep gel) on the area of interest.
Trunk kinetics
Kinematics were captured by Physilog 5 sensors (Mindmaze, Lausanne, Switzerland) that were placed on the vertebral bodies of C3, T1, T6, T10, L1, and S1 of the participants. Each sensor is a stand-alone 10 degree of freedom Micro-Electro-Mechanical System (MEMS) inertial measurement unit includes 3D accelerometer, 3D gyroscope, 3D magnetometer, and a barometric pressure sensor at a sampling frequency of 256 Hz. These inertial measurement units were synchronized via a proprietary radiofrequency protocol.
Clinical evaluations
Orthostatic challenge with tilt-table test
Participants were transferred to a supine position on a table capable of head-up tilt. We applied restraint straps to secure the patient on the knees, on the hips, and on the chest, with the feet stabilized. Resting supine blood pressure was recorded continuously for approximately 5-10 minutes to establish baseline values. Following this period, we tilted the participants upright up to a maximum of 70 degrees while recording hemodynamic values and symptoms of orthostatic tolerance. The time to reach the desired tilt angle from supine was achieved in less than 45 seconds. Participants were tilted until reaching their tolerance threshold or for a maximum duration of 10 minutes. They were asked not to talk during the test except to inform of and grade symptoms. The participant was asked to report any symptoms every 1-3 minutes. The participant was asked to rank the intensity of their symptoms between 1-10, 1 being no symptoms at all, and 10 being maximum feelings of dizziness, lightheadedness58, or nausea19,58. The participants were instructed to notify the research team if they needed to be returned to the supine position.
Blood draws and circulating norepinephrine
Blood samples were withdrawn at baseline after a period of rest in the supine position for 15 minutes and at 2 minutes (T2) after the onset of the tilt. During baseline testing, a total of 3 samples of blood were with-drawn. When testing with EES turned on, an additional sample of blood was withdrawn 5 minutes after the onset of EES in supine position. Blood samples allowed quantification of plasma cate-cholamines, including plasma norepinephrine, plasma epinephrine, and plasma dopamine. To avoid false positives, the participants were asked to abstain from any beta-blocker, dihydropyridine calcium channel blocker, phenoxybenzamine, anxiolytics, or decongestant medication. To ensure sample quality, prior to and immediately after, collection tubes were kept cold via immersion in crushed ice and kept in foil to reduce exposure to light. Samples were analyzed by the local Hospital Lab services at the University Hospital of Lausanne (Switzerland) and the Foothills Medical Center (Canada). Neuropathic pain prevented blood sampling in P8. Data and tests were not available at the time for P13 and were not part of the clinical protocol for P14.
Cerebral Blood Flow
Blood flow was measured in the seated position using Doppler ultrasound imaging equipment (Epiq Elite, Philips, United States). Measurements were conducted both with-out and with EES. The left and right proximal internal carotid arteries (ICA) and vertebral arteries (VA) were assessed using a 12-3 MHz broadband linear array probe (L12-3, Philips). Additionally, the right and left middle cerebral arteries (MCA) were examined using a transcranial 5-1MHz-pulsed TCD probe.
Postprandial blood pressure measurements
On 2 separate days at the same study time point, participants received an identical meal with equal caloric content and their blood pressure was monitored via an arm cuff prior to, during, and for 30 minutes to 1 hour following food intake. This assessment was performed 1 day without EES and 1 day with EESBP. For the assessment with EESBP, the stimulation was turned on during and after the meal. All points before starting the meal were averaged to determine a baseline value. Available cuff measures between 30 minutes and 60 minutes were averaged to quantify postprandial measures. The delta between baseline and postprandial averages was reported for the conditions without and with EESBP. Only participants who reported experiencing post-prandial hypotension were considered in the analysis (P3, P4, P5, P7, P10).
Questionnaires
Patient-Reported-Outcome-Measure questionnaires were sent to study participants via the electronic data capture system. Baseline questionnaires were sent once before surgery and then repeated every month after surgery. The following questionnaires are reported: Autonomic Dysreflexia following SCI questionnaire (ADFSCI), Study Sleep Scale (MOS-S) 53, Neurogenic Bowel Dysfunction Score (NBDS), and Usability SUS41.
Spirometer measurements
The study participants were asked to perform a systematic set of breathing and cough tasks. Volume, flow, pressure, and time data were recorded using the Spirometer (Spirobank II Smart Spirometer). During the coughing maneuver, EESTRUNK was precisely timed with the expiratory phase after releasing the glottis blockage. All participants from the HemON study performed the respiratory and coughing assessments at 3 and 6 months without stimulation and with EESTRUNK.
Canadian Occupational Performance Measure (COPM)
The COPM was employed to assess the functional performance and perceived satisfaction with activities of daily living of each participant. The COPM was administered by trained occupational therapists to follow standard procedures. At baseline, participants identified up to five activities they perceived as important but difficult to perform and rated their perceived performance and satisfaction on a 10-point scale. The general scores were calculated by simply averaging the performance and satisfaction ratings across activities. The same activities were then rated again by the participants 3 and 6 months after implantation.
Spasticity assessment
Participants’ upper and lower limb spasticity was assessed using the Modified Ashworth Scale (MAS) in 2 conditions: without and with EESSPASTICITY. Participants were either asked to use their standard-of-care medication or to not use their standard-of-care medication the morning prior to the assessment. Participants were asked to maintain the stimulation off in the mornings prior to the assessment, which were conducted in the late morning.
Rehabilitation outcomes
In the HemON trial, physiotherapists and occupational therapists kept track of the 1-month rehabilitation participants’ sessions and documented whether the different activities were performed without or with EES through a standardized form (Mapping of Rehab Training - MART).
Muscle sympathetic nerve activity
In the HEMO trial, efferent postganglionic Muscle sympathetic nerve activity (MSNA) was acquired from the left common peroneal (fibular) nerve using microneurography1,59,60. The location of the common fibular nerve was located using palpation. A 2-MΩ tungsten microelectrode (Frederick Haer and Co, Bowdoin) was then inserted percutaneously alongside a subdermal low impedance reference electrode 2cm adjacent. The amplifier and head stage (NeuroAmpEx, AD Instruments) were also grounded to the participant with a surface electrode on the patella. The recording electrode was adjusted until entering the nerve. Crossing into the nerve and into different nerve bundles was identified audibly by the microneurographer. Once in the nerve, the recording microelectrode was confirmed to be near nerve fibers directed towards skeletal muscle by auditory feedback during tapping/palpation of the tibialis anterior/peroneal muscles, which has been shown previously to be indicative of microelectrode proximity to efferent postganglionic muscle sympathetic nerves1. If auditory feedback was present during light stroking of skin on the dorsal foot/lower leg, this indicated proximity to post-ganglionic skin sympathetic fibers. The search was stopped when the electrode was in proximity to muscle sympathetic fibers.
Data processing
Intraoperative blood pressure data
During surgery (see Neurosurgical intervention), changes in blood pressure were recorded in response to different locations of EES using the arterial line (see Post-operative hemodynamic monitoring). Change in blood pressure or heart rate was defined as the difference in the average of 30-second windows before the start of the stimulation and at the end of the stimulation. For cross-subject analysis, the lower thoracic and lumbosacral stimulation locations in each subject were calculated by their relative distance to the center of the hemodynamic hotspot (T11) normalized by their respective spinal length (measured from T9 to the conus). Due to surgery time limits, P3 blood pressure mapping was not performed.
Postoperative blood pressure data
During a tilt-table test (See Orthostatic challenge), changes in blood pressure were recorded without stimulation or in response to different types of stimulation (continuous or closed-loop stimulation) using the Finometer (see Post-operative hemodynamic monitoring). Change in blood pressure or heart rate was defined as the difference in the average of a 60-second window before the start of the tilt and a 20-second window at 3 minutes of the challenge. If the participant could not tolerate at least 3 minutes of the test due to low blood pressure or other symptoms, an average of a 20-second window before the end of the tilt was used. All measurements in seated position were measured with stimulation on for 3 – 5 minutes. Change in blood pressure, or heart rate, was defined as the difference in the average of a 20-second window before the start of EES and the average of a 20-second window at 3 minutes, prior to stop stimulation. All signals were smoothed over a 10-second window for illustration. The same processing was used for post-operative, day 1 quantification. Data from continuous blood pressure monitoring was excluded if technical issues with the Finometer were observed including loss of the plethysmographic signal, incorrect calibration as measured by comparing continuous measuring to standard measurements from a brachial cuff, spasms, or due to study protocol deviations or adverse events.
Stimulation paradigms
The efficacy of continuous, closed-loop, and inertial measurement unit-based control of EES was measured by: 1) target error, defined as the error of systolic blood pressure with respect to a user-defined target on the systolic blood pressure, 2) error variability, defined as the standard deviation of the target error, to measure stability of the stimulation during dynamic orthostatic challenges. Both metrics were calculated from 1 minute before the dynamic challenges to 1 minute after the end of the challenge.
Cerebral Blood Flow (CBF) Calculations
Peak systolic and end-diastolic flow velocities were obtained from the Doppler spectrum analysis over three consecutive cycles in all arteries. Mean velocity was derived through the integration of the Doppler curve for each corresponding pair, with the mean value calculated across the three cycles. Arterial diameters were determined in B-mode for the left and right ICAs and VAs. The flow rate (Q) in each vessel, measured in milliliters per minute, was calculated using Equation 1, assuming a parabolic velocity profile. The total experimental cerebral blood flow (CBF_exp) was determined by summing the flow of the left and right ICAs and VAs, according to Equation 261.
Autonomic Dysfunction Following Spinal Cord Injury (ADFSCI)
The ADFSCI is a 24-item questionnaire divided into four sections: demographics, medication, Autonomic Dysreflexia (AD), and hypotension. The hypotension section consists of 7 items. Each item employs a 5-point scale to measure the frequency and severity of symptoms related to hypotension, including headaches, goose-bumps, dizziness, and light-headedness, across different situational contexts. The hypotension score corresponds to the sum of all the hypotension items. Participants were categorized as experiencing symptoms if the item score was higher or equal to 2, else they were categorized as not experiencing symptoms.
Analysis of kinematics
The acceleration and angular velocity collected by the Physilog 5 sensors were passed through a datafusion algorithm62 to extract sensor orientation with respect to the Earth in quaternion representation. Roll, pitch, and yaw angles were then calculated by using one of two representations (zyx or zyz), depending on the sensor placement with respect to its axis with the greatest motion, such that there would be no discontinuities in the pitch. The ratio of change in back curvature when turning the stimulation on is approximated as the ratio of the angles between the T6 and T10 vertebrae in the sagittal plane with and without stimulation.
Muscle sympathetic nerve activity
The MSNA signal was sampled at 20kHz, amplified (20,000×), and band-pass filtered (0.3–2.0 kHz) to obtain muscle sympathetic spike activity, and then rectified and integrated (0.1-s time constant) to obtain a multi-unit (mean voltage) neurogram (Labchart version 8, AD Instruments). After acquiring a stable recording site, a 3-min baseline commenced followed by step changes in epidural stimulation intensity for 1-2 minutes at each stage. Muscle sympathetic action potential spikes were identified as waveforms that matched a triphasic morphology with the main phase being negative. The negative deflections were only assessed if they were larger than 4.5x the standard deviation of the neurogram assessed during baseline (stimulator off).
Splanchnic impedance methods
In the HEMO trial, Bioelectrical Impedance (BEI) of the lower torso (abdominal cavity) was used to estimate fluid shifts within the splanchnic vascular regions during tilt-table testing63. BEI is measured by driving a small electrical current (0.7 mA with a frequency of 37 kHz) between electrodes and measuring the impedance with voltage-sensing electrodes on each end of the body segment of interest (Model 2994D/4 THRIM, UFI). For the lower torso, surface electrodes were placed at approximately the left 6th rib and around the left inguinal ligament. The impedance, measured in ohms, was low pass filtered at 0.5 Hz. Since body segment volume is inversely proportional to electrical resistance (V ∝1/R) such that an increase in resistance reflects a loss of body fluid, the impedance values were used to estimate body segment volume using the following equation: , where L is the segment length in meters, reff is the effective resistance 1.0 Wm, and R is the segment impedance in W.
Statistics
All data are reported as mean values and individual data points. No statistical methods were used to predetermine sample sizes, but our sample sizes are similar to those reported in previous publications64. All statistical analysis was performed in R using the base package ‘stats’, with primary implementation through the ‘tidyverse’ and ‘broom’ packages. Tests used included one or two-tailed paired or independent samples Student’s t-tests, one-way ANOVA for evaluations with more than two conditions, and one- or two-way repeated-measures ANOVA for assessments, when data were distributed normally, tested using a Shapiro–Wilk test. Post hoc Tukey tests were applied when appropriate. For regressions, mixed model linear regression was used in cases of multiple observations, or else standard linear modeling. In cases where group size was equal to or less than three null hypothesis testing was not completed. The significance level was set as P < 0.05. Exclusions of data are noted in the relevant methods sections. Unless stated otherwise, experiments were not randomized, and the investigators were not blinded to allocation during experiments and outcome assessment.
Funding
This work was financially supported by: Eurostars project E!113969 PREP2GO, DARPA subaward A21-0795-S001 (P.O. 1083297), Eurostars E!1748 IMPULSE, Medtronic (ERP-2020-12543), Onward Medical, PHRT-279, the Defitech foundation, CIHR, Cumming Medical Research Fund, Edith Rodie Estate, Hotchkiss Brain Institute, Libin Cardiovascular Institute.
Author contributions
Writing of the manuscript
G.C. and A.Phillips., wrote the manuscript with J.W.S., L.Asboth., J.Bloch., N.H., R.D., A.P.G., and S.H., and all the authors contributed to editing the manuscript.
Conducting therapy sessions or clinical evaluations
A.P.G., A.I., K.K., A. Phillips, R.M., J. Lee, J. Rimok, L. Asboth, M.B., J.B, L.B., K.D., N.H., N.I., J. Ledoux, E. Pralong, C.V., G.W., G.C., L.D.H., R.D.” S.H., N.M., N.L., A.N., F.A., M.D., K.G., C.Picq, J.Pradarelli, M.Rieger, R.V., P.D., D.S, E.A., L.Aviolat, O.B., T.B., F.B., M.B., J.Bloch, J.Brancato, K.B., S. Carda, K.D., N.Hermann, J.H., N.I., J.Ledoux, V.L., J.Miroz, A.Paley, J.Pillonel, M.Ramirez, M.V., X.V., L.W., G.W., J.W.S., P.K. and E.B.
Data analysis
A.P.G., A.I., K.K., R.M., P.D., D.S., J. Lee, J. Rimok, A. Phillips, L..Asboth, T.B., F.B., N.H., J.Ledoux, T.M., L.D.H., R.D., G.D., S.H., N.M., J.W.S., A.N., F.A., M.D., K.G., A.K., S.Morand, F.S., R.V. and E.B.
Hardware/software development
A.P.G., L.Asboth, N.H., C.H., T.M., C.S., G.C., R.D., G.D., S.H., N.M., A.R., J. Garcia, N. Kuster, E.N., T.N., M.C., M.D., V.D., J.D., L.D., Y.D., D.G., K.G., H.L., J.Murphy, E.Paoles, J.F., T.V., E.M., S.Mandija and C.B.
Regulatory affairs management
K.K., A.Phillips, L.Asboth, J.Bloch, L.B., N.H., G.C., L.D.H., R.D., A.P.G., I.V., P.B., M.D., L.D., K.G., C.J., H.L., C.Picq, M.Rieger, E.R., R.V. and A.W.
Illustration preparation
L.Asboth, N.H., G.C., R.D., A. Phillips, K.K., A.P.G., M.G., S.H., F.M., J.Ravier and J.W.S.
Funding acquisition and supervision
G.C., J.B. and A.Phillips.
Team Supervision
K.K., A.Phillips, L.Asboth, J.Bloch, G.E., L.H., T.M., P.S., K.V., G.W., G.C., R.D., J.W.S., N. Kuster, E.N., N.Keijsers, N.L., I.V., M.D., H.L., J.Murphy, E.R., A.W., S.Mandija and C.B.
Surgery
J.Bloch, S.Casha, F.G., E.K. and J.W.S.
Data availability
Software routines developed for the data analysis will be made available upon reasonable request to the corresponding authors.
Code availability
Codes are available upon reasonable requests to G.C. (gregoire.courtine{at}epfl.ch)
Competing interests
The authors declare competing financial interests: G.C., A.Phillips., J.W.S, J.Bloch., R.D., L.Asboth, T.M. and A.R. hold various patents in relation with the present work. G.C., A.Philips., J.Bloch, V.D., J.Murphy and H.L. are minority shareholders of ONWARD medical, a company with direct relationships with the presented work. G.C., J.Bloch., A.Phillips. and R.D. are consultants of ONWARD medical. P.B., M.D., L.D., K.G., C.J., H.L., C.Picq, M.Rieger, E.R., R.V., A.W., F.A., S.Morand, J.Pradarelli, E.Paoles, J.Murphy, Y.D., M.C., D.G., A.K, J.D., F.S. and V.D. are employees of ONWARD medical. N.Kuster and E.N. are shareholders of ZMT Zurich MedTech AG, which produces the Sim4Life software. The other authors declare no competing interests.
Additional information
Correspondence and requests for materials should be addressed to G.C., J.Bloch and A.Phillips.
Acknowledgments
The authors would like to thank the study participants for their commitment and trust. All participants gave their informed consent for publication of their images. The authors would also like to thank M. Boulingre, J. Bourgeot, K. Charlong, M. Kameni-Hekam, C. Magrini, L. Mathews, L. McCracken, C. Parcheminey, M. Pedone, C. Perrin, M. Picek, K. Plesniar, C. Schiff, M. Smith and K. Tourigny for various contributions including data collection, support in patient logistics, patient care and regulatory affairs; I. Benrikane, A. Di Russo, and Y. Qin for their contributions in data analysis.
Footnotes
↵# These authors jointly supervised this work
↵& Clinical Study Team Lorraine Aviolat7, Nadia Bérard23, Julie Brancato23, Krystel Bruyère7, Rebecca Charbonneau1,2, Chris Drummond-Main2, Sean Dukelow1,2, Grégoire Eberle7, John Gaudet25, Natacha Herrmann7, Julie Hervé5,6,7, Jamie Johnston2, Valentin Lupi7, John-Paul Miroz5,6,7, Aurélie Paley5,6,7, Julie Pillonel7, Carole Poulin-Kella7, Mélanie Ramirez7, Katrien Van Den Keybus7, Camille Varescon5,6,7, Molywan Vat5,6,7, Xavier Vo Pham7, Laurence Wenger7 and Patrick Koomen9
↵&& Onward Team Francesco Acquati6,11, Pierre Bessot11, Julien Dedelley11, Vincent Delattre11, Anahita Kyani11, Hendrik Lambert11, Sebastien Morand11, Chloé Picq11, Jared Pradarelli11, Francesca Stradolini11 and Rosanne van Dijsseldonk11
Missing Table 1 and 2 added.