Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Phase matters: A role for the subthalamic network during gait

  • Gabriele Arnulfo ,

    Contributed equally to this work with: Gabriele Arnulfo, Nicolò Gabriele Pozzi

    Roles Conceptualization, Formal analysis, Software, Writing – original draft

    Affiliations Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany, Department of Informatics, Bioengineering, Robotics and System Engineering, University of Genoa, Genoa, Italy

  • Nicolò Gabriele Pozzi ,

    Contributed equally to this work with: Gabriele Arnulfo, Nicolò Gabriele Pozzi

    Roles Conceptualization, Data curation, Investigation, Writing – original draft

    Affiliation Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

  • Chiara Palmisano,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft

    Affiliations Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany, Department of Electronics, Information and Bioengineering, MBMC Lab, Politecnico di Milano, Milan, Italy

  • Alice Leporini,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Software

    Affiliation Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

  • Andrea Canessa,

    Roles Conceptualization, Data curation, Formal analysis, Software

    Affiliations Department of Informatics, Bioengineering, Robotics and System Engineering, University of Genoa, Genoa, Italy, Fondazione Europea di Ricerca Biomedica (FERB Onlus), Cernusco s/N (Milan), Italy

  • Joachim Brumberg,

    Roles Conceptualization, Methodology, Software, Writing – original draft

    Affiliation Department of Nuclear Medicine, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

  • Gianni Pezzoli,

    Roles Funding acquisition, Project administration, Supervision, Writing – review & editing

    Affiliation Centro Parkinson ASST G. Pini-CTO, Milan, Italy

  • Cordula Matthies,

    Roles Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation Department of Neurosurgery, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

  • Jens Volkmann,

    Roles Conceptualization, Project administration, Supervision, Writing – review & editing

    Affiliation Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

  • Ioannis Ugo Isaias

    Roles Conceptualization, Investigation, Project administration, Supervision, Validation, Writing – review & editing

    Isaias_I@ukw.de

    Affiliation Department of Neurology, University Hospital and Julius-Maximillian-University, Wuerzburg, Germany

Abstract

The role of the subthalamic nucleus in human locomotion is unclear although relevant, given the troublesome management of gait disturbances with subthalamic deep brain stimulation in patients with Parkinson’s disease. We investigated the subthalamic activity and inter-hemispheric connectivity during walking in eight freely-moving subjects with Parkinson’s disease and bilateral deep brain stimulation. In particular, we compared the subthalamic power spectral densities and coherence, amplitude cross-correlation and phase locking value between resting state, upright standing, and steady forward walking. We observed a phase locking value drop in the β-frequency band (≈13-35Hz) during walking with respect to resting and standing. This modulation was not accompanied by specific changes in subthalamic power spectral densities, which was not related to gait phases or to striatal dopamine loss measured with [123I]N-ω-fluoropropyl-2β-carbomethoxy-3β-(4-iodophenyl)nortropane and single-photon computed tomography. We speculate that the subthalamic inter-hemispheric desynchronization in the β-frequency band reflects the information processing of each body side separately, which may support linear walking. This study also suggests that in some cases (i.e. gait) the brain signal, which could allow feedback-controlled stimulation, might derive from network activity.

Introduction

Human gait is a complex motor behavior and requires a constant and coordinated flow of information across functionally specialized brain areas [1]. The locomotor network can be schematically divided into a spinal and a supraspinal component comprising of different cortical areas (including the supplementary motor area, SMA), the basal ganglia, the subthalamic nucleus (STN), the mesencephalic locomotor region (MLR), and the cerebellum [14]. At a spinal level, a central pattern generator (CPG) can produce the basic motor pattern for stepping (rhythmic motor activity) autonomously, in the absence of supraspinal and sensory inputs[5]. The selection and modulation of gaits (e.g. initiation, termination, velocity, direction, and spatial orientation) are instead performed at the higher levels of the nervous system [6,7].

Locomotion is controlled by the tuned interplay of such neural circuits and can be studied by assessing their oscillatory activity [8,9]. In particular, the phase synchronization between oscillatory neural circuits is hypothesized to facilitate the efficacy of information exchange between different brain areas and could be particularly relevant for an efficient flow of information in large-scale networks, such as the locomotor network [913].

In this framework, the role of the basal ganglia is poorly defined but their involvement in locomotion is clearly suggested by the presence of gait impairment in neurological diseases with altered basal ganglia functioning, such as Parkinson’s disease (PD) [14]. In particular, the STN is an essential node of the supraspinal locomotor network, being connected to the SMA and projecting to the MLR directly and through the basal ganglia output nuclei (i.e. the globus pallidus pars interna [GPi] and the substantia nigra pars reticulate [SNr]) [15]. The STN is active during gait [4], and the modulation of its activity with deep brain stimulation (STN-DBS) impacts gait performances [16].

In PD, the loss of striatal dopaminergic innervation determines an excessive synchronization of subcortical neuronal oscillations to cortical rhythm, thus hindering subcortical motor control processing [1721]. These functional derangements are usually observed in the STN. Local field potentials (LFPs) recorded in unmedicated PD patients showed an excessive synchronization of the STN neuronal oscillations in the β-frequency band (≈13–35 Hz), which extends to connected nuclei (e.g. the GPi and motor cortices) and to the opposite hemisphere [19,2225]. Such synchronization correlates with the severity of parkinsonian symptoms, and diminishes with dopaminergic medications, deep brain stimulation (DBS) and during movement preparation and execution [26].

Published works have focused on the use of β-band activity as a control signal for closed-loop DBS in patients with PD [27,28]. Pursuant to such study, extensive neurophysiological work will be required to describe not only symptom-specific physiomarkers [18,20,22], but especially physiological markers of particular brain states (e.g. gait and sleep) [29]. To this end, under the hypothesis of a direct involvement of the STN during fine-tuning of complex and synchronized bilateral movements [19,30,31], we investigated the subthalamic local neuronal oscillations and inter-hemispheric coupling during gait. We also established whether STN oscillatory changes were related to specific phases of gait, assessed with a kinematic analysis, and influenced by striatal dopaminergic loss, measured with [123I]N-ω-fluoropropyl-2β-carbomethoxy-3β-(4-iodophenyl)nortropane (FP-CIT) and single-photon computed tomography (SPECT).

Methods

Patients and surgery

We enrolled eight subjects with PD implanted at Julius-Maximilians-University-Hospital Wuerzburg between December 2013 and May 2014 with the Activa PC+S® neurostimulation system (Medtronic, PLC). This system allows therapeutic DBS as well as on-demand recording of the LFPs from the implanted electrodes [32,33]. The Activa PC+S® system and the related hard- and software for programming and readout were provided under a request for application agreement with Medtronic, PLC. The company was not involved in study design, patient selection, data analysis, or reporting of the results.

