ORIGINAL RESEARCH article

Front. Hum. Neurosci., 08 September 2022

Sec. Brain Imaging and Stimulation

Volume 16 - 2022 | https://doi.org/10.3389/fnhum.2022.958521

Subthalamic high-beta oscillation informs the outcome of deep brain stimulation in patients with Parkinson's disease

  • 1. Division of Movement Disorders, Department of Neurology, Chang Gung Memorial Hospital, Taoyuan, Taiwan

  • 2. Neuroscience Research Center, Chang Gung Memorial Hospital, Taoyuan, Taiwan

  • 3. College of Medicine, Chang Gung University, Taoyuan, Taiwan

  • 4. Department of Neurosurgery, Chang Gung Memorial Hospital, Taoyuan, Taiwan

  • 5. Department of Mathematics, National Taiwan University, Taipei, Taiwan

  • 6. Department of Public Health, Biostatistics Consulting Center, College of Medicine, Chang Gung University, Taoyuan, Taiwan

  • 7. Department of Obstetrics and Gynecology, Chang Gung Memorial Hospital, Chiayi, Taiwan

  • 8. Department of Mathematics, Duke University, Durham, NC, United States

  • 9. Department of Statistical Science, Duke University, Durham, NC, United States

  • 10. Department of Neuroradiology, Chang Gung Memorial Hospital, Taoyuan, Taiwan

  • 11. Professor Lu Neurological Clinic, Taoyuan, Taiwan

Article metrics

View details

26

Citations

6k

Views

1,9k

Downloads

Abstract

Background:

The therapeutic effect of deep brain stimulation (DBS) of the subthalamic nucleus (STN) for Parkinson's disease (PD) is related to the modulation of pathological neural activities, particularly the synchronization in the β band (13–35 Hz). However, whether the local β activity in the STN region can directly predict the stimulation outcome remains unclear.

Objective:

We tested the hypothesis that low-β (13–20 Hz) and/or high-β (20–35 Hz) band activities recorded from the STN region can predict DBS efficacy.

Methods:

Local field potentials (LFPs) were recorded in 26 patients undergoing deep brain stimulation surgery in the subthalamic nucleus area. Recordings were made after the implantation of the DBS electrode prior to its connection to a stimulator. The maximum normalized powers in the theta (4–7 Hz), alpha (7–13 Hz), low-β (13–20 Hz), high-β (20–35 Hz), and low-γ (40–55 Hz) subbands in the postoperatively recorded LFP were correlated with the stimulation-induced improvement in contralateral tremor or bradykinesia–rigidity. The distance between the contact selected for stimulation and the contact with the maximum subband power was correlated with the stimulation efficacy. Following the identification of the potential predictors by the significant correlations, a multiple regression analysis was performed to evaluate their effect on the outcome.

Results:

The maximum high-β power was positively correlated with bradykinesia–rigidity improvement (rs = 0.549, p < 0.0001). The distance to the contact with maximum high-β power was negatively correlated with bradykinesia–rigidity improvement (rs = −0.452, p < 0.001). No significant correlation was observed with low-β power. The maximum high-β power and the distance to the contact with maximum high-β power were both significant predictors for bradykinesia–rigidity improvement in the multiple regression analysis, explaining 37.4% of the variance altogether. Tremor improvement was not significantly correlated with any frequency.

Conclusion:

High-β oscillations, but not low-β oscillations, recorded from the STN region with the DBS lead can inform stimulation-induced improvement in contralateral bradykinesia–rigidity in patients with PD. High-β oscillations can help refine electrode targeting and inform contact selection for DBS therapy.

Introduction

High-frequency stimulation directed at the subthalamic nucleus (STN) is an effective therapy for patients with Parkinson's disease (PD) (Krack et al., 2003; Weaver et al., 2009; Follett et al., 2010). The selection of the active stimulation contact in postoperative programming is typically based on empirical trials by experienced neurologists (Chen et al., 2003; Volkmann et al., 2006). With the advancements in imaging and connectomics, enhanced therapeutic benefits have been achieved through the selection of the optimal contact location (Vanegas-Arroyave et al., 2016; Akram et al., 2017; Horn et al., 2017b; Dembek et al., 2019; Boutet et al., 2021). However, whether the local field potentials (LFPs) in the STN recorded during the perioperative period predict the outcome of deep brain stimulation (DBS) remains inconclusive (Ray et al., 2008; Zaidel et al., 2010; Boex et al., 2018).

The efficacy of STN DBS may rely on the modulation of the pathological neural activity in the STN area (Zaidel et al., 2010; Accolla et al., 2016; Horn et al., 2017a; Milosevic et al., 2020; Kehnemouyi et al., 2021). Exaggerated β (13–35 Hz) oscillations in the “motor domain” of the cortical–subcortical network have been suggested as the key biomarker of motor impairments in patients with PD (Brown, 2003; Hammond et al., 2007; Eusebio et al., 2009; Little and Brown, 2014). Suppression of such β activity by levodopa or DBS is associated with the improvement of parkinsonism (Brown et al., 2001; Levy et al., 2002; Priori et al., 2004; Weinberger et al., 2006; Ray et al., 2008; Bronte-Stewart et al., 2009; Kühn et al., 2009; Eusebio et al., 2011). Electrophysiological and volume of tissue-activated modeling studies have indicated that the proximity of the chronic stimulation contact to the sensorimotor STN β oscillations predicts a favorable therapeutic outcome with STN DBS (Butson et al., 2007; Horn et al., 2017a).

Previous studies have also suggested that there were different pathophysiology and clinical relevance between low-β (13–20 Hz) and high-β frequency (20–35 Hz) activities in the STN (Priori et al., 2004; Marceglia et al., 2006; Oswal et al., 2016; Godinho et al., 2021; Neuville et al., 2021). However, whether and how the subthalamic β oscillations at different subbands predict the therapeutic outcome of DBS in PD remains unclear.

In the current study, we investigated the association between the therapeutic efficacy of STN DBS and the LFPs recorded by the implanted DBS electrode during the perioperative period. First, different frequency bands were correlated with the stimulation efficacy to evaluate whether there was a frequency-specific association. Second, we explored whether the proximity of the active contact to the maximum β oscillations was also correlated with the stimulation efficacy. Lastly, a multiple regression analysis was performed to evaluate the combined predictive value of the factors that were significantly associated with the stimulation efficacy.

Materials and methods

Patients and surgery

