Electrical Stimulation Artifacts Research Articles

Single-pulse electrical stimulation artifact removal using the novel matching pursuit-based artifact reconstruction and removal method (MPARRM)

Objective. Single-pulse electrical stimulation (SPES) has been widely used to probe effective connectivity. However, analysis of the neural response is often confounded by stimulation artifacts. We developed a novel matching pursuit-based artifact reconstruction and removal method (MPARRM) capable of removing artifacts from stimulation-artifact-affected electrophysiological signals. Approach. To validate MPARRM across a wide range of potential stimulation artifact types, we performed a bench-top experiment in which we suspended electrodes in a saline solution to generate 110 types of real-world stimulation artifacts. We then added the generated stimulation artifacts to ground truth signals (stereoelectroencephalography signals from nine human subjects recorded during a receptive speech task), applied MPARRM to the combined signal, and compared the resultant denoised signal with the ground truth signal. We further applied MPARRM to artifact-affected neural signals recorded from the hippocampus while performing SPES on the ipsilateral basolateral amygdala in nine human subjects. Main results. MPARRM could remove stimulation artifacts without introducing spectral leakage or temporal spread. It accommodated variable stimulation parameters and recovered the early response to SPES within a wide range of frequency bands. Specifically, in the early response period (5–10 ms following stimulation onset), we found that the broadband gamma power (70–170 Hz) of the denoised signal was highly correlated with the ground truth signal ( R=0.98±0.02 , Pearson), and the broadband gamma activity of the denoised signal faithfully revealed the responses to the auditory stimuli within the ground truth signal with 94%±1.47% sensitivity and 99%±1.01% specificity. We further found that MPARRM could reveal the expected temporal progression of broadband gamma activity along the anterior-posterior axis of the hippocampus in response to the ipsilateral amygdala stimulation. Significance. MPARRM could faithfully remove SPES artifacts without confounding the electrophysiological signal components, especially during the early-response period. This method can facilitate the understanding of the neural response mechanisms of SPES.

Design and Implementation ofReal-time Electrical Stimulation Artifact Suppression based on STM32

A biosignal is used as a control signal for electrical stimulation to restore weakened muscle function due to damage to the central nervous system. In patients with central nervous system damage, sufficient muscle contraction does not occur spontaneously. In this case, applying electrical stimulation can cause normal muscle contraction. However, it is necessary to remove the electrical stimulation artifact caused by the electrical stimulation. This paper describes a system design that removes electrical stimulation artifact in real time using a Cortex-M4-based STM32F processor. The STM32F is a very advantageous MCU for such DSPs, especially because it has a built-in floating point operator. Using STM32F's various high-performance peripherals (12-bit parallel ADC and 12-bit DAC, UART, Timer), an optimized embedded system was implemented.In this paper, the simulated and real-time results were compared and evaluated with the designed fir filter. In addition, the performance of the filter was evaluated through frequency analysis. As a result, it was verified that a high-performance 32-bit STM32F with floating point calculator and various peripherals is suitable for real-time signal processing

A Method for Suppressing Electrical Stimulation Artifacts from Electromyography.

When surface electromyography (EMG) signal is used in a real-time functional electrical stimulation (FES) system for feedback control, the artifact from electrical stimulation is a key challenge for EMG signal processing. To address this challenge, this study proposes a novel method to suppress stimulation artifacts in the EMG-driven closed-loop FES system. The proposed method is inspired by an experimental study that compares artifacts generated by electrical stimulations with different current intensities. It is found that (1) spikes of stimulation artifacts are susceptible to the current intensity and (2) tailing components are similar under different current intensities. Based on these observations, the proposed method combines the blanking and template subtracting strategies for suppressing stimulation artifact. The length of blanking window for suppressing the stimulation spike is adaptively determined by a spike detection algorithm and the first-order derivative analysis of signal. An autoregressive model is used to estimate the tailing part of stimulation artifact, which is an adaptive template for subtracting the artifact. The proposed method is evaluated on both semi-synthetic and experimental datasets. Verified on the semi-synthetic dataset, the proposed method achieves better performance than the classic blanking method. Validated on the experimental dataset, the proposed method substantially decreases the power of stimulation artifact in the EMG. These results indicate that the proposed method can effectively suppress the stimulation artifact while retains the useful EMG signal for an EMG-driven FES system.

Towards Safe Infrared Nerve Stimulation: A Systematic Experimental Approach.

Neural activation by infrared nerve stimulation (INS) gains growing interest as a potential alternative to conventional electric nerve stimulation, since unambiguous advantages like contact-free operation, enhanced spatial selectivity and lack of (electrical) stimulation artifacts are promising for both future electrophysiological research and clinical application. For the systematic investigation of laser nerve activation, we recently introduced a novel experimental approach. Comprising a defined focused beam profile, it enables remote controlled, contact-free pulsed laser stimulation of the rat sciatic nerve, simultaneous to high-speed temperature measurement in vivo. Up to now, successful neural activation with single laser pulses (2 - 6 mJ) was observed in all performed experiments, however, it strongly depended on the particular nerve location. Hence, we depict the investigation of spatial dependency of the nerve response and identify `regions of excitability' on the nerve surface, that are highly susceptible to INS. By means of thermal imaging, we simultaneously monitored the nerve surface temperature, where we observed progressing temperature build-up during single pulse stimulation with repetition rates above 4 Hz. In this work, we present current results of our ongoing research.

Diffusion geometry approach to efficiently remove electrical stimulation artifacts in intracranial electroencephalography

Objective. Cortical oscillations, electrophysiological activity patterns, associated with cognitive functions and impaired in many psychiatric disorders can be observed in intracranial electroencephalography (iEEG). Direct cortical stimulation (DCS) may directly target these oscillations and may serve as therapeutic approaches to restore functional impairments. However, the presence of electrical stimulation artifacts in neurophysiological data limits the analysis of the effects of stimulation. Currently available methods suffer in performance in the presence of nonstationarity inherent in biological data. Approach. Our algorithm, shape adaptive nonlocal artifact removal (SANAR) is based on unsupervised manifold learning. By estimating the Euclidean median of k-nearest neighbors of each artifact in a nonlocal fashion, we obtain a faithful representation of the artifact which is then subtracted. This approach overcomes the challenges presented by nonstationarity. Main results. SANAR is effective in removing stimulation artifacts in the time domain while preserving the spectral content of the endogenous neurophysiological signal. We demonstrate the performance in a simulated dataset as well as in human iEEG data. Using two quantitative measures, that capture how much of information from endogenous activity is retained, we demonstrate that SANAR’s performance exceeds that of one of the widely used approaches, independent component analysis, in the time domain as well as the frequency domain. Significance. This approach allows for the analysis of iEEG data, single channel or multiple channels, during DCS, a crucial step in advancing our understanding of the effects of periodic stimulation and developing new therapies.

Electrical stimulus artifact cancellation and neural spike detection on large multi-electrode arrays.

Simultaneous electrical stimulation and recording using multi-electrode arrays can provide a valuable technique for studying circuit connectivity and engineering neural interfaces. However, interpreting these measurements is challenging because the spike sorting process (identifying and segregating action potentials arising from different neurons) is greatly complicated by electrical stimulation artifacts across the array, which can exhibit complex and nonlinear waveforms, and overlap temporarily with evoked spikes. Here we develop a scalable algorithm based on a structured Gaussian Process model to estimate the artifact and identify evoked spikes. The effectiveness of our methods is demonstrated in both real and simulated 512-electrode recordings in the peripheral primate retina with single-electrode and several types of multi-electrode stimulation. We establish small error rates in the identification of evoked spikes, with a computational complexity that is compatible with real-time data analysis. This technology may be helpful in the design of future high-resolution sensory prostheses based on tailored stimulation (e.g., retinal prostheses), and for closed-loop neural stimulation at a much larger scale than currently possible.

Mechanomyography-Based Wearable Monitor of Quasi-Isometric Muscle Fatigue for Motor Neural Prostheses.

A motor neural prosthesis based on surface functional electrical stimulation (sFES) can restore functional movement (e.g., standing, walking) in patients with a spinal cord injury (SCI). sFES generates muscle contractions in antigravity muscles and allows balance-assisted standing. This induced standing has several benefits, such as improved cardiovascular function, decreased incidence of urinary infections, reduced joint contractures, and muscle atrophy. The duration of sFES assisted standing is limited due to the quick onset of muscle fatigue. Currently, there is no method available to reliably estimate real-time muscle fatigue during sFES. Simply monitoring the M-wave changes is not suitable due to the high signal disturbances that arise during multi-channel electrical stimulation. Mechanomyography (MMG) is immune to electrical stimulation artifacts and can be used to detect subtle vibrations on the surface of the skin related to activation of the underlying muscle's motor units (MU). The aim of this study was to develop a method for detecting muscle fatigue brought on by sFES. The method was tested in three different heads of the quadriceps muscle in SCI patients during electrically elicited quasi-isometric contraction. Six spinal cord-injured male volunteers, with no voluntary control of the quadriceps muscle participated in the study. Electrical bursts of voltage-controlled monophasic square pulses at frequencies of 1 kHz (50% duty cycle) at 50 Hz (15% duty cycle) were used to generate thigh muscle contractions that controlled the knee joint in the sagittal plane. The pulse amplitudes were set to position the knee joint at a 5° angle from the horizontal plane and when the knee angle dropped to 20° (e.g., the quadriceps were unable to hold the lower leg in the desired position), the test was terminated. Two data segments lasting 10 s each, at the beginning and end of each test, were analyzed. The muscle contraction was assessed by MMG sensors positioned on the rectus femoris, vastus lateralis, and vastus medialis muscles. Data segments were decomposed into 11 frequency bands using a Cauchy wavelet transform. In the initial time interval (non-fatigued muscle), the power peak was concentrated in the 11.31 Hz frequency band. In the final interval (muscle fatigued) this peak shifted to lower frequencies (2 and 6 Hz frequency bands). The decreased frequency was most prominent during the last 4 s of the recordings. It was shown that MMG could be used as a real-time indicator of muscle fatigue during FES-induced isometric contraction of quadriceps; hence, MMG could be used in closed-loop control as a fatigue detector. Subsequent studies for non-isometric contractions could possibly lead to prediction of muscle fatigue before contractile failure during functional use of the muscle.

Phase properties of transcranial electrical stimulation artifacts in electrophysiological recordings

Monitoring brain activity during transcranial electric stimulation (tES) is an attractive approach for causally studying healthy and diseased brain activity. Yet, stimulation artifacts complicate electrophysiological recordings during tES. Design and evaluation of artifact removal methods require a through characterization of artifact features, i.e. characterization of the transfer function that defines the relationship between the tES stimulation current and tES artifacts. Here we characterize the phase relationship between stimulation current and tES artifacts in EEG and MEG. We show that stimulation artifacts are not pure in-phase or anti-phase signals, but that non-linear mechanisms induce steady phase deflections relative to the stimulation current. Furthermore, phase deflections of stimulation artifacts are slightly modulated by each heartbeat and respiration. For commonly used stimulation amplitudes, artifact phase deflections correspond to signals several times bigger than normal brain signal. Moreover, the strength of phase deflections varies with stimulation frequency. These phase effects should be accounted for during artifact removal and when comparing recordings with different stimulation frequencies. We summarize our findings in a mathematical model of tES artifacts and discuss how this model can be used in simulations to design and evaluate artifact rejection techniques. To facilitate this research, all raw data of this study is made freely available.

Auditory steady state responses and cochlear implants: Modeling the artifact-response mixture in the perspective of denoising.

Auditory steady state responses (ASSRs) in cochlear implant (CI) patients are contaminated by the spread of a continuous CI electrical stimulation artifact. The aim of this work was to model the electrophysiological mixture of the CI artifact and the corresponding evoked potentials on scalp electrodes in order to evaluate the performance of denoising algorithms in eliminating the CI artifact in a controlled environment. The basis of the proposed computational framework is a neural mass model representing the nodes of the auditory pathways. Six main contributors to auditory evoked potentials from the cochlear level and up to the auditory cortex were taken into consideration. The simulated dynamics were then projected into a 3-layer realistic head model. 32-channel scalp recordings of the CI artifact-response were then generated by solving the electromagnetic forward problem. As an application, the framework’s simulated 32-channel datasets were used to compare the performance of 4 commonly used Independent Component Analysis (ICA) algorithms: infomax, extended infomax, jade and fastICA in eliminating the CI artifact. As expected, two major components were detectable in the simulated datasets, a low frequency component at the modulation frequency and a pulsatile high frequency component related to the stimulation frequency. The first can be attributed to the phase-locked ASSR and the second to the stimulation artifact. Among the ICA algorithms tested, simulations showed that infomax was the most efficient and reliable in denoising the CI artifact-response mixture. Denoising algorithms can induce undesirable deformation of the signal of interest in real CI patient recordings. The proposed framework is a valuable tool for evaluating these algorithms in a controllable environment ahead of experimental or clinical applications.

Feasibility of magnetoencephalographic source imaging in patients with thalamic deep brain stimulation for epilepsy

SummarySource localization of interictal spikes in patients with medically refractory epilepsy is the most common clinical application of magnetoencephalography (MEG). In recent decades, many patients with intractable epilepsy have been treated with various forms of neurostimulation, including thalamic deep brain stimulation (DBS). Patients with suboptimal seizure control after DBS might in some cases benefit from further investigations for resective epilepsy surgery, including MEG source imaging (MSI). We sought to determine the feasibility and accuracy of MSI in the setting of active thalamic DBS. Simultaneous EEG/MEG was obtained in a patient using an Elekta 306‐channel MEG system, with high‐frequency (100 Hz) DBS of the thalamic anterior nuclei cycling between on and off states. Magnetic artifacts associated with the DBS apparatus were successfully suppressed using the spatiotemporal signal space separation (tSSS) method. Electrical stimulation artifact was removed by standard digital low‐pass filtering. Dipole source modeling results for spike foci in frontal and posterior temporal regions were comparable between stimulation on and stimulation off states, and the source solutions corresponded well to the localization of spikes documented by intracranial EEG. MSI is thus feasible and source solutions can be accurate when performed in patients with active thalamic DBS for epilepsy.

Electrically and Hybrid-Induced Muscle Activations: Effects of Muscle Size and Fiber Type.

The effect of three electrical stimulation (ES) frequencies (10, 35, and 50 Hz) on two muscle groups with different proportions of fast and slow twitch fibers (abductor pollicis brevis (APB) and vastus lateralis (VL)) was explored. We evaluated the acute muscles’ responses individually and during hybrid activations (ES superimposed by voluntary activations). Surface electromyography (sEMG) and force measurements were evaluated as outcomes. Ten healthy adults (mean age: 24.4 ± 2.5 years) participated after signing an informed consent form approved by the university Institutional Review Board. Protocols were developed to: 1) compare EMG activities during each frequency for each muscle when generating 25% Maximum Voluntary Contraction (MVC) force, and 2) compare EMG activities during each frequency when additional voluntary activation was superimposed over ES-induced 25% MVC to reach 50% and 75% MVC. Empirical mode decomposition (EMD) was utilized to separate ES artifacts from voluntary muscle activation. For both muscles, higher stimulation frequency (35 and 50Hz) induced higher electrical output detected at 25% of MVC, suggesting more recruitment with higher frequencies. Hybrid activation generated proportionally less electrical activity than ES alone. ES and voluntary activations appear to generate two different modes of muscle recruitment. ES may provoke muscle strength by activating more fatiguing fast acting fibers, but voluntary activation elicits more muscle coordination. Therefore, during the hybrid activation, less electrical activity may be detected due to recruitment of more fatigue-resistant deeper muscle fibers, not reachable by surface EMG.

A wavelet based method for electrical stimulation artifacts removal in electromyogram

A technique for artifact removal based on the continuous wavelet transform is presented. It uses common mother wavelets to find the temporal localization of stimulation artifacts on electromyogram signal recording during an electrical surface evoked contraction of a muscle. This method is applied with different kinds of stimulation pulse parameters including shape and duration changes. This method is used with standard stimulation pulse waveforms such as monophasic or biphasic ones. It can also be applied when the artifacts and M waves are in the same range of amplitude where threshold techniques are inefficient. Lastly, a method to determine which mother wavelet efficiently removed artifacts is proposed, results indicate the Haar wavelet performs best among fourteen tested wavelets.

Infrared Neural Stimulation of Thalamocortical Brain Slices

Infrared neural stimulation (INS) has been well characterized in the peripheral nervous system, and has been shown to enable stimulation with high spatial precision and without causing the typical electrical stimulation artifact on the recording electrode. The next step in the development of INS is to demonstrate feasibility to stimulate neurons located in the central nervous system (CNS). Thalamocortical brain slices were used to establish feasibility of INS in the CNS and to optimize laser parameters. Infrared light was used to evoke action potentials in the brain slice with no electrical stimulation artifact. This response was blocked by the application of tetrodotoxin demonstrating neurological origin of the recorded signal. Threshold radiant exposure decreased as the absorption coefficient of the wavelength of light increased. Higher repetition rates lead to a decrease in threshold radiant exposure, and threshold radiant exposure was found to decrease as the spot size diameter increased. Additionally, neuronal responses to INS were intracellularly recorded demonstrating artifact free electrical recordings. The results from this paper lay the foundation for future in vivo studies to develop INS for CNS stimulation.

The Cortical Chronometry of Electrogustatory Event-related Potentials

Electrogustometry (EGM) is the standard tool to assess gustatory functions in clinical environments. The stimulation elicits a percept often described as metallic, sour or salty, also referred to as electric taste. To date, the neuronal mechanisms that underlie electric taste perception are not yet fully understood. Electroencephalographic (EEG) approaches will certainly complement behavioral procedures and, furthermore, extend the understanding of gustatory processing in general and disturbances of gustatory functions in particular. We used anodal pulses applied to the tip of the participants' tongue while EEG was recorded. The major disadvantage of combining EEG and EGM, namely the electrical stimulation artifact, was overcome by means of Independent Component Analysis (ICA), which separated the EGM artifact from the neural portion of the EEG. After artifact correction, we found a largely uncontaminated electrogustatory event-related potential (eGERP) at both individual and group level. Furthermore, source analysis revealed an early involvement of bilateral insular cortices and the adjacent operculi, the areas comprising the primary taste cortex. The procedures, described in detail, pave the way for the eGERP to become an affordable and objective tool for the assessment of taste function, and thus to complement behavioral measures (i.e. EGM detection thresholds). Furthermore, they render the access to different levels of the electrogustatory processing pathway possible and by doing so they may aid the identification and localisation of lesions that cause taste disturbances.

Noncontact Stimulation of the Cavernous Nerves in the Rat Prostate Using a Tunable-Wavelength Thulium Fiber Laser

Laser nerve stimulation has recently been studied in neuroscience as an alternative to electrical stimulation. Its advantages include noncontact stimulation, better spatial selectivity, and elimination of electrical stimulation artifacts. This study explored laser stimulation of the rat cavernous nerves as a potential alternative to electrical nerve mapping during nerve-sparing radical prostatectomy. The cavernous nerves were surgically exposed in 10 male rats. A thulium fiber laser stimulated the nerves, with a wavelength of 1870 nm, pulse energy of 7.5 mJ, radiant exposure of 1 J/cm2, pulse duration of 2.5 msec, pulse rate of 10 Hz, and 1-mm laser spot diameter, for a stimulation time of 60 sec. A significant increase in the intracavernosal pressure was detected on laser stimulation, with pressure returning to baseline values after stimulation ceased. This study demonstrates the feasibility of noncontact stimulation of the cavernous nerves using near-infrared laser radiation.

The removal of electrical stimulation artifact from surface EMG

The removal of electrical stimulation artifact from surface EMG

Spatial localization of saccade targets. I. Compensation for stimulation-induced perturbations in eye position.

1. Monkeys were trained to look to brief visual targets presented in a completely darkened room. On some trials, after the visual target disappeared but before a saccade to the target could be initiated, the eyes were driven to another position in the orbit by electrical stimulation of the superior colliculus. Retinocentric models of the saccadic system predict that a saccade with a predetermined distance and direction based entirely on retinal error will occur. If this were the case, gaze would miss the target location by a distance and direction equal to the vector of the stimulation-induced movement. Spatial models assume that the retinal error signal will be combined with information about the change in eye position produced by stimulation and predict that the animal will look to the position of the target in space. 2. Results confirm the predictions of spatial models. Animals compensated for the stimulation-induced perturbation by looking to the position of the target in space. The result predicted by retinocentric models-a saccade with a direction and amplitude based on retinal error alone-was never observed. 3. The eye movement that compensated for the change in eye position produced by stimulation was a saccade, not a passive, lowvelocity movement to an orbital position of mechanical equilibrium established by a tonic pattern of motoneuron activation specified by the visual target. This indicates that a new saccade command, based on stored information about the location of the retinal image and information about the new position of the eyes, had been issued. Computation of the vector of the compensatory saccade does not necessarily increase the latency of target acquisition. The interval between the end of a stimulation-induced saccade and the beginning of the compensatory saccade was frequently 20 ms or less, permitting the animal to acquire the target with a normal latency.