Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Neurobiology of Disease

Subthalamic–Cortical Network Reorganization during Parkinson's Tremor

Peter M. Lauro, Shane Lee, Umer Akbar and Wael F. Asaad
Journal of Neuroscience 24 November 2021, 41 (47) 9844-9858; DOI: https://doi.org/10.1523/JNEUROSCI.0854-21.2021
Peter M. Lauro
1Department of Neuroscience, Brown University, Providence, Rhode Island 02912
2Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island 02912
3The Warren Alpert Medical School, Brown University, Providence, Rhode Island 02903
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shane Lee
1Department of Neuroscience, Brown University, Providence, Rhode Island 02912
2Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island 02912
4Norman Prince Neurosciences Institute, Rhode Island Hospital, Providence, Rhode Island 02903
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Shane Lee
Umer Akbar
1Department of Neuroscience, Brown University, Providence, Rhode Island 02912
2Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island 02912
3The Warren Alpert Medical School, Brown University, Providence, Rhode Island 02903
4Norman Prince Neurosciences Institute, Rhode Island Hospital, Providence, Rhode Island 02903
5Department of Neurology, Rhode Island Hospital, Providence, Rhode Island 02903
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wael F. Asaad
1Department of Neuroscience, Brown University, Providence, Rhode Island 02912
2Robert J. and Nancy D. Carney Institute for Brain Science, Brown University, Providence, Rhode Island 02912
3The Warren Alpert Medical School, Brown University, Providence, Rhode Island 02903
4Norman Prince Neurosciences Institute, Rhode Island Hospital, Providence, Rhode Island 02903
6Department of Neurosurgery, Rhode Island Hospital, Providence, Rhode Island 02903
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Wael F. Asaad
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Tremor, a common and often primary symptom of Parkinson's disease, has been modeled with distinct onset and maintenance dynamics. To identify the neurophysiologic correlates of each state, we acquired intraoperative cortical and subthalamic nucleus recordings from 10 patients (9 male, 1 female) performing a naturalistic visual–motor task. From this task, we isolated short epochs of tremor onset and sustained tremor. Comparing these epochs, we found that the subthalamic nucleus was central to tremor onset, as it drove both motor cortical activity and tremor output. Once tremor became sustained, control of tremor shifted to cortex. At the same time, changes in directed functional connectivity across sensorimotor cortex further distinguished the sustained tremor state.

SIGNIFICANCE STATEMENT Tremor is a common symptom of Parkinson's disease (PD). While tremor pathophysiology is thought to involve both basal ganglia and cerebello–thalamic–cortical circuits, it is unknown how these structures functionally interact to produce tremor. In this article, we analyzed intracranial recordings from the subthalamic nucleus and sensorimotor cortex in patients with PD undergoing deep brain stimulation surgery. Using an intraoperative task, we examined tremor in two separate dynamic contexts: when tremor first emerged, and when tremor was sustained. We believe that these findings reconcile several models of Parkinson's tremor, while describing the short-timescale dynamics of subcortical–cortical interactions during tremor for the first time. These findings may describe a framework for developing proactive and responsive neurostimulation models for specifically treating tremor.

  • deep brain stimulation
  • directed functional connectivity
  • electrocorticography
  • Parkinson's disease
  • subthalamic nucleus
  • tremor

Introduction

Tremor, a cardinal symptom of Parkinson's disease (PD), typically manifests as a 4–6 Hz oscillatory movement of the distal limbs during rest or sustained posture (Lance et al., 1963). While often the presenting motor symptom of PD, tremor (and its response to dopamine replacement therapy) is highly variable across patients (Koller, 1984, 1986; Zach et al., 2015; Dirkx et al., 2017, 2019). PD tremor neurophysiology has been described by the “dimmer switch” model where an “on-off” mechanism is separable from a magnitude controller (Helmich et al., 2012). Specifically, functional magnetic resonance imaging (fMRI) BOLD activity from basal ganglia nuclei such as the globus pallidus pars interna (GPi) correlates with the presence or absence of tremor, whereas immediate tremor amplitude better correlates with BOLD signal from structures in cerebello–thalamo–cortical circuits such as motor cortex (Helmich et al., 2011; Helmich, 2018). The GPi and the monosynaptically connected subthalamic nucleus (STN; Albin et al., 1989) are common therapeutic targets for deep brain stimulation (DBS). Indeed, DBS in each nucleus is equally effective in reducing tremor (Wong et al., 2020). However, the precise role of the STN and its interactions with cortex in these tremor dynamics is unknown.

Low-frequency (4–8 Hz) oscillatory bursting has been observed in both the STN and GPi MPTP primate models of PD (Bergman et al., 1994; Raz et al., 2000). This bursting, although present in the absence of tremor, becomes highly synchronized with tremor once it emerges. STN recordings from patients with PD have similarly revealed theta/tremor frequency (3–8 Hz) activity that is coherent with electromyography (EMG) recordings of tremulous limbs (Levy et al., 2000; Reck et al., 2009, 2010). Accordingly, STN tremor frequency oscillations (along with higher-frequency oscillations) have been used to predict clinical measures of tremor (Hirschmann et al., 2016; Telkes et al., 2018; Asch et al., 2020). Further, studies applying STN DBS at tremor frequencies entrained tremor to the phase of the stimulation, consistent with a direct modulatory role of STN on tremor (Cagnan et al., 2014).

At the same time, tremor reorganizes cortical activity. Magnetoencephalography studies of patients with PD identified a broad cortical tremor network comprising “intrinsic” (ventrolateral anterior thalamus, premotor and motor cortex) and “extrinsic” (cerebellum, ventrolateral intermedius (VIM), somatosensory cortex) loops hypothesized to initialize and stabilize tremor, respectively (Volkmann et al., 1996; Timmermann et al., 2003). This corticocortical synchronization at single and double tremor frequencies extends to STN local field potential (LFP) and EMG recordings as well (Hirschmann et al., 2013). Meanwhile, intraoperative studies combining electrocorticography (ECoG) and STN LFP recordings found decreases in alpha (8–13 Hz) and beta (13–30 Hz) coherence during tremor (Qasim et al., 2016). Despite this broad corticocortical synchronization at tremor frequencies, it remains unclear whether these neurophysiological changes are specific to tremor onset or maintenance. In addition, although STN and sensorimotor cortex become coherent during tremor, the manner in which tremor-related activity is coordinated across structures and how different networks of activity may reflect the different stages of tremor production and maintenance are unknown.

Thus, to understand whether there are indeed distinct neurophysiological mechanisms of tremor initiation and maintenance, and to better understand what neurophysiological interactions characterize these states, we recorded local field potential activity from the STN along with ECoG from sensorimotor cortices while subjects with PD engaged in a task that elicited initiation and persistence of tremor. Specifically, we tested whether the STN (like the GPi) drove tremor specifically during onset, while cortical structures drove sustained tremor.

Materials and Methods

Participants.