Twenty-six patients with advanced PD who underwent bilateral STN DBS from 2009 to 2016 at Chang Gung Memorial Hospital (CGMH), Taiwan, were recruited (nine women; age, 61.0 ± 7.6 years; disease duration, 14.4 ± 5.8 years; preoperative baseline OFF-medication Unified Parkinson's Disease Rating Scale [UPDRS] III score, 45.0 ± 15.0). The Ethics Review Board in CGMH approved the study (CGMH-IRB No. 97-2345B), and informed consent was obtained from all patients. The patients' clinical characteristics are summarized in Table 1.

Table 1

No.SexAge (y)Disease duration (y)Pre-op UPDRS III ON/OFF MEDPre-op LEDDMotor subtypeMain disabling symptomsChronic stimulation contactsa
1F601019/371,358ARMotor fluctuationC3, C10
2M711122.5/38.51,836ARMotor fluctuation, DyskinesiaC2, C10
3F55149/48.51,524TMotor fluctuation, DyskinesiaC3, C11
4F70621/42900ARMotor fluctuationC1, C10
5M561617/441,214TMotor fluctuationC1, C11
6F602324.5/27.51,248TMotor fluctuationC3, C11
7M473016/29980ARMotor fluctuation, DyskinesiaC3, C10
8M651834/381,400ARMotor fluctuation, DyskinesiaC2, C9
9M711215/27832TMotor fluctuation, DyskinesiaC3, C11
10F691128/611,364ARMotor fluctuationC2, C11
11M48547/741,836ARTruncal dystoniaC1, C10
12M531320/311,820ARAxial rigidity, DyskinesiaC3, C10
13M652419.5/31.51,696ARMotor fluctuation, FOGC2, C11
14F691223/41850ARMotor fluctuationC2, C11
15M491523/481,530TDyskinesiaC3, C11
16F68734/761,200TMotor fluctuation, DyskinesiaC1, C9
17M701613/291,150ARMotor fluctuation, FOG, VH to DAC3, C11
18M57923/47799ARMotor fluctuationC2, C9
19F541718/34.51,790ARMotor fluctuation, DyskinesiaC2, C8
20M532428/44710ARBiphasic dyskinesiaC0, C10
21M681573/783,502.5ARMotor fluctuation, DyskinesiaC3, C10
22M601737/40940ARAxial rigidity, FOGC1, C9
23M561326/55892TMotor fluctuation, Axial rigidityC3, C11
24F71816/301,100ARMotor fluctuation, FOGC3, C10
25M621427/50948ARMotor fluctuationC3, C10
26M581422/67.51,750ARMotor fluctuationC1, C11

Patient characteristics.

a

Lead contacts C0–3 are on the left; C8–C11 are on the right.

F, female; M, male; T, tremor; AR, akinesia-rigidity; FOG, freezing of gait; VH, visual hallucinations; DA, dopamine agonist; UPDRS, Unified Parkinson's Disease Rating Scale; LEDD, levodopa equivalent daily dose.

The surgical procedures have been described in detail elsewhere (Chen et al., 2010, 2021). Preoperative magnetic resonance imaging (MRI) was done for the 26 patients, with 5 of them having high-resolution 1.5 tesla MRI records. Whole-brain stereotactic unenhanced computed tomography was obtained after the application of the Cosman-Roberts-Wells frame (Integra Radionics, Burlington, MA, USA) with a slice thickness of 1 mm. The images were transferred to the StealthStation S7 navigation system (Medtronic, Minneapolis, MN, USA) and superimposed to define the location corresponding to the STN in the atlas of Schaltenbrand and Wahren (Schaltenbrand and Wahren, 1977). The trajectories were then aimed at the center of STN under direct visualization on T2-weighted axial, coronal, and sagittal MRI.

To aid targeting, microelectrode recording (MER) was performed by an experienced neurophysiologist (C.C.C.) for all patients to ensure that the trajectories passed through the STN and to confirm the entry and exit levels of the STN, as well as the level of substantia nigra reticulata (SNr), based on the recognition of the typical transition of multi-unit firing patterns. The DBS lead (model 3389, Medtronic) was then implanted, with the lower border of contact C0 aimed at the level of the lower border of the STN or SNr, so that the contact C2 or C3 would lie in the sensorimotor domain of the STN. The correct placement of DBS leads in the region of the STN was verified by: (1) effective intraoperative macrostimulation, (2) postoperative T2-weighted MRI compatible with the placement of at least one lead contact in the STN region, and (3) when assessed 6 months postoperatively, a significant improvement in the UPDRS motor score during chronic DBS off medication compared to UPDRS off medication with stimulator switched off (44.3 ± 10.9 %, mean reduction ± SD; p < 10−14 with paired t-test).

LFP recordings and signal processing

Resting-state LFPs were recorded within 5 days postoperatively, before the externalized leads were connected to the DBS pulse generator. Unipolar LFPs were recorded for ~200 s (213.0 ± 20.7 s) at a sampling rate of 2048 Hz via the TMSi-Porti amplifier (Twente Medical Systems International, Oldenzaal, Netherlands), using a common average reference as the built-in montage. LFPs were loaded onto a computer using custom software. For offline analysis, the LFP data were converted to an analyzable format using Spike2 software (Cambridge Electronics Design, Cambridge, UK).

The unipolar LFPs were analyzed using MATLAB 2019b (MathWorks, Natick, MA, Massachusetts). All LFPs were visually examined, and segments with artifacts were rejected. Two STNs (the right STN of patient 11 and the left STN of patient 25) were excluded due to poor recording quality, and the remaining 50 STNs were selected for further analysis. Locally weighted scatterplot smoothing with a span of 512 data points (0.25 s) was used for detrending. Spectral powers in different frequency bands (θ, 4–7 Hz; α, 7–13 Hz; low-β, 13–20 Hz; high-β, 20–35 Hz; and low-γ, 40–55 Hz) were calculated using a discrete Fourier transform. All powers were then normalized to the total power of 5–55 and 65–95 Hz frequency bands, avoiding contamination by movement artifacts (<5 Hz) and the power mains artifact (60 Hz in Taiwan). The maximum normalized power of each frequency band among four contacts was selected for further analysis of the correlation with the stimulation outcomes and with the clinical variables.

Clinical evaluation

Initial programming of the stimulation parameters was done at 1 month postoperatively. During the initial programming session, a systematic evaluation of the contacts was performed. The stimulation was usually set at a frequency of 130 Hz and a pulse width of 60 μs, and the voltage was progressively increased at each contact to evaluate the stimulation effect on contralateral rigidity, bradykinesia, and tremor, as well as side effects. The contact that achieved the greatest improvement of motor symptoms with the least side effects was selected for chronic stimulation. The programming process was blind to any electrophysiological results. The stimulation intensity and dopaminergic medicine were adjusted gradually to avoid the occurrence of dyskinesia during the follow-up stage.

For clinical assessment, UPDRS-III ON/OFF stimulation after overnight withdrawal of dopaminergic medication was evaluated at about 1 year (0.8 ± 0.3 years, mean ± SD), with all subitems registered. The UPDRS hemibody scores were separated into the tremor (items 20 and 21), the bradykinesia–rigidity (items 22–26), and the axial (items 27–31) scores. The stimulation efficacy used in the main analysis was defined as follows:

Statistical analysis

Statistical analysis was performed using SPSS for Windows version 25 (IBM, Chicago, IL, USA). Due to the presence of extreme values, nonparametric Spearman's correlation was used for assessing the correlations between normalized power and stimulation efficacy, as well as for the correlation between distance and stimulation efficacy. Paired t-test was used to compare the maximum power locations in different frequency bands relative to the stimulation contacts. Spearman's correlation was used to explore the potential relationship between the baseline clinical variables (age, sex, disease duration, presence of levodopa-induced dyskinesia, daily levodopa equivalent dose, UPDRS-III during ON/OFF medication, and levodopa response) and the stimulation efficacy or the β powers; the Mann–Whitney U-test was used if dichotomized groups were compared (sex and presence of levodopa-induced dyskinesia). Lastly, a multiple linear regression model was computed to estimate the combined effects of the variables that were significantly correlated with the stimulation efficacy. Significance was indicated by p < 0.05.

Results

Figure 1A illustrates an example of raw unipolar LFPs and power spectra from four contacts of the left STN in patient 14. The power spectra indicate a clear increase in power in the β frequency range (13–35 Hz). In this example of STN, two subpeaks centered at ~13 Hz and 24 Hz are evident. Figure 1B shows the reconstruction of the lead locations in five patients, and Figure 1C displays all power spectra from every contact of all patients. Figures 1D,E show the PSDs and the peak frequency distribution of the contacts with the maximum broad-band β power, respectively.

Figure 1

Figure 1

Example of DBS lead location and LFP. (A) Left: Lead-DBS reconstruction of an example lead located on the left STN of patient 14. Middle: The common average referenced unipolar LFPs from the same macroelectrode. Right: The normalized power spectra derived from the LFPs. Spectral peaks were observed at ~13 Hz and 24 Hz. The power spectra were generated using Welch's method and normalized so that the Y scale represents the spectral power relative to the total power (5–55 Hz and 65–95 Hz) in each 1-Hz bin. (B) Lead-DBS reconstruction of five patients' DBS leads, revealing that the electrodes were appropriately located at the sensorimotor STN. (C) Power spectra from every contact of all analyzed electrodes (N = 50 sides). Each thin line is an individual power spectrum. The thick solid line in each plot represents the averaged PSD, and the colored shadow indicates the 95% confidence interval of mean. (D) Power spectra from the contacts with the maximum normalized broad-band β power (13–35 Hz). The transparent rectangles in cyan and green mark the low-β (13–20 Hz) and high-β (20–35 Hz) frequencies, respectively. (E) Distribution of the peak frequencies (the local maximum within the 10–35 Hz range) of the power spectra from the contacts with the maximum broad-band β. LFP, local field potential; STN, subthalamic nucleus; SNr, substantia nigra pars reticulata; PSD, power spectral density.

Correlation between DBS efficacy and β oscillations

We explored whether the maximum oscillatory power in different frequency bands might be associated with the improvement of motor impairments due to STN DBS in PD. We correlated the maximum LFP power among four lead contacts at five frequencies and the improvement of motor impairments in response to DBS. The results of the correlation analysis are summarized in Table 2.

Table 2

UPDRS improvementNormalized power
θαlow-βhigh-βlow-γ
(4–7 Hz)(7–13 Hz)(13–20 Hz)(20–35 Hz)(40–55 Hz)
Bradykinesia + rigidity−0.039−0.131−0.0550.549***−0.110
(n = 50)p = 0.786p = 0.366p = 0.704p < 0.0001p = 0.446
Tremor0.1000.069−0.3020.150−0.109
(n = 39)p = 0.544p = 0.896p = 0.061p = 0.362p = 0.507
Axialb0.1280.135−0.2460.0260.289
(n = 26)p = 0.534p = 0.510p = 0.226p = 0.900p = 0.152

Correlations between oscillation frequency and stimulation efficacya.

a

Spearman's correlation coefficients between the postoperative local field potential (LFP) powers and the improvement in contralateral motor symptoms with STN stimulation.

b

Average power over bilateral STN was used for the correlation with the axial scores.

***

p < 0.001.

UPDRS, Unified Parkinson's Disease Rating Scale.

The maximum LFP in the high-β frequency range was positively correlated with the stimulation-induced improvement in contralateral bradykinesia–rigidity (rs = 0.549, p < 0.0001). This correlation was frequency-specific, as no such association was noted between stimulation-related improvement in bradykinesia–rigidity and the maximum power in the theta (rs = −0.039, p = 0.786), alpha (rs = −0.131, p = 0.366), low-β (rs = −0.055, p = 0.704), or low-gamma frequency range (rs = −0.110, p = 0.446). In other words, the higher the high-β LFP power, the greater the improvement in contralateral bradykinesia–rigidity with stimulation (R2 = 0.267 with simple linear regression). However, improvement in the tremor on contralateral limbs was not significantly correlated with the maximum LFP power in any frequency range. The improvement in axial symptoms was not correlated with the LFP power in any frequency.

Correlation between DBS efficacy and the distance between the contact with the maximum β oscillations and that used for chronic stimulation

Figure 2 shows the relationship between the distance to the maximum low-/high-β and the DBS efficacy for bradykinesia–rigidity. The DBS efficacy had a negative correlation with the center-to-center distance between the chronic stimulation contact and the contact with the maximum high-β power (rs = −0.452, p < 0.001; Figure 2A). Stimulation at contact with the maximum high-β power was associated with greater therapeutic efficacy. This effect was frequency-specific because the distance to the contact with the maximum low-β power did not affect the stimulation outcome (rs = −0.067, p = 0.645; Figure 2B).

Figure 2

Figure 2

Association between the DBS efficacy and the distance between the contact with the maximum oscillations and the active contact. (A) Bradykinesia-rigidity improvement was not correlated with distance to max low-β. (B) Bradykinesia-rigidity improvement was negatively correlated with distance to max high-β. The distance was measured center-to-center between contacts. The fitted line was used to mark a trend. B, bradykinesia; R, rigidity; rs, Spearman's correlation coefficient.

The depths of the contacts with the maximum low-β and high-β power were compared. No significant difference was observed between the average depth of maximum low-β and maximum high-β (p = 0.830 with paired t-test).

Results of the multiple regression of the high-β variables on the DBS efficacy

In the aforementioned correlation analyses, only two variables were identified to be significantly correlated with the stimulation efficacy on bradykinesia–rigidity, namely, the maximum high-β power and the distance to the maximum high-β. Thus, the two variables were put into a multiple linear regression as predictors to examine their combined effects on the contralateral bradykinesia–rigidity outcome. However, a moderate correlation between the two factors was found (rs = −0.395, p = 0.005). To examine whether this collinearity is problematic, we performed a partial correlation analysis, which showed that both variables were still significantly correlated with the outcome when the effect of the other factor was partialized out (data not shown), and therefore they were adequately independent to each other in the current model. The regression model is shown in Figure 3. In this regression model, the slope coefficient for the maximum high-β power was 0.425 (95% CI: [0.139, 0.711], p = 0.0045), and the coefficient for the distance effect was −0.026 (95% CI: [−0.044, −0.008], p = 0.0058), which were both significant. The overall regression model was significant [R2 = 0.374, F(2,47) = 14.061, p = 1.6 × 10−5]. Together, the power of the maximum high-β and the distance to the maximum high-β accounted for 37.4% of the variance in the therapeutic outcome in contralateral bradykinesia–rigidity.

Figure 3

Figure 3

The maximum high-β power in the STN and the distance between the contact with maximum high-β oscillation and the active contact are significant predictors for the improvement in bradykinesia and rigidity with DBS. The 3D graph demonstrates the results of the multiple linear regression model. The higher the high-β power and the shorter the distance between the active contact and the contact with the maximum high-β oscillations, the greater the improvement in B + R due to DBS on contralateral limbs. The fitting formula is inset. The overall regression is significant [R2 = 0.374, F(2,47) = 14.061, p =1.6 × 10−5]. The regression plane is displayed. B, bradykinesia; R, rigidity.

Relationship between low-/high-β oscillations and clinical characteristics

In order to identify a possible association between the baseline clinical characteristics and β oscillations, we correlated the patient's age, disease duration, motor scores, and levodopa responsiveness with maximum LFP power in the low-β and high-β frequency ranges. No significant association was found between the power in the two β frequency ranges and any of the clinical features (Table 3). The normalized low-β and high-β power were not correlated with the bradykinesia–rigidity score assessed postoperatively in the stimulation-off and medication-off state (rs = −0.228, p = 0.112, and rs = −0.025, p = 0.864, respectively). The occurrence of levodopa-induced dyskinesia (18 of 26 patients) was not associated with DBS outcome or LFP power in either β band (Table 3).

Table 3

Stim efficacy: B + RMax LB powerMax HB power
rsprsprsp
Pre-operative
Agea−0.0270.8950.2160.3110.0740.731
Sexa,b0.6270.6760.858
Disease durationa0.0550.789−0.1890.3750.1280.551
Dyskinesiaa,b0.3970.3780.974
LEDDa0.1220.5530.441*0.0310.0001.000
OFF UPDRS-IIIa−0.0250.9040.1810.3960.0100.961
ON UPDRS-IIIa−0.0780.7050.2830.181−0.1820.396
Levodopa response: UPDRS-IIIa−0.0470.821−0.1640.4430.2200.302
OFF UPDRS: B + R−0.2280.1040.1800.2100.0260.856
ON UPDRS: B + R−0.1610.2550.2300.109−0.1280.375
Levodopa response: B + R−0.1000.482−0.0870.5500.1190.412
Post-operative
OFF UPDRS: B + R−0.1620.252−0.2280.112−0.0250.864

Association of clinical characteristics with DBS efficacy and β frequency power.

a

Average power over bilateral STN was used for the correlation with person-bound variables: age, sex, duration, dyskinesia, LEDD, and UPDRS-III.

bMann–Whitney U-test was used for binary comparisons.

*

p < 0.05.

LB, low-β; HB, high-β; B, bradykinesia; R, rigidity; UPDRS, Unified Parkinson's Disease Rating Scale; LEDD, levodopa equivalent daily dose; rs, Spearman's correlation coefficient.

Discussion

We demonstrated that the high-β oscillations, as revealed by using subthalamic LFPs, helped in the estimation of the outcome of DBS therapy in patients with PD. Higher high-β power in the STN region and the proximity of the active contact to the maximum high-β oscillations were associated with a greater response to the DBS of the STN. A multiple linear regression model indicated that the spectral power at the high-β frequency band and the distance between the depth of the contact with the maximum high-β power and the active contact accounted for a 37.4% variance of the improvement in the contralateral bradykinesia–rigidity score on stimulation. This effect was frequency-specific, as this relationship was not observed with the low-β frequency band.

Pathophysiology of low-β and high-β oscillations in the STN

The therapeutic outcome of DBS in PD is sought to be related to the modulation of the pathological oscillations in the STN area (Godinho et al., 2006; Accolla et al., 2016; Milosevic et al., 2020; Tamir et al., 2020; Kehnemouyi et al., 2021). Exaggerated synchronization in the β frequency range (13–35 Hz) is considered as an important biomarker of Parkinson's disease and linked to motor impairment (Brown, 2003; Hammond et al., 2007; Eusebio et al., 2009; Little and Brown, 2014). The suppression of β activities was evident with stimulation (Bronte-Stewart et al., 2009; Eusebio et al., 2011) and correlates with motor improvement (Kehnemouyi et al., 2021). Although studies have suggested different pathophysiological mechanisms and clinical relevance for the two β bands (Priori et al., 2004; Marceglia et al., 2006; Oswal et al., 2016; Godinho et al., 2021; Neuville et al., 2021), it remains unclear whether low-β (13–20 Hz) or high-β (20–35 Hz) oscillations recorded from the STN area in patients with PD would correlate with the DBS efficacy differently. Previous reports have suggested that low-β activity is suppressed with levodopa (Priori et al., 2004; Marceglia et al., 2006; Little et al., 2013), correlates with baseline symptoms (Neumann et al., 2016; van Wijk et al., 2016; West et al., 2016), relates to movement slowing (Lofredi et al., 2019), and is considered to be more directly pathological. One recent report demonstrated that low-β oscillation power was higher at the clinically chosen stimulation contacts than the non-active ones, suggesting that the therapeutic outcome of DBS might be through the modulation of the pathological low-β oscillations (Horn et al., 2017a).

However, it has also been proposed that stimulation at the surrounding tracts in the STN region is involved in the improvement of parkinsonian symptoms, particularly the hyperdirect pathway. The high-β oscillations were found to reflect the synchrony between the STN and the cortex and may indicate the hyperdirect pathway (Hirschmann et al., 2011; Litvak et al., 2011; Whitmer et al., 2012; van Wijk et al., 2016; Tinkhauser et al., 2018b). The cortical–STN coupling in the high-β range was reduced by DBS (Whitmer et al., 2012), and selective stimulation at the hyperdirect pathway effectively ameliorated parkinsonian symptoms (Gradinaru et al., 2009; Li et al., 2012; Miocinovic et al., 2018).

Though exaggerated low-β activities in the STN area were considered pathological and relevant to motor impairment in PD (Hammond et al., 2007), recent computational modeling studies suggested that the pathological synchrony in the low-β frequencies was induced by the hyperdirect pathway (Holgado et al., 2010; Pavlides et al., 2015; Oswal et al., 2021), an important functional connection between the motor cortex and the STN. Furthermore, previous work also demonstrated that high-β activities were an electrophysiological biomarker of the hyperdirect pathway (Oswal et al., 2016). Therefore, modulation of the high-β oscillations by DBS is likely to suppress the generation of pathological low-β oscillations and lead to the improvement of motor impairment in PD.

The significant association between the maximum high-β power and the improvement in bradykinesia–rigidity due to DBS in our study was consistent with the previous works using microelectrode recordings (Zaidel et al., 2010) and using LFP with a machine learning approach (Hirschmann et al., 2022). Both results suggested the high-β activities recorded in the STN region predicted the therapeutic benefit of DBS.

The maximum high-β power correlated with the outcome of DBS of the STN

The significant correlation between high-β power at the STN and clinical improvement might have important clinical implications. The sampling of the maximum high-β power was at least partially relevant to the surgical trajectory of DBS lead implantation, given that the amplitude and power of LFP are likely to be affected by the deviation of the electrode from the oscillatory generator even for only a few millimeters (Chen et al., 2006b; Zaidel et al., 2010; Telkes et al., 2016). Therefore, accurate placement of DBS electrode in the STN region with greater high-β power may lead to a favorable outcome for DBS. However, high-β oscillations may also be related to the pathophysiology of PD. In particular, patients' demographic features, such as age, disease duration, the occurrence of dyskinesia, baseline motor impairment, and preoperative levodopa responsiveness, have been reported to be associated with DBS outcomes for PD (Charles et al., 2002; Welter et al., 2002; Jaggi et al., 2004; Pahwa et al., 2005; Kleiner-Fisman et al., 2006). We accordingly tested the relationship between these clinical variables and the oscillatory low-β and high-β powers, but no significant association was found.

The correlation between high-β power and improvement in bradykinesia–rigidity might be driven by the correlation between high-β and baseline motor impairment because patients with severe motor impairments tended to benefit from DBS more (Chen et al., 2006a; Schuepbach et al., 2019). We tested this hypothesis by correlating the motor scores off-stimulation with the two β subbands. However, no significant correlation was observed. This result is consistent with that of a recent study demonstrating that high-β oscillations in the STN particularly predicted the reduction of bradykinesia–rigidity due to stimulation but were not relevant to the patient's motor symptoms off-stimulation (Hirschmann et al., 2022).

However, it is important to stress that the possible relevance between the subband β activities and the clinical features still could not be excluded entirely. Studies using different approaches to analyze LFPs in the STN region have shown a direct correlation between the parkinsonian symptoms and the β stability index (Little et al., 2012), the β burst duration (Neuville et al., 2021), or the spatially extended coherence (Pogosyan et al., 2010) in the high-β range. It remains to be explored whether other more sophisticated approaches, which take into account the nonlinear relationship of LFP signals, can further identify the association between the β subbands activities and the clinical features of PD.

On the contrary to our data that demonstrated a significant positive correlation between high-β power in the STN and improvement in bradykinesia–rigidity due to DBS, a previous study reported a negative correlation between the β band oscillations and the therapeutic outcome of DBS (Ray et al., 2008). These divergent results might be explained by the different approaches employed in the two studies. Ray and colleagues measured the chronic DBS efficacy by contrasting the preoperative and postoperative motor scores. In contrast, the improvement in motor symptoms in our study was assessed when stimulation was switched on and off after surgery. Therefore, the DBS efficacy in our study could be more directly attributed to the stimulation itself and the potential confounding factors, such as disease progression or interaction with medication adjustment, could be avoided.

Another factor that may affect the direction of correlation is the “stun effect” (Chen et al., 2006b) in which pathological signals in the STN is temporarily disrupted during the penetration of the DBS electrode, and patients' symptoms improved even without stimulation. It has been argued that the better-localized electrodes might result in more disruption of β oscillations and ensure the benefits of STN stimulation. In our group, we routinely descended the DBS electrode slowly in steps of 1 mm during the electrode implantation to minimize the stun effect, allowing more residual pathological oscillatory activity to be recorded. Ultimately, by recording STN LFP from patients with chronically implanted bidirectional devices that allow sustained LFP data retrieval after surgery (Cummins et al., 2021), it might be possible to disentangle the confound of a stun effect on the correlation between the β power in the STN and the improvement in bradykinesia–rigidity.

The association between LFP power and clinical outcome was not only frequency-specific but also symptom-specific. No significant correlation was observed between the improvements in tremor and the low-β or high-β power in the present study. This is in line with the previous evidence showing that β oscillations were only correlated with bradykinesia–rigidity, but not tremor (Ray et al., 2008; Kühn et al., 2009; Zaidel et al., 2010). These observations imply that the pathology of tremors might involve different networks (Helmich et al., 2012). However, we did see a trend, although not significant, toward a negative correlation between the low-β power and the improvement in tremor. Recent studies have also shown that tremor is associated with the reduction of β oscillations (Qasim et al., 2016; Asch et al., 2020). Further study that takes into account the relationships between the oscillatory activities in the β and other frequency bands might provide more insights into the efficacy of DBS on tremors. The lack of association between β oscillations and improvement in tremor may also be explained by the ceiling effect of DBS efficacy for tremor, as the tremor was completely suppressed by stimulation in more than half of our patient cohort (20 out of 39 sides).

DBS efficacy was correlated with the distance between the contact with the maximum high-β and that used for chronic stimulation

In addition to the accuracy of electrode implantation, the selection of the active contact may affect the therapeutic outcome. Previous studies suggest that stimulation at the DBS contact closest to the site with the maximum β power in the STN produces the greatest improvement in parkinsonian symptoms (Ince et al., 2010; Yoshida et al., 2010; Tinkhauser et al., 2018a; Milosevic et al., 2020). In the current study, we demonstrated that the distance between the depth of active contact and the maximum high-β was negatively correlated with the improvement in bradykinesia–rigidity in response to DBS. The greatest improvement was observed when the depth of active contact coincided with that of maximum high-β power. These results may help to identify the optimal therapeutic target in the STN area for DBS. Note that though this correlation was specific to the high-β frequency and not seen in low-β frequency, no significant difference was found between the depths of maximum low-β and high-β in the present study. It remains unclear whether the spatial distribution of low-β and high-β activities is identical in the STN area or not. One report suggested that high-β oscillations lie dorsally to low-β oscillations (Miyagi et al., 2009), but this has not been confirmed by later studies using multi-unit recording or the probabilistic LFP mapping method (Zaidel et al., 2010; Horn et al., 2017a). This ambiguity may be circumvented by intraoperative LFP recordings, which allow neuronal activities in the STN area to be recorded in 1-mm or smaller descending steps using the DBS electrode during electrode implantation; as such, the spatial resolution will be higher than the fixed bipolar recordings (Chen et al., 2006b).

The utility of an LFP-based approach to guide DBS programming is supported by a recent study combining multiple LFP spectral features to inform the selection of stimulation contact in PD (Shah et al., 2022). Compared with using β-range oscillations as the single feature, adding LFP features recorded from resting or movement states improved the accuracy of predicting the most satisfactory stimulating contact. This algorithm was particularly helpful for programming with multicontact directional DBS leads. In contrast to our results indicating that high-β activities in the STN were a good indicator for the outcome of DBS, Shah et al. reported that low-β activities correlated with clinical efficacy. However, they limited the clinical assessment only to upper-limb rigidity, whereas we examined tremor, axial symptoms, and combined rigidity and bradykinesia in the upper and lower limbs. STN activities in high-β were more connected to the motor cortex and considerably disrupted motor function. Modulation of high-β oscillations is likely to be related to different symptoms, such as bradykinesia–rigidity in four limbs (Little et al., 2012) or bradykinesia in the lower limbs (Tinkhauser et al., 2019). We determined the correlation of electrophysiological features with various parkinsonian symptoms, and the results might have crucial clinical applications to inform more personalized treatment strategies for patients with PD.

High-β power and proximity of the contact for stimulation to the contact with maximum high-β inform the DBS efficacy

We identified at least two nonredundant, if not complementary, factors that were associated with the outcome of DBS therapy: the maximum high-β power and the distance between the contact with chronic stimulation and that of the maximum high-β activities. We combined these factors in a multiple regression model, as this method has been used in several electrophysiological studies to improve the predictive value (Ozkurt et al., 2011; Little et al., 2012; Kehnemouyi et al., 2021). We demonstrated that the combination of these factors predicted ~37.4% of the variation of improvement in the contralateral bradykinesia–rigidity scores in response to DBS. These results might inform both the therapeutic potential of an implanted electrode and the optimal active contact for chronic stimulation. This model may be used to predict more accurately the improvements in bradykinesia–rigidity of individual patients in response to DBS in the STN. In addition, our findings are also applicable to the development of closed-loop DBS. In contrast to the conventional DBS that suppressed broad spectral activities in the STN, the closed-loop DBS is more selective in its suppression of LFP. It remains to be seen whether the closed-loop DBS that is selectively triggered by the high-β oscillations is more effective than that triggered by the broad-band β oscillations.

The results of our hypothesis-driven study are consistent with those of a recent study that used machine learning techniques to analyze multiple features of LFP. Hirschmann et al. demonstrated that local high-β power was among the most crucial local features to predict DBS outcomes (Hirschmann et al., 2022). The prediction was not driven by symptom severity. Despite differences in methods and scope, the two studies reached the same conclusion that high-β power in the STN estimated improvement in bradykinesia–rigidity.

In contrast to their finding that the prediction of the DBS outcome was not due to the distance to the anatomical “sweet spot” in their study, we observed disparities in the depth selected for chronic stimulation and that the depth of the maximum high-β, which were purely based on electrophysiology rather than on imaging reconstruction, had a significant effect on DBS efficacy. Therefore, the present study indicated that high-β oscillations recorded from the DBS electrode can help evaluate the optimal stimulation contact for use in chronic stimulation.

Limitations

There were several limitations in our study. First, one limitation inherent in this study is that the postoperative MR reconstruction was only available in five patients. It was not possible to compare the predicting value of electrophysiology to anatomy. Future work that correlated distance between anatomical sweet spot deriving from electrode reconstruction and stimulating contact with the clinical outcome will help to answer this question. This limitation is due to the fact that preoperative high-resolution MRI was only available in five of the patients. High-resolution MRI and tissue activation modeling approaches have suggested that stimulating the dorsolateral, sensorimotor STN produced the best outcome of DBS (Dembek et al., 2019; Kehnemouyi et al., 2021). Though the postoperative MRI attested to the accuracy of the electrode placement in our study, the locations of the contacts were presumptive and have not been confirmed by postoperative MR reconstruction. We elected to use an electrophysiological approach based on the DBS macroelectrode LFP recordings to test how well the high-β oscillations could predict the DBS efficacy, which proved to be valid.

Second, the present study only identified high-β oscillations as a single feature to estimate the therapeutic outcome. Combining multiple electrophysiological features and imaging markers resulted in higher prediction accuracy than that obtained by using beta activity alone (Hirschmann et al., 2022; Shah et al., 2022). However, the use of the dynamic high-β oscillations as a single feature is more intuitive and computation saving and can be a rational target for closed-loop DBS. The real-time rapid estimation of the LFP amplitude to deliver stimulation and the reduction of consumed electric energy are crucial for the closed-loop DBS regime.

Third, normalization by total power can be misleading if a strong tremor peak at 5 Hz exists. However, only 7 of the 26 patients in this study had tremor-dominant PD. A tremor in patients often subsided during LFP recording due to the stun effect occurring a few days after electrode implantation. Moreover, if a strong tremor peak causes a medium beta peak to appear small in the normalized version, it would lead to the underestimation of significance between beta oscillations and clinical improvement.

It should also be noted that our results came from a dataset with a limited number of patients recorded in the same center. Whether the result is generalizable requires further investigation.

Conclusion

The high-β oscillations were associated with the bradykinesia–rigidity improvement from subthalamic DBS for PD patients. Both the maximum high-β activity recorded by the DBS macroelectrode and the proximity of the contact of stimulation to the contact with the max high-β were significantly correlated with stimulation efficacy. Our findings are important in informing the electrode implantation and selecting the optimal stimulation contact for chronic DBS. Combing these two factors may provide a more accurate estimation of a patient's response to DBS therapy.

Funding

This work was supported by the Ministry of Science and Technology, Taiwan (MOST108-2314-B-182-014-MY3 and MOST111-2321-B-A49-002), National Health Research Institutes, Taiwan (NHRI-EX111-11104NI), and the Chang Gung Memorial Hospital, Taiwan (CMRPG3B1432).

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors upon request, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Institutional Review Board and Ethics Committee of Chang Gung Memorial Hospital: CGMH-IRB No. 97-2345B. The patients/participants provided their written informed consent to participate in this study.

Author contributions

C-CC gave the conception and established the team. Y-CC, P-LC, C-SL, and C-CC interviewed the patients and collected the data. P-LC and Y-CC organized the database. P-LC, H-TW, M-CC, and T-CL performed the signal analysis and the statistical analysis. P-HT and M-CY operated on the patients. C-HY managed brain images. P-LC and C-CC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Acknowledgments

The authors are grateful to Prof. Hagai Bergman for his critical review of the manuscript and helpful comments. The authors also thank Prof. Shih-Tseng Lee for operating on some patients. We would like to thank Miss Chih-Ting Lin and Mr. Ming-Kai Yang for assisting with the experiments.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AccollaE. A.Herrojo RuizM.HornA.SchneiderG. H.Schmitz-HubschT.DraganskiB.et al. (2016). Brain networks modulated by subthalamic nucleus deep brain stimulation. Brain139, 25032515. 10.1093/brain/aww182

  • 2

    AkramH.SotiropoulosS. N.JbabdiS.GeorgievD.MahlknechtP.HyamJ.et al. (2017). Subthalamic deep brain stimulation sweet spots and hyperdirect cortical connectivity in Parkinson's disease. Neuroimage158, 332345. 10.1016/j.neuroimage.2017.07.012

  • 3

    AschN.HerschmanY.MaozR.Auerbach-AschC. R.ValskyD.Abu-SninehM.et al. (2020). Independently together: subthalamic theta and beta opposite roles in predicting Parkinson's tremor. Brain Commun. 2, fcaa074. 10.1093/braincomms/fcaa074

  • 4

    BoexC.TyrandR.HorvathJ.FleuryV.SadriS.CorniolaM.et al. (2018). What is the best electrophysiologic marker of the outcome of subthalamic nucleus stimulation in parkinson disease?World Neurosurg. 120, e1217e1224. 10.1016/j.wneu.2018.09.047

  • 5

    BoutetA.GermannJ.GwunD.LohA.EliasG. J. B.NeudorferC.et al. (2021). Sign-specific stimulation 'hot' and 'cold' spots in Parkinson's disease validated with machine learning. Brain Commun. 3, fcab027. 10.1093/braincomms/fcab027

  • 6

    Bronte-StewartH.BarberiniC.KoopM. M.HillB. C.HendersonJ. M.WingeierB.et al. (2009). The STN beta-band profile in Parkinson's disease is stationary and shows prolonged attenuation after deep brain stimulation. Exp. Neurol.215, 2028. 10.1016/j.expNeurol2008.09.008

  • 7

    BrownP. (2003). Oscillatory nature of human basal ganglia activity: relationship to the pathophysiology of Parkinson's disease. Mov. Disord.18, 357363. 10.1002/mds.10358

  • 8

    BrownP.OlivieroA.MazzoneP.InsolaA.TonaliP.Di LazzaroV.et al. (2001). Dopamine dependency of oscillations between subthalamic nucleus and pallidum in Parkinson's disease. J. Neurosci.21, 10331038. 10.1523/JNEUROSCI.21-03-01033.2001

  • 9

    ButsonC. R.CooperS. E.HendersonJ. M.McIntyreC. C. (2007). Patient-specific analysis of the volume of tissue activated during deep brain stimulation. Neuroimage34, 661670. 10.1016/j.neuroimage.2006.09.034

  • 10

    CharlesP. D.Van BlercomN.KrackP.LeeS. L.XieJ.BessonG.et al. (2002). Predictors of effective bilateral subthalamic nucleus stimulation for PD. Neurology59, 932934. 10.1212/WNL.59.6.932

  • 11

    ChenC. C.BruckeC.KempfF.KupschA.LuC. S.LeeS. T.et al. (2006a). Deep brain stimulation of the subthalamic nucleus: a two-edged sword. Curr. Biol.16, R952953. 10.1016/j.cub.2006.10.013

  • 12

    ChenC. C.HsuY. T.ChanH. L.ChiouS. M.TuP. H.LeeS. T.et al. (2010). Complexity of subthalamic 13-35 Hz oscillatory activity directly correlates with clinical impairment in patients with Parkinson's disease. Exp. Neurol.224, 234240. 10.1016/j.expNeurol2010.03.015

  • 13

    ChenC. C.LeeS. T.WuT.ChenC. J.ChenM. C.LuC. S.et al. (2003). Short-term effect of bilateral subthalamic stimulation for advanced Parkinson's disease. Chang. Gung. Med. J.26, 344351.

  • 14

    ChenC. C.PogosyanA.ZrinzoL. U.TischS.LimousinP.AshkanK.et al. (2006b). Intra-operative recordings of local field potentials can help localize the subthalamic nucleus in Parkinson's disease surgery. Exp. Neurol.198, 214221. 10.1016/j.expNeurol2005.11.019

  • 15

    ChenY. C.WuH. T.TuP. H.YehC. H.LiuT. C.YeapM. C.et al. (2021). Theta oscillations at subthalamic region predicts hypomania state after deep brain stimulation in Parkinson's disease. Front Hum Neurosci.15, 797314. 10.3389/fnhum.2021.797314

  • 16

    CumminsD. D.KochanskiR. B.GilronR.SwannN. C.LittleS.HammerL. H.et al. (2021). Chronic sensing of subthalamic local field potentials: comparison of first and second generation implantable bidirectional systems within a single subject. Front. Neurosci.15, 725797. 10.3389/fnins.2021.725797

  • 17

    DembekT. A.RoedigerJ.HornA.RekerP.OehrnC.DafsariH. S.et al. (2019). Probabilistic sweet spots predict motor outcome for deep brain stimulation in Parkinson disease. Annals of Neurology86, 527538. 10.1002/ana.25567

  • 18

    EusebioA.PogosyanA.WangS. Y.AverbeckB.DoyleL.CantiniauxS.et al. (2009). Resonance in subthalamo-cortical circuits in Parkinson's disease. Neurology72, A222A222. 10.1093/brain/awp079

  • 19

    EusebioA.ThevathasanW.Doyle GaynorL.PogosyanA.ByeE.FoltynieT.et al. (2011). Deep brain stimulation can suppress pathological synchronisation in parkinsonian patients. J. Neurol. Neurosurg. Psychiatry82, 569573. 10.1136/jnnp.2010.217489

  • 20

    FollettK. A.WeaverF. M.SternM.HurK.HarrisC. L.LuoP.et al. (2010). Pallidal versus subthalamic deep-brain stimulation for Parkinson's disease. N. Engl. J. Med.362, 20772091. 10.1056/NEJMoa0907083

  • 21

    GodinhoF.Fim NetoA.BianquetiB. L.de LuccasJ. B.VarjaoE.Terzian FilhoP. R.et al. (2021). Spectral characteristics of subthalamic nucleus local field potentials in Parkinson's disease: Phenotype and movement matter. Eur J. Neurosci.53, 28042818. 10.1111/ejn.15103

  • 22

    GodinhoF.ThoboisS.MagninM.GuenotM.PoloG.BenatruI.et al. (2006). Subthalamic nucleus stimulation in Parkinson's disease : anatomical and electrophysiological localization of active contacts. J. Neurol.253, 13471355. 10.1007/s00415-006-0222-z

  • 23

    GradinaruV.MogriM.ThompsonK. R.HendersonJ. M.DeisserothK. (2009). Optical deconstruction of parkinsonian neural circuitry. Science324, 354359. 10.1126/science.1167093

  • 24

    HammondC.BergmanH.BrownP. (2007). Pathological synchronization in Parkinson's disease: networks, models and treatments. Trends Neurosci.30, 357364. 10.1016/j.tins.2007.05.004

  • 25

    HelmichR. C.HallettM.DeuschlG.ToniI.BloemB. R. (2012). Cerebral causes and consequences of parkinsonian resting tremor: a tale of two circuits?Brain135, 32063226. 10.1093/brain/aws023

  • 26

    HirschmannJ.OzkurtT. E.ButzM.HomburgerM.ElbenS.HartmannC. J.et al. (2011). Distinct oscillatory STN-cortical loops revealed by simultaneous MEG and local field potential recordings in patients with Parkinson's disease. Neuroimage55, 11591168. 10.1016/j.neuroimage.2010.11.063

  • 27

    HirschmannJ.SteinaA.VesperJ.FlorinE.SchnitzlerA. (2022). Neuronal oscillations predict deep brain stimulation outcome in Parkinson's disease. Brain Stimul15, 792802. 10.1016/j.brs.2022.05.008

  • 28

    HolgadoA. J.TerryJ. R.BogaczR. (2010). Conditions for the generation of beta oscillations in the subthalamic nucleus-globus pallidus network. J. Neurosci.30, 1234012352. 10.1523/JNEUROSCI.0817-10.2010

  • 29

    HornA.NeumannW-. J.DegenK.SchneiderG-. H.KühnA. A. (2017a). Toward an electrophysiological “sweet spot” for deep brain stimulation in the subthalamic nucleus. Human Brain Mapping38, 33773390. 10.1002/hbm.23594

  • 30

    HornA.ReichM.VorwerkJ.LiN.WenzelG.FangQ.et al. (2017b). Connectivity Predicts deep brain stimulation outcome in Parkinson disease. Ann. Neurol.82, 6778. 10.1002/ana.24974

  • 31

    InceN. F.GupteA.WichmannT.AsheJ.HenryT.BeblerM.et al. (2010). Selection of optimal programming contacts based on local field potential recordings from subthalamic nucleus in patients with Parkinson's disease. Neurosurgery67, 390397. 10.1227/01.NEU.0000372091.64824.63

  • 32

    JaggiJ. L.UmemuraA.HurtigH. I.SiderowfA. D.ColcherA.SternM. B.et al. (2004). Bilateral stimulation of the subthalamic nucleus in Parkinson's disease: surgical efficacy and prediction of outcome. Stereotact. Funct. Neurosurg.82, 104114. 10.1159/000078145

  • 33

    KehnemouyiY. M.WilkinsK. B.AnidiC. M.AndersonR. W.AfzalM. F.Bronte-StewartH. M.et al. (2021). Modulation of beta bursts in subthalamic sensorimotor circuits predicts improvement in bradykinesia. Brain144, 473486. 10.1093/brain/awaa394

  • 34

    Kleiner-FismanG.HerzogJ.FismanD. N.TammaF.LyonsK. E.PahwaR.et al. (2006). Subthalamic nucleus deep brain stimulation: summary and meta-analysis of outcomes. Mov. Disord. 21 Suppl 14, S290304. 10.1002/mds.20962

  • 35

    KrackP.BatirA.Van BlercomN.ChabardesS.FraixV.ArdouinC.et al. (2003). Five-year follow-up of bilateral stimulation of the subthalamic nucleus in advanced Parkinson's disease. N. Engl. J. Med.349, 19251934. 10.1056/NEJMoa035275

  • 36

    KühnA. A.TsuiA.AzizT.RayN.BrückeC.KupschA.et al. (2009). Pathological synchronisation in the subthalamic nucleus of patients with Parkinson's disease relates to both bradykinesia and rigidity. Exp. Neurol.215, 380387. 10.1016/j.expNeurol2008.11.008

  • 37

    LevyR.AshbyP.HutchisonW. D.LangA. E.LozanoA. M.DostrovskyJ. O.et al. (2002). Dependence of subthalamic nucleus oscillations on movement and dopamine in Parkinson's disease. Brain 125(Pt 6), 1196-1209. 10.1093/brain/awf128

  • 38

    LiQ.KeY.ChanD. C.QianZ. M.YungK. K.KoH.et al. (2012). Therapeutic deep brain stimulation in Parkinsonian rats directly influences motor cortex. Neuron76, 10301041. 10.1016/j.neuron.2012.09.032

  • 39

    LittleS.BrownP. (2014). The functional role of beta oscillations in Parkinson's disease. Parkinsonism and Relat. Disord.20, S4448. 10.1016/S1353-8020(13)70013-0

  • 40

    LittleS.PogosyanA.KuhnA. A.BrownP. (2012). β band stability over time correlates with Parkinsonian rigidity and bradykinesia. Exp. Neurol.236, 383388. 10.1016/j.expNeurol2012.04.024

  • 41

    LittleS.TanH.AnzakA.PogosyanA.KuhnA.BrownP.et al. (2013). Bilateral functional connectivity of the basal ganglia in patients with Parkinson's disease and its modulation by dopaminergic treatment. PLoS ONE8, e82762. 10.1371/journal.pone.0082762

  • 42

    LitvakV.JhaA.EusebioA.OostenveldR.FoltynieT.LimousinP.et al. (2011). Resting oscillatory cortico-subthalamic connectivity in patients with Parkinson's disease. Brain134, 359-374. 10.1093/brain/awq332

  • 43

    LofrediR.TanH.NeumannW-. J.YehC-. H.SchneiderG-. H.KühnA. A.et al. (2019). Beta bursts during continuous movements accompany the velocity decrement in Parkinson's disease patients. Neurobiol. Dis.127, 462471. 10.1016/j.nbd.2019.03.013

  • 44

    MarcegliaS.FoffaniG.BianchiA. M.BaselliG.TammaF.EgidiM.et al. (2006). Dopamine-dependent non-linear correlation between subthalamic rhythms in Parkinson's disease. J. Physiol.571, 579-591. 10.1113/jphysiol.2005.100271

  • 45

    MilosevicL.SchererM.CebiI.GuggenbergerR.MachetanzK.NarosG.et al. (2020). Online mapping with the deep brain stimulation lead: a novel targeting tool in Parkinson's disease. Mov. Disord.35, 15741586. 10.1002/mds.28093

  • 46

    MiocinovicS.de HemptinneC.ChenW.IsbaineF.WillieJ. T.OstremJ. L.et al. (2018). Cortical potentials evoked by subthalamic stimulation demonstrate a short latency hyperdirect pathway in humans. J. Neurosci.38, 91299141. 10.1523/JNEUROSCI.1327-18.2018

  • 47

    MiyagiY.OkamotoT.MoriokaT.TobimatsuS.NakanishiY.AiharaK.et al. (2009). Spectral analysis of field potential recordings by deep brain stimulation electrode for localization of subthalamic nucleus in patients with Parkinson's disease. Stereotact. Funct. Neurosurg. 87, 211218. 10.1159/000225974

  • 48

    NeumannW-. J.DegenK.SchneiderG-. H.BrückeC.HueblJ.BrownP.et al. (2016). Subthalamic synchronized oscillatory activity correlates with motor impairment in patients with Parkinson's disease. Mov. Disord.31, 17481751. 10.1002/mds.26759

  • 49

    NeuvilleR. S.PetrucciM. N.WilkinsK. B.AndersonR. W.HoffmanS. L.ParkerJ. E.et al. (2021). Differential effects of pathological beta burst dynamics between Parkinson's disease phenotypes across different movements. Front Neurosci.15, 733203. 10.3389/fnins.2021.733203

  • 50

    OswalA.BeudelM.ZrinzoL.LimousinP.HarizM.FoltynieT.et al. (2016). Deep brain stimulation modulates synchrony within spatially and spectrally distinct resting state networks in Parkinson's disease. Brain139, 14821496. 10.1093/brain/aww048

  • 51

    OswalA.CaoC.YehC. H.NeumannW. J.GratwickeJ.AkramH.et al. (2021). Neural signatures of hyperdirect pathway activity in Parkinson's disease. Nat. Commun.12, 5185. 10.1038/s41467-021-25366-0

  • 52

    OzkurtT. E.ButzM.HomburgerM.ElbenS.VesperJ.WojteckiL.et al. (2011). High frequency oscillations in the subthalamic nucleus: a neurophysiological marker of the motor state in Parkinson's disease. Exp. Neurol.229, 324331. 10.1016/j.expNeurol2011.02.015

  • 53

    PahwaR.WilkinsonS. B.OvermanJ.LyonsK. E. (2005). Preoperative clinical predictors of response to bilateral subthalamic stimulation in patients with Parkinson's disease. Stereotact. Funct. Neurosurg.83, 8083. 10.1159/000086866

  • 54

    PavlidesA.HoganS. J.BogaczR. (2015). Computational models describing possible mechanisms for generation of excessive beta oscillations in Parkinson's disease. PLoS Comput Biol11, e1004609. 10.1371/journal.pcbi.1004609

  • 55

    PogosyanA.YoshidaF.ChenC. C.Martinez-TorresI.FoltynieT.LimousinP.et al. (2010). Parkinsonian impairment correlates with spatially extensive subthalamic oscillatory synchronization. Neuroscience171, 245257. 10.1016/j.Neurosci.ence.2010.08.068

  • 56

    PrioriA.FoffaniG.PesentiA.TammaF.BianchiA. M.PellegriniM.et al. (2004). Rhythm-specific pharmacological modulation of subthalamic activity in Parkinson's disease. Exp. Neurol.189, 369379. 10.1016/j.expNeurol2004.06.001

  • 57

    QasimS. E.de HemptinneC.SwannN. C.MiocinovicS.OstremJ. L.StarrP.A. (2016). Electrocorticography reveals beta desynchronization in the basal ganglia-cortical loop during rest tremor in Parkinson's disease. Neurobiol Dis86, 177186. 10.1016/j.nbd.2015.11.023

  • 58

    RayN. J.JenkinsonN.WangS.HollandP.BrittainJ. S.JointC.et al. (2008). Local field potential beta activity in the subthalamic nucleus of patients with Parkinson's disease is associated with improvements in bradykinesia after dopamine and deep brain stimulation. Exp. Neurol.213, 108113. 10.1016/j.expNeurol2008.05.008

  • 59

    SchaltenbrandG.WahrenW. (1977). Atlas of Stereotaxy of the Human Brain.Stuttgart: Thieme.

  • 60

    SchuepbachW. M. M.TonderL.SchnitzlerA.KrackP.RauJ.HartmannA.et al. (2019). Quality of life predicts outcome of deep brain stimulation in early Parkinson disease. Neurology92, e1109e1120. 10.1212/WNL.0000000000007037

  • 61

    ShahA.NguyenT. K.PetermanK.KhawaldehS.DeboveI.ShahS. A.et al. (2022). Combining multimodal biomarkers to guide deep brain stimulation programming in parkinson disease. Neuromodulation23:S1094-7159(22)00038-1. 10.1016/j.neurom.2022.01.017

  • 62

    TamirI.WangD.ChenW.OstremJ. L.StarrP. A.de HemptinneC.et al. (2020). Eight cylindrical contact lead recordings in the subthalamic region localize beta oscillations source to the dorsal STN. Neurobiol. Dis. 146, 105090. 10.1016/j.nbd.2020.105090

  • 63

    TelkesI.Jimenez-ShahedJ.ViswanathanA.AboschA.InceN. F. (2016). Prediction of STN-DBS electrode implantation track in parkinson's disease by using local field potentials. Front. Neurosci.10, 198. 10.3389/fnins.2016.00198

  • 64

    TinkhauserG.PogosyanA.DeboveI.NowackiA.ShahS. A.SeidelK.et al. (2018a). Directional local field potentials: A tool to optimize deep brain stimulation. Mov. Disord.33, 159164. 10.1002/mds.27215

  • 65

    TinkhauserG.ShahS. A.FischerP.PetermanK.DeboveI.NygyuenK.et al. (2019). Electrophysiological differences between upper and lower limb movements in the human subthalamic nucleus. Clin. Neurophysiol.130, 727738. 10.1016/j.clinph.2019.02.011

  • 66

    TinkhauserG.TorrecillosF.DuclosY.TanH.PogosyanA.FischerP.et al. (2018b). Beta burst coupling across the motor circuit in Parkinson's disease. Neurobiol. Dis. 117, 217225. 10.1016/j.nbd.2018.06.007

  • 67

    van WijkB. C.BeudelM.JhaA.OswalA.FoltynieT.HarizM. I.et al. (2016). Subthalamic nucleus phase-amplitude coupling correlates with motor impairment in Parkinson's disease. Clin. Neurophysiol. 127, 20102019. 10.1016/j.clinph.2016.01.015

  • 68

    Vanegas-ArroyaveN.LauroP. M.HuangL.HallettM.HorovitzS. G.ZaghloulK. A.et al. (2016). Tractography patterns of subthalamic nucleus deep brain stimulation. Brain139, 12001210. 10.1093/brain/aww020

  • 69

    VolkmannJ.MoroE.PahwaR. (2006). Basic algorithms for the programming of deep brain stimulation in Parkinson's disease. Mov. Disord. 21 Suppl 14, S284289. 10.1002/mds.20961

  • 70

    WeaverF. M.FollettK.SternM.HurK.HarrisC.MarksW. J.et al. (2009). Bilateral deep brain stimulation vs best medical therapy for patients with advanced Parkinson disease: a randomized controlled trial. JAMA301, 6373. 10.1001/jama.2008.929

  • 71

    WeinbergerM.MahantN.HutchisonW. D.LozanoA. M.MoroE.HodaieM.et al. (2006). Beta oscillatory activity in the subthalamic nucleus and its relation to dopaminergic response in Parkinson's disease. J. Neurophysiol.96, 32483256. 10.1152/jn.00697.2006

  • 72

    WelterM. L.HouetoJ. L.Tezenas du MontcelS.MesnageV.BonnetA.M.PillonB.et al. (2002). Clinical predictive factors of subthalamic stimulation in Parkinson's disease. Brain125, 575583. 10.1093/brain/awf050

  • 73

    WestT.FarmerS.BerthouzeL.JhaA.BeudelM.FoltynieT.et al. (2016). The Parkinsonian subthalamic network: measures of power, linear, and non-linear synchronization and their relationship to L-DOPA treatment and OFF state motor severity. Front. Hum. Neurosci.10, 517. 10.3389/fnhum.2016.00517

  • 74

    WhitmerD.de SolagesC.HillB.YuH.HendersonJ.M.Bronte-StewartH. (2012). High frequency deep brain stimulation attenuates subthalamic and cortical rhythms in Parkinson's disease. Front. Hum. Neurosci.6, 155. 10.3389/fnhum.2012.00155

  • 75

    YoshidaF.Martinez-TorresI.PogosyanA.HollE.PetersenE.ChenC. C.et al. (2010). Value of subthalamic nucleus local field potentials recordings in predicting stimulation parameters for deep brain stimulation in Parkinson's disease. J. Neurol. Neurosurg. Psychiatry81, 885889. 10.1136/jnnp.2009.190918

  • 76

    ZaidelA.SpivakA.GriebB.BergmanH.IsraelZ. (2010). Subthalamic span of beta oscillations predicts deep brain stimulation efficacy for patients with Parkinson's disease. Brain J. Neurol.133, 20072021. 10.1093/brain/awq144

Summary

Keywords

deep brain stimulation, Parkinson's disease, beta oscillations, subthalamic nucleus (STN), stimulation efficacy

Citation

Chen P-L, Chen Y-C, Tu P-H, Liu T-C, Chen M-C, Wu H-T, Yeap M-C, Yeh C-H, Lu C-S and Chen C-C (2022) Subthalamic high-beta oscillation informs the outcome of deep brain stimulation in patients with Parkinson's disease. Front. Hum. Neurosci. 16:958521. doi: 10.3389/fnhum.2022.958521

Received

31 May 2022

Accepted

09 August 2022

Published

08 September 2022

Volume

16 - 2022

Edited by

Vladimir Litvak, University College London, United Kingdom

Reviewed by

Alberto Averna, Bern University Hospital, Switzerland; Jan Hirschmann, Heinrich Heine University of Düsseldorf, Germany; Chunyan Cao, Shanghai Jiao Tong University School of Medicine, China

Updates

Copyright

*Correspondence: Chiung-Chu Chen

This article was submitted to Brain Imaging and Stimulation, a section of the journal Frontiers in Human Neuroscience

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics