Abstract
Dysfunction and diseases of the gastrointestinal (GI) tract are a major driver of medical care. The vagus nerve innervates and controls multiple organs of the GI tract and vagus nerve stimulation (VNS) could provide a means for affecting GI function and treating disease. However, the vagus nerve also innervates many other organs throughout the body, and off-target effects of VNS could cause major side effects such as changes in blood pressure. In this study, we aimed to achieve selective stimulation of populations of vagal afferents using a multi-contact cuff electrode wrapped around the abdominal trunks of the vagus nerve. Four-contact nerve cuff electrodes were implanted around the dorsal (N = 3) or ventral (N = 3) abdominal vagus nerve in six ferrets, and the response to stimulation was measured via a 32-channel microelectrode array (MEA) inserted into the left or right nodose ganglion. Selectivity was characterized by the ability to evoke responses in MEA channels through one bipolar pair of cuff contacts but not through the other bipolar pair. We demonstrated that it was possible to selectively activate subpopulations of vagal neurons using abdominal VNS. Additionally, we quantified the conduction velocity of evoked responses to determine what types of nerve fibers (i.e., Aδ vs. C) responded to stimulation. We also quantified the spatial organization of evoked responses in the nodose MEA to determine if there is somatotopic organization of the neurons in that ganglion. Finally, we demonstrated in a separate set of three ferrets that stimulation of the abdominal vagus via a four-contact cuff could selectively alter gastric myoelectric activity, suggesting that abdominal VNS can potentially be used to control GI function.
Similar content being viewed by others
Introduction
The abdominal vagus nerve is involved in the control of gastrointestinal (GI) tract motility and inflammation, as well as pancreatic function, and consequently plays critical roles in GI disease, obesity, and diabetes1. The mechanisms through which vagal pathways mediate these effects, however, have been difficult to determine due to complex organ innervation, and also because the abdominal vagus nerve contains a variety of fiber types, of which approximately 80% are afferent fibers and 20% are efferent fibers2. Strategies for understanding vagal function have relied primarily on pharmaceutical and ablation approaches but with limited success3,4. Vagus nerve stimulation (VNS) with electrical current has been used to interrogate vagal function, as well as for therapeutic applications, including control of epileptic seizures and obesity5,6. However, the therapeutic effects of existing VNS approaches have been modest, with significant off-target effects, including hoarseness, cough, dyspnea, pain, paresthesia, nausea, and headache7. Importantly, available clinical VNS systems use electrodes with contacts that fully wrap around the nerve and simultaneously engage many functional pathways, thus providing little insight into mechanisms of action of nerve stimulation and limiting the potential to tune VNS to limit side-effects8.
Many studies have demonstrated that multi-contact nerve cuff electrodes, with contacts spaced around the circumference of the nerve, can achieve a selective interface with peripheral nerves, allowing for targeting of specific functions while avoiding off-target effects9. For example, studies have shown that selective stimulation of the median, radial, and ulnar nerves can activate particular muscles of the hand and arm10,11. Often, selectivity of stimulation of efferent neurons is measured by recording evoked EMG activity in skeletal muscles11,12 or by measuring joint torques in response to muscle contractions13,14. For studies focused on stimulation of sensory afferents, selectivity can be measured by psychophysical measurement of perception of those sensations15, EMG recording of reflexive muscle responses to stimulation16, or by recording compound action potentials via nerve cuff recordings17. Measuring the selectivity of VNS is challenging because the physiological responses to stimulation are often not well understood or may only occur following minutes or hours of stimulation; for example, changes in immune responses or gastrointestinal rhythms18. Techniques that quantify selectivity by measuring effects on end-organ function or by recording evoked responses in nerve branches are also highly challenging because the vagus nerve branches are small, variable, and often inaccessible; they also innervate organs throughout the abdomen19.
The nodose ganglion is a compelling location to measure the selectivity of VNS, particularly for stimulation of the abdominal vagus nerve. The nodose ganglia are enlargements of the vagus nerve containing the cell bodies of autonomic sensory neurons projecting from organs throughout the body including the heart, lungs, and alimentary tract20,21. Action potentials can be recorded from these cell bodies using microelectrode arrays (MEA)22.
In this study, we sought to assess if specific populations of axons in the abdominal vagus nerve could be selectively activated using VNS. Experiments were conducted in ferrets, which have several advantages as a model for assessing abdominal vagal function because of similarity to humans in gastric anatomy, vagal regulation of gastric motility, and emetic responses23. In six animals, a four-contact cuff electrode was wrapped around the abdominal vagus nerve and compound action potential (CAP) signals elicited by stimulation through different contact pairs were recorded using a 32-channel MEA inserted in the nodose ganglion. We varied stimulation parameters (pulse amplitude and pulse width) to ascertain the maximal number of MEA channels that responded to stimulation through one cuff electrode pair or the other. We also measured the conduction velocity of the nerve fibers that were activated by stimulation, and quantified the somatotopic organization of the responses in the nodose ganglia through a nearest-neighbor analysis. Finally, we demonstrated that our abdominal VNS approach can selectively drive changes in gastric myoelectric activity, demonstrating the promise of this technique for treating GI diseases.
Results
Overview
The primary goal of this study was to determine whether subpopulations of axons in the abdominal vagus nerve could be selectively activated as a means to control GI function while avoiding off-target side effects. In six ferrets, a cuff electrode with two bipolar pairs of contacts (Fig. 1) was wrapped around the dorsal (N = 3) or ventral (N = 3) abdominal vagus nerve, and the response to stimulation was measured by recording evoked CAPs through a 32-channel MEA implanted into the nodose ganglion (left for ventral vagus nerve and right for dorsal vagus nerve, as these are the predominant pathways for neurons innervating the ventral and dorsal regions of the stomach, respectively)24. Stimulation amplitude (0–3000 µA) and pulse width (0.1, 0.5, and 1.0 ms) were varied using a binary search algorithm to quantify threshold (i.e., the minimum stimulation amplitude that evoked a response in the nodose for each pulse width) and their effect on selectivity of VNS. For each unique amplitude/pulse width combination, a train of 120 pulses was delivered through a bipolar cuff contact pair at 2 Hz. A stimulus-triggered average of the response recorded from each MEA channel was calculated to reduce noise and facilitate detection of CAPs (Fig. 2), which often had a peak-to-peak amplitude of less than 10 µV. For each pulse width, the binary search varied stimulation amplitude values with a 20 μA target resolution. If a response was detected on any MEA channel at a stimulation amplitude, amplitude would be decreased by half the difference from the next lowest amplitude tested, and if no response was detected, amplitude would be increased by half the distance from the next highest amplitude tested. Selectivity was determined based on the number of MEA channels that recorded a CAP in response to abdominal VNS through only one cuff contact pair. Further, the somatotopic organization of the responses in the nodose ganglia was determined. Finally, we demonstrated that low amplitude stimulation through different cuff contact pairs can evoke unique effects on gastric myoelectric activity.
Abdominal VNS evokes CAP responses in the nodose ganglion
In all six animals, abdominal VNS through each cuff pair evoked CAP responses in multiple MEA channels. For pulse widths of 0.1, 0.5, and 1.0 ms, Table 1 shows the threshold stimulation amplitude that evoked a response in at least one MEA channel. It is important to note that we did not control the rotational orientation of the cuff on the vagus nerve, so the designation of cuff contact pairs 1:2 and 3:4 is arbitrary. All animals exhibited responses to stimulation with a pulse width of 0.1 ms, except one cuff pair in animal F22-19 and both cuff pairs in animal F25-19. In both of those animals, stimulation with a pulse width of 0.4 ms evoked responses in multiple MEA channels. Mean thresholds were 1027, 307, and 248 µA for pulse widths of 0.1, 0.5, and 1.0 ms, respectively, in accordance with the expected strength-duration relationship in which threshold amplitude decreases exponentially as pulse width increases25.
Selectivity of abdominal VNS
To quantify selectivity of abdominal VNS, we characterized recruitment within the MEA (Fig. 3) as stimulation amplitude increased from threshold. We also calculated a selectivity index (SI), which, when maximized, selects stimulation parameters for each pair of cuff contacts that maximize the number of MEA channels with evoked responses while minimizing the number of overlapping MEA channels that respond to stimulation through both cuff contact pairs. Table 2 shows the stimulation amplitudes that maximized SI for each cuff contact pair at each pulse width in all animals. For all animals (Fig. 4), as amplitude or pulse width increased, stimulation evoked responses on multiple additional MEA channels across the entire array. Maximizing the selectivity index while limiting overlap to no more than 3 MEA channels (i.e. 10% of 32 total channels) demonstrated that stimulation through pairs of cuff contacts could selectively activate an average of 2.5, 8.25, and 7 MEA channels with pulse widths of 0.1, 0.5, and 1 ms, respectively. It is important to note, that even with 1 ms pulses at 3 mA, stimulation through an individual electrode typically evoked responses in ~ 20 MEA channels (gray bars in Fig. 4), with some MEA channels never detecting any CAP response to stimulation.
Conduction velocities of VNS evoked responses
To determine the types of nerve fibers (e.g. Aδ vs C) activated by abdominal VNS, the conduction velocities of the evoked CAPs in the nodose ganglion were calculated (Fig. 5). At both threshold and the maximum stimulation amplitude tested (i.e. 3 mA), the majority of evoked responses (70.1% and 91.7%, respectively) had conduction velocities between 0 and 3 m/s (i.e. C fibers)26. At threshold, evoked responses recorded by an individual MEA channel (Figure S1) were almost always either exclusively C (0.1 ms: 25%; 0.5 ms: 72.37%; 1.0 ms: 72.97%) or Aδ (0.1 ms: 68.75%; 0.5 ms: 16.36%; 1.0 ms: 21.62%), with only a small portion of electrodes (0.1 ms: 6.25%; 0.5 ms: 10.91%; 1.0 ms: 5.41%) recording both types of responses on the same MEA channel at threshold. While we did not quantify the receptive fields or modalities of these neurons, their conduction velocities suggest C-fibers were primarily recruited, along with a smaller population of Aδ fibers.
Patterns of recruitment across the nodose ganglion
To understand if there is somatotopic organization of the nodose ganglion with respect to the abdominal vagus nerve, we performed a nearest neighbor analysis for the MEA channels that responded to stimulation for each animal. Trials where at least one of the cuff pairs elicited no selective responses at the stimulation amplitude that maximized SI were excluded from this analysis. MEA channels that recorded non-selective responses from both cuff pairs were excluded from calculations for the opposing bipolar cuff pair because the distance to the nearest responding channel from the opposing bipolar cuff pair would always be zero but the distance to the nearest responding channel from the same bipolar cuff pair would always be greater than zero for these electrodes, resulting in a non-informative skew in the result. Histograms of this analysis across 6 ferrets are shown in Fig. 6a,c,e.
To further understand how recruitment of neurons at the nodose ganglion changed as a function of stimulation amplitude, we calculated the centroid location of all responding MEA channels for each bipolar cuff pair. We calculated these centroids at the stimulation amplitude that maximized SI and at the maximum stimulation amplitude tested (i.e., 3 mA). We then calculated the distance between the centroid of each cuff pair and compared these distances to determine if there was a change in the overall pattern of activation across the MEA as a function of stimulation amplitude (Fig. 6b,d,f). For all animals and pulse widths, a Wilcoxon signed rank test shows that there was no statistical difference between conditions (p = 0.25, 0.125, and 0.0625 for pulse widths of 0.1, 0.5, and 1.0 ms, respectively).
Functional effects of selective abdominal VNS
In a separate set of experiments in three additional ferrets, we examined the effects of abdominal VNS on GI myoelectric activity. The goal of these experiments was to determine if low amplitude VNS through bipolar pairs of cuff contacts on the abdominal vagus nerve could selectively drive functional changes in the stomach. For these experiments, a four-contact cuff was implanted around the ventral abdominal vagus nerve and four four-contact planar electrodes were sutured to the serosal surface of the stomach to record GI myoelectric activity. Across the three ferrets, 10 GI myoelectric signals displayed a statistically significant dominant frequency peak at baseline (9.62 ± 0.63 counts per minute; cpm). In all ferrets, VNS through at least one of the bipolar cuff pairs selectively evoked a distinct effect on normal GI myoelectric activity. VNS was delivered at a frequency of 15 Hz and a pulse width of 0.1 ms in all ferrets. For Ferret 15–18, when stimulation amplitude was 0.1 mA, there was no change in the GI myoelectric activity for either bipolar pair. At 0.2 mA, stimulation on one pair of bipolar contacts (1–2) did not result in any change in the GI myoelectric activity, while stimulation on the other pair of bipolar contacts (3–4) resulted in a decrease in the power in the normogastric range (8–11 cpm) to nearly 0 (Fig. 7a). This change in GI myoelectric activity was not accompanied by any overt behavioral response (retching or emesis). The waterfall plots in Fig. 7b display the change in the power spectral density of the signal in response to stimulation on bipolar contacts 1–2 and 3–4. For stimulation on contacts 3–4 we observed an immediate suppression of GI myoelectric activity upon stimulation onset.
For Ferret 13–18, there was no change in the GI myoelectric activity for either bipolar pair when the stimulation amplitude was below 1 mA. At 1 mA, stimulation on one pair of bipolar contacts (1–2) resulted in a decrease in the power in the normogastric range (8–11 cpm) to 20% of the total power in the recorded signal (Fig. 7e). Approximately 1.5 min after the onset of stimulation we observed retching followed by emesis. The waterfall plot in Fig. 7f displays the change in the power spectral density of the signal leading up to the first retch, in response to stimulation on bipolar contacts 1–2. Stimulation on the other pair of bipolar contacts (3–4) did not result in any change in the GI myoelectric activity and was not accompanied by any overt behavioral response.
Similarly, for Ferret 14–18, there was no change in the GI myoelectric activity for either bipolar pair when the stimulation amplitude was below 1 mA. At 1 mA, stimulation on one pair of bipolar contacts (3–4) resulted in a decrease in the power in the normogastric range (8–11 cpm) to 20% of the total power in the recorded signal (Fig. 7i). This decrease in power was accompanied by retching and emesis approximately 1.1 min after the onset of stimulation. The waterfall plot in Fig. 7j displays the change in the power spectral density of the signal leading up to the first retch, in response to stimulation on bipolar contacts 3–4. Stimulation at 1 mA on the other pair of bipolar contacts (1–2) resulted in a decrease in the power in the normogastric range (8–11 cpm) to 38% of the total power in the recorded signal. However, this was not accompanied by any overt behavioral response (retching or emesis). For ferrets 15–18 and 13–18, the decrease in power in the normogastric range, roughly coincided with the onset of stimulation whereas for Ferret 14–18, the suppression of GI myoelectric activity occurred 20–30 s after the onset of stimulation. For all three animals, stimulation did not drive a change in breathing rate (Fig. 7c,g,k), heart rate (Fig. 7d,h,l), body temperature, or systolic blood pressure (Figure S2) at amplitudes that drove changes in GI myoelectric activity.
Discussion
The goal of this study was to selectively stimulate sub-populations of fibers within the abdominal vagus nerve to drive changes in GI activity that could be used for treating diseases such as gastroparesis and obesity. In six ferrets, we demonstrated that a multi-contact cuff electrode wrapped around the abdominal vagus nerve can selectively activate independent populations of nerve fibers, and in three other ferrets, we showed that selective VNS can drive changes in GI myoelectric activity. We assessed the selectivity of stimulation by recording evoked CAPs through a MEA inserted into the nodose ganglion of the vagus and quantifed the number of microelectrode channels activated by one or both bipolar pairs of cuff contacts. In all animals, stimulation through both cuff pairs with at least one pulse width drove selective activation of fibers recorded by at least one MEA channel. Pulse width had a critical effect on the ability to selectively activate axons in the abdominal vagus nerve. At a short pulse width (i.e. 0.1 ms), evoked CAPs were recorded from relatively few MEA channels, with no responses in two animals at maximal stimulation amplitude. Conversely, at a long pulse width (i.e. 1 ms), evoked CAPs were recorded from most MEA channels, although very few were selectively activated by only one of the bipolar cuff pairs. A pulse width of 0.5 ms produced a balance between these two extremes and the most selective activation of MEA channels, with 1–7 MEA channels selectively activated by each bipolar cuff pair in each animal.
Abdominal VNS primarily activated fibers with conduction velocities below 3 m/s, although there were also responses from a smaller set of fibers with faster conduction velocities. The abdominal vagus nerve contains both afferent and efferent pathways that convey a variety of information to and from the abdominal organs27,28,29, as well as polysynaptic vago-vagal reflex pathways that project to the hindbrain and return to the periphery30. Previous studies have shown that the nerve is comprised primarily of Aδ and C fibers, which have conduction velocities of 3–30 and 0–3 m/s, respectively26. The larger diameter Aδ fibers are expected to have a lower threshold for extracellular electrical stimulation, although the much higher prevalence of C fibers was likely the primary reason we observed so many more responses from this fiber type27,31. It is important to consider, however, that the polysynaptic vago-vagal reflex pathways may have introduced additional conduction delays in our recordings which would have the effect of causing faster conducting Aδ fiber responses to appear as slower C fiber responses. Future studies should address this important limitation with methods that can distinguish direct activation from multi-synaptic responses, such as vagal rhizotomy and distal vagotomy32,33.
We also examined the patterns of activation across the MEA to determine if there was a structural relationship between the location of stimulation at the abdominal vagus nerve and the location of activation in the nodose ganglion. If there was a strong somatotopic relationship between the abdominal vagus nerve and the nodose ganglion, we would expect that the MEA channels with CAPs evoked by one bipolar stimulation pair would be spatially separate from those channels with CAPs evoked by the other bipolar stimulation pair. Instead, we found that there was no difference between the distances of MEA channels with responses evoked by the same bipolar pair and those evoked by different bipolar pairs. Further, high amplitude stimulation at the abdominal vagus nerve evoked responses on channels throughout the MEA, and there was no statistical difference in the location of centroids of MEA channels recording responses from the two bipolar cuff pairs, either for stimulation that maximized SI or at the maximum tested amplitude. These results suggest that it may be challenging to achieve selective recruitment of abdominal vagus nerve fibers with stimulation at the nodose ganglion because there is not substantial somatotopic organization within that structure. It also reveals a potential limitation of stimulation at the abdominal vagus, in that it may be difficult to predict a priori what population of vagal afferents will be activated by stimulation through a particular pair of cuff contacts. As such, clinical devices may require post-implantation tuning to achieve desired effects and tune out unwanted side-effects.
In addition to quantifying the selectivity of abdominal VNS, we also measured its effect on GI function through changes in GI myoelectric activity. In three animals, we demonstrated that stimulation through one bipolar pair of cuff contacts had little or no effect on myoelectric activity, while stimulation through the other pair of contacts strongly suppressed that activity. In all animals, stimulation through both pairs of cuff contacts drove little or no change in breathing rate, heart rate, blood pressure, or body temperature. These preliminary results support that VNS can have selective and differential effects on the GI system. It is unclear whether the effects of abdominal VNS on gastric myoelectric responses is due to stimulation of afferent or efferent vagal pathways.
While these results demonstrate the promise of abdominal VNS for selectively stimulating the GI system, several important limitations should be addressed. Even when we maximized SI, for most animals selective responses were still limited to fewer than 5 MEA channels per bipolar cuff pair (i.e. < 33% of all MEA channels). Still, this level of selective activation may be sufficient to achieve functional effects, as was evident in three animals in which stimulation selectively evoked a substantial change in GI myoelectric activity. Multiple factors may have limited the selectivity of abdominal VNS in these experiments. The binary search method used here is an efficient approach for quickly finding threshold, but it results in irregular sampling of the parameter space that likely limited our ability to optimize stimulation parameters post hoc. Future work should use parameter space sampling methods that are specifically designed for optimizing stimulation parameters for selectivity9. Additionally, while the ferret is a good model for GI function (e.g. because of the intact emetic response), the ferret abdominal vagus is small (500–750 µm diameter) and monofascicular, which likely limited our ability to achieve selective stimulation. Future work should focus on animal models with vagus nerve anatomy similar to humans, such as the pig34.
In addition to issues of VNS selectivity, there are broader limitations of the current study, including effects of anesthesia and physiological variability across animals. Our studies were performed using inhalational isoflurane. Isoflurane is known to dampen activity of neuronal cell bodies in peripheral ganglia35, which could have biased our data. Although this is likely not a problem in the context of VNS trials because even quiescent neurons can be stimulated, we did use basal neural activity to determine if the MEA insertion into the nodose ganglia had adequate numbers of active channels for VNS testing. Moreover, it is not well documented but isoflurane is reported to affect gastric myoelectric activity36. There are also substantial physiological differences between animals, which we have observed in in gastric myoelectric activity in awake and anesthetized ferrets37. In this study, we quantified selectivity via CAP recordings in one set of ferrets and measured changes in GI myoelectric activity with similar stimulation parameters in a different set of animals. This approach was necessary because implantation of the MEA in the nodose may have disrupted GI myoelectric activity, but it required us to quantify selectivity and measure functional changes in separate animals. Future studies should demonstrate selectivity, optimize VNS, and quantify functional effects in each subject.
Our focus in the current study, and the primary potential clinical impact, is in the optimization of VNS selectivity. Current commercial cuff electrode geometry is limited in the number of circumferential contacts, and thus selectivity. In the ferret, we have tested 4, 6, and 8 contact vagus nerve cuff electrodes, and 6 contacts (two rows of 3 circumferential contacts) appear to be the current limit for consistent manufacturing quality for abdominal vagus nerve of the ferret. This limitation further justifies that future experiments should be conducted in species with larger vagus nerve diameters, such as the cat or pig. In addition, there is potential to increase the efficiency of algorithmic optimization of VNS. We used recordings of the nodose ganglia to assess selectivity but it is unlikely this approach could be applied clinically; therefore, it will be important to use more accessible targets, such as serosal gastric myoelectric activity or cuff electrode recordings from the vagus nerve to confirm VNS selectivity [e.g.37,38]. Such an approach may lead to the application of closed-loop techniques to record and apply VNS to treat diseases of the abdominal cavity, such as GI disease, obesity, diabetes, and inflammatory disorders.
Methods
Animals
Nine adult male ferrets (weight: 1–1.7 kg; Marshall BioResources, North Rose, NY, USA) were used in this study. All experimental procedures were approved by the University of Pittsburgh Institutional Animal Care and Use Committee and were performed in accordance with ARRIVE and other relevant guidelines and regulations. Animals were housed in wire cages (62 × 74 × 46 cm) under a 12-h standard light cycle (lights on at 0700 h), in a temperature (20–24 °C) and humidity (30–70%) controlled environment. Food (ferret kibble: Mazuri Exotic Animal Nutrition, St. Louis, MI) and drinking water were freely available. Food was removed 3 h before induction of anesthesia. After each experiment, animals were euthanized with an injection of a 5 ml solution of SomnaSol (390 mg/ml pentobarbital sodium; 5 mg/ml phenytoin sodium; SomnaSol EUTHANASIA-III Solution, Henry Schein Animal Health, Dublin, Ohio, USA).
Surgical procedure
Anesthesia was induced and maintained with inhaled isoflurane (5% induction, 1–3% maintenance), and a tracheotomy was performed followed by insertion of an intratracheal tube to monitor respiration and deliver the anesthetic agent. Vital signs were monitored throughout the experiment, including blood pressure, heart rate, body temperature, and respiration rate, and isoflurane was adjusted to maintain a surgical plane of anesthesia (i.e., non-responsive to toe pinch). In six of the ferrets, a laparotomy was performed, followed by implantation of a four-contact nerve-cuff electrode (Micro Leads, Inc. Somerville, MA) around the ventral (n = 3) or the dorsal (n = 3) abdominal vagus trunk (Fig. 1). The cuff electrode contacts (area: 0.5 mm2, spaced 1 mm × 0.6 mm, 600 µm inner diameter) were arranged in two bipolar pairs with electrodes in each pair offset from each other by 90°, and the pairs spaced equally around the circumference of the nerve to provide current steering for targeted stimulation. A 32-channel MEA (4-by-8 electrode arrangement with 400 µm inter-electrode pitch, 1 mm long shanks; Black-Rock Microsystems, Salt Lake City, UT) was implanted into the nodose ganglion with a pneumatic inserter for rapid insertion through the epineurium. Recordings of spontaneous single-unit activity were used to verify insertion of the MEA into the nodose, and if necessary, additional impacts were applied to insert the device further. A platinum wire was placed near the nodose to act as a reference and another platinum wire was inserted under the skin as the recording ground. In the additional three ferrets, an identical procedure was used to implant a four-contact nerve-cuff electrode around the ventral abdominal vagus nerve and four planar electrodes, each with four contacts (Micro Leads, Inc. Somerville, MA), were sutured to the ventral gastric surface. The locations of these planar electrodes on the ventral gastric surface was similar to our prior study37.
Stimulation and data acquisition
A Grapevine Neural Interface Processor (Ripple, Salt Lake City, UT) and stimulation headstage (Nano2 + Stim) were used to deliver stimulation to the pairs of electrodes on the abdominal cuff, while a recording headstage (Nano2) was used to record evoked CAP signals from the nodose MEA and GI myoelectric activity from the serosal surface of the stomach. Nodose recordings were sampled at 30 kHz and filtered with a high pass filter at 150 Hz and low pass filter at 7500 Hz. A notch filter from 50 to 70 Hz was also applied to eliminate line noise. GI myoelectric signals were sampled at 30 kHz with a high pass filter at 0.1 Hz and a low pass filter at 7500 Hz. An adapter was inserted between the cuff and stimulation headstage to connect four output channels in parallel to increase the maximum stimulation amplitude per channel from 1.5 to 6 mA. For all experiment sessions, the impedance of electrodes was measured in response to a 1 kHz sinusoidal stimulus pre- and post-implantation as well as after each recording session to track changes in the electrode-tissue interface and ensure the connectivity of the system.
Stimulus-triggered averaging of compound action potentials
MEA recordings were analyzed using an automated algorithm written in Matlab (version 2017a, Mathworks, Natick, MA) to detect stimulation-evoked CAPs and minimize the time required to determine the parameters for threshold (amplitude and pulse width) for each bipolar cuff pair. For a given pulse width and amplitude combination, 120 pulses were delivered at 2 Hz through a longitudinal bipolar pair of electrodes in the cuff. The data recorded from each electrode within the MEA was segmented into a window around each stimulation event starting at the time point corresponding to a conduction velocity of 30 m/s (adjusted for each animal based on distance between the cuff and MEA) and ending at 498 ms post-stimulation. This interval removed any artifact generated by stimulation and allowed for detection of action potentials from neurons with slow conduction velocities and corresponding long latencies. The ensemble average of these windows was calculated to improve the signal-to-noise ratio. A sliding 1 ms moving root mean squared (RMS; Fig. 2) window with a step size of 0.1 ms was used to smooth this stimulus-triggered averaged (STA) signal. To remove artifacts in the signal caused by EMG activity or missing data packets, blanking was implemented over portions of the signal crossing 8 mV with a linearly interpolated signal between the start and end of the artifact periods. Baseline noise levels were measured for each MEA channel from stimulation recordings by first blanking a window centered at each stimulation event with a duration of 3 ms longer than the pulse width and secondly calculating the ensemble average of 120 randomly selected windows from this signal. The threshold for detecting a CAP response was set to between 2.4 and 2.6 standard deviations above the mean of the ensemble average of 120 segments of baseline recording. This multiplier of standard deviation was determined by creating a subset of ground truth data from expert evaluation of MEA recordings and generating ROC curves from a range of gain values, and selecting a multiplier that achieved a false positive rate of less than 10%. For each bipolar pair of cuff electrode contacts, stimulation threshold was defined for each pulse width as the minimum amplitude required to activate at least one MEA channel. In one animal (21–19), background noise in the MEA recordings prevented automated detection of responses to stimulation at 3000 µA and 1000 µs for one bipolar cuff pair, so responses were detected manually by an expert observer.
Threshold amplitude was determined for three different pulse widths (0.1, 0.5, and 1 ms). In one animal with a ventral cuff and one animal with a dorsal cuff, we also tested 0.4 ms pulses because no response was measured at 0.1 ms with amplitudes up to 3 mA. To determine if supra-threshold stimulation consistently and reliably evoked responses in nodose ganglion neurons, we randomly selected nine MEA channels across the six animals in which a response was detected, and examined the individual response to each of the 120 stimuli. At peri-threshold amplitudes, for 1080 total stimulus pulses, we found that 1058 evoked an MEA response. For the two randomly selected MEA channels that had fewer than 120 responses, the evoked response was only slightly larger than the background noise, which may have made the responses difficult to detect. Together, these results confirm that each supra-threshold stimulus pulse reliably evoked a response in the MEA and that we did not see effects such as conduction failure in our recordings.
The conduction velocity (CV) of nerve fibers responding to stimulation was calculated using the distance between the cuff and the electrode on the ventral or dorsal trunk, divided by the time between the stimulation event and when the RMS signal crossed the detection threshold. Signals were divided into time windows corresponding to 0.5 m/s increments, and only a single response could be detected within each time window, although multiple responses could be detected across different conduction velocities.
Optimization of selectivity stimulation parameters
In order to quantify our ability to selectively stimulate subpopulations of neurons in the abdominal vagus with a multi-contact cuff electrode, we performed a post-hoc optimization of stimulation parameters to maximize the number of nodose ganglion MEA channels with a CAP while minimizing the number of overlapping MEA channels with CAPs driven by both pairs of cuff electrodes. For a given pulse width, we varied pulse amplitude and calculated a selectivity index function of the form:
where SI is the selectivity index, TR is the total number of nodose electrodes with an evoked CAP response to stimulation with either cuff contact pair, TS is the total number of nodose electrode channels that recorded a response from stimulation through both cuff contact pairs (i.e., overlap in stimulation), T1:2 is the total number of nodose MEA channels with an evoked CAP response from stimulation through cuff contact pair 1:2, and T3:4 is the total number of nodose MEA electrodes that detected a response from stimulation through cuff contact pair 3–4. The first term in this equation calculates the number of MEA channels with non-overlapping responses to only one of the two cuff pairs, penalizing large overlaps in the number of MEA channels responding to stimulation through both cuff contact pairs. The second term penalizes imbalance between the responses from the two cuff contact pairs to prevent the condition in which stimulation through one cuff contact pair is maximized and no stimulation is delivered through the other cuff contact pair (i.e., it ensures that the responses are balanced between the two cuff contact pairs). We also excluded responses with more than three overlapping MEA channels (e.g. more than 10% overlap), to avoid producing results with excessive overlap that would not be functionally useful in tuning the effects of stimulation. Each term is multiplied by 0.5 so that SI ranges between 0 and 1. By maximizing this SI equation, we determined the stimulation amplitude that maximized the number of selectively responding MEA channels.
Nearest-neighbor analysis to quantify somatopic organization of nodose ganglion
To quantify the effect of stimulation amplitude on the location of activated MEA channels in the nodose ganglion, a nearest-neighbor analysis was performed for each pulse width at the stimulation amplitude that maximized SI and at the maximum amplitude that was tested (i.e. 3 mA). Animals where the maximum SI resulted in no selective responses for at least one cuff pair were excluded. For each MEA channel that recorded a response to abdominal VNS, we calculated the Euclidian distance to the nearest channel that recorded a response due to (1) stimulation from the same bipolar cuff pair with the same stimulation parameters and (2) from the opposing bipolar cuff pair with the stimulation parameters that elicited the maximum SI response. MEA channels that recorded non-selective responses from both cuff pairs were excluded from calculations for the opposing bipolar cuff pair because the distance to the nearest responding channel from the opposing bipolar cuff pair would always be zero but the distance to the nearest responding channel from the same bipolar cuff pair would always be greater than zero for these electrodes, resulting in a non-informative skew in the result.
To further quantify changes in the location of activation across the nodose ganglion in response to abdominal VNS, the centroid of all MEA channels recording responses from VNS were calculated for each cuff pair at the stimulation amplitude that maximized SI and at the maximum stimulation amplitude tested (i.e. 3 mA). The distance from one corner of the MEA was calculated for each channel, and the average distances of all channels that recorded a response from abdominal VNS was calculated to determine the centroid of activation for each set of stimulation parameters. The Euclidian distance between these centroids for each cuff pair was then calculated to quantify the relative distance across the array between locations of activation for the two bipolar cuff pairs. Violin plots were created using the custom MATLAB function violin.m created by Holger Hoffman39.
Gastrointestinal myoelectric activity as a functional measure of selectivity
We have previously shown that GI myoelectric activity can be used to identify the physiological state (normal, distended, pre-retch) of the stomach37. For this study, we analyzed data from a subset of those experiments to demonstrate the selective effect of VNS on GI myoelectric activity. Stimulation was delivered at the abdominal vagus through each bipolar cuff pair at a rate of 15 Hz with either 0.1 or 0.5 ms/phase symmetric pulses. Stimulation amplitude was increased on consecutive trials (0.2, 0.4, 0.6, 0.8, 1, and 5 mA) until it resulted in retching. For each trial, GI myoelectric activity, heart rate, breathing rate, body temperature, and blood pressure were recorded at baseline for 5 min followed by 2 min of stimulation.
GI myoelectric recordings were analyzed post-hoc using MATLAB (Mathworks, Natick, MA). For every planar electrode, the waveform recorded on each of the four contacts was averaged to generate a single GI myoelectric waveform for that electrode. Analysis methods for GI myoelectric activity were adopted from our prior study37. Briefly, each planar-averaged GI myoelectric signal was filtered using a low-pass Butterworth filter with a 2.5 Hz (150 cpm, 4th order) cutoff. The filtered signal was then downsampled to 10 Hz and a second low-pass Butterworth filter with a cut-off frequency of 0.3 Hz (18 cpm, 2nd order) was applied. GI signals that displayed a statistically significant dominant frequency in the 0 to 15 counts per minute (cpm) range were retained for all analysis. Each GI myoelectric signal was partitioned into 60-s segments with a 54-s (90%) overlap between consecutive segments and the power spectrum for each segment was computed using the fast Fourier transform (fft, bin size: 0.1 cpm). Additionally, the fraction of power in the normogastric range (8–11 cpm) during baseline and stimulation was compared to identify differential effects of vagus stimulation.
Data availability
All study data will be made available through the NIH SPARC data portal (https://sparc.science). Code that generates figures from the paper can be downloaded from https://github.com/cchorn/SPARC_vagus_stim_optimization.
References
Berthoud, H. R. & Neuhuber, W. L. Functional and chemical anatomy of the afferent vagal system. Auton. Neurosci. 85, 1–17 (2000).
Andrews, P. L. Vagal afferent innervation of the gastrointestinal tract. Prog. Brain Res. 67, 65–86 (1986).
Browning, K. N. Role of central vagal 5-HT3 receptors in gastrointestinal physiology and pathophysiology. Front. Neurosci. 9, 413 (2015).
Hillsley, K. & Grundy, D. Plasticity in the mesenteric afferent response to cisplatin following vagotomy in the rat. J. Auton. Nerv. Syst. 76, 93–98 (1999).
Kwon, C.-S. et al. Epilepsy and neuromodulation-randomized controlled trials. Brain Sci. 8, 69 (2018).
Apovian, C. M. et al. Two-year outcomes of vagal nerve blocking (vBloc) for the treatment of obesity in the ReCharge trial. Obes. Surg. 27, 169–176 (2017).
Ben-Menachem, E. Vagus nerve stimulation, side effects, and long-term safety. J. Clin. Neurophysiol. 18, 415–418 (2001).
Horn, C. C., Ardell, J. L. & Fisher, L. E. Electroceutical targeting of the autonomic nervous system. Physiology 34, 150–162 (2019).
Fisher, L. E., Tyler, D. J. & Triolo, R. J. Optimization of selective stimulation parameters for multi-contact electrodes. J. Neuroeng. Rehabil. 10, 25 (2013).
Brill, N. A. et al. Evaluation of high-density, multi-contact nerve cuffs for activation of grasp muscles in monkeys. J. Neural Eng. 15, 036003 (2018).
Polasek, K. H., Hoyen, H. A., Keith, M. W., Kirsch, R. F. & Tyler, D. J. Stimulation stability and selectivity of chronically implanted multicontact nerve cuff electrodes in the human upper extremity. IEEE Trans. Neural Syst. Rehabil. Eng. 17, 428–437 (2009).
Tigra, W. et al. Exploring selective neural electrical stimulation for upper limb function restoration. Eur. J. Transl. Myol. 26, 6035 (2016).
Fisher, L. E., Tyler, D. J., Anderson, J. S. & Triolo, R. J. Chronic stability and selectivity of four-contact spiral nerve-cuff electrodes in stimulating the human femoral nerve. J. Neural Eng. 6, 046010 (2009).
Tarler, M. D. & Mortimer, J. T. Comparison of joint torque evoked with monopolar and tripolar-cuff electrodes. IEEE Trans. Neural Syst. Rehabil. Eng. 11, 227–235 (2003).
Tan, D. W., Schiefer, M. A., Keith, M. W., Anderson, J. R. & Tyler, D. J. Stability and selectivity of a chronic, multi-contact cuff electrode for sensory stimulation in human amputees. J. Neural Eng. 12, 026002 (2015).
Wagner, F. B. et al. Targeted neurotechnology restores walking in humans with spinal cord injury. Nature 563, 65–71 (2018).
Nanivadekar, A. C., Ayers, C. A., Gaunt, R. A., Weber, D. J. & Fisher, L. E. Selectivity of afferent microstimulation at the DRG using epineural and penetrating electrode arrays. J. Neural Eng. 17, 016011 (2019).
Koopman, F. A. et al. Vagus nerve stimulation inhibits cytokine production and attenuates disease severity in rheumatoid arthritis. Proc. Natl. Acad. Sci. USA 113, 8284–8289 (2016).
McCrea, E. D. The abdominal distribution of the Vagus. J. Anat. 59, 18–40 (1924).
Harper, A. A., McSwiney, B. A. & Suffolk, S. F. Afferent fibres from the abdomen in the vagus nerves. J. Physiol. 85, 267–276 (1935).
Kupari, J., Häring, M., Agirre, E., Castelo-Branco, G. & Ernfors, P. An Atlas of vagal sensory neurons and their molecular specialization. Cell Rep. 27, 2508–2523 (2019).
Vallone, F. et al. Simultaneous decoding of cardiovascular and respiratory functional changes from pig intraneural vagus nerve signals. bioRxiv 2020.06.01.127050. https://doi.org/10.1101/2020.06.01.127050 (2020).
Pfeiffer, C. J. Surface topology of the stomach in man and the laboratory ferret. J. Ultrastruct. Res. 33, 252–262 (1970).
Fitzakerley, J. L. & Lucier, G. E. Connections of a vagal communicating branch in the ferret I. Pathways and cell body location. Brain Res. Bull. 20, 189–196 (1988).
Geddes, L. A. & Bourland, J. D. The strength-duration curve. IEEE Trans. Biomed. Eng. BME 32, 458–459 (1985).
Paintal, A. S. Vagal sensory receptors and their reflex effects. Physiol. Rev. 53, 159–227 (1973).
Page, A. J. & Blackshaw, L. A. An in vitro study of the properties of vagal afferent fibres innervating the ferret oesophagus and stomach. J. Physiol. 512, 907–916 (1998).
Duclaux, R., Mei, N. & Ranieri, F. Conduction velocity along the afferent vagal dendrites: a new type of fibre. J. Physiol. 260, 487–495 (1976).
Andrews, P. L. R. & Taylor, T. V. An electrophysiological study of the posterior abdominal vagus nerve in man. Clin Sci (Lond). 63(2), 169–73 (1982).
Andrews, P. L. & Wood, K. L. Vagally mediated gastric motor and emetic reflexes evoked by stimulation of the antral mucosa in anaesthetized ferrets. J. Physiol. 395, 1–16 (1988).
Mackay, T. W. & Andrews, P. L. A comparative study of the vagal innervation of the stomach in man and the ferret. J. Anat. 136, 449–481 (1983).
Norgren, R. & Smith, G.P. A method for selective section of vagal afferent or efferent axons in the rat. Am J Physiol. 267, R1136–41 (1994).
Andrews, P. L., & Scratcherd, T. The gastric motility patterns induced by direct and reflex excitation of the vagus nerves in the anaesthetized ferret. J Physiol. 302(1), 363–378. https://doi.org/10.1113/jphysiol.1980.sp013248 (1980).
Settell, M. L. et al. Functional vagotopy in the cervical vagus nerve of the domestic pig: implications for the study of vagus nerve stimulation. J. Neural Eng. 17, 026022 (2020).
Puil, E. & Gimbarzevsky, B. Modifications in membrane properties of trigeminal sensory neurons during general anesthesia. J. Neurophysiol. 58, 87–104 (1987).
Cheng, W., Chow, B. & Tam, P. K. Electrogastrographic changes in children who undergo day-surgery anesthesia. J. Pediatr. Surg. 34, 1336–1338 (1999).
Nanivadekar, A. C. et al. Machine learning prediction of emesis and gastrointestinal state in ferrets. PLoS ONE 14, e0223279 (2019).
Matthew P. et al. Flexible platform for biofeedback-driven control and personalization of electrical nerve stimulation therapy. IEEE Trans. Neural Syst. Rehabilitation Eng. 23(3), 475–484. https://doi.org/10.1109/TNSRE.2014.2351271 (2015).
Hoffman, H. violin.m—Simple Violin Plot Using Matlab Default Kernel Density Estimation. INRES (University of Bonn) (2015).
Acknowledgements
This work was supported by National Institutes of Health funding (Common fund SPARC Program award U18TR002205).
Author information
Authors and Affiliations
Contributions
L.E.F., C.C.H., J.A.S., and A.C.N. designed the experiments. C.C.H. and D.M.M. conducted the surgery. B.J.Y., D.M.M., S.F., and M.S. maintained anesthesia and monitored physiology. J.O., L.W., and B.M. designed and manufactured the multi-contact cuff electrodes. J.A.S., D.W.B., and A.C.N. performed data analysis. J.A.S., D.W.B., A.C.N., C.C.H. and L.E.F. drafted the manuscript. All authors read and supplied edits for the final manuscript.
Corresponding author
Ethics declarations
Competing interests
JO, LW, and BM are employees of Micro-Leads Inc. All other authors have no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Shulgach, J.A., Beam, D.W., Nanivadekar, A.C. et al. Selective stimulation of the ferret abdominal vagus nerve with multi-contact nerve cuff electrodes. Sci Rep 11, 12925 (2021). https://doi.org/10.1038/s41598-021-91900-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-91900-1
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.