Automatic analysis and visualization of microelectrode recording trajectories to the subthalamic nucleus: preliminary results

A Systematic Review of Neurophysiology-Based Localization Techniques Used in Deep Brain Stimulation Surgery of the Subthalamic Nucleus

OBJECTIVE

This systematic review is conducted to identify, compare, and analyze neurophysiological feature selection, extraction, and classification to provide a comprehensive reference on neurophysiology-based subthalamic nucleus (STN) localization.

MATERIALS AND METHODS

The review was carried out using the methods and guidelines of the Kitchenham systematic review and provides an in-depth analysis on methods proposed on STN localization discussed in the literature between 2000 and 2021. Three research questions were formulated, and 115 publications were identified to answer the questions.

RESULTS

The three research questions formulated are answered using the literature found on the respective topics. This review discussed the technologies used in past research, and the performance of the state-of-the-art techniques is also reviewed.

CONCLUSION

This systematic review provides a comprehensive reference on neurophysiology-based STN localization by reviewing the research questions other new researchers may also have.

Correlating Beta Oscillations from Intraoperative Microelectrode and Postoperative Implanted Electrode in Patients Undergoing Subthalamic Nucleus Deep Brain Stimulation for Parkinson Disease; A Feasibility Study

OBJECTIVE

We sought to investigate the feasibility of intraoperative local field potential (LFP) recording from the microelectrode during deep brain stimulation surgery for patients with Parkinson disease.

METHODS

Sixteen subthalamic nucleus recordings from 10 Parkinson disease patients who underwent deep brain stimulation surgery were included in this study. Signals from microelectrodes were amplified and differently filtered to display real-time single-unit neuronal activity and LFP simultaneously during surgery. LFP recordings were also recorded postoperatively from the implanted macroelectrodes and, power spectral density and peak frequency of beta oscillation of LFP (beta LFP) between 2 conditions were compared.

RESULTS

Stable intraoperative beta LFP were observed in 68.75% (11 of 16) cases. There was no significant difference of peak frequency between intraoperative and postoperative beta-LFP but significant difference of mean percentage of beta LFP was noted between 2 conditions.

CONCLUSIONS

Despite low signal-to-noise ratio and susceptibility to noises from external sources, this study shows that intraoperative recording of beta LFP using microelectrode is feasible. And, given that no significant difference in peak frequency of beta LFP between intraoperative and postoperative LFP was found, we suggest that not only intraoperative beta LFP can be used as a reliable surrogate for postoperative beta LFP, but it can also provide us an information for estimating the location with maximal power of beta oscillation within the subthalamic nucleus.

Deep brain stimulation: a review of the open neural engineering challenges

OBJECTIVE

Deep brain stimulation (DBS) is an established and valid therapy for a variety of pathological conditions ranging from motor to cognitive disorders. Still, much of the DBS-related mechanism of action is far from being understood, and there are several side effects of DBS whose origin is unclear. In the last years DBS limitations have been tackled by a variety of approaches, including adaptive deep brain stimulation (aDBS), a technique that relies on using chronically implanted electrodes on 'sensing mode' to detect the neural markers of specific motor symptoms and to deliver on-demand or modulate the stimulation parameters accordingly. Here we will review the state of the art of the several approaches to improve DBS and summarize the main challenges toward the development of an effective aDBS therapy.

APPROACH

We discuss models of basal ganglia disorders pathogenesis, hardware and software improvements for conventional DBS, and candidate neural and non-neural features and related control strategies for aDBS.

MAIN RESULTS

We identify then the main operative challenges toward optimal DBS such as (i) accurate target localization, (ii) increased spatial resolution of stimulation, (iii) development of in silico tests for DBS, (iv) identification of specific motor symptoms biomarkers, in particular (v) assessing how LFP oscillations relate to behavioral disfunctions, and (vi) clarify how stimulation affects the cortico-basal-ganglia-thalamic network to (vii) design optimal stimulation patterns.

SIGNIFICANCE

This roadmap will lead neural engineers novel to the field toward the most relevant open issues of DBS, while the in-depth readers might find a careful comparison of advantages and drawbacks of the most recent attempts to improve DBS-related neuromodulatory strategies.