All patients were diagnosed with PD according to the UK Parkinson Disease Brain Bank criteria [34] and met the inclusion criteria for bilateral STN-DBS [35]. None of them had cognitive decline or mood disturbances as evaluated with standardized rating scales (i.e. Parkinson neuropsychometric dementia assessment, Mattis dementia rating scale, Hamilton depression rating scale and the non-motor symptoms scale). For clinical assessment, patients were evaluated within one month before surgery (pre-DBS) with the United Parkinson’s Disease Rating Scales motor score (UPDRS-III) after overnight (>12 h) withdrawal of all dopaminergic medications (meds-off) and upon receiving 1 to 1.5-times (range 200/50-300/75 mg of levodopa/benserazide) the levodopa-equivalent of the morning dose (meds-on) [36]. After surgery (post-DBS), patients were assessed with the UPDRS-III in four conditions: (i) stimulation off for at least two hours (stim-off); (ii) with bilateral STN stimulation (stim-on) [monopolar stimulation (6/7 subjects); mean amplitude: 3.1 mA (SD: 0.76 mA); pulse width 60 μs; mean frequency 138 Hz (SD 31 Hz)]; (iii) meds-on (as before surgery); (iv) meds-on and stim-on. The therapeutic response to DBS or levodopa was expressed as percentage of improvement according to the formula, adapted from Isaias and coll. [37]: (1) where (a) is meds-off UPDRS-III score and (b) is meds-on UPDRS-III score at pre-DBS or meds-off, stim-on UPDRS-III score at the time of the study. In this study, all subjects were tested at least six months after surgery. At the time of the experiment, all patients were on stable dopaminergic treatment (for at least two months) and chronically stimulated (at least one month with unchanged DBS parameters). The Institutional Review Board of the University Hospital Wuerzburg approved the study and all patients gave written informed consent. The applied methods were conducted in accordance with the Declaration of Helsinki.

The surgical procedure has been previously described [38]. In brief, implantation was performed under local anesthesia using Leksell’s Frame (Elekta, Leksell Stereotaxy System, Stockholm, Sweden). The DBS electrode used was model 3389 (Medtronic, PLC) with four platinum–iridium cylindrical contacts of 1.5 mm each and a contact-to-contact separation of 0.5 mm. Contacts 0 and 8 were the lowermost and contacts 3 and 11 the uppermost (E0-3 refers to right- and E8-11 to the left-hemisphere). The intended coordinates for STN were 12 mm lateral, 2 mm posterior, 4 mm ventral to the mid-commissural point and were adjusted according to individual STN delineation on T2-weighted and SWI images (MAGNETOM Trio, SIEMENS Healthcare, Erlangen, Germany) and with intraoperative microelectrode recordings. Micro- and macro-electrode stimulation and intraoperative CT scan also served to confirm targeting. Postoperative scanning (1 mm slice-thickness CT scan fusion with the pre-operative MRI) confirmed the electrode location. The precise localization within the STN of the active contacts used for chronic stimulation was further confirmed by means of SureTuneTM (Medtronic, PLC) [39].

Protocol, set-up and biomechanical data processing

Recordings were performed at the gait laboratory of the Department of Neurology, University Hospital Wuerzburg. Patients were investigated in the morning at least 12 h after their last dose of antiparkinsonian medication and at two hours after pausing the stimulation. Patients performed at least three trials (range 3–7) according to their clinical conditions for an assessment program involving three conditions: (i) resting (i.e. quietly sitting for 60 s with eyes open during random thinking); (ii) upright standing (for 60 s); (iii) over-ground bare-foot walking at self-selected speed.

For the biomechanical analysis, we considered the steps performed at steady state velocity (steady state linear walking) in a calibrated volume of 8 x 2 x 5 m (total number of 227 strides, range: 19–37 strides per patient). Kinematics of body segments was measured using an optoelectronic system (SMART-DX400, BTS, Milano, Italy), which computed the 3D coordinates of 29 spherical retro-reflective markers (15 mm diameter) fixed to anatomical landmarks [40,41]. The marker coordinates were low-pass filtered (cut-off frequency of 7 Hz) and interpolated. Kinematic parameters were automatically extracted by ad hoc algorithms developed in Matlab® ambient (Matlab 2015b, The MathWorks, Inc., Natick, Massachusetts, USA) and then checked by visual inspection. Steady state linear walking was described by spatio-temporal gait parameters: stride length, duration, and velocity (normalized to subject’s height), stance and double-support duration. Temporal parameters (i.e. stance and double-support) were time-normalized as a percentage of the stride duration. For each subject and condition, all variables have been averaged over the trials. With the same experimental setup, we also evaluated 11 healthy controls (9 males, 2 females; median age 60 years, range 50–66 years) matched for age and anthropometric measurements.

Differences in spatio-temporal gait parameters between controls and PD groups were analyzed with Steel-Dwass all pairs test. We used JMP software version 13 (SAS Institute Inc., Cary, North Carolina, USA).

SPECT imaging

SPECT data acquisition and reconstruction has been previously described [42,43]. All patients performed a SPECT with FP-CIT to measure the dopamine reuptake transporter (DAT) density within three months before surgery (in the best medical and DBS treatment state). Scans were started 180 min after injection of 182.4±3.6 MBq of FP-CIT on a dual-headed integrated SPECT/CT system (Symbia T2; Siemens, Erlangen, Germany). SPECT data were spatially normalized onto a FP-CIT MNI-based template and volumes of interest (i.e. left and right caudate nucleus, left and right putamen, left and right striatum) as well as a reference region in the occipital cortex were defined using the volume of interest (VOI) of the automated anatomical labeling method [44]. The non-displaceable binding potential (BPND) for caudate nucleus, putamen and whole striatum (for both hemispheres) was then assessed using average regional uptake values from VOI analysis and the occipital cortex as reference region [45]. Based on molecular imaging data, we identified the hemisphere with less (STN–) or more (STN+) dopaminergic innervation. We used the whole striatum, rather than its motor part (i.e. the putamen), as the boundaries between the putamen and the caudate nucleus are uncertain in SPECT images. Asymmetry indexes (AI) for striatal BPND were calculated using the formula: (2)

Contralateral refers in this case to the striatum opposite to the clinically most affected side. For healthy subjects, we conventionally refer to the right side as ipsilateral [46,47]. Striatal DAT binding measurements of the patient group were compared with normal values of 15 healthy subjects (4 males, 11 females; median age 67 years, range 44–74). The Institutional Review Board of the Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico di Milano approved the investigation, which was conducted in accordance with the Declaration of Helsinki, and all patients gave written informed consent.

Differences in BPND of DAT between controls and PD groups were analyzed with the Steel-Dwass all pairs test. We used JMP software version 13 (SAS Institute Inc., Car, North Carolina, USA).

LFPs recordings and synchronization

We recorded subthalamic LFPs with a single bipolar derivation for each STN, amplified by 2000 and sampled at 422 Hz. The recording contacts were selected accordingly to clinical effect of chronic stimulation: we selected a bipolar montage crossing the chronically active electrode [48].

A limitation of the Activa PC+S® system is the recording of the ECG artefact by the ventral contacts [49]. We removed ECG waves in LFP recordings of three out of 14 STNs by detecting the signal peaks representing QRS complexes. Details about the ECG correction algorithms and relevant discussion can be found in [19].

Synchronicity across STN recordings and kinematic measures has been previously described [19]. In brief, synchronization was achieved by means of a transcutaneous electrical nerve stimulation (TENS) artefact that was introduced on-demand into the acquisition systems. The electrodes for the TENS were placed in the supraclavicular fossa, directly on the cable that connects the impulse generator with the electrodes, and over the burr-hole site. The same TENS artefact fed into the Activa PC+S® and the SMART-DX400. We introduced two TENS artefacts: one at the beginning and one at the end of each recording session, thus accounting for the different nominal and real sampling frequencies of the acquisition systems. We detected the TENS artefact in all recordings and used the sharp drop-off as synchronization instant across modalities. After synchronization, the analyses were further performed with the dataset containing LFPs resampled at 400 Hz, excluding the time windows with the TENS artefacts.

Power spectral analyses

We characterized the STN power modulation of resting, upright standing and walking by means of the multi tapers method (five tapers) [50]. Given the relative novelty of the Activa PC+S® system, we first defined the quality of the acquired signal with respect to the internal noise. Manufacturer datasheet as well as previous publications reported a nominal noise floor (NF) for the beta band of 150nV/rtHz [32,33]. We normalized the PSD dividing by the square of NF and transformed it into dB/Hz (10Log10(PSD/NF^2)). After this conversion, 0dB/Hz corresponds to the minimum threshold for reliability (Fig 1). We then estimated group level variance using the bootstrapping (20 repetitions; resampling with replacement) technique and estimated the confidence intervals between 5th and 95th percentiles of the bootstrap distributions.

thumbnail
Fig 1. Spectral profiles (average) during resting, upright standing, and gait.

Average spectral power of the STN local field potential during resting (blue line), standing (pink line) and gait (green line) for the two hemispheres, with less (–) and more (+) striatal dopamine innervation. Results are corrected for the nominal noise floor level of the device (150nV/rtHz, [32,33]), all values above 0dB/Hz are reliable. Shaded areas represent the confidence intervals (5–95%) of the group mean. The peak at 32 Hz is a known artefact of the Activa PC+S® system tied to clock settings or due to a triggered check of the battery status.

https://doi.org/10.1371/journal.pone.0198691.g001

We then analyzed the STN activity within the gait phases. We used the wavelet transform with 60 Morlet wavelets with 10 cycles/s (from 6 to 90 Hz) [51]. Frequencies below 6 Hz were discarded as the selected time resolution of 10 cycles/s exceeds the stride duration. This allowed exclusion of low-frequency biological artefacts. We calculated the event-related power relative changes (event related synchronization [ERS] and desynchronization [ERD]), normalizing the mean low-β, high-β and γ power by subtracting and dividing for the average of the entire recording (baseline), to reduce the effects of single-trial normalization: (3)

Since stride length varied across subjects, we normalized the ERS and ERD of each stride to a reference stride duration by means of piece-wise polynomial time-warp algorithm, similarly to that proposed by Sadeghi and coll. [52]. Briefly, the algorithm performs a polynomial interpolation of the time axis against a reference period (the stride in our case). We used a third order polynomial interpolation with a 2 s time-window centered at each velocity peak (VP) of the swing foot. This time-window exceeded the duration of the stride and ensured that the border-effect did not influence the analysis window. We analyzed the STN activity (either STN–or STN+) contralateral to each foot VP (VPcontra). A time-warp algorithm was applied for each single stride and stride phases (i.e. single- and double-support phases of the stance, acceleration and deceleration phases of the swing) defined according to the instants: heel contact, toe off and velocity peak of the swing foot.

Network analyses

We investigated the inter-hemispheric connectivity of the two STNs by computing the spectral coherency (Coh), the phase locking value (PLV) [53] and the amplitude cross-correlation (CC). We selected this metric as they estimate the synchronized rhythmic excitability of neural ensembles, which is considered useful for the flow of information between nodes of a functional network as it might strengthen the information transmission [9,53,54]. The Coh represents one of the most widely accepted metric of functional connectivity [9]. The PLV and the CC disentangle the two main components of the Coh, respectively the phase and the amplitude correlation of the investigated signals [53,55].

We divided resting, standing and walking recordings in multiple epochs of 3 s and computed all the following interaction metrics in each epoch and averaged across epochs for each subject separately. For spectral coherence, we used multi-tapers spectral profiles computed over 13 frequency bins (2 up to 90 Hz) defined to minimize the overlaps and with constant relative bandwidth equal to 0.5.

Coherency is defined as: (4) where Sxx, Syy represent the auto-spectrum of signal X and Y, respectively and Sxy represents the cross-spectrum. This measure (i.e. Coh) is dependent from both phase and amplitude, which were further separately measured. We band-pass filtered the (artefact free) LFPs with 13 finite impulse response (FIR) filters with central frequencies 2 up to 90 Hz and with constant relative bandwidth of 0.5. The band-pass filtered data were then Hilbert transformed to extract the analytic signal (5)

The PLV is defined as: (6) where Δφi represents the phase difference between X and Y at the i-th time sample. The PLV describes the phase of the signals, being independent from their amplitude, and it helps understanding the Coh, especially when the signal-to-noise ratio of the recording is questioned [55].

We then computed the correlation of the signals amplitude with the CC analysis. To quantify CC between amplitudes (7) of signal X and Y we used the Pearson correlation coefficient (r) defined as: (8) where (μ) represent the sample mean and (s) the sample standard deviation.

All three connectivity measures range from 0 to 1, meaning “no” or “complete” coherence, phase synchrony or CC, respectively. For these measures, we quantified surrogate noise levels by time-rotation of the original time series. Specifically, under the null-hypothesis of no correlation between signals x and y, we split-and-rotate the time course of y at random onsets and quantified all three connectivity metrics on the time-rotated y against the original signal of x. We repeated this procedure 100 times, assessing the distribution of the null-hypothesis. For each of the three metrics separately, we computed the p-value as the number of elements from the surrogate distribution that exceeded the observed measure. We considered only significant values (p<0.05) for further statistical analyses. As for power modulation analyses, we quantified the group level variance using similar bootstrapping approach. We compared all possible combinations (i.e. rest vs. walking, standing vs. walking, rest vs. standing) for all the computed indices by means of the Wilcoxon test matched pairs. Of relevant note, we were not able to compare inter-hemispheric connectivity by gait phases because of the overlay between the two hemibodies (i.e. the single-support phase of one leg is the swing phase of the other leg, and vice versa).

We instead investigated the correlation between inter-hemispheric connectivity and striatal dopaminergic denervation. To this end, we estimated the overall striatal dopaminergic loss, which was calculated as the mean of the dopaminergic loss of left and right striatum. However, this measure prevented the distinction between hemispheres (as for PSD), only illustrating the overall dopaminergic loss of each subject. To overcome this limitation, we also correlated the PLV with the AI (Eq 2) as a measure for the dopaminergic degeneration imbalance between the two hemispheres.

All methods and analyses were implemented in Matlab (Matlab 2013a, The MathWorks, Inc., Natick, Massachusetts, USA) using ad hoc algorithms, fieldtrip [56] and integrated in Brainstorm [57]. Due to strict privacy laws, the datasets and analyses of this study shall be made available only upon personal request.

Results

Clinical data

Before surgery (pre-DBS), the median UPDRS-III score in meds-off condition was 44 points (range: 26–55) and 11 points (range: 4–24) in meds-on condition. After surgery (post-DBS), the median UPDRS-III score in meds-off and stim-off condition was 42 points (range: 27–65), in stim-on (and meds-off) condition was 15 points (range: 9–31), in meds-on (and stim-off) was 12 points (range: 8–23), and in meds-on and stim-on condition was 10 points (range: 5–16). The median percentage of improvement (see formula in Patients and surgery section) due to dopaminergic medication was 76.79% (range: 42.50–92.73%) and the one related to STN-DBS was 66.56% (range: 52.50–83.64%). Complete demographic and clinical information of each subject are listed in Table 1.

Biomechanical data

In agreement with previous studies, subjects with PD showed shortened stride length and lower stride velocities (both average and maximum velocity) with respect to matched healthy subjects [41,58]. Of note, the duration of the stance and double-support phase (as percentage of the stride) was within the normal range, thus suggesting that despite the disease-related biomechanical alterations, the patients preserved a normal timing of the gait cycle (i.e. no shuffling gait or festination) (Table 2).

SPECT data

As expected, all PD patients showed significantly reduced non-displaceable binding potential (BPND) of DAT in the caudate nucleus and the putamen. The striatal dopaminergic loss exceeded 50% bilaterally in six out of eight subjects (Table 3). The clinically most affected side (higher UPDRS-III scores) always corresponded to the striatum with less nigro-striatal dopaminergic innervation. The most dopamine-depleted hemisphere was the right one in three subjects.

thumbnail
Table 3. Non-displaceable binding potentials of dopamine reuptake transporters.

https://doi.org/10.1371/journal.pone.0198691.t003

Local subthalamic oscillations

We did not find any clear modulation of local STN activity when comparing the power profiles of walking and standing with resting state (Fig 1). The STN activity was also not modulated during specific gait phases (Fig 2). The STN neuronal oscillations were not related to the striatal dopamine content, despite the wide range of DAT binding loss (46–86%) in our patient collective (Fig 3). Still, the STN+ showed an increase of the average β-band power during walking, which did not reach a statistical significance (Fig 1). The lack of a clear pattern was particularly evident when looking at each STN independently. We showed a bilateral reduction of the β-activity during walking in five STNs (of four subjects: wue09 STN–, wue10 STN–, wue07 STN–and STN+, wue03 STN+), whereas five STNs (of four subjects: wue04 STN–and STN+, wue10 STN+, wue6 STN+ and wue11 STN+) showed an increase of the β-power (Fig 3).

thumbnail
Fig 2. Modulation of the spectral power during the gait cycle.

Event related synchronization (ERS) and desynchronization (ERD) in low β- (top) high β- (middle) and γ-frequency band (bottom). Subthalamic power changes of the phases of gait are shown as the average relative change of the whole stride of all subjects. Shaded areas represent the confidence intervals (5–95%) of the group mean. We analyzed the power changes of STN–and STN+ during the gait cycle of the contralateral foot (but they could be also referred to the matched gait phases of the ipsilateral one). Stance is the period during which the foot is on the ground (dark and light orange bars). The stance phase includes a period of bilateral foot contact with the floor (double-support phases [dark orange bars]), and a period of unilateral foot contact (single-support phase [light orange bar]). The swing phase (light green and dark green bars) is the interval in which the foot is lifted from the floor. Thanks to the velocity peak (VP) of the marker placed on the lateral malleolus, we identified an acceleration (light green) a deceleration (dark green) sub-phase of the swing phase. HS = heel strike; TO = toe off; VP = velocity peak; lower case subscript indicates the foot contralateral (contra) or ipsilateral (ipsi) to STN–or STN+.

https://doi.org/10.1371/journal.pone.0198691.g002

thumbnail
Fig 3. Spectral profiles (single subject) during resting, upright standing and gait.

Single subject spectral power of the STN local field potential during resting (blue line), standing (pink line) and gait for the two hemispheres, with less (–) and more (+) striatal dopamine innervation. Axial slices are left-right flipped to match the corresponding STN. The peak at 32 Hz is a known artefact of the Activa PC+S® system tied to clock settings or due to a triggered check of the battery status. SPECT scans (central column) show striatal dopaminergic loss as percentage decline with respect to healthy subjects (calculated from BPND of DAT, Table 3).

https://doi.org/10.1371/journal.pone.0198691.g003

Inter-hemispheric subthalamic coupling

We showed a reduction in the inter-hemispheric subthalamic coupling during walking with respect to the upright standing and resting state condition. In particular, the phase synchronization (PLV) of the two STNs diminished significantly in the β-frequency band (Fig 4). We did not find any correlation between the PLV values and the overall striatal dopaminergic loss (resting: ρ: 0.47, p = 0.24; standing: ρ: 0.19, p = 0.65; walking: ρ: 0.05, p = 0.90) or the AI (resting: ρ: -0.61, p = 0.1; standing: ρ: 0.32, p = 0.44; walking: ρ: 0.35, p = 0.40).

thumbnail
Fig 4. Inter-hemispheric coherence.

Inter-hemispheric coherence (Coh, top panel), phase locking value (PLV, bottom left plot) and amplitude cross-correlation (CC, bottom right plot) during resting state (blue line), upright standing (pink line) and walking (green line). Statistical significance (red bar, paired Wilcoxon test, p<0.05 uncorrected) was reached for the PLV selectively in the β-frequency band between resting state and walking. Shaded areas represent the confidence intervals (5–95%) of the group mean.

https://doi.org/10.1371/journal.pone.0198691.g004

Discussion

In this study, we showed in PD patients an inter-hemispheric subthalamic coupling in the β-frequency band during resting and standing, and a desynchronization during gait (Fig 4).

The framework of human locomotion control starts in the SMA, reaches the basal ganglia for refinement and then the MLR, where the cerebellar inputs converge, to descend to the medullary and pontine reticular formations and the spinal cord [59]. The STN is a cornerstone of this network, receiving direct afferences from the SMA and projecting to both the MLR and the basal ganglia output nuclei (i.e. the GPi and the SNr), which also project to the MLR [15]. The STN modulates the integration of cortical and cerebellar information by activating or inhibiting the MLR via direct glutamatergic projection or basal ganglia GABAergic output [59]. These feed-forward pathways can support the information processing in the MLR and the medial reticular formation (i.e. rhythm generating system)[59], which projects to the CPGs in the spinal cord.

We suggest that subthalamic inter-hemispheric desynchronization may reflect the downstream conveying of locomotor information for each body side separately, which may facilitate the spontaneous alternating activity of the spinal CPGs for steady linear walking. Of relevance, the movement-related modulation of the subthalamic coupling was narrowed to the β-frequency band (Fig 4). This result further supports the relevance of β oscillations in motor control. The functional circuits of the brain are multiple and segregated by frequency, so that the precise tuning to distinct frequencies could provide a way of marking and segregating different motor networks, over and above any anatomical segregation of processing streams [60].