All patients undergoing routine, awake placement of deep brain stimulating electrodes for intractable, idiopathic PD between November 2015 and September 2017 were invited to participate in this study. Patients with PD were selected and offered the surgery by a multidisciplinary team based solely on clinical criteria, and the choice of the target (STN vs GPi) was made according to each patient's particular circumstance (disease manifestations, cognitive status, and goals; Akbar and Asaad, 2017). In this report, we focused on 10 patients [9 male (M), 1 female (F)] undergoing STN DBS with intraoperative ECoG recordings. Patients were off all anti-parkinsonian medications for at least 12 h in advance of the surgical procedure (Unified Parkinson's Disease Rating Scale (UPDRS) part III score, 48.2 ± 15.6). Four patients were considered tremor dominant, and six patients had average tremor UPDRS III scores >2 in their right hand (Jankovic et al., 1990). Approximately age-matched control subjects [3 M, 11 F (often patients' partners)] also participated in this study (n = 14 subjects); patients were 55.6–78.5 years of age (mean age, 65.2 ± 7.4 years), and control subjects were 48.3–79.2 years of age (mean age, 62.4 ± 10.0 years) at the time of testing (Mann–Whitney U test comparing ages, p > 0.05). Control subjects were required simply to be free of any diagnosed or suspected movement disorder and to have no physical limitation preventing them from seeing the display or manipulating the joystick. There was a strong male bias in the patient population (9 M, 1 F) and a female preponderance in the control population (3 M, 11 F), reflecting weaker overall biases in the prevalence of PD and the clinical use of DBS therapy (Accolla et al., 2007; Hariz et al., 2011; Rumalla et al., 2018). All subjects were right handed. Patients and other subjects agreeing to participate in this study signed informed consent forms, and experimental procedures were undertaken in accordance with an approved Rhode Island Hospital human research protocol (Lifespan Internal Review Board Protocol #263157) and the Declaration of Helsinki.

Surgical procedure.

Microelectrode recordings (MERs) from the region of the STN of awake patients are routinely obtained to map the target area and guide DBS electrode implantation. A single dose of short-acting sedative medication (typically, propofol) was administered before the start of each procedure, at least 60–90 min before MERs. The initial trajectory was determined on high-resolution (typically, 3 T) magnetic resonance (MR) images coregistered with computed tomography (CT) images demonstrating previously implanted skull-anchor fiducial markers (version 3.0; FHC). Localization of the target relied on a combination of direct and indirect targeting, using the visualized STN as well as standard stereotactic coordinates relative to the anterior and posterior commissures. Appropriate trajectories to the target were then selected to avoid critical structures and to maximize the length of intersection with the STN. A 3-D printed stereotactic platform (STarFix Micro-Targeting System, FHC) was then created such that it could be affixed to these anchors, providing a precise trajectory to each target (Konrad et al., 2011). Microdrives were attached to the platform and then loaded with microelectrodes. Recordings were typically conducted along the anterior, center, and posterior trajectories (with respect to the initial MRI-determined trajectory) separated by 2 mm, corresponding to the axis of highest anatomic uncertainty based on the limited visualization of the STN on MRI. Bilateral ECoG strips were placed posteriorly along sensorimotor cortices through the same burr hole used for MER insertion for temporary recordings. MER began ∼10–12 mm above the MRI-estimated target, which was chosen to lie near the inferior margin of the STN, approximately two-thirds of the distance laterally from its medial border. The STN was identified electrophysiologically as a hyperactive region typically first encountered ∼3–6 mm above estimated target (Gross et al., 2006). At variable intervals, when at least one electrode was judged to be within the STN, electrode movement was paused to assess neural activity and determine somatotopic correspondence, as per routine clinical practice. At these times, if patients were willing and able, additional recordings were obtained in conjunction with patient performance of the visual–motor task.

Neurophysiological signals and analysis.

Microelectrode signals were recorded using “NeuroProbe” tungsten electrodes (Alpha Omega). ECoG signals were acquired using 8-contact subdural strips with 10 mm contact-to-contact spacing (Ad-Tech Medical). All signals were acquired at 22–44 kHz and synchronized using Neuro Omega data acquisition systems (Alpha Omega). Microelectrode impedances were typically 400–700 kΩ, while ECoG contact impedances were typically 10–30 kΩ. Patients performed up to four sessions of the task, with microelectrodes positioned at different depths for each session. As microelectrodes were not independently positionable, some signals may have necessarily been acquired outside of the STN. All recorded signals were nevertheless considered and analyzed.

Neural data were analyzed using the “numpy/scipy” Python 3 environment (Harris et al., 2020; Virtanen et al., 2020; https://numpy.org/, https://www.scipy.org/). Offline, ECoG contacts were rereferenced to a common median reference within a strip (Liu et al., 2015). All resulting signals were bandpass filtered between 2 and 600 Hz, and notch filtered at 60 Hz and its harmonics. Time series were z-scored and artifacts above 4 SDs were removed. These resulting time series were then downsampled to 1 kHz. Time series were bandpass filtered using a Morlet wavelet convolution (wave number 7) at 1 Hz intervals, covering 3–400 Hz. The instantaneous power and phase at each frequency was then acquired by the Hilbert transform. To analyze broad frequency bands, we grouped frequencies as follows: theta, 3–8 Hz; alpha, 8–12 Hz; low beta, 12–20 Hz; high beta, 20–30 Hz; low gamma, 30–60 Hz; mid-gamma, 60–100 Hz; high gamma, 100–200 Hz; and hfo; 200–400 Hz. For interregional analyses [phase-locking value (PLV), phase slope index (PSI), and Granger prediction], we focused on frequencies up to 100 Hz; spectral or time series data were subsequently downsampled to 250 Hz.

Anatomical reconstruction of recording sites.

Patients underwent preoperative, intraoperative, and postoperative imaging per routine clinical care. Preoperatively, stereotactic protocol MR images were obtained (Vario 3.0 T scanner, Siemens) that included T1- and T2-weighted sequences [T1, MPRAGE sequence: repetition time (TR), 2530 ms; echo time (TE), 2.85 ms; matrix size, 512 × 512 voxels; in-plane resolution, 0.5 × 0.5 mm2; 224 sagittal slices; 1 mm slice thickness; T2, SPACE sequence: TR, 3200 ms; TE, 409 ms; matrix size, 512 × 512 voxels; in-plane resolution, 0.5 × 0.5 mm2; 224 sagittal slices; 1 mm slice thickness]. Preoperative, intraoperative, and postoperative (in some cases) CT scans were also acquired [Lightspeed VCT Scanner for extraoperative imaging, GE Healthcare (tube voltage, 120 kV; tube current, 186 mA; data acquisition diameter, 320 mm; reconstruction diameter, 250 mm; matrix size, 512 × 512 voxels; in-plane resolution, 0.488 × 0.488 mm2; 267 axial slices; slice thickness, 0.625 mm); AIRO Mobile CT System Scanner for intraoperative scanning, Mobius Imaging (tube voltage, 120 kV; tube current, 240 mA; data acquisition diameter, 1331 mm; reconstruction diameter, 337 mm; matrix size, 512 × 512 voxels; in-plane resolution, 0.658 × 0.658 mm2; 182 axial slices; slice thickness, 1 mm)]. Postoperative MR images [Aera 1.5 T scanner, Siemens (T1, MPRAGE sequence: TR, 2300 ms; TE, 4.3 ms; matrix size, 256 × 256 voxels; in-plane resolution, 1.0 × 1.0 mm2; 183 axial slices; slice thickness, 1 mm; specific absorption rate, <0.1 W/g)] were typically obtained 1–2 d after the operation to confirm proper final electrode location.

To reconstruct recording locations, MR and CT images were coregistered using the WayPoint Navigator software (FHC). The raw DICOM images and the linear transform matrices were exported and applied to reconstructed image volumes using the AFNI command “3dAllineate,” bringing them into a common coordinate space (Cox, 1996; Li et al., 2016). Microelectrode depths were calculated by combining intraoperative recording depth information with electrode reconstructions obtained from postoperative images using methods described previously (Lauro et al., 2016, 2018). To determine the anatomic distribution of microelectrode recording sites across patients, preoperative T1-weighted MR images were registered to a T1-weighted MNI reference volume (MNI152_T1_2009c) using the AFNI command “3dQwarp” (Fonov et al., 2009). The resulting patient-specific transformation was then applied to recording site coordinates. MNI-warped recording coordinates were then assessed for proximity to structures such as the STN as delineated on the MNI PD25 atlas (Xiao et al., 2012, 2015, 2017). ECoG contacts were segmented from intraoperative CT volumes using the same DBStar processing as microelectrodes. Contacts were then projected onto individual cortical surface reconstructions generated from preoperative T1 volumes (Dale et al., 1999; Fischl et al., 2002; Saad and Reynolds, 2012; Trotta et al., 2018). Individual cortical surface reconstructions were coregistered to a standard Desikan–Destrieux surface parcellation (Argall et al., 2006; Desikan et al., 2006; Destrieux et al., 2010). Contacts were labeled and grouped as “premotor cortex” (PMC), “motor cortex” (MC), “somatosensory cortex” (SC), or “parietal cortex” if they contained the following anatomic parcellation labels: PMC, ctx_lh_G_front_sup, ctx_lh_G_front_middle; MC, ctx_lh_G_precentral; SC, ctx_lh_G_postcentral; and posterior parietal cortex (PPC), ctx_lh_G_parietal_sup, ctx_lh_G_pariet_inf-Supramar.

If a contact had more than one label (8 of 80 contacts), they were removed from further analysis.

Experimental design.

We used a visual–motor target tracking task to estimate the degree of motor dysfunction in a continuous fashion. Specifically, while patients with PD reclined on the operating bed in a “lawn-chair” position, a joystick was positioned within their dominant hand, and a boom-mounted display was positioned within their direct line of sight at a distance of ∼1 m. The task was implemented in MonkeyLogic (Asaad and Eskandar, 2008a,b; Asaad et al., 2013) and required subjects to follow a green target circle that moved smoothly around the screen by manipulating the joystick with the goal of keeping the white cursor within the circle (Fig. 1A). The target circle followed one of several possible paths (invisible to the subject), with each trial lasting 10–30 s. Each session consisted of up to 36 trials (∼13 min of tracking data), and subjects performed one to four sessions during the operation. Control subjects performed this task in an extraoperative setting.

Speed quantification.

To calculate movement speed, x- and y-joystick traces were 3 Hz low-pass filtered, and the euclidean change of cursor position was calculated over time. To standardize movement speed within patients, movement speed values within a session were minimum–maximum normalized into a measure of “slowness,” where 0 = highest speed and 1 = lowest speed.

Tremor amplitude quantification.

To calculate tremor, x- and y-joystick traces were 3–8 Hz bandpass filtered, and a one-dimensional linear projection of the filtered traces was calculated. Tremor amplitude and phase were calculated using the Hilbert transform of the resulting one-dimensional time series.

Tremor epoch design.

To standardize tremor amplitude across patients, tremor amplitude values from control subjects and patients were averaged into 4 s contiguous, nonoverlapping epochs. We chose our 4 s window size based in part on fMRI studies (Helmich et al., 2011), which based estimates of tremor amplitude on the timescale of echoplanar imaging TRs, which were 1–2 s. In addition, we calculated the autocorrelation of tremor amplitude within individual patient sessions and found that across all patients the central peak (>3 SDs above the mean) spanned 2 s. To capture the transition from one discrete state (no tremor) to another (sustained tremor), we chose a window size of 4 s to capture both states within one “tremor onset” epoch.

The resulting average and SD of the control tremor amplitude distribution were used to z-transform control subject and PD patient tremor amplitude epochs (Fig. 1B). To determine a cutoff to optimally differentiate control and PD population tremor data, receiver operating characteristic (ROC) tests were performed between supracutoff population data for cutoff values ranging from −2 (the lowest observed in both populations) to 10. The maximum area under the curve (AUC) value was observed for z = 3 (ROC AUC = 0.85), which was used for subsequent analyses.

“No tremor” and “sustained tremor” epochs were identified by 4 s epochs where the average tremor amplitude was subthreshold or suprathreshold. Potential tremor onset epochs were detected by taking the continuous tremor amplitude (1 ms samples) and identifying those time points where tremor amplitude crossed from subthreshold to suprathreshold levels. Epochs were then classified as “onset” if the mean of tremor amplitude samples in the subsequent 2 s were >3 SDs, and if the mean of tremor amplitude samples in the preceding 2 s was <3 SDs.

While there was no within-condition epoch overlap (i.e., each sustained tremor epoch was nonoverlapping), there was slight overlap between no tremor epochs with the pretrigger segment of tremor onset epochs (12 of 575 epochs across all subjects; mean ± SD of overlap, 1.275 ± 0.582 s). For sustained tremor, there were 18 of 171 epochs with some overlap with the post-trigger segment of tremor onset (mean ± SD of overlap, 1.008 ± 0.824 s).

Tremor frequency calculation and UPDRS correlation.

To calculate each patient's dominant tremor frequency (i.e., the frequency with the largest amplitude), a distribution of tremor amplitude was created by aggregating each patient's tremor amplitude epochs. In parallel, a frequency distribution was created by calculating the dominant (highest-power) tremor frequency within each epoch. A patient-specific dominant tremor frequency was then calculated as the frequency containing the highest aggregate tremor amplitude.

Correlations between task-derived tremor amplitude and UPDRS were conducted with subscores pertaining to the upper extremity relevant to the patient's task performance (rest tremor, postural tremor, finger taps, hand opening/closing, rapid alternating movements, rigidity). Each patient UPDRS subscore was Spearman correlated with the median of each patient's tremor amplitude distribution and was assessed for significance using a bootstrap null distribution (1000 iterations) where tremor medians were randomly shuffled with respect to UPDRS subscores.

Tremor/speed-spectral power correlation.

To determine whether spectral power across frequencies correlated with changes in tremor amplitude or slowness, linear mixed models were fit to 4 s epochs of averaged tremor/slowness and spectral magnitude of canonical frequency bands (theta, alpha, low beta, high beta, low gamma, high gamma, hfo). Models were fit within entire task sessions for each band and were specified as follows: Tremor/Slowness∼Powerband+(1|Subject).

Tremor epoch spectral power modulation.

To determine whether the spectral band power at each structure differed by tremor epoch type, linear mixed models were used to compare spectral band power across epoch types by the following model: Powerband∼C(TremorEpochType)+(1|Subject).

Tremor-Neural theta phase-locking value.

To determine whether theta (3–8 Hz) in tremor and neural recordings were synchronized, the PLV was calculated with tremor and neural theta phase per trial (Lachaux et al., 1999). Theta-phase estimates for neural spectral data were calculated by taking the circular/angular mean for narrow-band phase estimates between 3 and 8 Hz at each time point (t), as follows: PLVTremor-Neuralθ=1T|∑t=1Tei(θTremor(t)−θNeural(t))|.(1)

To determine whether tremor-neural theta phase synchrony at each structure differed by tremor epoch type, linear mixed models were used to compare PLV values across epoch types by the following model: PLVTremor-Neuralθ∼C(TremorEpochType)+(1|Subject). All PLV-related analyses were also calculated with the pairwise phase consistency (PPhC) measure to control for differences in number of trials across conditions (Vinck et al., 2010; Aydore et al., 2013). PPhC=NtrialsNtrials−1(PLV2−1Ntrials).(2)

As PLV and PPhC results were qualitatively similar, we reported PLV results.

Tremor–neural theta phase slope index.

To understand the lag–lead relationship between tremor (a bandpassed signal) and neural theta-phase locking, the PSI was calculated for the theta band (3–8 Hz) with 1 Hz frequency resolution (Nolte et al., 2008) using the “spectral_connectivity” Python toolbox (https://github.com/Eden-Kramer-Lab/spectral_connectivity, https://doi.org/10.5281/zenodo.4088934).

As the spectral_connectivity toolbox uses the multitaper transform for spectral analysis, the number of necessary tapers (L) was calculated by first calculating the time–half-bandwidth product (TW) using the desired frequency resolution (Δf; 1 Hz for parity with wavelet spectral analyses) and the time window of the entire trial (N, 4 s; Prerau et al., 2017), as follows: TW=NΔf2.(3)

We subsequently used TW to calculate L using the floor function (⌊⌋), as follows: L=⌊2TW−1⌋.(4)

With our parameters, three Slepian tapers were used for whole-trial single-window PSI estimates, as follows: PSITremor,Neural=ℑ(∑fϵFCTremor,Neural*(f)·CTremor,Neural(f+Δf)).(5)

PSI was then estimated from the imaginary (ℑ) component of the complex coherency (C) between tremor and neural theta, where the complex coherency was calculated as follows from the cross-spectral density matrix (S) between the two signals: CTremor,Neural(f)=STremor,Neural(f)STremor,Tremor(f)·SNeural,Neural(f).(6)

Phase offsets between 1 Hz frequency bands (Δf) within theta (F) were used to calculate the phase slope. Because of our short-timescale windowed application of PSI, we did not normalize values of PSI by their SD (Young et al., 2017). To determine whether tremor or neural recordings exhibited directional theta influence, the empirical PSI was compared with a null distribution of 1000 PSI values generated from shuffling time series of one signal across trials. The p values were calculated empirically from the resulting distribution and corrected for multiple comparisons with the Benjamini–Hochberg method at q = 0.05.

Tremor epoch interregional phase-locking value.

To compare time-varying phase synchrony across structures, the PLV was calculated across each structure pair (j, k) per 1 Hz frequency band (f) from 1 to 100 Hz using wavelet-derived spectral data, as follows: PLVf(t)=1Ntrials|∑n=1Ntrialsei(θj(f,t,n)−θk(f,t,n))|.(7)

To determine whether pairwise frequency band PLVs differed by tremor epoch type, linear mixed models were used to compare PLV values across epoch types by the following model: PLVband∼C(TremorEpochType)+(1|Subject).

Tremor epoch interregional Granger prediction.

To understand whether tremor epoch-related dynamic changes in spectral power or synchrony were driven by dynamic directional influences of one structure onto another, nonparametric spectral Granger prediction (GP) was calculated between each structure pair using the “spectral_connectivity” python toolbox. Specifically, frequency information (1 Hz frequency resolution) for each structure–time series pair was calculated using a single 4000 ms multitaper window (three tapers). From there, a frequency-based estimation of information flow between structures was calculated using a cross-density spectral matrix (Dhamala et al., 2008). Subsequently, frequency-specific (f) GP (i.e., the log-ratio of total frequency power over nonpredicted frequency power) was calculated between structure pairs (j, k) for each epoch type using S, the spectral transfer matrix (H), and the noise covariance matrix (∑), as follows: GPj→k(f)=ln(Skk(f)Skk(f)−(∑jj−∑jk2∑kk)|Hjk(f)|2).(8)

To determine whether one structure exhibited frequency-specific Granger prediction on another, the empirical GP was compared with a null distribution of 1000 GP values generated from shuffling time series of one structure across trials. The p values for each frequency were calculated empirically from the resulting distribution and corrected for multiple comparisons with the Benjamini–Hochberg method at q = 0.05.

To understand how GP varied as a function of time, frequency information for each structure–time series pair were calculated in 2000 ms windows with 100 ms overlap using the multitaper transform for each event trial. To maintain the same number of tapers (three tapers) between static and dynamic GP analyses, frequency resolution was increased to 2 Hz for dynamic GP calculation. To determine whether one structure exhibited time-varying directional influence on another, the empirical GP was compared with a null distribution of 1000 GP values generated from shuffling time series of one structure across trials. The p values for each time and frequency point were calculated empirically from the resulting distribution and corrected for multiple comparisons with the Benjamini–Hochberg method at q = 0.05. Resulting significant time–frequency clusters were additionally filtered by considering only clusters whose area was greater than the 95th percentile of all Benjamini–Hochberg method-corrected significant clusters.

Tremor epoch interregional phase slope index.

To calculate theta directed connectivity across structures, the PSI was used for the theta band (3–8 Hz) with 1 Hz frequency resolution across structures. Frequency information (1 Hz frequency resolution) for each structure–time series pair was calculated in a single 4000 ms window using the multitaper transform (three tapers) for each event trial. To determine whether one structure exhibited PSI influence on another, the empirical PSI was compared with a null distribution of 1000 PSI values generated from shuffling one structure's time series across trials. The p values were calculated empirically from the resulting distribution and corrected for multiple comparisons with the Benjamini–Hochberg method at q = 0.05.

To calculate the time-varying PSI between broad frequency bands, PSI was calculated using a 2000 ms window sliding by 100 ms (three tapers with 2 Hz frequency resolution). A bootstrap was then performed, and empirical p values for each time point were corrected for multiple comparisons with the Benjamini–Hochberg method at q = 0.05.

Statistical analysis.

Data in text are represented as mean ± SD. Because data were aggregated across multiple subjects, linear mixed models performed with the “statsmodels” Python toolbox were used to disentangle the fixed effects of subject population, event condition, or spectral band power from the random effects of each subject's dataset (Lindstrom and Bates, 1988; Seabold and Perktold, 2010). All linear mixed models were random intercept models, where each subject's dataset was assigned a random intercept 1|Subject. Once a model was fit, p values were calculated from z-scored parameter estimates (parameters estimates divided by their SEs) against the normal distribution. Because directed connectivity measures (PSI, GP) use multiple epochs for a single estimate of directed connectivity, linear mixed models were not able to account for individual patient variability in these results. Instead, we used bootstrapping where recordings were shuffled across all epochs aggregated across all subjects, and p values were calculated empirically from the resulting distribution. All other statistical tests, unless otherwise specified, were conducted in the “scipy” Python environment. The p values were controlled for multiple comparisons by using the Benjamini–Hochberg procedure at q = 0.05 (Benjamini and Hochberg, 1995).

Data availability.

The datasets supporting the current study have not been deposited in a public repository because they contain patient information but are available along with analysis code on request.

Results

Intraoperative behavioral and neural data acquisition

Ten patients with PD undergoing DBS implantation and 14 age-matched control subjects (see Materials and Methods) performed a simple visual–motor task where they followed an onscreen target using a joystick-controlled cursor with their dominant hand (Fig. 1A). Each patient performed one to four sessions of this target-tracking task during the procedure for a total of 27 sessions, while control subjects performed 1 session each for a total of 14 sessions. Tremor amplitude and cursor speed were quantified continuously from the x- and y-joystick traces (control subjects, n = 1856 epochs; PD, n = 3400 epochs). As patients were not all clinically tremor dominant, and because all subjects contributed variable amounts of data, linear mixed models were used to quantify the difference of tremor/speed across subject populations while accounting for individual subject variability. While the resulting PD and control speed distributions were distinct (linear mixed-model coefficient = 0.884, z = – 3.138, p = 0.002), PD distributions trended toward having increased tremor relative to control subjects (linear mixed-model coefficient = 0.180, z = 1.859, p = 0.063; Fig. 1B,C). Nevertheless, the partial overlap of the PD and control tremor distributions (indicative of periods without tremor in PD patients), along with the long right tail of the PD distribution, gave us a large dynamic range of tremor to analyze with respect to neural data. The dominant tremor frequency across patients was 4.48 ± 0.57 Hz. While tremor amplitude correlated with the resting tremor UPDRS subscore across patients (Spearman ρ = 0.92, p < 0.001, bootstrap test), it did not with the postural tremor subscore (ρ = 0.54, p = 0.065, bootstrap test). Based on the distinct tremor frequency peak and its correlation with clinical measures of resting tremor, we interpreted our task-derived tremor as reflective of resting tremor (Dirkx et al., 2018).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Tremor and movement speed calculated from the intraoperative visual–motor task. A, Left, Schematic of task target (green) and joystick (gray) traces from a single trial. Middle top, Bandpass-filtered x- and y-joystick traces from the task trial. Middle bottom, Low-pass-filtered x- and y-joystick traces from the task trial. Right top, One-dimensional projection of bandpass-filtered traces (black), with tremor amplitude measured from the envelope (orange). Right bottom, Cursor speed measured from low-pass-filtered traces (black). B, Distribution of 4 s tremor amplitude epochs for control subject and PD patient populations. ∘, Degrees of visual angle. Vertical dashed line indicates ROC-derived cutoff value between control and PD populations. While there is overlap on the left side of the distribution (patients with PD can exhibit control-like performance), the PD distribution is highly skewed on the right side of the distribution, allowing a large range of tremor expression. C, Distribution of 4 s speed epochs for control subject and PD patient populations. The bimodality of the control distribution corresponded to the preprogrammed speed of the onscreen target. Despite this, note that the PD distribution is shifted toward lower speed values. D, Coronal view of microelectrode recording density on an MNI reference volume. The inset panel displays a close-up view of the subthalamic nucleus (outlined in black). L, Left. E, Recording density of ECoG contacts on an MNI reference surface.

Across the 10 patients with PD, we obtained 81 microelectrode recordings within the STN (peak recording density: MNI coordinates, x=−13,y=−11,z=−5; Fig. 1D) as well as 72 ECoG recordings from cortex, including PMC (n = 27 recordings), MC (n = 16), SC (n = 15), and PPC (n = 14; Fig. 1E). As all patients were right handed, all STN and cortical recordings were obtained from the left hemisphere.

Tremor is a neurophysiologically distinct motor feature of Parkinson's disease

To understand the relationship of broadband neural activity to tremor expression, we examined the correlation between tremor amplitude and spectral power in neurophysiological recordings. Sorting session-wide spectral data by tremor epochs (rather than according to time) revealed informative band-specific patterns of activity (STN, n = 81 session recordings; PMC, n = 75; MC, n = 49; SC, n = 51; PPC, n = 41; Fig. 2A). Specifically, across cortical structures with the exception of PMC, spectral power in low and high beta frequency ranges (12–30 Hz) were found to negatively correlate with tremor amplitude (linear mixed-model coefficients = −0.325 to −0.902, z = −5.000 to −18.931, p ≤ 5.77 * 10–7; Fig. 2B). Interestingly, beta power appeared to drop off fairly quickly with even low levels of tremor becoming evident (SC: power curve fit, r2 = 0.77; linear fit, r2 = 0.54). Meanwhile, theta power positively correlated with tremor amplitude in the STN, MC, and SC (linear mixed-model coefficients = 0.076–0.732, z = 3.875–6.569, p ≤ 1.07 * 10–4).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Tremor and slowness exhibit distinct spectral power correlations with intracranial recordings. A, Population-averaged task session spectral power, sorted by the tremor amplitude of each epoch (left) or slowness (right). For ease of visualization, frequency power was z-scored within frequencies across epochs. B, Average session-wide narrowband (1 Hz) spectral Spearman correlation (ρ) with tremor amplitude and slowness. Note that while beta frequencies exhibited an opposing relationship with tremor and slowness, theta frequencies exhibited a distinct positive correlation with tremor.

To compare tremor-related neural activity with a distinct PD motor feature (specifically bradykinesia), neural data were also analyzed with respect to movement “slowness” during the same target-tracking task. Note that PD subjects appeared to lack a higher mode of movement velocity that was clearly present in control subjects, reflecting an inability to move the cursor consistently as quickly as the target (Fig. 1C). We calculated a minimum-maximum normalized measure of inverse cursor speed (0 = highest speed, 1 = lowest speed) to capture this effect as a positive pathologic sign, parallel to the sign of tremor. In contrast to tremor, we observed positive correlations between slowness and alpha/beta (8–30 Hz) power in all cortical structures (linear mixed-model coefficients = 0.159–1.141, z = 6.937–20.587, p ≤ 4.01 * 10–12; Fig. 2B). However, theta did not show a significant correlation with slowness in any structure (p > 0.05). Thus, theta appeared to relate specifically to tremor, whereas the relationship to beta activity was generally reversed between these PD-related motor manifestations. So while there was broadly the appearance of a symmetric opposition between tremor and slowness in terms of their correlations with neural activity across frequencies (Fig. 2B), this difference in the theta frequency relationship, as well as perhaps a consistent difference in mid-gamma (in which the correlation with tremor was typically close to 0, but the correlation with slowness was typically greater in magnitude and negative in direction), suggest that these motor features are not simply opposite ends of a single spectrum but rather have distinct fingerprints in neural activity.

Subthalamic theta preceded tremor at onset

Because lower-frequency oscillations, particularly theta, were most consistently and strongly positively associated with tremor across structures, and because they encompassed the range of observed tremor frequencies from a behavioral perspective (4–6 Hz), we next turned our attention to understanding the relationship of theta band activity within each structure to tremor-defined epochs. Using a control versus PD subject ROC-derived tremor threshold (see Materials and Methods), behavioral and spectral data were organized into 4 s epochs and categorized as follows: no tremor epochs (n = 575 epochs, 2300 s), tremor-onset epochs (n = 406 epochs, 1624 s), and sustained tremor epochs (n = 171 epochs, 684 s; Fig. 3A). All 10 patients contributed at least one epoch to the No Tremor and Tremor Onset conditions, with 6 patients contributing at least one epoch to the Sustained Tremor condition. The resulting behavioral and spectral data were aggregated across subjects (No Tremor: n = 1725 tremor–STN paired recording epochs, Tremor Onset: n = 1218, Sustained Tremor: n = 513).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Spectral power during different tremor dynamic states. A, Tremor event design. Based on a population-based tremor ROC threshold, epochs representing different states of tremor dynamics were isolated. For each event type, the average tremor amplitude (±SE) in patients with PD relative to control subjects is displayed over time. Horizontal dashed line denotes the tremor threshold (3 SDs relative to control subjects). Vertical dashed line (t = 0) in tremor-onset events represents the “trigger” where tremor amplitude crossed the tremor threshold. B, Average spectral power (±SE) across frequencies for each tremor event type, by recording site. Vertical dashed lines represent frequency band borders. While theta oscillations increased in power across STN, MC, and SC, increased tremor was associated with increased alpha/low beta power in PMC and PPC.

STN theta power was indeed significantly elevated during tremor onset (linear mixed-model coefficient = 0.070, z = 8.039, p=9.047*10−16) and sustained tremor (linear mixed-model coefficient = 0.129, z = 9.729, p=2.264*10−22) relative to no tremor (Fig. 3B). Likewise, phase synchrony (measured as PLV) between STN theta and tremor was increased during tremor onset (linear mixed-model coefficient = 0.080, z = 13.331, p=1.54*10−40) and sustained tremor (linear mixed-model coefficient = 0.126, z = 13.738, p=5.99*10−43; Fig. 4A).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Neural theta exhibited structure-specific temporal relationships with tremor. A, Histograms of per-trial PLVs between tremor and neural theta by tremor state. Solid lines indicate normal distribution fit to each tremor state PLV histogram, while vertical dashed lines indicate the median of each tremor state PLV histogram. The y-axis indicates the proportion of trials within each PLV histogram bin. Note that STN histograms for tremor onset and sustained tremor are highly overlapping. B, PSI between tremor and neural theta by tremor state. Positive values indicated that tremor phase preceded neural phase, while negative values indicated that neural phase preceded tremor. Magenta asterisks indicate significant (p < 0.05, bootstrap test) PSI effects.

In light of this close relationship between STN theta and tremor, we next examined the temporal relationship between STN theta and tremor phase. Specifically, we calculated the PSI between tremor and STN theta phase. Because the PSI considers multiple phase relationships within a range of frequencies, it can succeed in determining the net leading or lagging oscillation in a manner that avoids the circularity problem inherent in methods such as the PLV (Nolte et al., 2008). Here, the PSI revealed STN theta led tremor exclusively during tremor onset (p = 0.011, bootstrap test; Fig. 4B), consistent with a causal role for the STN in the initiation but not necessarily the maintenance of tremor.

Somatosensory cortex theta consistently followed tremor

Like the STN, SC theta power positively correlated with tremor amplitude. Therefore, we investigated whether this spectral–tremor relationship varied similarly with tremor state (No Tremor, n = 1256 tremor–SC paired recording epochs; Tremor Onset, n = 746; Sustained Tremor, n = 150). SC theta power was indeed significantly elevated during tremor onset (linear mixed-model coefficient = 0.010, z = 4.831, p=1.36*10−6) and sustained tremor (linear mixed-model coefficient = 0.020, z = 5.475, p=4.38*10−8), relative to no tremor (Fig. 3B). SC-tremor theta PLV also was increased during tremor onset (linear mixed-model coefficient = 0.039, z = 5.967, p=2.42*10−9) and sustained tremor (linear mixed-model coefficient = 0.180, z = 15.793, p=7.50*10−56; Fig. 4A).

However, in contrast to the STN, phase-slope analysis of tremor and SC theta phase revealed that SC theta phase followed tremor phase during both tremor onset and sustained tremor (p ≤ 0.002, bootstrap test; Fig. 4B). Therefore, the strong tremor-related theta oscillation seen in SC was reflective rather than causal of tremor.

Motor cortex theta consistently preceded tremor

Like the STN and SC, MC theta power showed a clear graded relationship with tremor magnitude (Fig. 2). Examining MC theta power across tremor states (No Tremor, n = 1066 tremor–MC paired recording epochs; Tremor Onset, n = 692; Sustained Tremor, n = 312) revealed that it was relatively increased during tremor onset (linear mixed-model coefficient = 0.006, z = 2.701, p = 0.007) but not sustained tremor (linear mixed-model coefficient = 0.002, z = 0.429, p = 0.668; Fig. 3B). Furthermore, MC-tremor theta PLV increased from no-tremor to tremor-onset (linear mixed-model coefficient = 0.018, z = 2.646, p = 0.008) to sustained tremor (linear mixed-model coefficient = 0.105, z = 10.184, p=2.34*10−24; Fig. 4A). Interestingly, examining the PSI for MC theta and tremor revealed that MC theta led tremor during both tremor onset and sustained tremor (p ≤ 0.014, bootstrap test; Fig. 4B). Thus, in contrast to SC, MC theta preceded tremor output.

Tremor-related theta transitioned from STN to cortex during tremor onset

Because both STN and MC theta power were elevated during tremor onset, and STN and MC theta phase led tremor phase during tremor onset, we investigated the dynamics of STN–MC coupling during the dynamics of tremor initiation (No Tremor, n = 3198 STN–MC paired recording epochs; Tremor Onset, n = 2076; Sustained Tremor, n = 936). Static phase slope analysis of STN and MC revealed that STN theta led MC theta during tremor onset (p < 0.001, bootstrap test; Fig. 5A). To understand whether this phase relationship was time locked to increasing tremor, we calculated STN-MC theta PSI as a function of time within the tremor onset window. Within this epoch, STN theta preceded MC theta most consistently ∼0.5 s after tremor detection (t = 0) to the end of the tremor onset epoch (t = 0.5–1.0 s; p < 0.05, bootstrap test; Fig. 5B). At no point in this window did MC theta appear to precede STN theta.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Tremor initiation was driven by the subthalamic nucleus. A, Static PSI between STN and MC recordings during tremor states. Magenta asterisks indicate significant (p < 0.05, bootstrap test) PSI effects. B, Dynamic PSI between STN and MC theta during tremor onset. Highlighted regions indicate significant PSI (p < 0.05, bootstrap test). Vertical dashed line (t = 0) indicates tremor-onset trigger. C, Directed GP between STN and MC theta during tremor onset. Vertical dashed line (t = 0) indicates tremor-onset trigger. Highlighted regions indicate significant Granger prediction (p < 0.001, bootstrap test). D, Static PSI between STN and SC recordings during tremor states. Magenta asterisks indicate significant (p < 0.05, bootstrap test) PSI effects. E, Dynamic PSI between STN and SC theta during tremor onset. Highlighted regions indicate significant PSI (p < 0.05, bootstrap test). Vertical dashed line (t = 0) indicates tremor-onset trigger. F, Directed GP between STN and SC theta during tremor onset. Vertical dashed line (t = 0) indicates tremor-onset trigger. Highlighted regions indicate significant Granger prediction (p < 0.001, bootstrap test). G, Summary of theta PSI results. Solid lines represent directed functional connectivity between neural regions and tremor.

We also investigated whether STN theta and MC theta power influenced each other by calculating time-varying nonparametric spectral GP (see Materials and Methods). Briefly, a nonzero GP at a particular frequency indicated that spectral power in one structure was predictive of spectral power in another. Unlike the PSI, GP allows the disentangling of asymmetric, bidirectional influences across two signals (Dhamala et al., 2008). As with PSI, STN theta power predicted MC theta power from 200 ms after the tremor onset trigger to the end of epoch (t = 0.2–1.0 s; p < 0.05, bootstrap test; Fig. 5C). Again, MC theta did not predict STN theta at any point in the epoch. Together, these results converged to suggest STN theta drove MC theta during tremor onset.

Once tremor was established, however, the theta phase slope relationship flipped, with MC theta phase preceding STN theta phase (Fig. 5A), revealing a dynamic transition with increasing tremor. Together with the loss of STN theta influence over tremor during sustained tremor (Fig. 4B), tremor output appeared to become cortically driven rather than STN driven as tremor became established.

Because the STN and SC both exhibited positive correlations between theta power and increasing tremor, we also investigated whether STN/SC dynamics varied during tremor onset (No Tremor, n = 3768 STN–SC paired recording epochs; Tremor Onset, n = 2238; Sustained Tremor, n = 450). Like MC, static phase slope analysis of STN and SC theta revealed that STN theta led SC during tremor onset (p < 0.001, bootstrap test; Fig. 5D). Dynamic STN–SC PSI additionally revealed that STN theta led SC theta between 200 ms after the tremor onset trigger to the end of the epoch (t = 0.2–1.0 s; p < 0.001, bootstrap test; Fig. 5E). Simultaneously, STN theta power predicted SC theta power from 400 ms before the tremor-onset trigger to end of the tremor-onset epoch (t = – 0.4 to +1.0 s; p < 0.001, bootstrap test; Fig. 5F). During sustained tremor epochs, however, the theta-phase slope relationship between STN and SC became ambiguous (p = 0.091, bootstrap test), again representing a loss of STN influence over cortical theta activity (Fig. 5D). Altogether, although the STN drove both tremor and cortical theta as tremor emerged, the transition to sustained tremor was accompanied by a decoupling of the STN from cortex in the theta band (Fig. 5G).

Motor cortex lost influence over posterior cortices with increasing tremor

As STN–MC theta phase influence flipped from tremor onset to sustained tremor, we investigated whether the functional connectivity of MC extended to other cortical regions with increasing tremor. To understand whether tremor-mediated corticocortical interactions occurred in frequency bands other than theta, we calculated both nondirected (PLV) and directed (GP) functional connectivity between the MC and other cortical regions across the 3–100 Hz spectrum (paired recording epochs from No Tremor, Tremor Onset, and Sustained Tremor conditions: MC–PMC: n = 2692, 2064, 757; MC–SC: n = 2190, 1130, 210; MC–PPC: n = 1458, 1074, 935). MC–SC PLV across any frequency band did not modulate by tremor state (PLV; linear mixed model, p > 0.05; Fig. 6A). To identify whether synchrony detected by the PLV was driven by one structure in the pair, broad-spectrum GP was calculated. In the absence of tremor, we found that MC predicted SC high beta/low gamma power (p < 0.001, bootstrap test) and SC predicted MC theta/alpha/low beta/mid-gamma power (p < 0.001, bootstrap test; Fig. 6B). During sustained tremor, however, MC theta now predicted SC theta (GP, 2.34-fold difference, p < 0.001, bootstrap test), and SC low beta/low gamma predicted MC low beta/low gamma (GP, 1.99- to 2.18-fold difference; p < 0.001, bootstrap test).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

During sustained tremor, gamma coupling between premotor/motor and somatosensory/parietal cortices decreased. A, PLV between MC and other cortical regions. Lines ± shaded borders represent average ± SE PLV. Highlighted frequency ranges indicate increased (orange) or decreased (blue) PLV with increasing tremor. B, Pairwise GP between MC and other cortical regions. The title of each subpanel indicates the directionality of the structure pair GP. Highlighted frequency ranges indicate increased (orange) or decreased (blue) GP with increasing tremor. Note that MC broad-spectrum coupling with SC and PPC generally decreased with increasing tremor. C, PLV between SC and other cortical regions. Lines ± shaded borders represent average ± SE PLV. Highlighted frequency ranges indicate increased (orange) or decreased (blue) PLV with increasing tremor. D, Pairwise GP between SC and other cortical regions. The title of each subpanel indicates the directionality of the structure pair GP. Highlighted frequency ranges indicate increased (orange) or decreased (blue) GP with increasing tremor. Note that tremor generally shifted the frequency of coupling among SC, PPC, and PMC from gamma to alpha/low beta with increasing tremor. E, PLV between the STN and cortical regions. Lines ± shaded borders represent average ± SE PLV. Highlighted frequency ranges indicate increased (orange) PLV with increasing tremor. F, Pairwise GP between the STN and cortical regions. The title of each subpanel indicates the directionality of the structure pair GP. Highlighted frequency ranges indicate increased (orange) GP with increasing tremor, and increased GP specific to Tremor Onset (gray). Note that while broad-spectrum STN–cortical PLV generally increased with increasing tremor, directional changes were less distinct. For ease of visualization, GP curves were low-pass filtered and frequencies within 58–62 Hz were masked. Vertical dashed lines represent frequency band borders.

MC–PPC PLV similarly did not modulate as tremor increased (PLV linear mixed model, p > 0.05; Fig. 6A). However, Granger analysis revealed that PPC theta/alpha/low beta predicted MC theta/alpha/low beta regardless of tremor state (GP, 1.02–4.20-fold difference; p < 0.001, bootstrap test in all tremor states; Fig. 6B). In contrast, while MC high beta/low gamma/mid-gamma predicted PPC high beta/low gamma/mid-gamma in the absence of tremor (GP, 1.07–1.48-fold difference; p < 0.001, bootstrap test), this relationship flipped during sustained tremor, with PPC high beta/low gamma/mid-gamma predicting MC high beta/low gamma/mid-gamma (GP, 1.50–1.99-fold difference; p < 0.001, bootstrap test).

In sum, MC exerted less influence over posterior (SC, PPC) and anterior (PMC) cortical regions with increasing tremor. Specifically, MC–PMC PLV decreased within low beta/low gamma (12–20/30–60 Hz) specifically during sustained tremor (PLV linear mixed-model coefficients, −0.010 to −0.017; z = −2.100 to −2.323, p ≤ 0.035; Fig. 6A). While not within the same frequency range, PMC theta/alpha also appeared to predict MC theta/alpha during sustained tremor (GP, 2.68- to 2.71-fold increase; p < 0.001, bootstrap test; Fig. 6B).

Premotor cortex coupled with posterior cortices during tremor

Because SC decoupled from the STN during sustained tremor while still reflecting tremor output, we investigated whether SC instead coupled with other cortical regions as tremor increased (paired recording epochs from No Tremor, Tremor Onset, and Sustained Tremor conditions: SC–PMC: n = 3442, 2292, 459; SC–PPC: n = 1284, 808, 165; PMC–PPC: n = 1780, 1462, 1029). SC–PPC PLV similarly did not modulate as tremor increased (PLV, linear mixed model, p > 0.05; Fig. 6C). Although PLV did not significantly modulate with tremor, PPC theta/alpha/low beta (8–20 Hz) predicted SC theta/alpha/low beta during sustained tremor (GP, 1.61- to 8.55-fold difference; p < 0.001, bootstrap test; Fig. 6D). While SC low gamma/mid-gamma predicted PPC low gamma/mid-gamma during the absence of tremor (GP, 1.23- to 1.32-fold difference; p < 0.001, bootstrap test), this did not hold for sustained tremor. Thus, SC–PPC connectivity shifted to a distinct state during sustained tremor, with PPC predicting lower frequencies (theta, alpha, low beta) in SC. At the same time, higher frequency (gamma) directed connectivity between SC and PPC decreased as tremor increased.

SC and PMC interactions exhibited decreases in functional connectivity, with decreased high beta–low gamma (20–60 Hz) PLV (PLV linear mixed-model coefficients, −0.022 to −0.007; z = −1.975 to −4.488; p < = 0.048) with increasing tremor (Fig. 6C). Like PPC, SC alpha/low beta was driven by PMC alpha/low beta specifically during sustained tremor (GP, 4.41- to 13.10-fold difference; p < 0.001, bootstrap test; Fig. 6D). Thus, in contrast to MC, which lost influence over posterior cortical regions, SC became increasingly influenced by both posterior (PPC) and anterior (PMC) cortices with increasing tremor. However, this increase in connectivity was specific to alpha/low beta frequencies while gamma coupling decreased between SC and PMC/PPC.

To follow the spread of tremor-related cortical coupling, we investigated whether PMC and PPC interacted during sustained tremor. Here, we observed an exaggerated version of the same tremor-induced frequency shift (gamma to beta) of power and phase synchrony. When analyzing tremor epoch-related spectral power in PMC and PPC in Figure 3B, both regions demonstrated tremor-related decreases in low beta/high beta frequencies (linear mixed-model coefficients, −0.007 to −0.019; z = −2.358 to −5.256; p ≤ 0.018). At the same time, PMC exhibited decreases in low gamma power during sustained tremor relative to no tremor (linear mixed-model coefficients, −0.003; z = – 2.895; p = 0.004).

These similar changes in power were mirrored by changes in PMC–PPC PLV synchrony (Fig. 6C). PMC–PPC low to mid-gamma PLV decreased as tremor increased (PLV linear mixed-model coefficients, −0.016 to−0.020; z = −3.021 to −3.367; p ≤ 0.003), while PMC–PPC alpha/low beta PLV increased with tremor (PLV linear mixed-model coefficients, 0.018–0.020; z = 2.348–3.253; p ≤ 0.018). Regardless of tremor state, PMC–PPC phase synchrony was driven by PMC onto PPC. When tremor was absent, PMC low gamma/mid-gamma predicted PPC low gamma/mid gamma (GP, 1.98- to 2.12-fold difference; p < 0.001, bootstrap test; Fig. 6D). During sustained tremor, PMC low beta/high beta power predicted PPC low beta/high beta power (GP, 2.87- to 6.26-fold difference; p < 0.001, bootstrap test).

Overall, tremor was associated with a frequency shift (gamma to beta) of power and phase synchrony among PMC, PPC, and SC. Specifically, PMC exerted increasing influence over posterior regions (SC, PPC) in lower frequencies (alpha, low beta) with increasing tremor. However, this increase in lower-frequency coupling coincided with decreases in higher-frequency coupling (gamma). In addition, directional gamma influence between MC and PPC flipped with increasing tremor (MC → PPC in the absence of tremor, PPC → MC during sustained tremor), revealing that sustained tremor is a state of altered gamma synchrony across sensorimotor cortex.

STN broadly synchronized with, but selectively influenced, sensorimotor cortex during sustained tremor

Finally, to understand whether STN influence over cortex extended beyond theta, functional and directed connectivity were calculated between the STN and sensorimotor cortex (paired recording epochs from No Tremor, Tremor Onset, and Sustained Tremor conditions: STN–PMC: n = 4680, 3732, 1281; STN–MC: n = 3198, 2076, 936; STN–SC: n = 3768, 2238, 450; STN–PPC: n = 2154, 1698, 1437). STN–cortical theta PLV synchrony (with the exception of SC) increased as a function of tremor (PLV linear mixed-model coefficients, 0.012–0.025; z = 2.873–6.827; p ≤ 0.004; Fig. 6E). However, directed connectivity between the STN and cortex was specific to tremor state (Fig. 6F). STN theta power (4–6 Hz) predicted both MC theta power (GP, 2.29-fold difference; p < 0.001, bootstrap test) and SC theta power (GP, 1.79-fold difference; p < 0.001, bootstrap test) exclusively during tremor onset. In contrast, STN theta power predicted PPC theta power (GP, 1.86-fold difference; p < 0.001, bootstrap test) and PMC theta power (GP, 1.59-fold difference; p < 0.001, bootstrap test) only during sustained tremor. Thus, consistent with the PSI results, the STN shifted its influence over cortex in the theta band (STN → MC/SC during tremor onset; STN → PMC/PPC during sustained tremor) across dynamic tremor states.

Discussion

Using a naturalistic behavioral task, we were able to characterize tremor dynamics and to isolate specific tremor states, particularly tremor onset and maintenance. Across structures, we found that theta power positively and beta power negatively correlated with tremor, as has been found in previous reports (Hirschmann et al., 2013; Qasim et al., 2016; Asch et al., 2020). However, our study is the first to dissect electrophysiological correlates of tremor onset and sustained tremor. During the emergence of tremor, not only did STN and motor cortical theta power increase, but STN and motor cortical theta phase preceded the phase of tremor. Moreover, STN theta activity drove motor cortical theta during tremor onset, suggesting a direct role of the STN in initiating tremor output.

Once tremor emerged, however, motor cortex appeared to sustain tremor. At the same time, motor cortex became less coupled with somatosensory and parietal cortices, despite the presence of prominent somatosensory cortex theta power, which closely followed tremor. Instead, premotor cortex synchronized via low beta frequencies with posterior cortices (somatosensory, parietal) at the expense of gamma frequency synchronization observed in the absence of tremor. This low beta synchrony was notably asymmetric across these structures, with premotor cortex exerting influence over posterior cortices.

Together, although tremor amplitude corresponded to global changes in theta and beta power, the relationship between these frequency bands to tremor output was highly structure specific. While STN–motor cortical interactions appeared to initiate tremor, premotor cortex-driven network effects may help sustain tremor. This STN-mediated dynamic reorganization of cortical connectivity is consistent with both the “dimmer switch” model and the intrinsic and extrinsic cortical loops of Parkinson's tremor (Helmich et al., 2011; Volkmann et al., 1996; Fig. 7). Like the GPi, we revealed that the STN acted as a “switch” to mediate the onset of tremor by influencing motor cortex (Dirkx et al., 2016). While these STN–motor cortical interactions formed the intrinsic loop of tremor output, we expanded this model to reveal that shifts from gamma to beta synchrony across premotor–parietal cortices reflected the extrinsic loop in the stable tremor state.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Synthetic model of subcortical–cortical interactions during tremor. Solid lines represent directed functional connectivity between neural regions and tremor. Dashed lines during Sustained Tremor represent interactions from the No Tremor state that are no longer present.

Tremor onset was mediated by subthalamic theta driving motor cortex

STN theta amplitude positively correlated with tremor amplitude regardless of tremor dynamic states. While the phase of STN theta consistently preceded tremor phase during tremor onset, it did not during sustained tremor. However, STN theta activity was still significantly phase locked to tremor during sustained tremor. This mixed relationship to tremor may reflect several roles of STN: interconnections with GPi contribute to tremor initiation, while disynaptic connections with cerebellum may influence ongoing monitoring of tremor output (Bostan et al., 2010; Helmich et al., 2011). Indeed, STN projections to cerebellar cortex may perhaps propagate tremor–frequency oscillations within the basal ganglia to motor cortical–thalamo–cerebellar loops (Wu and Hallett, 2013; Bostan and Strick, 2018). Future experiments combining STN and cerebellar recordings could describe this tremor-onset mechanism, while trying to disentangle the neural control of tremor amplitude and phase (Cagnan et al., 2014; Helmich et al., 2021).

Regardless, STN theta drove motor cortex activity during tremor onset. Although tremor has previously been found to decrease beta coherence between STN and motor cortex (Qasim et al., 2016), while increasing theta coherence (Hirschmann et al., 2013), we demonstrated directed theta-phase interactions from STN to motor cortex specifically during tremor onset. While a previous case study of tremor onset displayed local STN and cortical alpha/beta power changes with tremor onset (Hirschmann et al., 2019), we show here that STN and motor cortical theta activity are directionally linked. We also demonstrated that during sustained tremor, the STN–motor cortex theta phase slope relationship reversed, suggesting the theta influence over sustained tremor shifted source from STN to cortex.

Motor cortex desynchronized with posterior cortices while sustaining tremor

As tremor progressed, motor cortex theta increasingly drove tremor. While previous studies have correlated motor cortical activity to tremor (Timmermann et al., 2003; Helmich et al.,2011), this is the first study, to our knowledge, that has demonstrated a directed relationship between ECoG recordings and tremor. Although motor cortex was synchronized to tremor, motor cortex appeared to desynchronize with other cortical structures with the exception of premotor cortex, as has been found previously (Timmermann et al., 2003; Qasim et al., 2016). While other studies have found that motor cortex increased its synchrony with premotor and parietal cortices during tremor (Hirschmann et al., 2013), this was calculated only at tremor and double-tremor frequencies.

Tremor reorganized premotor and parietal cortical coupling

Although premotor and parietal cortices did not exhibit a direct theta relationship to tremor, changes in tremor initiated a frequency shift in premotor–parietal coupling dynamics. In the absence of tremor, these regions were functionally coupled at higher frequencies (high beta, low to mid-gamma). fMRI studies in patients with PD have found that these regions exhibit overactive BOLD activity during self-initiated sequential hand movements (Samuel et al., 1997), which is hypothesized to compensate for decreased BOLD activity in frontostriatal circuits in the dopamine-depleted state (Wu et al., 2011). Furthermore, cortical gamma frequency power and synchrony are associated specifically with voluntary movement (Crone et al., 1998; Miller et al., 2007). In our study, this bidirectional premotor–parietal gamma activity may have reflected task monitoring and spatial tracking (motor output) using sensory information.

During sustained tremor, however, parietal and premotor cortices both exhibited increases in low beta power. This low beta activity was also functionally coupled, with premotor driving parietal cortex. Elevated low beta oscillations have been observed in premotor cortex recordings in MPTP nonhuman primates with predominantly akinetic/rigid symptoms (Wang et al., 2017). While not observed in our study, increased premotor high beta influence over the STN has also been found to correlate with akinetic/rigid symptoms (Sharott et al., 2018). Premotor low beta oscillations may function here in a similar antikinetic fashion with other cortical structures during tremor.

In any case, with increasing tremor premotor–parietal gamma activity diminished while premotor low beta activity drove parietal activity. These frequency shifts may be best understood in the framework of communication-through-coherence theory (Fries, 2015). Specifically, while symmetric or bottom-up gamma oscillations permit effective and precise transmission of motor-related information across structures, lower-frequency oscillations such as alpha/beta act as top-down feedback. Here, task-related gamma synchrony observed across sensorimotor cortex decreased with tremor. In contrast, lower-frequency oscillations such as low beta increased in synchrony, perhaps reflecting an absence of voluntary movement, which normally acts to suppress tremor (Naros et al., 2018).

Implications for closed-loop deep brain stimulation

Because of the clinical interest in developing adaptive closed-loop DBS to more precisely treat PD symptoms, various electrophysiological observations have been investigated as potential tremor biomarkers to inform stimulation (Hirschmann et al., 2017; Shah et al., 2018; Yao et al., 2020). While promising, the features used for tremor detection do not take into account the dynamic nature of tremor—namely, the distinct neurophysiological signature of tremor onset. Because of the breadth of STN beta-frequency oscillation research in PD, initial closed-loop DBS efforts have focused on using beta oscillations as a proxy for bradykinesia symptoms (Little et al., 2013, 2016a,b; Velisar et al., 2019). However, beta-driven DBS has been shown to worsen tremor in some patients (He et al., 2020; Piña-Fuentes et al., 2020).

Here, we demonstrated that subthalamic theta was present whether tremor was emerging or sustained. The addition of STN theta-based biomarkers to closed-loop DBS could help treat the separate symptom axis of tremor. Further, we have provided the best evidence to date that cortical ECoG theta is a robust marker for tremor. Specifically, we found that motor cortical theta was synchronized to STN theta during tremor states, and that somatosensory theta was a reliable indicator of immediate tremor amplitude.

These results overall argue for a combined subcortical–cortical stimulation/recording paradigm not unlike cortical–thalamic closed-loop DBS for ET (Opri et al., 2020). By combining recordings from the STN and sensorimotor cortex, an algorithm could infer whether tremor was about to emerge (STN and MC theta) or was already present (SC theta). In particular, somatosensory cortical recordings could allow for continuous monitoring of tremor despite any stimulus artifact or competing oscillations in the STN. Ideally, DBS for a patient with a mixed-motor phenotype could be optimized between STN beta for bradykinesia symptoms and SC theta oscillations for tremor.

Limitations and conclusions

Because all tremor data were quantified from patients as they were moving their upper limb during our tracking task, our tremor conditions do not reflect a pure “rest” tremor. However, as parkinsonian tremor can often emerge as patients maintain a posture or perform a task, our approach still captured meaningful aspects of tremor. Because of our PD population receiving mostly STN DBS for clinical reasons, we were unable to assess the role of the GPi and cerebellar thalamus (VIM) neurophysiology to tremor onset and/or maintenance. In addition, as increased cognitive load has been found to exacerbate tremor, our observed tremor-related changes in nontremor frequencies within cortex may have reflected cognitive or visuomotor processes (e.g., eye movements) not directly related to tremor (Dirkx et al., 2020). Although we attempted to overcome the influence of individual subjects in our tremor epoch datasets by using linear mixed models, we were unable to apply linear mixed models to our directed connectivity analyses (PSI, GP) and thus may be susceptible to individual subject influence. However, our directed connectivity results were often reinforced by nondirected measures of functional connectivity (PLV), suggesting that directed results reflected the same underlying phenomena. Nevertheless, our awake behaving intraoperative recordings revealed that the STN and motor cortex work together to initiate tremor, and tremor is in part sustained by premotor–parietal synchrony.

Footnotes

  • This work was supported by a National Institutes of Health (NIH) National Institute of Neurological Disorders and Stroke Training Grant (T32-MH-020068) to P.M.L.; a Doris Duke Clinical Scientist Development Award (#2014101) to W.F.A.; an NIH COBRE Award from the National Institute of General Medical Sciences Grant P20-GM-103645 (Principal Investigator, Jerome Sanes) supporting W.F.A.; a Neurosurgery Research and Education Foundation grant to W.F.A.; the Lifespan Norman Prince Neurosciences Institute; and the Brown University Robert J. and Nancy D. Carney Institute for Brain Science. Part of this research was conducted using computational resources and services at the Center for Computation and Visualization at Brown University, with funding provided by NIH Office of the Director Grant S10-OD-025181. W.F.A. has received proprietary equipment and technical support for unrelated research through the Medtronic external research program. We thank the patients in this study for participation in the study. We also thank Kelsea Laubenstein-Parker for technical assistance, Karina Bertsch for administrative support, and Ann Duggan-Winkle for clinical support. In addition, we also thank Minkyu Ahn, David Segar, Tina Sankhla, and Daniel Shiebler for helping in the development of the motor task experiment.

  • The authors have patents and patent applications broadly relevant to Parkinson's disease (but not directly based upon this work). W.F.A. has received proprietary equipment and technical support for unrelated research through the Medtronic external research program. The authors declare no other competing financial interests.

  • Correspondence should be addressed to Peter M. Lauro at peter_lauro{at}brown.edu or Wael F. Asaad at wfasaad{at}alum.mit.edu

SfN exclusive license.

References

  1. ↵
    1. Accolla E,
    2. Caputo E,
    3. Cogiamanian F,
    4. Tamma F,
    5. Mrakic-Sposta S,
    6. Marceglia S,
    7. Egidi M,
    8. Rampini P,
    9. Locatelli M,
    10. Priori A
    (2007) Gender differences in patients with Parkinson's disease treated with subthalamic deep brain stimulation. Mov Disord 22:1150–1156. doi:10.1002/mds.21520 pmid:17469208
    OpenUrlCrossRefPubMed
  2. ↵
    1. Akbar U,
    2. Asaad WF
    (2017) A comprehensive approach to deep brain stimulation for movement disorders. R I Med J 100:30–33. pmid:28564666
    OpenUrlPubMed
  3. ↵
    1. Albin RL,
    2. Young AB,
    3. Penney JB
    (1989) The functional anatomy of basal ganglia disorders. Trends Neurosci 12:366–375. doi:10.1016/0166-2236(89)90074-x pmid:2479133
    OpenUrlCrossRefPubMed
  4. ↵
    1. Argall BD,
    2. Saad ZS,
    3. Beauchamp MS
    (2006) Simplified intersubject averaging on the cortical surface using SUMA. Hum Brain Mapp 27:14–27. doi:10.1002/hbm.20158 pmid:16035046
    OpenUrlCrossRefPubMed
  5. ↵
    1. Asaad WF,
    2. Eskandar EN
    (2008a) Achieving behavioral control with millisecond resolution in a high-level programming environment. J Neurosci Methods 173:235–240. doi:10.1016/j.jneumeth.2008.06.003 pmid:18606188
    OpenUrlCrossRefPubMed
  6. ↵
    1. Asaad WF,
    2. Eskandar EN
    (2008b) A flexible software tool for temporally-precise behavioral control in Matlab. J Neurosci Methods 174:245–258. doi:10.1016/j.jneumeth.2008.07.014 pmid:18706928
    OpenUrlCrossRefPubMed
  7. ↵
    1. Asaad WF,
    2. Santhanam N,
    3. McClellan S,
    4. Freedman DJ
    (2013) High-performance execution of psychophysical tasks with complex visual stimuli in MATLAB. J Neurophysiol 109:249–260. doi:10.1152/jn.00527.2012 pmid:23034363
    OpenUrlCrossRefPubMed
  8. ↵
    1. Asch N,
    2. Herschman Y,
    3. Maoz R,
    4. Auerbach-Asch CR,
    5. Valsky D,
    6. Abu-Snineh M,
    7. Arkadir D,
    8. Linetsky E,
    9. Eitan R,
    10. Marmor O,
    11. Bergman H,
    12. Israel Z
    (2020) Independently together: subthalamic theta and beta opposite roles in predicting Parkinson's tremor. Brain Commun 2:fcaa074. doi:10.1093/braincomms/fcaa074 pmid:33585815
    OpenUrlCrossRefPubMed
  9. ↵
    1. Aydore S,
    2. Pantazis D,
    3. Leahy RM
    (2013) A note on the phase locking value and its properties. Neuroimage 74:231–244. doi:10.1016/j.neuroimage.2013.02.008 pmid:23435210
    OpenUrlCrossRefPubMed
  10. ↵
    1. Benjamini Y,
    2. Hochberg Y
    (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Stat Methodol 57:289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bergman H,
    2. Wichmann T,
    3. Karmon B,
    4. DeLong MR
    (1994) The primate subthalamic nucleus. II. Neuronal activity in the MPTP model of parkinsonism. J Neurophysiol 72:507–520. doi:10.1152/jn.1994.72.2.507 pmid:7983515
    OpenUrlCrossRefPubMed
  12. ↵
    1. Bostan AC,
    2. Dum RP,
    3. Strick PL
    (2010) The basal ganglia communicate with the cerebellum. Proc Natl Acad Sci U S A 107:8452–8456. doi:10.1073/pnas.1000496107 pmid:20404184
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Bostan AC,
    2. Strick PL
    (2018) The basal ganglia and the cerebellum: nodes in an integrated network. Nat Rev Neurosci 19:338–350. doi:10.1038/s41583-018-0002-7 pmid:29643480
    OpenUrlCrossRefPubMed
  14. ↵
    1. Cagnan H,
    2. Little S,
    3. Foltynie T,
    4. Limousin P,
    5. Zrinzo L,
    6. Hariz M,
    7. Cheeran B,
    8. Fitzgerald J,
    9. Green AL,
    10. Aziz T,
    11. Brown P
    (2014) The nature of tremor circuits in parkinsonian and essential tremor. Brain 137:3223–3234. doi:10.1093/brain/awu250 pmid:25200741
    OpenUrlCrossRefPubMed
  15. ↵
    1. Cox RW
    (1996) AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res 29:162–173. doi:10.1006/cbmr.1996.0014 pmid:8812068
    OpenUrlCrossRefPubMed
  16. ↵
    1. Crone NE,
    2. Miglioretti DL,
    3. Gordon B,
    4. Lesser RP
    (1998) Functional mapping of human sensorimotor cortex with electrocorticographic spectral analysis. II. Event-related synchronization in the gamma band. Brain 121:2301–2315. doi:10.1093/brain/121.12.2301
    OpenUrlCrossRefPubMed
  17. ↵
    1. Dale AM,
    2. Fischl B,
    3. Sereno MI
    (1999) Cortical surface-based analysis: I. Segmentation and surface reconstruction. Neuroimage 9:179–194. doi:10.1006/nimg.1998.0395 pmid:9931268
    OpenUrlCrossRefPubMed
  18. ↵
    1. Desikan RS,
    2. Ségonne F,
    3. Fischl B,
    4. Quinn BT,
    5. Dickerson BC,
    6. Blacker D,
    7. Buckner RL,
    8. Dale AM,
    9. Maguire RP,
    10. Hyman BT,
    11. Albert MS,
    12. Killiany RJ
    (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31:968–980. doi:10.1016/j.neuroimage.2006.01.021 pmid:16530430
    OpenUrlCrossRefPubMed
  19. ↵
    1. Destrieux C,
    2. Fischl B,
    3. Dale A,
    4. Halgren E
    (2010) Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature. Neuroimage 53:1–15. doi:10.1016/j.neuroimage.2010.06.010 pmid:20547229
    OpenUrlCrossRefPubMed
  20. ↵
    1. Dhamala M,
    2. Rangarajan G,
    3. Ding M
    (2008) Analyzing information flow in brain networks with nonparametric granger causality. Neuroimage 41:354–362. doi:10.1016/j.neuroimage.2008.02.020 pmid:18394927
    OpenUrlCrossRefPubMed
  21. ↵
    1. Dirkx MF,
    2. den Ouden H,
    3. Aarts E,
    4. Timmer M,
    5. Bloem BR,
    6. Toni I,
    7. Helmich RC
    (2016) The cerebral network of Parkinson's tremor: an effective connectivity fMRI study. J Neurosci 36:5362–5372. doi:10.1523/JNEUROSCI.3634-15.2016 pmid:27170132
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Dirkx MF,
    2. den Ouden HEM,
    3. Aarts E,
    4. Timmer MHM,
    5. Bloem BR,
    6. Toni I,
    7. Helmich RC
    (2017) Dopamine controls Parkinson's tremor by inhibiting the cerebellar thalamus. Brain 140:721–734. doi:10.1093/brain/aww331 pmid:28073788
    OpenUrlCrossRefPubMed
  23. ↵
    1. Dirkx MF,
    2. Zach H,
    3. Bloem BR,
    4. Hallett M,
    5. Helmich RC
    (2018) The nature of postural tremor in Parkinson disease. Neurology 90:e1095–e1103. doi:10.1212/WNL.0000000000005215 pmid:29476038
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Dirkx MF,
    2. Zach H,
    3. van Nuland A,
    4. Bloem BR,
    5. Toni I,
    6. Helmich RC
    (2019) Cerebral differences between dopamine-resistant and dopamine-responsive Parkinson's tremor. Brain 142:3144–3157. doi:10.1093/brain/awz261 pmid:31509182
    OpenUrlCrossRefPubMed
  25. ↵
    1. Dirkx MF,
    2. Zach H,
    3. van Nuland AJ,
    4. Bloem BR,
    5. Toni I,
    6. Helmich RC
    (2020) Cognitive load amplifies Parkinson's tremor through excitatory network influences onto the thalamus. Brain 143:1498–1511. doi:10.1093/brain/awaa083 pmid:32355951
    OpenUrlCrossRefPubMed
  26. ↵
    1. Fischl B,
    2. Salat DH,
    3. Busa E,
    4. Albert M,
    5. Dieterich M,
    6. Haselgrove C,
    7. van der Kouwe A,
    8. Killiany R,
    9. Kennedy D,
    10. Klaveness S,
    11. Montillo A,
    12. Makris N,
    13. Rosen B,
    14. Dale AM
    (2002) Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron 33:341–355. doi:10.1016/s0896-6273(02)00569-x pmid:11832223
    OpenUrlCrossRefPubMed
  27. ↵
    1. Fonov V,
    2. Evans A,
    3. McKinstry R,
    4. Almli C,
    5. Collins D
    (2009) Unbiased nonlinear average age-appropriate brain templates from birth to adulthood. Neuroimage 47:S102. doi:10.1016/S1053-8119(09)70884-5
    OpenUrlCrossRef
  28. ↵
    1. Fries P
    (2015) Rhythms for cognition: communication through coherence. Neuron 88:220–235. doi:10.1016/j.neuron.2015.09.034 pmid:26447583
    OpenUrlCrossRefPubMed
  29. ↵
    1. Gross RE,
    2. Krack P,
    3. Rodriguez-Oroz MC,
    4. Rezai AR,
    5. Benabid AL
    (2006) Electrophysiological mapping for the implantation of deep brain stimulators for Parkinson's disease and tremor. Mov Disord 21 [Suppl 14]:S259–S283. doi:10.1002/mds.20960 pmid:16810720
    OpenUrlCrossRefPubMed
  30. ↵
    1. Hariz GM,
    2. Nakajima T,
    3. Limousin P,
    4. Foltynie T,
    5. Zrinzo L,
    6. Jahanshahi M,
    7. Hamberg K
    (2011) Gender distribution of patients with Parkinson's disease treated with subthalamic deep brain stimulation; a review of the 2000–2009 literature. Parkinsonism Relat Disord 17:146–149. doi:10.1016/j.parkreldis.2010.12.002 pmid:21195012
    OpenUrlCrossRefPubMed
  31. ↵
    1. Harris CR,
    2. Millman KJ,
    3. van der Walt SJ,
    4. Gommers R,
    5. Virtanen P,
    6. Cournapeau D,
    7. Wieser E,
    8. Taylor J,
    9. Berg S,
    10. Smith NJ,
    11. Kern R,
    12. Picus M,
    13. Hoyer S,
    14. van Kerkwijk MH,
    15. Brett M,
    16. Haldane A,
    17. del Río JF,
    18. Wiebe M,
    19. Peterson P,
    20. Gérard-Marchant P, et al.
    (2020) Array programming with NumPy. Nature 585:357–362. doi:10.1038/s41586-020-2649-2 pmid:32939066
    OpenUrlCrossRefPubMed
  32. ↵
    1. He S,
    2. Mostofi A,
    3. Syed E,
    4. Torrecillos F,
    5. Tinkhauser G,
    6. Fischer P,
    7. Pogosyan A,
    8. Hasegawa H,
    9. Li Y,
    10. Ashkan K,
    11. Pereira E,
    12. Brown P,
    13. Tan H
    (2020) Subthalamic beta targeted neurofeedback speeds up movement initiation but increases tremor in Parkinsonian patients. eLife 9:e60979. doi:10.7554/eLife.60979
    OpenUrlCrossRef
  33. ↵
    1. Helmich RC
    (2018) The cerebral basis of Parkinsonian tremor: a network perspective. Mov Disord 33:219–231. doi:10.1002/mds.27224 pmid:29119634
    OpenUrlCrossRefPubMed
  34. ↵
    1. Helmich RC,
    2. Janssen MJR,
    3. Oyen WJG,
    4. Bloem BR,
    5. Toni I
    (2011) Pallidal dysfunction drives a cerebellothalamic circuit into Parkinson tremor. Ann Neurol 69:269–281. doi:10.1002/ana.22361 pmid:21387372
    OpenUrlCrossRefPubMed
  35. ↵
    1. Helmich RC,
    2. Hallett M,
    3. Deuschl G,
    4. Toni I,
    5. Bloem BR
    (2012) Cerebral causes and consequences of parkinsonian resting tremor: a tale of two circuits? Brain 135:3206–3226. doi:10.1093/brain/aws023 pmid:22382359
    OpenUrlCrossRefPubMed
  36. ↵
    1. Helmich RC,
    2. Van den Berg KRE,
    3. Panyakaew P,
    4. Cho HJ,
    5. Osterholt T,
    6. McGurrin P,
    7. Shamim EA,
    8. Popa T,
    9. Haubenberger D,
    10. Hallett M
    (2021) Cerebello-cortical control of tremor rhythm and amplitude in Parkinson's disease. Mov Disord 36:1727–1729. doi:10.1002/mds.28603
    OpenUrlCrossRef
  37. ↵
    1. Hirschmann J,
    2. Hartmann CJ,
    3. Butz M,
    4. Hoogenboom N,
    5. Özkurt TE,
    6. Elben S,
    7. Vesper J,
    8. Wojtecki L,
    9. Schnitzler A
    (2013) A direct relationship between oscillatory subthalamic nucleus–cortex coupling and rest tremor in Parkinson's disease. Brain 136:3659–3670. doi:10.1093/brain/awt271 pmid:24154618
    OpenUrlCrossRefPubMed
  38. ↵
    1. Hirschmann J,
    2. Butz M,
    3. Hartmann CJ,
    4. Hoogenboom N,
    5. Özkurt TE,
    6. Vesper J,
    7. Wojtecki L,
    8. Schnitzler A
    (2016) Parkinsonian rest tremor is associated with modulations of subthalamic high-frequency oscillations. Mov Disord 31:1551–1559. doi:10.1002/mds.26663 pmid:27214766
    OpenUrlCrossRefPubMed
  39. ↵
    1. Hirschmann J,
    2. Schoffelen JM,
    3. Schnitzler A,
    4. van Gerven MAJ
    (2017) Parkinsonian rest tremor can be detected accurately based on neuronal oscillations recorded from the subthalamic nucleus. Clin Neurophysiol 128:2029–2036. doi:10.1016/j.clinph.2017.07.419 pmid:28841506
    OpenUrlCrossRefPubMed
  40. ↵
    1. Hirschmann J,
    2. Abbasi O,
    3. Storzer L,
    4. Butz M,
    5. Hartmann CJ,
    6. Wojtecki L,
    7. Schnitzler A
    (2019) Longitudinal recordings reveal transient increase of alpha/low-beta power in the subthalamic nucleus associated with the onset of parkinsonian rest tremor. Front Neurol 10:145. doi:10.3389/fneur.2019.00145
    OpenUrlCrossRef
  41. ↵
    1. Jankovic J,
    2. McDermott M,
    3. Carter J,
    4. Gauthier S,
    5. Goetz C,
    6. Golbe L,
    7. Huber S,
    8. Koller W,
    9. Olanow C,
    10. Shoulson I
    (1990) Variable expression of Parkinson's disease: a base-line analysis of the DATATOP cohort. The Parkinson Study Group. Neurology 40:1529–1534. doi:10.1212/wnl.40.10.1529 pmid:2215943
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Koller WC
    (1984) The diagnosis of Parkinson's disease. Arch Intern Med 144:2146–2147. doi:10.1001/archinte.1984.04400020048005
    OpenUrlCrossRefPubMed
  43. ↵
    1. Koller WC
    (1986) Pharmacologic treatment of Parkinsonian tremor. Arch Neurol 43:126–127. doi:10.1001/archneur.1986.00520020020009 pmid:3947248
    OpenUrlCrossRefPubMed
  44. ↵
    1. Konrad PE,
    2. Neimat JS,
    3. Yu H,
    4. Kao CC,
    5. Remple MS,
    6. D'Haese P-F,
    7. Dawant BM
    (2011) Customized, miniature rapid-prototype stereotactic frames for use in deep brain stimulator surgery: initial clinical methodology and experience from 263 patients from 2002 to 2008. Stereotact Funct Neurosurg 89:34–41. doi:10.1159/000322276 pmid:21160241
    OpenUrlCrossRefPubMed
  45. ↵
    1. Lachaux JP,
    2. Rodriguez E,
    3. Martinerie J,
    4. Varela FJ
    (1999) Measuring phase synchrony in brain signals. Hum Brain Mapp 8:194–208. doi:10.1002/(SICI)1097-0193(1999)8:4<194::AID-HBM4>3.0.CO;2-C
    OpenUrlCrossRefPubMed
  46. ↵
    1. Lance JW,
    2. Schwab RS,
    3. Peterson EA
    (1963) Action tremor and the cogwheel phenomenon in Parkinson's disease. Brain 86:95–110. doi:10.1093/brain/86.1.95 pmid:13928399
    OpenUrlCrossRefPubMed
  47. ↵
    1. Lauro PM,
    2. Vanegas-Arroyave N,
    3. Huang L,
    4. Taylor PA,
    5. Zaghloul KA,
    6. Lungu C,
    7. Saad ZS,
    8. Horovitz SG
    (2016) DBSproc: an open source process for DBS electrode localization and tractographic analysis. Hum Brain Mapp 37:422–433. doi:10.1002/hbm.23039
    OpenUrlCrossRefPubMed
  48. ↵
    1. Lauro PM,
    2. Lee S,
    3. Ahn M,
    4. Barborica A,
    5. Asaad WF
    (2018) DBStar: an open-source tool kit for imaging analysis with patient-customized deep brain stimulation platforms. Stereotact Funct Neurosurg 96:13–21. doi:10.1159/000486645
    OpenUrlCrossRef
  49. ↵
    1. Levy R,
    2. Hutchison WD,
    3. Lozano AM,
    4. Dostrovsky JO
    (2000) High-frequency synchronization of neuronal activity in the subthalamic nucleus of parkinsonian patients with limb tremor. J Neurosci 20:7766–7775. doi:10.1523/JNEUROSCI.20-20-07766.2000
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Li X,
    2. Morgan PS,
    3. Ashburner J,
    4. Smith J,
    5. Rorden C
    (2016) The first step for neuroimaging data analysis: DICOM to NIfTI conversion. J Neurosci Methods 264:47–56. doi:10.1016/j.jneumeth.2016.03.001 pmid:26945974
    OpenUrlCrossRefPubMed
  51. ↵
    1. Lindstrom MJ,
    2. Bates DM
    (1988) Newton-Raphson and EM algorithms for linear mixed-effects models for repeated-measures data. J Am Stat Assoc 83:1014–1022. doi:10.2307/2290128
    OpenUrlCrossRef
  52. ↵
    1. Little S,
    2. Pogosyan A,
    3. Neal S,
    4. Zavala B,
    5. Zrinzo L,
    6. Hariz M,
    7. Foltynie T,
    8. Limousin P,
    9. Ashkan K,
    10. FitzGerald J,
    11. Green AL,
    12. Aziz TZ,
    13. Brown P
    (2013) Adaptive deep brain stimulation in advanced Parkinson disease. Ann Neurol 74:449–457. doi:10.1002/ana.23951
    OpenUrlCrossRefPubMed
  53. ↵
    1. Little S,
    2. Beudel M,
    3. Zrinzo L,
    4. Foltynie T,
    5. Limousin P,
    6. Hariz M,
    7. Neal S,
    8. Cheeran B,
    9. Cagnan H,
    10. Gratwicke J,
    11. Aziz TZ,
    12. Pogosyan A,
    13. Brown P
    (2016a) Bilateral adaptive deep brain stimulation is effective in Parkinson's disease. J Neurol Neurosurg Psychiatry 87:717–721. doi:10.1136/jnnp-2015-310972 pmid:26424898
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Little S,
    2. Tripoliti E,
    3. Beudel M,
    4. Pogosyan A,
    5. Cagnan H,
    6. Herz D,
    7. Bestmann S,
    8. Aziz T,
    9. Cheeran B,
    10. Zrinzo L,
    11. Hariz M,
    12. Hyam J,
    13. Limousin P,
    14. Foltynie T,
    15. Brown P
    (2016b) Adaptive deep brain stimulation for Parkinson's disease demonstrates reduced speech side effects compared to conventional stimulation in the acute setting. J Neurol Neurosurg Psychiatry 87:1388–1389. doi:10.1136/jnnp-2016-313518 pmid:27530809
    OpenUrlFREE Full Text
  55. ↵
    1. Liu Y,
    2. Coon WG,
    3. de Pesters A,
    4. Brunner P,
    5. Schalk G
    (2015) The effects of spatial filtering and artifacts on electrocorticographic signals. J Neural Eng 12:056008. doi:10.1088/1741-2560/12/5/056008 pmid:26268446
    OpenUrlCrossRefPubMed
  56. ↵
    1. Miller KJ,
    2. Leuthardt EC,
    3. Schalk G,
    4. Rao RPN,
    5. Anderson NR,
    6. Moran DW,
    7. Miller JW,
    8. Ojemann JG
    (2007) Spectral changes in cortical surface potentials during motor movement. J Neurosci 27:2424–2432. doi:10.1523/JNEUROSCI.3886-06.2007 pmid:17329441
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Naros G,
    2. Grimm F,
    3. Weiss D,
    4. Gharabaghi A
    (2018) Directional communication during movement execution interferes with tremor in Parkinson's disease. Mov Disord 33:251–261. doi:10.1002/mds.27221 pmid:29427344
    OpenUrlCrossRefPubMed
  58. ↵
    1. Nolte G,
    2. Ziehe A,
    3. Nikulin VV,
    4. Schlögl A,
    5. Krämer N,
    6. Brismar T,
    7. Müller KR
    (2008) Robustly estimating the flow direction of information in complex physical systems. Phys Rev Lett 100:234101. doi:10.1103/PhysRevLett.100.234101 pmid:18643502
    OpenUrlCrossRefPubMed
  59. ↵
    1. Opri E,
    2. Cernera S,
    3. Molina R,
    4. Eisinger RS,
    5. Cagle JN,
    6. Almeida L,
    7. Denison T,
    8. Okun MS,
    9. Foote KD,
    10. Gunduz A
    (2020) Chronic embedded cortico-thalamic closed-loop deep brain stimulation for the treatment of essential tremor. Sci Transl Med 12:eaay7680. doi:10.1126/scitranslmed.aay7680
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Piña-Fuentes D,
    2. van Dijk JMC,
    3. van Zijl JC,
    4. Moes HR,
    5. van Laar T,
    6. Oterdoom DLM,
    7. Little S,
    8. Brown P,
    9. Beudel M
    (2020) Acute effects of adaptive deep brain stimulation in Parkinson's disease. Brain Stimul 13:1507–1516. doi:10.1016/j.brs.2020.07.016
    OpenUrlCrossRef
  61. ↵
    1. Prerau MJ,
    2. Brown RE,
    3. Bianchi MT,
    4. Ellenbogen JM,
    5. Purdon PL
    (2017) Sleep neurophysiological dynamics through the lens of multitaper spectral analysis. Physiology (Bethesda) 32:60–92. doi:10.1152/physiol.00062.2015 pmid:27927806
    OpenUrlCrossRefPubMed
  62. ↵
    1. Qasim SE,
    2. de Hemptinne C,
    3. Swann NC,
    4. Miocinovic S,
    5. Ostrem JL,
    6. Starr PA
    (2016) Electrocorticography reveals beta desynchronization in the basal ganglia-cortical loop during rest tremor in Parkinson's disease. Neurobiol Dis 86:177–186. doi:10.1016/j.nbd.2015.11.023 pmid:26639855
    OpenUrlCrossRefPubMed
  63. ↵
    1. Raz A,
    2. Vaadia E,
    3. Bergman H
    (2000) Firing patterns and correlations of spontaneous discharge of pallidal neurons in the normal and the tremulous 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine vervet model of parkinsonism. J Neurosci 20:8559–8571. doi:10.1523/JNEUROSCI.20-22-08559.2000
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Reck C,
    2. Florin E,
    3. Wojtecki L,
    4. Krause H,
    5. Groiss S,
    6. Voges J,
    7. Maarouf M,
    8. Sturm V,
    9. Schnitzler A,
    10. Timmermann L
    (2009) Characterisation of tremor-associated local field potentials in the subthalamic nucleus in Parkinson's disease. Eur J Neurosci 29:599–612. doi:10.1111/j.1460-9568.2008.06597.x pmid:19187268
    OpenUrlCrossRefPubMed
  65. ↵
    1. Reck C,
    2. Himmel M,
    3. Florin E,
    4. Maarouf M,
    5. Sturm V,
    6. Wojtecki L,
    7. Schnitzler A,
    8. Fink GR,
    9. Timmermann L
    (2010) Coherence analysis of local field potentials in the subthalamic nucleus: differences in parkinsonian rest and postural tremor. Eur J Neurosci 32:1202–1214. doi:10.1111/j.1460-9568.2010.07362.x pmid:21039956
    OpenUrlCrossRefPubMed
  66. ↵
    1. Rumalla K,
    2. Smith KA,
    3. Follett KA,
    4. Nazzaro JM,
    5. Arnold PM
    (2018) Rates, causes, risk factors, and outcomes of readmission following deep brain stimulation for movement disorders: analysis of the U.S. Nationwide Readmissions Database. Clin Neurol Neurosurg 171:129–134. doi:10.1016/j.clineuro.2018.06.013 pmid:29909184
    OpenUrlCrossRefPubMed
  67. ↵
    1. Saad ZS,
    2. Reynolds RC
    (2012) SUMA. Neuroimage 62:768–773. doi:10.1016/j.neuroimage.2011.09.016 pmid:21945692
    OpenUrlCrossRefPubMed
  68. ↵
    1. Samuel M,
    2. Ceballos-Baumann AO,
    3. Blin J,
    4. Uema T,
    5. Boecker H,
    6. Passingham RE,
    7. Brooks DJ
    (1997) Evidence for lateral premotor and parietal overactivity in Parkinson's disease during sequential and bimanual movements. A PET study. Brain 120:963–976. doi:10.1093/brain/120.6.963
    OpenUrlCrossRefPubMed
  69. ↵
    1. Seabold S,
    2. Perktold J
    (2010) Statsmodels: econometric and statistical modeling with Python. Paper presented at SciPy 2010: the 9th Python in Science Conference, Austin, TX, June. doi:10.25080/Majora-92bf1922-011
    OpenUrlCrossRef
  70. ↵
    1. Shah SA,
    2. Tinkhauser G,
    3. Chen CC,
    4. Little S,
    5. Brown P
    (2018) Parkinsonian tremor detection from subthalamic nucleus local field potentials for closed-loop deep brain stimulation. In: 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society: Learning from the Past, Looking to the Future, pp 2320–2324. Piscataway, NJ: IEEE.
  71. ↵
    1. Sharott A,
    2. Gulberti A,
    3. Hamel W,
    4. Köppen JA,
    5. Münchau A,
    6. Buhmann C,
    7. Pötter-Nerger M,
    8. Westphal M,
    9. Gerloff C,
    10. Moll CKE,
    11. Engel AK
    (2018) Spatio-temporal dynamics of cortical drive to human subthalamic nucleus neurons in Parkinson's disease. Neurobiol Dis 112:49–62. doi:10.1016/j.nbd.2018.01.001 pmid:29307661
    OpenUrlCrossRefPubMed
  72. ↵
    1. Telkes I,
    2. Viswanathan A,
    3. Jimenez-Shahed J,
    4. Abosch A,
    5. Ozturk M,
    6. Gupte A,
    7. Jankovic J,
    8. Ince NF
    (2018) Local field potentials of subthalamic nucleus contain electrophysiological footprints of motor subtypes of Parkinson's disease. Proc Natl Acad Sci U S A 115:E8567–E8576. doi:10.1073/pnas.1810589115 pmid:30131429
    OpenUrlAbstract/FREE Full Text
  73. ↵
    1. Timmermann L,
    2. Gross J,
    3. Dirks M,
    4. Volkmann J,
    5. Freund HJ,
    6. Schnitzler A
    (2003) The cerebral oscillatory network of parkinsonian resting tremor. Brain 126:199–212. doi:10.1093/brain/awg022 pmid:12477707
    OpenUrlCrossRefPubMed
  74. ↵
    1. Trotta MS,
    2. Cocjin J,
    3. Whitehead E,
    4. Damera S,
    5. Wittig JH,
    6. Saad ZS,
    7. Inati SK,
    8. Zaghloul KA
    (2018) Surface based electrode localization and standardized regions of interest for intracranial EEG. Hum Brain Mapp 39:709–721. doi:10.1002/hbm.23876 pmid:29094783
    OpenUrlCrossRefPubMed
  75. ↵
    1. Velisar A,
    2. Syrkin-Nikolau J,
    3. Blumenfeld Z,
    4. Trager MH,
    5. Afzal MF,
    6. Prabhakar V,
    7. Bronte-Stewart H
    (2019) Dual threshold neural closed loop deep brain stimulation in Parkinson disease patients. Brain Stimulation 12:868–876. doi:10.1016/j.brs.2019.02.020
    OpenUrlCrossRefPubMed
  76. ↵
    1. Vinck M,
    2. van Wingerden M,
    3. Womelsdorf T,
    4. Fries P,
    5. Pennartz CMA
    (2010) The pairwise phase consistency: a bias-free measure of rhythmic neuronal synchronization. Neuroimage 51:112–122. doi:10.1016/j.neuroimage.2010.01.073 pmid:20114076
    OpenUrlCrossRefPubMed
  77. ↵
    1. Virtanen P,
    2. Gommers R,
    3. Oliphant TE,
    4. Haberland M,
    5. Reddy T,
    6. Cournapeau D,
    7. Burovski E,
    8. Peterson P,
    9. Weckesser W,
    10. Bright J,
    11. van der Walt SJ,
    12. Brett M,
    13. Wilson J,
    14. Millman KJ,
    15. Mayorov N,
    16. Nelson ARJ,
    17. Jones E,
    18. Kern R,
    19. Larson E,
    20. Carey CJ, et al.
    (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods 17:352. doi:10.1038/s41592-020-0772-5
    OpenUrlCrossRef
  78. ↵
    1. Volkmann J,
    2. Joliot M,
    3. Mogilner A,
    4. Ioannides AA,
    5. Lado F,
    6. Fazzini E,
    7. Ribary U,
    8. Llinás R
    (1996) Central motor loop oscillations in parkinsonian resting tremor revealed magnetoencephalography. Neurology 46:1359–1370. doi:10.1212/wnl.46.5.1359 pmid:8628483
    OpenUrlAbstract/FREE Full Text
  79. ↵
    1. Wang J,
    2. Johnson LA,
    3. Jensen AL,
    4. Baker KB,
    5. Molnar GF,
    6. Johnson MD,
    7. Vitek JL
    (2017) Network-wide oscillations in the parkinsonian state: alterations in neuronal activities occur in the premotor cortex in parkinsonian nonhuman primates. J Neurophysiol 117:2242–2249. doi:10.1152/jn.00011.2017
    OpenUrlCrossRefPubMed
  80. ↵
    1. Wong JK,
    2. Viswanathan VT,
    3. Nozile-Firth KS,
    4. Eisinger RS,
    5. Leone EL,
    6. Desai AM,
    7. Foote KD,
    8. Ramirez-Zamora A,
    9. Okun MS,
    10. Wagle Shukla A
    (2020) STN versus GPi deep brain stimulation for action and rest tremor in Parkinson's disease. Front Hum Neurosci 14:578615. doi:10.3389/fnhum.2020.578615
    OpenUrlCrossRef
  81. ↵
    1. Wu T,
    2. Hallett M
    (2013) The cerebellum in Parkinson's disease. Brain 136:696–709. doi:10.1093/brain/aws360 pmid:23404337
    OpenUrlCrossRefPubMed
  82. ↵
    1. Wu T,
    2. Wang L,
    3. Hallett M,
    4. Chen Y,
    5. Li K,
    6. Chan P
    (2011) Effective connectivity of brain networks during self-initiated movement in Parkinson's disease. Neuroimage 55:204–215. doi:10.1016/j.neuroimage.2010.11.074 pmid:21126588
    OpenUrlCrossRefPubMed
  83. ↵
    1. Xiao Y,
    2. Beriault S,
    3. Pike GB,
    4. Collins DL
    (2012) Multicontrast multiecho FLASH MRI for targeting the subthalamic nucleus. Magn Reson Imaging 30:627–640. doi:10.1016/j.mri.2012.02.006 pmid:22503090
    OpenUrlCrossRefPubMed
  84. ↵
    1. Xiao Y,
    2. Fonov V,
    3. Bériault S,
    4. Al Subaie F,
    5. Chakravarty MM,
    6. Sadikot AF,
    7. Pike GB,
    8. Collins DL
    (2015) Multi-contrast unbiased MRI atlas of a Parkinson's disease population. Int J Comput Assist Radiol Surg 10:329–341. doi:10.1007/s11548-014-1068-y pmid:24841147
    OpenUrlCrossRefPubMed
  85. ↵
    1. Xiao Y,
    2. Fonov V,
    3. Chakravarty MM,
    4. Beriault S,
    5. Al Subaie F,
    6. Sadikot A,
    7. Pike GB,
    8. Bertrand G,
    9. Collins DL
    (2017) A dataset of multi-contrast population-averaged brain MRI atlases of a Parkinson's disease cohort. Data Brief 12:370–379. doi:10.1016/j.dib.2017.04.013 pmid:28491942
    OpenUrlCrossRefPubMed
  86. ↵
    1. Yao L,
    2. Brown P,
    3. Shoaran M
    (2020) Improved detection of Parkinsonian resting tremor with feature engineering and Kalman filtering. Clin Neurophysiol 131:274–284. doi:10.1016/j.clinph.2019.09.021 pmid:31744673
    OpenUrlCrossRefPubMed
  87. ↵
    1. Young CK,
    2. Ruan M,
    3. McNaughton N
    (2017) A critical assessment of directed connectivity estimates with artificially imposed causality in the supramammillary-septo-hippocampal circuit. Front Syst Neurosci 11:72. doi:10.3389/fnsys.2017.00072 pmid:29033799
    OpenUrlCrossRefPubMed
  88. ↵
    1. Zach H,
    2. Dirkx M,
    3. Bloem BR,
    4. Helmich RC
    (2015) The clinical evaluation of Parkinson's tremor. J Parkinsons Dis 5:471–474. doi:10.3233/JPD-150650 pmid:26406126
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 41 (47)
Journal of Neuroscience
Vol. 41, Issue 47
24 Nov 2021
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Subthalamic–Cortical Network Reorganization during Parkinson's Tremor
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Subthalamic–Cortical Network Reorganization during Parkinson's Tremor
Peter M. Lauro, Shane Lee, Umer Akbar, Wael F. Asaad
Journal of Neuroscience 24 November 2021, 41 (47) 9844-9858; DOI: 10.1523/JNEUROSCI.0854-21.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Subthalamic–Cortical Network Reorganization during Parkinson's Tremor
Peter M. Lauro, Shane Lee, Umer Akbar, Wael F. Asaad
Journal of Neuroscience 24 November 2021, 41 (47) 9844-9858; DOI: 10.1523/JNEUROSCI.0854-21.2021
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • deep brain stimulation
  • directed functional connectivity
  • electrocorticography
  • Parkinson's disease
  • subthalamic nucleus
  • tremor

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Identification of the Acid-Sensitive Site Critical for Chloral Hydrate (CH) Activation of the Proton-Activated Chloride Channel
  • Disentangling Object Category Representations Driven by Dynamic and Static Visual Input
  • Irrelevant Threats Linger and Affect Behavior in High Anxiety
Show more Research Articles

Neurobiology of Disease

  • Early TNF-Dependent Regulation of Excitatory and Inhibitory Synapses on Striatal Direct Pathway Medium Spiny Neurons in the YAC128 Mouse Model of Huntington's Disease
  • Activation of PPARα exhibits therapeutic efficacy in a mouse model of juvenile neuronal ceroid lipofuscinosis
  • Peroxynitrite contributes to behavioral responses, increased trigeminal excitability, and changes in mitochondrial function in a preclinical model of migraine
Show more Neurobiology of Disease
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.