Stop! border ahead: Automatic detection of subthalamic exit during deep brain stimulation surgery

BACKGROUND

Microelectrode recordings along preplanned trajectories are often used for accurate definition of the subthalamic nucleus (STN) borders during deep brain stimulation (DBS) surgery for Parkinson's disease. Usually, the demarcation of the STN borders is performed manually by a neurophysiologist. The exact detection of the borders is difficult, especially detecting the transition between the STN and the substantia nigra pars reticulata. Consequently, demarcation may be inaccurate, leading to suboptimal location of the DBS lead and inadequate clinical outcomes.

METHODS

We present machine-learning classification procedures that use microelectrode recording power spectra and allow for real-time, high-accuracy discrimination between the STN and substantia nigra pars reticulata.

RESULTS

A support vector machine procedure was tested on microelectrode recordings from 58 trajectories that included both STN and substantia nigra pars reticulata that achieved a 97.6% consistency with human expert classification (evaluated by 10-fold cross-validation). We used the same data set as a training set to find the optimal parameters for a hidden Markov model using both microelectrode recording features and trajectory history to enable real-time classification of the ventral STN border (STN exit). Seventy-three additional trajectories were used to test the reliability of the learned statistical model in identifying the exit from the STN. The hidden Markov model procedure identified the STN exit with an error of 0.04 ± 0.18 mm and detection reliability (error < 1 mm) of 94%.

CONCLUSIONS

The results indicate that robust, accurate, and automatic real-time electrophysiological detection of the ventral STN border is feasible. © 2016 International Parkinson and Movement Disorder Society.

Prediction of STN-DBS Electrode Implantation Track in Parkinson's Disease by Using Local Field Potentials

Optimal electrophysiological placement of the DBS electrode may lead to better long term clinical outcomes. Inter-subject anatomical variability and limitations in stereotaxic neuroimaging increase the complexity of physiological mapping performed in the operating room. Microelectrode single unit neuronal recording remains the most common intraoperative mapping technique, but requires significant expertise and is fraught by potential technical difficulties including robust measurement of the signal. In contrast, local field potentials (LFPs), owing to their oscillatory and robust nature and being more correlated with the disease symptoms, can overcome these technical issues. Therefore, we hypothesized that multiple spectral features extracted from microelectrode-recorded LFPs could be used to automate the identification of the optimal track and the STN localization. In this regard, we recorded LFPs from microelectrodes in three tracks from 22 patients during DBS electrode implantation surgery at different depths and aimed to predict the track selected by the neurosurgeon based on the interpretation of single unit recordings. A least mean square (LMS) algorithm was used to de-correlate LFPs in each track, in order to remove common activity between channels and increase their spatial specificity. Subband power in the beta band (11–32 Hz) and high frequency range (200–450 Hz) were extracted from the de-correlated LFP data and used as features. A linear discriminant analysis (LDA) method was applied both for the localization of the dorsal border of STN and the prediction of the optimal track. By fusing the information from these low and high frequency bands, the dorsal border of STN was localized with a root mean square (RMS) error of 1.22 mm. The prediction accuracy for the optimal track was 80%. Individual beta band (11–32 Hz) and the range of high frequency oscillations (200–450 Hz) provided prediction accuracies of 72 and 68% respectively. The best prediction result obtained with monopolar LFP data was 68%. These results establish the initial evidence that LFPs can be strategically fused with computational intelligence in the operating room for STN localization and the selection of the track for chronic DBS electrode implantation.

[Deep brain recording and length of surgery in stereotactic and functional neurosurgery for movement disorders]

OBJECTIVE

Our objectives were to study the length of multi-unit recordings (MURs) of brain activity in 20 years of movement disorder neurosurgeries and to determine the number of times in which it was necessary for the teams using single-unit recording (SUR) to explore all the electrode tracks in the simultaneously recorded sites (SRS).

MATERIAL AND METHOD

This was a retrospective descriptive statistical analysis of MUR length on 4,296 tracks in 952 surgeries. The exclusion criteria were: tracks with fewer than 5 recorded signals, tracks that had a signal length different from the habitual 2s, or there being unusual situations not related to the MUR, as well as the first 20 surgeries of each surgical target. This yielded a total of 3,448 tracks in 805 surgeries. We also determined the number of the total 952 surgeries in which all the tracks in the SURs of the SRS were explored.