Besides this physiological interpretation, we cannot exclude that the changes in subthalamic dynamics during gait in PD may be compensatory or maladaptive. The lack of correlation between the PLV and the overall striatal dopaminergic loss and the AI argues against this interpretation. However, it is well known that in the untreated state, PD is dominated by pathologically exaggerated synchronization and coherence in the basal ganglia-cortical circuit, which implies a shift from a dynamic system to stability [17,38,6167]. Therefore, the subthalamic inter-hemispheric desynchronization during gait, with respect to resting and standing, might also reflect a need to compensate for an excessive static connectivity (persistent β-coupling) of the motor network.

Overall, these findings are of particular interest considering the evolution of stimulation therapy towards (adaptive) closed-loop systems directly triggered by biological signals [29]. By timing the stimulation bursts to the phase of tremor oscillations, Cagnan and coll. showed that phase-specific stimulation could control tremor severity in subjects with essential tremor and dystonia [68]. In this study, the stimulation was triggered by kinematic measurements of tremor, but the trigger could be easily extended to recordings of the LFP, using the strength of different oscillators to determine the stimulation phase in real-time [68]. During the development of such an applicable closed-loop stimulation, besides identifying reliable markers of a pathological brain activity it would be crucial to characterize the activity patterns of common and frequent motor behaviors (e.g. gait) [29]. Our study represents the first step in this direction.

Much of the interpretation of our findings presently remains speculative, given also the absolute novelty of our recordings. A subthalamic inter-hemispheric network has been consistently described and investigated in subjects with PD [19,2225], but only during resting state or simple unilateral movements. In particular, voluntary movements modulate the subthalamic coupling in the θ- and α-frequency bands, but not in the β-frequency band [19,30]. This led to the hypothesis that the subthalamic coherence in the β-frequency band in PD could directly reflect a striatal dopamine loss [19], but its inconsistent suppression by levodopa challenges this idea [23]. Alternatively, inter-hemispheric synchronization and desynchronization in the β-frequency band may serve selectively complex and bilateral motor control processing [60], such as locomotion.

A secondary finding of our study is the inter-hemispheric subthalamic coupling during resting, which likely reflect polysynaptic inter-hemispheric neural interactions [22]. This result is in line with previous studies in PD [19,2225] and complements the inter-hemispheric cortical coupling [64,69] and cortico-subthalamic coupling [64,70] in defying a long-distance synchronization network (β-network) in PD [22]. A number of anatomical connections may contribute to such a network, including (i) bilateral cortical inputs and synchronization via the corpus callosum [64,6971] of the hyperdirect [72] or the direct pathways [15], (ii) bilateral projections of the striatum to both STNs [73], and (iii) subthalamic connections through the GPi/SNr and the MLR [23].

In contrast to network activity, we did not find a clear modulation of local subthalamic activity (i.e. PSD) during walking or standing with respect to resting state (Fig 1). This (unexpected) result is, however, in line with two previous studies that also failed to show a clear suppression of the β-band power in the STN of 15 PD patients during walking [49,74]. The lack of β-suppression in the STN during walking may be due to the small β-activity recorded in the recruited patients (Fig 3) and to the technical limitations of the recording device, which showed a poor signal-to-noise ratio. However, we have corrected our analysis for the nominal noise level of the device by normalizing each PSD to 150 muV/Hz [32,33] and we were able to capture with the same device a clear subthalamic β-suppression during upper-arm movements in three patients of the current sample (i.e. wue02, wue09, and wue11) [19]. It seems therefore unlikely that the lack of consistent β-suppression relies on device limitations. Instead, it should be considered that in this study we did not assess the subthalamic β-activity during transition phases (i.e. when changing an ongoing motor action, such as at gait initiation or termination), and the subthalamic β-activity signaled the likelihood of a forthcoming motor or cognitive action [7577]. Therefore, the sensorimotor systems in PD could maintain a status quo (lack of β-suppression) during standing and unperturbed walking as steady state motor actions.

The mismatch between local and network activity, and the lack of a correlation with striatal dopamine loss (Fig 3), might explain the difficulties in ameliorating gait disturbances with STN-DBS as well as the poor levodopa responsiveness of some of these symptoms [78].

Despite these negative findings, we would like to report an anecdotal subthalamic β-suppression during the swing phase of the gait cycle (Fig 2, epochs between TO and HS [ipsi- and contralateral]). This β-suppression was limited to the STN, thus possibly reflecting the striatal dopamine influence of movement control [19].

This study has several limitations. Firstly, without histological verification the placement of the electrodes in STN remains presumptive. Secondly, although our analytical techniques are biased against the detection of stochastic, non-oscillatory activity, we cannot exclude having captured such oscillations (that characterize the ventro-medial area of the STN) in our analysis. Still, we verified the electrode location by fusing the pre- and post-operative imaging of each subject and we recorded from the chronically-active electrodes, which were effective in improving motor symptoms (Table 1). Furthermore, none of the patients reported stimulation-induced behavioral or mood changes related to ventral STN stimulation, thus indirectly confirming that we recorded the activity of the dorsolateral part of the STN, the area with predominant oscillatory activity [7982]. Thirdly, the subthalamic coherence can be due to the volume conduction of a synchronous activity from a third source (e.g. the cerebral cortex). Measurements of coherence are indeed susceptible to zero-phase volume conduction effects; however, this should not be the case for inter-hemispheric interactions. Previous studies demonstrated that the extent of the volume conduction effect depends on several issues such as the source density, the orientation, and the conducting media [83]. This led to a wide range of possible distances (from 0.6 to 5 mm) at which volume conduction may occur [83], but all below the distance between the two STNs. The distance across STNs saved our recordings from the common input of the background magnetic field (i.e. volume conduction), but did not prevent the possible interference from a third (unknown) source of oscillations, which might influence the synchronization analysis. Although the STN afferents are mainly ipsilateral, studies in rats have identified a few thalamic and peduncolopontine axons that also projected to the contralateral STN [84,85].

Despite these limitations, it is worth mentioning that we recorded the LFPs months after surgery and with a fully implanted device, thus limiting the influence of the high impedance variability and the microlesioning effect that bias the immediate post-operative recordings [70].

In summary, this study highlights the inter-hemispheric subthalamic desynchronization in the β-frequency band during walking in subjects with PD. This result has a twofold bearing. First, it expands our understanding of locomotor control in subjects with PD. This topic is of particular relevance given the high incidence of gait disturbances in these patients and the related consequences (i.e. falls, fractures, institutionalization, loss of independence, and increased mortality) [86]. Second, this study suggests that in some cases (i.e. gait), the brain signals, which could allow feedback-controlled stimulation techniques, might derive from an inter-hemispheric network activity rather than local oscillations (e.g. single STN beta-band activity) [74].

Acknowledgments

The authors would like to thank Uri Ramirez Pasos, for critically reviewing the manuscript, Dr. Frank Steigerwald and Dr. Martin M. Reich for patient referral. The manuscript was proofread by Deborah Nock (Medical WriteAway, Norwich, UK). The data here presented are available upon request, inquires can be sent to the data access point manager (Annegret Hauck; University Hospital Würzburg, Department of Neurology, Direktionssekretariat; Josef-Schneider-Straße 11, 97080 Würzburg; phone +49 (0)931/201 2375; fax fax +49 (0)931/201 23946; hauck_a3@ukw.de) or to the corresponding author (Ioannis U. Isaias, MD, PhD; University Hospital Würzburg, Department of Neurology; Josef-Schneider-Straße 11, 97080 Würzburg; phone +49 (0)931/201 23605; fax +49 (0)931/201 23946; E-Mail: Isaias_I@ukw.de).

References

  1. 1. Takakusaki K. Functional Neuroanatomy for Posture and Gait Control. J Mov Disord. Korean Movement Disorders Society; 2017;10: 1–17. pmid:28122432
  2. 2. Kiehn O. Decoding the organization of spinal circuits that control locomotion. Nat Rev Neurosci. Nature Publishing Group; 2016;17: 224–238. pmid:26935168
  3. 3. la Fougère C, Zwergal A, Rominger A, Förster S, Fesl G, Dieterich M, et al. Real versus imagined locomotion: A [18F]-FDG PET-fMRI comparison. Neuroimage. 2010; pmid:20034578
  4. 4. Tard C, Delval A, Devos D, Lopes R, Lenfant P, Dujardin K, et al. Brain Metabolic Abnormalities During Gait With Freezing in Parkinson’S Disease. Neuroscience. 2015;307: 281–301. pmid:26341909
  5. 5. Clark DJ. Automaticity of walking: functional significance, mechanisms, measurement and rehabilitation strategies. Front Hum Neurosci. Frontiers Media SA; 2015;9: 246. pmid:25999838
  6. 6. Rossignol S, Dubuc R, Gossard J-P. Dynamic Sensorimotor Interactions in Locomotion. Physiol Rev. 2005;86.
  7. 7. Hooper SL, Büschges A. Neurobiology of motor control: fundamental concepts and new directions. ISBN: 978-1-118-87340-3.
  8. 8. Singer W. Neuronal synchrony: a versatile code for the definition of relations? Neuron. 1999;24: 49–65, 111–25. pmid:10677026
  9. 9. Fries P. Rhythms for Cognition: Communication through Coherence. Neuron. NIH Public Access; 2015;88: 220–35. pmid:26447583
  10. 10. Salinas E, Sejnowski TJ. Correlated neuronal activity and the flow of neural information. Nat Rev Neurosci. 2001;2: 539–550. pmid:11483997
  11. 11. Fries P. A mechanism for cognitive dynamics: Neuronal communication through neuronal coherence. Trends Cogn Sci. 2005;9: 474–480. pmid:16150631
  12. 12. Bastos AM, Vezoli J, Fries P. Communication through coherence with inter-areal delays. Curr Opin Neurobiol. 2015;31: 173–180. pmid:25460074
  13. 13. Cagnan H, Duff EP, Brown P. The relative phases of basal ganglia activities dynamically shape effective connectivity in Parkinson’s disease. Brain. 2015; pmid:25888552
  14. 14. Nutt JG, Marsden CD, Thompson PD. Human walking and higher-level gait disorders, particularly in the elderly. Neurology. 1993;43: 268–79. pmid:8437689
  15. 15. Parent a., Hazrati LN. Functional anatomy of the basal ganglia. II. The place of subthalamic nucleus and external pallidum in basal ganglia circuitry. Brain Res Rev. 1995;20: 128–154. pmid:7711765
  16. 16. Pötter-Nerger M, Volkmann J. Deep brain stimulation for gait and postural symptoms in Parkinson’s disease. Mov Disord. 2013;28: 1609–1615. pmid:24132849
  17. 17. Brown P, Oliviero A, Mazzone P, Insola A, Tonali P, Di Lazzaro V. Dopamine dependency of oscillations between subthalamic nucleus and pallidum in Parkinson’s disease. J Neurosci. 2001;21: 1033–8. pmid:11157088
  18. 18. de Hemptinne C, Swann NC, Ostrem JL, Ryapolova-Webb ES, San Luciano M, Galifianakis NB, et al. Therapeutic deep brain stimulation reduces cortical phase-amplitude coupling in Parkinson’s disease. Nat Neurosci. 2015;18: 779–786. pmid:25867121
  19. 19. Canessa A, Pozzi NG, Arnulfo G, Brumberg J, Reich MM, Pezzoli G, et al. Striatal Dopaminergic Innervation Regulates Subthalamic Beta-Oscillations and Cortical-Subcortical Coupling during Movements: Preliminary Evidence in Subjects with Parkinson’s Disease. Front Hum Neurosci. 2016;10: 611. pmid:27999534
  20. 20. de Hemptinne C, Ryapolova-Webb ES, Air EL, Garcia PA, Miller KJ, Ojemann JG, et al. Exaggerated phase-amplitude coupling in the primary motor cortex in Parkinson disease. Proc Natl Acad Sci U S A. National Academy of Sciences; 2013;110: 4780–5. pmid:23471992
  21. 21. Litvak V, Jha A, Eusebio A, Oostenveld R, Foltynie T, Limousin P, et al. Resting oscillatory cortico-subthalamic connectivity in patients with Parkinson’s disease. Brain. 2011; pmid:21147836
  22. 22. Little S, Tan H, Anzak A, Pogosyan A, Kühn A, Brown P. Bilateral functional connectivity of the basal ganglia in patients with Parkinson’s disease and its modulation by dopaminergic treatment. PLoS One. 2013; pmid:24376574
  23. 23. Hohlefeld FU, Huchzermeyer C, Huebl J, Schneider GH, Br??cke C, Sch??necker T, et al. Interhemispheric functional interactions between the subthalamic nuclei of patients with Parkinson’s disease. Eur J Neurosci. 2014;40: 3273–3283. pmid:25195608
  24. 24. Kato K, Yokochi F, Taniguchi M, Okiyama R, Kawasaki T, Kimura K, et al. Bilateral coherence between motor cortices and subthalamic nuclei in patients with Parkinson’s disease. Clin Neurophysiol. 2015; pmid:25591829
  25. 25. de Solages C, Hill BC, Koop MM, Henderson JM, Bronte-Stewart H. Bilateral symmetry and coherence of subthalamic nuclei beta band activity in Parkinson’s disease. Exp Neurol. Elsevier Inc.; 2010;221: 260–266. pmid:19944098
  26. 26. Brittain JS, Brown P. Oscillations and the basal ganglia: Motor control and beyond. Neuroimage. Elsevier Inc.; 2014;85: 637–647. pmid:23711535
  27. 27. Little S, Pogosyan A, Neal S, Zavala B, Zrinzo L, Hariz M, et al. Adaptive deep brain stimulation in advanced Parkinson disease. Ann Neurol. 2013; 449–457. pmid:23852650
  28. 28. Rosa M, Arlotti M, Marceglia S, Cogiamanian F, Ardolino G, Fonzo A Di, et al. Adaptive deep brain stimulation controls levodopa-induced side effects in Parkinsonian patients. Mov Disord. Wiley-Blackwell; 2017;32: 628–629. pmid:28211585
  29. 29. Kühn AA, Volkmann J. Innovations in deep brain stimulation methodology. Mov Disord. 2017;32: 11–19. pmid:27400763
  30. 30. Darvas F, Hebb AO. Task specific inter-hemispheric coupling in human subthalamic nuclei. Front Hum Neurosci. 2014; pmid:25249965
  31. 31. Kato K, Yokochi F, Iwamuro H, Kawasaki T, Hamada K, Isoo A, et al. Frequency-Specific Synchronization in the Bilateral Subthalamic Nuclei Depending on Voluntary Muscle Contraction and Relaxation in Patients with Parkinson’s Disease. Front Hum Neurosci. 2016; pmid:27064969
  32. 32. Rouse AG, Stanslaski SR, Cong P, Jensen RM, Afshar P, Ullestad D, et al. A chronic generalized bi-directional brain-machine interface. J Neural Eng. IOP Publishing; 2011;8: 36018. pmid:21543839
  33. 33. Stanslaski S, Afshar P, Cong P, Giftakis J, Stypulkowski P, Carlson D, et al. Design and validation of a fully implantable, chronic, closed-loop neuromodulation device with concurrent sensing and stimulation. IEEE Trans Neural Syst Rehabil Eng. 2012;20: 410–21. pmid:22275720
  34. 34. Hughes AJ, Daniel SE, Ben-Shlomo Y, Lees AJ. The accuracy of diagnosis of parkinsonian syndromes in a specialist movement disorder service. Brain. 2002;125: 861–870. pmid:11912118
  35. 35. Pollak P. Deep brain stimulation for Parkinson’s disease—patient selection. Handb Clin Neurol. 2013;116: 97–105. pmid:24112888
  36. 36. Allert N, Volkmann J, Dotse S, Hefter H, Sturm V, Freund HJ. Effects of bilateral pallidal or subthalamic stimulation on gait in advanced Parkinson’s disease. Mov Disord. 2001;16: 1076–1085. pmid:11748738
  37. 37. Isaias IU, Alterman RL, Tagliati M. Outcome predictors of pallidal stimulation in patients with primary dystonia: The role of disease duration. Brain. 2008;131: 1895–1902. pmid:18567622
  38. 38. Steigerwald F, Pötter M, Herzog J, Pinsker M, Kopper F, Mehdorn H, et al. Neuronal activity of the human subthalamic nucleus in the parkinsonian and nonparkinsonian state. J Neurophysiol. 2008;100: 2515–2524. pmid:18701754
  39. 39. Reich MM, Brumberg J, Pozzi NG, Marotta G, Roothans J, Åström M, et al. Progressive gait ataxia following deep brain stimulation for essential tremor: adverse effect or lack of efficacy? Brain. 2016;139: 2948–2956. pmid:27658421
  40. 40. Isaias IU, Volkmann J, Marzegan A, Marotta G, Cavallari P, Pezzoli G, et al. The Influence of Dopaminergic Striatal Innervation on Upper Limb Locomotor Synergies. Sgambato-Faure V, editor. PLoS One. Public Library of Science; 2012;7: e51464. pmid:23236504
  41. 41. Dipaola M, Pavan EE, Cattaneo A, Frazzitta G, Pezzoli G, Cavallari P, et al. Mechanical Energy Recovery during Walking in Patients with Parkinson Disease. Arias-Carrion O, editor. PLoS One. 2016;11: e0156420. pmid:27258183
  42. 42. Isaias IU, Marotta G, Hirano S, Canesi M, Benti R, Righini A, et al. Imaging essential tremor. Mov Disord. 2010;25: 679–686. pmid:20437537
  43. 43. Isaias IU, Moisello C, Marotta G, Schiavella M, Canesi M, Perfetti B, et al. Dopaminergic Striatal Innervation Predicts Interlimb Transfer of a Visuomotor Skill. J Neurosci. 2011;31: 14458–14462. pmid:21994362
  44. 44. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated Anatomical Labeling of Activations in SPM Using a Macroscopic Anatomical Parcellation of the MNI MRI Single-Subject Brain. NeuroImage. 2002. pmid:11771995
  45. 45. Innis RB, Cunningham VJ, Delforge J, Fujita M, Gjedde A, Gunn RN, et al. Consensus nomenclature for in vivo imaging of reversibly binding radioligands. J Cereb Blood Flow Metab. 2007;27: 1533–9. pmid:17519979
  46. 46. Isaias IU, Benti R, Goldwurm S, Zini M, Cilia R, Gerundini P, et al. Striatal dopamine transporter binding in Parkinson’s disease associated with theLRRK2 Gly2019Ser mutation. Mov Disord. 2006;21: 1144–1147. pmid:16671078
  47. 47. Isaias IU, Canesi M, Benti R, Gerundini P, Cilia R, Pezzoli G, et al. Striatal dopamine transporter abnormalities in patients with essential tremor. Nucl Med Commun. 2008;29: 349–353. pmid:18317299
  48. 48. Devos D, Szurhaj W, Reyns N, Labyt E, Houdayer E, Bourriez JL, et al. Predominance of the contralateral movement-related activity in the subthalamo-cortical loop. Clin Neurophysiol. 2006;117: 2315–2327. pmid:16926112
  49. 49. Quinn EJ, Blumenfeld Z, Velisar A, Koop MM, Shreve LA, Trager MH, et al. Beta oscillations in freely moving Parkinson’s subjects are attenuated during deep brain stimulation. Mov Disord. 2015;30: 1750–1758. pmid:26360123
  50. 50. Thomson DJ. Spectrum estimation and harmonic analysis. Proc IEEE. 1982;70: 1055–1096.
  51. 51. Torrence C, Compo GP, Torrence C, Compo GP. A Practical Guide to Wavelet Analysis. Bull Am Meteorol Soc. 1998;79: 61–78.
  52. 52. Sadeghi H, Allard P, Shafie K, Mathieu PA, Sadeghi S, Prince F, et al. Reduction of gait data variability using curve registration. Gait Posture. 2000;12: 257–64. pmid:11154937
  53. 53. Lachaux JP, Rodriguez E, Martinerie J, Varela FJ. Measuring phase synchrony in brain signals. Hum Brain Mapp. 1999;8: 194–208. pmid:10619414
  54. 54. Varela F, Lachaux J-P, Rodriguez E, Martinerie J. The brainweb: Phase synchronization and large-scale integration. Nat Rev Neurosci. Nature Publishing Group; 2001;2: 229–239. pmid:11283746
  55. 55. Bastos AM, Schoffelen J-M. A Tutorial Review of Functional Connectivity Analysis Methods and Their Interpretational Pitfalls. Front Syst Neurosci. Frontiers; 2016;9: 175. pmid:26778976
  56. 56. Oostenveld R, Fries P, Maris E, Schoffelen J-M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci. Hindawi; 2011;2011: 156869. pmid:21253357
  57. 57. Tadel F, Baillet S, Mosher JC, Pantazis D, Leahy RM. Brainstorm: a user-friendly application for MEG/EEG analysis. Comput Intell Neurosci. Hindawi Publishing Corp.; 2011;2011: 879716. pmid:21584256
  58. 58. Morris ME, Huxham F, McGinley J, Dodd K, Iansek R. The biomechanics and motor control of gait in Parkinson disease. Clin Biomech. 2001;16: 459–470.
  59. 59. Takakusaki K. Forebrain control of locomotor behaviors. Brain Res Rev. 2008;57: 192–198. pmid:17764749
  60. 60. Fogelson N, Pogosyan A, Kühn A a., Kupsch A, Van Bruggen G, Speelman H, et al. Reciprocal interactions between oscillatory activities of different frequencies in the subthalamic region of patients with Parkinson’s disease. Eur J Neurosci. 2005;22: 257–266. pmid:16029215
  61. 61. Priori A, Foffani G, Pesenti A, Tamma F, Bianchi AM, Pellegrini M, et al. Rhythm-specific pharmacological modulation of subthalamic activity in Parkinson’s disease. Exp Neurol. 2004;189: 369–79. pmid:15380487
  62. 62. Brown P, Williams D. Basal ganglia local field potential activity: Character and functional significance in the human. Clin Neurophysiol. 2005;116: 2510–2519. pmid:16029963
  63. 63. Kühn A a, Kempf F, Brücke C, Gaynor Doyle L, Martinez-Torres I, Pogosyan A, et al. High-frequency stimulation of the subthalamic nucleus suppresses oscillatory beta activity in patients with Parkinson’s disease in parallel with improvement in motor performance. J Neurosci. 2008;28: 6165–6173. pmid:18550758
  64. 64. Litvak V, Jha A, Eusebio A, Oostenveld R, Foltynie T, Limousin P, et al. Resting oscillatory cortico-subthalamic connectivity in patients with Parkinson’s disease. Brain. 2011;134: 359–74. pmid:21147836
  65. 65. Doyle Gaynor LMF, Kühn a. a., Dileone M, Litvak V, Eusebio a., Pogosyan a., et al. Suppression of beta oscillations in the subthalamic nucleus following cortical stimulation in humans. Eur J Neurosci. 2008;28: 1686–1695. pmid:18657185
  66. 66. Giannicola G, Marceglia S, Rossi L, Mrakic-Sposta S, Rampini P, Tamma F, et al. The effects of levodopa and ongoing deep brain stimulation on subthalamic beta oscillations in Parkinson’s disease. Exp Neurol. Elsevier Inc.; 2010;226: 120–127. pmid:20713047
  67. 67. Quinn EJ, Blumenfeld Z, Velisar A, Miller Koop M, Shreve LA, Trager MH, et al. Beta Oscillations in Freely Moving Parkinson’s Subjects Are Attenuated During Deep Brain Stimulation. pmid:26360123
  68. 68. Cagnan H, Pedrosa D, Little S, Pogosyan A, Cheeran B, Aziz T, et al. Stimulating at the right time: phase-specific deep brain stimulation. Brain. Oxford University Press; 2017;140: 132–145. pmid:28007997
  69. 69. Silberstein P, Pogosyan A, Kühn AA, Hotton G, Tisch S, Kupsch A, et al. Cortico-cortical coupling in Parkinson’s disease and its modulation by therapy. Brain. 2005;128: 1277–1291. pmid:15774503
  70. 70. Lalo E, Thobois S, Sharott A, Polo G, Mertens P, Pogosyan A, et al. Patterns of bidirectional communication between cortex and basal ganglia during movement in patients with Parkinson disease. J Neurosci. 2008;28: 3008–3016. pmid:18354004
  71. 71. Gooijers J, Swinnen SP. Interactions between brain structure and behavior: The corpus callosum and bimanual coordination. Neurosci Biobehav Rev. 2014;43: 1–19. pmid:24661987
  72. 72. Nambu A, Tokuno H, Takada M. Functional significance of the cortico-subthalamo-pallidal “hyperdirect” pathway. Neurosci Res. 2002;43: 111–117. pmid:12067746
  73. 73. Singh A, Mewes K, Gross RE, DeLong MR, Obeso JA, Papa SM. Human striatal recordings reveal abnormal discharge of projection neurons in Parkinson’s disease. Proc Natl Acad Sci U S A. National Academy of Sciences; 2016;113: 9629–34. pmid:27503874
  74. 74. Syrkin-Nikolau J, Koop MM, Prieto T, Anidi C, Afzal MF, Velisar A, et al. Subthalamic neural entropy is a feature of freezing of gait in freely moving people with Parkinson’s disease. Neurobiol Dis. 2017;108: 288–297. pmid:28890315
  75. 75. Oswal a., Litvak V, Sauleau P, Brown P. Beta Reactivity, Prospective Facilitation of Executive Processing, and Its Dependence on Dopaminergic Therapy in Parkinson’s Disease. J Neurosci. 2012;32: 9909–9916. pmid:22815506
  76. 76. Oswal A, Brown P, Litvak V. Synchronized neural oscillations and the pathophysiology of Parkinson’s disease. Curr Opin Neurol. 2013;26: 662–70. pmid:24150222
  77. 77. Engel AK, Fries P. Beta-band oscillations—signalling the status quo? Curr Opin Neurobiol. 2010;20: 156–165. pmid:20359884
  78. 78. Peterson DS, Horak FB. Neural Control of Walking in People with Parkinsonism. Physiology. 2016;31: 95–107. pmid:26889015
  79. 79. Moran a., Bergman H, Israel Z, Bar-Gad I. Subthalamic nucleus functional organization revealed by parkinsonian neuronal oscillations and synchrony. Brain. 2008;131: 3395–3409. pmid:18986993
  80. 80. Zaidel A, Spivak A, Shpigelman L, Bergman H, Israel Z. Delimiting subterritories of the human subthalamic nucleus by means of microelectrode recordings and a hidden Markov model. Mov Disord. 2009;24: 1785–1793. pmid:19533755
  81. 81. Seifried C, Weise L, Hartmann R, Gasser T, Baudrexel S, Szelényi A, et al. Intraoperative microelectrode recording for the delineation of subthalamic nucleus topography in Parkinson’s disease. Brain Stimul. 2012;5: 378–387. pmid:21782543
  82. 82. Moshel S, Shamir RR, Raz A, de Noriega FR, Eitan R, Bergman H, et al. Subthalamic nucleus long-range synchronization-an independent hallmark of human Parkinson’s disease. Front Syst Neurosci. Frontiers; 2013;7: 79. pmid:24312018
  83. 83. Kajikawa Y, Schroeder CE. How Local Is the Local Field Potential? Neuron. 2011;72: 847–858. pmid:22153379
  84. 84. Marani E, Heida T, Lakke EAJF, Usunoff KG. The Subthalamic Nucleus: Part I: Development, Cytology, Topography and Connections. Adv Anat Embryol Cell Biol. 2008;198:1–113, vii. pmid:18727483
  85. 85. Castle M, Aymerich MS, Sanchez-Escobar C, Gonzalo N, Obeso JA, Lanciego JL. Thalamic innervation of the direct and indirect basal ganglia pathways in the rat: Ipsi- and contralateral projections. J Comp Neurol. 2005;483: 143–153. pmid:15678473
  86. 86. Bloem BR, Grimbergen YA, Cramer M, Willemsen M, Zwinderman AH. Prospective assessment of falls in Parkinson’s disease. J Neurol. 2001;248: 950–8. pmid:11757958