RESULTS

The mean and its confidence interval (P=.05) of time per MUR track were 5.49±0.16min in subthalamic nucleus surgery, 8.82±0.24min in the medial or internal globus pallidus) and 18.51±1.31min in the ventral intermediate nucleus of the thalamus. For the total sum of tracks per surgery, in 75% of cases the total time was less than 39min in subthalamic nucleus, almost 42min in the medial or internal globus pallidus and less than 1h and 17min in ventral intermediate nucleus of the thalamus. All the tracks in the SUR SRS were explored in only 4.2% of the surgeries.

CONCLUSIONS

The impact of MUR on surgical time is acceptable for this guide in objective localization for surgical targets, without having to use several simultaneous electrodes (not all indispensable in most of the cases). Consequently, there is less risk for the patient.

Microelectrode recording duration and spatial density constraints for automatic targeting of the subthalamic nucleus

BACKGROUND

Accurate detection of the boundaries of the subthalamic nucleus (STN) in deep brain stimulation (DBS) surgery using microelectrode recording (MER) is considered to refine localization and may therefore improve clinical outcome. However, MER tends to extend operation time and its cost-utility balance has been debated.

OBJECTIVES

To quantify the tradeoff between accuracy of STN localization and the spatial and temporal parameters of MER that effect the operation time using an automated detection method.

METHODS

We retrospectively estimated the accuracy of STN detection on data from 100 microelectrode trajectories. Our dense (average step = 0.12 mm) and long (average duration = 22.5 s) MER data was downsampled in the spatial and temporal domains. Then, the STN borders were detected automatically on both the downsampled and original data and compared to each other.

RESULTS

With a recording duration of 16 s, average accuracy for detecting STN entry ranged from 0.06 mm for a 0.1-mm step to 0.51 mm for a 1.0-mm step. Smaller effects were found along the temporal axis. For example, a 0.1-mm recording step yielded an STN entry average accuracy ranging from 0.06 mm for a 16-second recording duration to 0.16 mm for 0.1 s.

CONCLUSIONS

STN entry detection error was about half of the step size. Sampling duration of STN activity can be minimized to 1 s/record without compromising accuracy. We conclude that bilateral DBS surgery time utilizing MER may be significantly shortened without compromising targeting accuracy.

Depth-time interpolation of feature trends extracted from mobile microelectrode data with kernel functions

BACKGROUND/AIMS

Microelectrode recording (MER) is necessary for precision localization of target structures such as the subthalamic nucleus during deep brain stimulation (DBS) surgery. Attempts to automate this process have produced quantitative temporal trends (feature activity vs. time) extracted from mobile MER data. Our goal was to evaluate computational methods of generating spatial profiles (feature activity vs. depth) from temporal trends that would decouple automated MER localization from the clinical procedure and enhance functional localization in DBS surgery.

METHODS

We evaluated two methods of interpolation (standard vs. kernel) that generated spatial profiles from temporal trends. We compared interpolated spatial profiles to true spatial profiles that were calculated with depth windows, using correlation coefficient analysis.

RESULTS

Excellent approximation of true spatial profiles is achieved by interpolation. Kernel-interpolated spatial profiles produced superior correlation coefficient values at optimal kernel widths (r = 0.932-0.940) compared to standard interpolation (r = 0.891). The choice of kernel function and kernel width resulted in trade-offs in smoothing and resolution.

CONCLUSIONS

Interpolation of feature activity to create spatial profiles from temporal trends is accurate and can standardize and facilitate MER functional localization of subcortical structures. The methods are computationally efficient, enhancing localization without imposing additional constraints on the MER clinical procedure during DBS surgery.

Novel signal-dependent filter bank method for identification of multiple basal ganglia nuclei in Parkinsonian patients

Microelectrode recordings are a valuable tool for assisting localization targets during deep brain stimulation procedures in Parkinson's disease neurosurgery. Attempts to automate and standardize this process have been limited by variability in patient neurophysiology and strong dynamics of microelectrode recordings. In this paper, a methodology for the identification of basal ganglia nuclei is presented that is based on a signal-dependent filter bank method using microelectrode recordings. The method is a customized realization of the discrete wavelet transform via the lifting scheme that is optimally tuned by genetic algorithms. Using this method, unique mother wavelet functions that exhibit an adaptable spectrum to the microelectrode recording dynamic are generated. Additionally, by extracting morphological features from the space-transformed microelectrode recording, it is possible to integrate them into three-dimensional (3D) feature spaces with maximum class separability. Finally, high discriminant feature spaces are fed into basic classifiers to recognize up to four basal nuclei. Comparison with several existing wavelets highlights the characteristics of new mother wavelets. Additionally, classification results show that identification of addressed nuclei in the basal ganglia can be performed with 95% confidence.

On detecting determinism and nonlinearity in microelectrode recording signals: approach based on non-stationary surrogate data methods

Two new surrogate methods, the Small Shuffle Surrogate (SSS) and the Truncated Fourier Transform Surrogate (TFTS), have been proposed to study whether there are some kind of dynamics in irregular fluctuations and if so whether these dynamics are linear or not, even if this fluctuations are modulated by long term trends. This situation is theoretically incompatible with the assumption underlying previously proposed surrogate methods. We apply the SSS and TFTS methods to microelectrode recording (MER) signals from different brain areas, in order to acquire a deeper understanding of them. Through our methodology we conclude that the irregular fluctuations in MER signals possess some determinism.

Towards high accuracy classification of MER signals for target localization in Parkinson's disease

In recent years Microelectrode recording (MER) analysis has proved to be a powerful localization tool of basal ganglia for Parkinson disease's treatment, especially the Subthalamic Nucleus (STN). In this paper, a signal-dependent method is presented for identification of the STN and other brain zones in Parkinsonian patients. The proposed method, refereed as optimal wavelet feature extraction method (OWFE), is constructed by lifting schemes (LS), which are a flexible and fast implementation of the wavelet transform (WT). The operators in the LS are optimized by means of Genetic Algorithms and Lagrange multipliers considering information contained in MER signals. Then a basic Bayesian classifier (LDC) is used to identify STN and other types of basal ganglia nuclei. The proposed method introduced several advantages from similar works reported in literature. First, the method is signal-dependent and non a priori information is required to decompose the MER signal. Second, the classification accuracy is mostly depended on the feature selection stage because it is not enhanced by elaborated classifiers such as support vector machines or hidden Markov models. Finally, the generalization property of the OWFE has been validated with two databases and different types of classifiers such as k-NN classifier and quadratic Bayesian classifier (QDC). Results have shown that proposed method is able to identify the STN with average accuracy superior than 97%.

Automated 3-dimensional brain atlas fitting to microelectrode recordings from deep brain stimulation surgeries

OBJECTIVE

Deep brain stimulation (DBS) surgeries commonly rely on brain atlases and microelectrode recordings (MER) to help identify the target location for electrode implantation. We present an automated method for optimally fitting a 3-dimensional brain atlas to intraoperative MER and predicting a target DBS electrode location in stereotactic coordinates for the patient.

METHODS

We retrospectively fit a 3-dimensional brain atlas to MER points from 10 DBS surgeries targeting the subthalamic nucleus (STN). We used a constrained optimization algorithm to maximize the MER points correctly fitted (i.e., contained) within the appropriate atlas nuclei. We compared our optimization approach to conventional anterior commissure-posterior commissure (AC/PC) scaling, and to manual fits performed by four experts. A theoretical DBS electrode target location in the dorsal STN was customized to each patient as part of the fitting process and compared to the location of the clinically defined therapeutic stimulation contact.

RESULTS

The human expert and computer optimization fits achieved significantly better fits than the AC/PC scaling (80, 81, and 41% of correctly fitted MER, respectively). However, the optimization fits were performed in less time than the expert fits and converged to a single solution for each patient, eliminating interexpert variance.

CONCLUSIONS AND SIGNIFICANCE

DBS therapeutic outcomes are directly related to electrode implantation accuracy. Our automated fitting techniques may aid in the surgical decision-making process by optimally integrating brain atlas and intraoperative neurophysiological data to provide a visual guide for target identification.