Abstract
Hyperkinetic states are common in human movement disorders, but their neural basis remains uncertain. One such condition is dyskinesia, a serious adverse effect of medical and surgical treatment for Parkinson's disease (PD). To study this, we used a novel, totally implanted, bidirectional neural interface to obtain multisite long-term recordings. We focus our analysis on two patients with PD who experienced frequent dyskinesia and studied them both at rest and during voluntary movement. We show that dyskinesia is associated with a narrowband gamma oscillation in motor cortex between 60 and 90 Hz, a similar, though weaker, oscillation in subthalamic nucleus, and strong phase coherence between the two. Dyskinesia-related oscillations are minimally affected by voluntary movement. When dyskinesia persists during therapeutic deep brain stimulation (DBS), the peak frequency of this signal shifts to half the stimulation frequency. These findings suggest a circuit-level mechanism for the generation of dyskinesia as well as a promising control signal for closed-loop DBS.
SIGNIFICANCE STATEMENT Oscillations in brain networks link functionally related brain areas to accomplish thought and action, but this mechanism may be altered or exaggerated by disease states. Invasive recording using implanted electrodes provides a degree of spatial and temporal resolution that is ideal for analysis of network oscillations. Here we used a novel, totally implanted, bidirectional neural interface for chronic multisite brain recordings in humans with Parkinson's disease. We characterized an oscillation between cortex and subcortical modulators that is associated with a serious adverse effect of therapy for Parkinson's disease: dyskinesia. The work shows how a perturbation in oscillatory dynamics might lead to a state of excessive movement and also suggests a possible biomarker for feedback-controlled neurostimulation to treat hyperkinetic disorders.
- deep brain stimulation
- dyskinesia
- electrocorticography
- local field potentials
- motor cortex
- Parkinson's disease
Introduction
Synchronized oscillatory activity is a property of healthy brain circuits, which may be dynamically modulated to coordinate long-distance communication across distributed brain networks (Fries, 2005; Siegel et al., 2012). However, pathological changes in oscillatory activity may contribute to brain disorders. In untreated Parkinson's disease (PD), excessive synchronized activity in the beta band (13–30 Hz) may constrain neural activity into an inflexible pattern, preventing dynamic changes necessary for normal movement generation (Brown and Williams, 2005; Hammond et al., 2007; Moran et al., 2008). This excessive synchronization may be responsible for some of the cardinal symptoms of PD, such as akinesia and rigidity (Kühn et al., 2006, 2008; de Hemptinne et al., 2015).
Network oscillations might also play a role in the adverse effects caused by PD therapies. One such effect is dyskinesia. Dyskinesia is characterized by involuntary choreiform movements associated with dopaminergic medication and/or deep brain stimulation (DBS; Bastide et al., 2015). These movements may preclude optimal therapeutic benefit from medication or DBS. Local field potential (LFP) recordings from cortex and basal ganglia in a rodent model of parkinsonism have shown that levodopa-induced dyskinesias are associated with an increase in gamma oscillatory power in both structures (Halje et al., 2012). Network changes associated with dyskinesia have not been established in humans.
Invasive human data have been obtained mainly from the basal ganglia, through temporarily externalized leads in the early postoperative period (Brown and Williams, 2005; Kühn et al., 2006, 2008; Hammond et al., 2007). Narrowband gamma oscillations have been reported in human basal ganglia recordings, but the relationship to dyskinesia and the consistency of the effect are not clear (Cassidy et al., 2002; Williams et al., 2002; Alonso-Frech et al., 2006; Alegre et al., 2012; Weinberger et al., 2012; Cagnan et al., 2014). Timing and duration of dyskinesia can be unpredictable, and they are often reduced immediately following DBS surgery, limiting the utility of short-term recording strategies.
To circumvent these shortcomings, we used a totally implantable bidirectional neural interface in humans requiring DBS implantation for the treatment of motor fluctuations and medication-induced dyskinesia. This investigational device, Activa PC+S (Medtronic), has the same stimulating capabilities as standard DBS devices, but also allows long-term recording and storage of electrocorticography (ECoG) potentials and LFPs. With the view that abnormal movement in PD arises from cortical–basal ganglia interactions (Silberstein et al., 2005; de Hemptinne et al., 2013, 2015), we simultaneously recorded potentials from a permanently implanted, four-contact ECoG strip placed over motor cortex, as well as from the therapeutic subthalamic nucleus (STN) four-contact lead. We focused our analysis on data from two patients collected over 12 months. We show that dyskinesia is associated with the emergence of a narrowband gamma oscillation throughout the basal ganglia–thalamocortical motor loop, which is modulated by DBS. This oscillatory activity suggests new strategies for feedback-controlled DBS that could limit hyperkinetic adverse effects.
Materials and Methods
Consent, regulatory approvals, and patient selection.
This protocol was approved by the University of California, San Francisco institutional review board (Protocol 13-10878) under a physician-sponsored investigational device exemption (G120283). The study was registered at ClinicalTrials.gov (Identifier NCT01934296). Informed consent was obtained under the Declaration of the Principles of Helsinki. During this study, five patients with PD (two female, three male) were implanted with the Activa PC+S at our center, but only two patients (one female, one male) had a large number of recordings with and without dyskinesia, both on and off medication and on and off DBS, as well as favorable signal-to-noise ratios at gamma band frequencies. Therefore, our statistical analysis of grouped data is restricted to these two patients. A much smaller data set is available from two additional patients, so for these, only individual examples are shown. One of these patients had only rare dyskinesia after DBS implantation, and the other had a low signal-to-noise ratio for cortical recordings such that frequencies >70 Hz were near or below the noise floor of the device and not reliably detected. The final patient did not experience dyskinesia after implantation, so no data from that subject were included. Baseline characteristics of the two patients included in the statistical analysis as well as the two additional patients for whom we present example data are provided in Table 1.
Clinical characterization.
Study patients were evaluated by a movement disorders neurologist and met criteria for a diagnosis of PD (i.e., presence of bradykinesia and at least one other parkinsonian cardinal sign and responsiveness to levodopa). Baseline motor function in the on and off medication states were characterized using the Unified Parkinson's Disease Rating Scale, motor subscale (UPDRS III). Patients were evaluated by a neuropsychologist to exclude significant cognitive impairment or untreated mood disorder.
Surgery.
A quadripolar subthalamic lead (model 3389, Medtronic) was placed using frame-based stereotaxy and confirmed by microelectrode recording in the awake state using standard methods (Starr et al., 2002). Proper location in the motor territory of the STN was verified by eliciting movement-related single-cell discharge patterns (Starr et al., 2002). A quadripolar cortical ECoG lead was placed in the subdural space through the same burr hole used for the subthalamic lead. For Patient 1, we used a cylindrical lead (model 3391, Medtronic) with 3 mm contact length and 4 mm intercontact spacing; for all other patients, we used a flat lead (model 3587A25, Medtronic), with 4 mm contact diameter and 10 mm intercontact spacing. At least one contact covered the posterior precentral gyrus (presumed primary motor cortex), ∼3 cm from the midline on the medial aspect of the hand knob (Yousry et al., 1997). Adequate localization of the ECoG strip was confirmed using intraoperative CT merged to the patient's preoperative MRI, as described previously (Shahlaie et al., 2011). Functional localization of the ECoG strip was verified with somatosensory-evoked potentials, as reported previously (Crowell et al., 2012). If time permitted, a movement task described previously (de Hemptinne et al., 2015) was also performed to ensure that canonical movement-related beta decreases and broadband gamma increases were produced (Crone et al., 1998a; Crone et al., 1998b). The exiting wire from the cortical contact array was secured to the skull with a titanium miniplate. The free ends of the cortical and subthalamic leads were coiled under the ipsilateral parietal scalp.
The remaining hardware was placed under general anesthesia. In the posterior parietal area, the free ends of the cortical and subthalamic leads were connected to 40 cm lead extenders (model 37087, Medtronic), which were tunneled down the neck to a Medtronic Activa PC+S placed in a pocket over the pectoralis muscle (Fig. 1A). The Activa PC+S is identical in shape and size to the standard Activa PC. For all patients except Patient 2, a contralateral STN electrode was implanted and attached to a separate Activa SC pulse generator, for clinical purposes only, so that bilateral therapy could be delivered. Patient 2 only received unilateral therapy, as clinically indicated.
Experimental design.
Following surgery, each patient participated in research visits at regular intervals during which ECoG potentials and STN LFPs were recorded and downloaded wirelessly. These visits included recordings on and off medication in the month before initial DBS programming, followed by visits on and off of therapeutic DBS. One month after implantation, programming of the STN lead(s) was performed to achieve the best clinical result. The cortical lead was not used for stimulation.
During formal study visits, STN and cortical field potentials were recorded at rest, during an iPad reaching task that has been described previously (de Hemptinne et al., 2015), and during walking. Rest and walking recordings were 1 to 2 min long, and iPad reaching task recordings were 3–5 min long. Some additional brain recordings were initiated in our clinic to document specific phenomena, such as briefly shifting DBS frequency during an episode of dyskinesia to document a shift in the gamma peak frequency (described in Results).
Patients were also given a home data-collection triggering device (Intercept Patient Programmer, model 37441), allowing the patient to initiate brain recordings at home for a prespecified duration (1 min). Patients were instructed to initiate such recordings if they were experiencing dyskinesia, particularly severe “off” periods, or if they were feeling especially asymptomatic. They were also instructed to take notes on how they felt during each recording, detailing any notable symptoms.
Recordings.
Investigator-initiated recordings were activated via the Sensing Programmer (model 8181) and Sensing Programmer software (model 8180), which are part of the Medtronic Activa PC+S system. The Activa PC+S device allows a maximum of two time domain channels (one from each lead) to be recorded simultaneously. Thus, we typically recorded from one bipolar contact pair in motor cortex and one from the STN, with a sampling rate of 800 Hz, unless otherwise specified. To capture data as rapidly as possible from multiple cortical electrodes, we occasionally performed brief “montage recordings,” during which the device sampled data from each cortical electrode pair sequentially. Data from each electrode pair were recorded for 30 s, before proceeding to the next pair. For these recordings, a sampling rate of 422 Hz was used.
In general, the cortical contact pair used for recordings was selected based on which pair showed the clearest somatosensory evoked potential, the strongest movement-related broadband gamma response, and/or strongest beta peak at rest, at the time of initial surgical insertion. Typically, there was good overlap for these three measures in terms of identifying one or two contact pairs. When ambiguity remained, the preoperative MRI merged to the postoperative CT was used to help select an optimal motor cortex contact pair. Additional contacts were sometimes used, particularly if more than one contact pair had shown strong movement-related signals as described above.
The selection criteria for the STN recording configuration differed before and after therapeutic DBS was activated at 1 month postimplantation, since optimal recording contact pairs depended on which contacts were used for stimulation [stimulation and recording could not be performed from the same contact(s) simultaneously]. Stimulation contacts were not known before the patient's initial DBS programming. Before stimulation, recording contacts were based on which contact pair had the largest beta peak the day after surgery (off PD medications). If selection was ambiguous, the center contacts (1 and 2) were used. After stimulation, recording contacts bordering stimulation contacts on either side were used to minimize artifact. Occasionally, no STN contacts were recorded either because an unusually long cortical recording was desired or more recordings than usual were performed, precluding recording from two sites due to the limited Activa PC+S memory storage capacity.
The Activa PC+S has several built-in filters. We avoided filters when possible, so our initial recordings used minimal filters (0.5 Hz high-pass filter and a built in 260 Hz antialiasing filter) and maximum gain (2000). However, after our early patient initiated therapeutic stimulation, we realized that saturation of the amplifier was sometimes occurring, so we began using a 100 Hz low pass filter, and, in cases of especially large stimulation artifact, gain was reduced to 1000 for the STN. The Activa PC+S also has a data compression feature that we did not use, to avoid compression-associated reduction in the signal-to-noise ratio.
Data analysis.
Analyses were performed in Matlab using a combination of built-in Matlab functions, EEGlab functions (Delorme and Makeig, 2004), and custom functions/scripts. Power spectral density (PSD) calculations used the Welch method (“pwelch” in Matlab, with a window length of 512 ms and a FFT length of 1024). The log of the PSD was then taken. To calculate coherence, we first filtered both the ECoG and the LFP signals at frequencies ranging from 2 to 50 Hz, with a 2 Hz bandwidth using a two-way finite impulse response 1 (FIR1) filter (“eegfilt” in EEGlab). Complex signals were then obtained for each filtered signal by taking the Hilbert transform of the filtered signal. Coherence was calculated as the cross-spectrum of the two signals, normalized by each signal's autospectrum. For phase coherence, a similar procedure was used, but instead of considering the entire complex signal, only the phase was extracted from each filtered signal using a Hilbert transform. Then the instantaneous phase difference between the two signals was taken (accounting for the fact that phase is a circular signal). Plots that show the instantaneous phase differences are derived from these phase difference values (i.e., instantaneous phase differences for the entire signal). To obtain phase coherence, the absolute value of the average of the instantaneous phase differences was calculated (i.e., the vector length, which signifies the consistency of the phase difference between regions).
Categorization of files for group analysis.
For each patient, we separated files into groups with and without dyskinesia based on the dyskinesia rating scales (Unified Dyskinesia Rating Scale Part III; Goetz et al., 2008) obtained just before or just after the relevant data recording session by a movement disorder neurologist, who did not analyze the electrophysiological data. If the dyskinesia rating scale score was greater than zero for the contralateral arm or trunk, the file was considered “with dyskinesia,” and if ratings were zero, the file was considered “without dyskinesia.” To ensure correct categorization, we also evaluated video collected during the data recordings (acquired during many, but not all, researcher-initiated data recordings) to confirm the presence or absence of dyskinesia. Recordings during rest, the iPad reaching task, and walking were included in the analysis. Additionally, for home recordings, if the patients noted that they were experiencing either “dyskinesia” or “involuntary movement” that was not tremor, the file was considered “with dyskinesia,” and if the patient's notes did not indicate dyskinesia, it was considered “without dyskinesia.”
Recordings where the DBS settings were changed during the recording (wash-in or washout or a change in DBS parameters) were excluded from the group analysis. Likewise, recordings from cortical contact pairs that never showed the narrowband gamma oscillation described in Results (regardless of dyskinesia status) were presumed to be insensitive detectors and were excluded. (However, we conducted a control analysis including these recordings and results did not differ.) Additionally, home recordings where the patient did not take notes or where notes were not clearly associated with a specific recording were excluded. Finally “montage recordings,” where data were sampled at 422 Hz (instead of 800 Hz), were excluded.
Measurement of gamma oscillation parameters.
For each recording, the maximum cortical or STN PSD peak between 62–83 Hz was extracted, and a peak “height” was calculating by subtracting the log PSD from the average log PSDs 5 Hz above and below the frequency of the peak. [To determine the frequency of the peak, a normalization procedure was performed. First, the spectrum was approximately flattened by fitting it to a fourth-order polynomial, excluding artifacts, to correct for the 1/f pattern (Pritchard, 1992). Then, the entire spectrum was normalized by subtracting the mean and dividing by the standard deviation across all frequencies. This method was used only to identify the frequency of the peak. The height of the peak at this frequency was calculated based on the original PSD as described above.] The method of quantifying gamma peak height is illustrated in Figure 2B. The width of the peak at half the height was also calculated.
Coherence and phase coherence did not require a correction for the 1/f decrease because coherence values did not exhibit the strong 1/f decrease in amplitude present in the PSD (Fig. 2A). Additionally, because coherence is internally normalized, raw coherence values were used for the statistical analyses rather than first calculating the peak height relative to neighboring frequencies.
Measurement of beta oscillation parameters.
Due to the importance of beta band synchronization in the motor system and in PD in particular (Crone et al., 1998b; Hammond et al., 2007), we also examined our four measures of interest (motor cortex PSD, STN PSD, coherence, and phase coherence) in the beta frequency range (13–30 Hz). Log PSD values between 13 and 30 Hz were averaged to derive the PSD measures. Coherence values were derived in the same manner as described for the gamma oscillation analysis. Note that because of DBS artifacts in the beta range (at folded subharmonics), we included only recordings made in the absence of stimulation for this analysis.
Statistics.
The values for all four variables of interest (motor cortex PSD, STN PSD, coherence, and phase coherence) were compared between groups (dyskinesia vs no dyskinesia) using a two-tailed, nonparametric Wilcoxon rank sum test. We tested this for each of the two patients separately (Fig. 3) as well grouped together (Fig. 2C). To determine whether the stimulation artifact may be driving any results, we also ran the same analyzes including only recordings with DBS off (Fig. 4). To determine how well the height of the gamma PSD peak or the coherence values would work as classifiers for dyskinesia and to characterize the specificity and sensitivity of the biomarkers, we derived receiver operating characteristic (ROC) curves (Zwieg and Campbell, 1993). These were derived using the “perfcurve” function in Matlab (Fig. 2D). Values were first fit to a logistical regression model (“fitglm”) using a binomial distribution.
Results
Recording locations, signal characteristics, and data overview
Recording sites in cortex and basal ganglia and example recordings from each site are shown in Figure 1, B and C. Average root mean square voltage for recordings with DBS off, were 6.1 μV for STN and 19 μV for motor cortex. Signal amplitudes for recordings from both cortex and STN were stable for the 12 months of the study (Fig. 1D). The number of recordings used for group analyses are reported in Table 2. For each recording, the longest artifact-free segment of data was used for analysis. The mean duration of recordings (excluding artifacts) was 96.5 s (SD, 72 s; range, 28–302 s). There was no difference in the recording length for recordings with and without dyskinesia (p > 0.88). One-hundred seven of the recordings were initiated by study personnel during formal study visits at regular intervals, while fifty-two were initiated at home by the patients. Dyskinesia occurred almost exclusively in the on medication state, but could be present on or off of DBS.
Gamma oscillations in cortex and STN are associated with dyskinesia
Visual inspection of PSD plots for each cortical recording revealed that many of those collected during episodes of dyskinesia showed a discrete, narrowband peak (local maximum at peak frequencies between 62 and 83 Hz). A similar but smaller peak was often present in the STN, which had strong coherence with motor cortex (both magnitude and phase) in the same frequency range (Fig. 2A). To characterize the relationship of these oscillatory phenomena to the presence of dyskinesia, we quantified the oscillation (Fig. 2B; see Materials and Methods) and segregated its amplitude by dyskinesia status. Results are presented in Figure 2C, with statistical details in Table 2. All four putative biomarkers (motor cortex and STN PSD, coherence, and phase coherence for the local maximum between 62 and 83 Hz) distinguished the dyskinetic state from the nondyskinetic state, with the strongest distinction for cortical PSD and cortex–STN phase coherence. Separation of biomarker amplitudes by dyskinesia status was also present when data from each patient were analyzed independently, with the sole exception of coherence in Patient 1 (which was at the trend level, p = 0.08; Fig. 3A, Table 2), showing the effects were not driven by a single patient. Separation of biomarker amplitudes also persisted when data were restricted to those with DBS off (Fig. 4, Table 3), demonstrating that the effects were not driven by stimulation artifacts, which could have interfered with accurate peak height analysis in the gamma frequency range.
ROC curves are provided for each measure in Figure 2D, and the area under the curve and sensitivities and specificities for detection of the dyskinetic state are listed in Table 2. Consistent with the statistical evaluation above, the gamma oscillation amplitude in motor cortex and gamma phase coherence between motor cortex and STN were better classifiers for the presence or absence of dyskinesia than the gamma oscillation amplitude in STN or coherence between motor cortex and STN. Since dyskinesia occurred almost exclusively in the “on medication state,” it is important to dissociate a medication effect from a dyskinesia effect. To show that the oscillatory biomarkers are more closely related to the dyskinetic state than the dopaminergic state (on vs off of dopaminergic medications), we also performed an ROC analysis for medication state (Fig. 5). This analysis included the same data shown in Figure 2D, but grouped files according to medication status. For recordings in-clinic, “off medication” recordings were after withholding medications for at least 12 h. Recordings that the patients initiated at home were considered “on medication” unless patients indicated that they had been off medications overnight, since patients were maintaining their regular medication regimen. This analysis showed that the area under the curve for the ROC analysis derived by medication status is lower for all four gamma oscillation derived biomarkers, than when using the biomarkers to predict the presence or absence of dyskinesia.
The fact that phase coherence distinguished the two clinical states to a greater degree than coherence suggests that the hyperkinetic state is more closely related to the phase relationship between motor cortex and basal ganglia gamma oscillations than to amplitude correlations between them. To visualize this effect, we pooled the instantaneous phase angles between motor cortex and STN for all time points, all recordings, and both patients. We showed a consistent phase difference during dyskinesia, but not without dyskinesia (Fig. 6A). The mean phase angle during dyskinesia was 71°. This phase relationship in the dyskinetic state persisted regardless of stimulation status (on, or off; Fig. 6B), demonstrating that the presence of stimulation artifact in those files recorded with stimulation on was not responsible for the observed consistency of the phase relationship between cortex and STN.
We also examined whether the height of the narrowband gamma oscillation correlated with dyskinesia severity. While we did observe a correlation (r = 0.63, p = 3.58 × 10−13; Fig. 7), this effect was driven most strongly by presence versus absence of dyskinesia, rather than by its severity. This is consistent with the observation in rodents that the relationship between gamma power and dyskinesia severity is sigmoidal (Halje et al., 2012).
Beta oscillations in cortex and STN and association with dyskinesia
In contrast to the gamma oscillation, in the beta range, neither motor cortex PSD, STN PSD, nor coherence distinguished the dyskinetic from the nondyskinetic states (Fig. 8). There was a significant difference in phase coherence (p = 7.24 × 10−4), with lower beta phase coherence between STN and motor cortex when dyskinesia was present compared to when it was absent (the opposite pattern that we observed for the gamma oscillation.) Of note we also repeated the beta range analysis with recordings categorized by medication (as performed in Fig. 5). For this analysis we also found a significant difference in phase coherence (p = 0.012) with lower phase coherence between motor cortex and STN when patients were on medications compared to off. We also observed lower STN PSD when patients were on medications compared to off (p = 0.0356), as has been observed in previous studies (Priori et al., 2004; Kühn et al., 2006).
Relationship between gamma oscillation characteristics in cortex vs STN during dyskinesia
The relationship between gamma oscillatory characteristics in cortex versus STN, for DBS off recordings, is described in Table 4. Peaks were of greater amplitude in cortex compared to STN. The peak frequency varied between recordings, but their distributions in the two patients studied overlapped (Fig. 9A). Of note, the exact frequency of the gamma peak may be related to the time of recording relative to levodopa dose (Halje et al., 2012). There was a significant correlation (Pearson's r = 0.597, p = 0.0016) between the frequency of the gamma peak in cortex and STN (Fig. 9B). Although only two subjects in the study (1 and 2) had a sufficient number of recordings in the dyskinetic state for statistical analyses (explained further in methods), very similar frequency and amplitude characteristics of the dyskinesia-associated cortical PSD peak were seen in two other subjects who had a small number of high signal-to-noise ratio recordings in the dyskinetic state (Fig. 10). The fifth study subject had no dyskinesia after DBS implantation and did not have a narrowband gamma oscillation in cortex or STN.
The gamma oscillation and voluntary movement
To be of optimal utility as a driver of closed loop deep brain stimulation, a biomarker of abnormal movement should not be strongly affected by normal movement. To assess the effect of voluntary movement on gamma oscillation derived biomarkers, we conducted an analysis similar to that in Figure 2B, but separated recordings that included voluntary movement [walking or performing an arm movement (iPad) task] from those without voluntary movement (Fig. 11). In the dyskinetic state there were no significant differences in the gamma oscillation amplitude for any measure between recordings obtained with and without voluntary movement. In contrast, the difference in the gamma oscillation for all comparisons with and without dyskinesia was highly significant.
In sensorimotor cortex, a broadband power increase is known to be associated with movement (Crone et al., 1998a; Miller et al., 2007); however, this gamma increase is over a much broader frequency range (usually 50–250 Hz) and likely has a different etiology than the narrowband peak we see here (Manning et al., 2009). Additionally, the movement-related broadband power change is small in amplitude compared to the height of the peak associated with dyskinesia (Fig. 12). This point is developed further in Discussion.
Localization of cortical gamma oscillation
Precise localization of cortical contact(s) over the region of the strongest gamma oscillatory signal will be important for clinical use of a cortical detector in feedback controlled DBS. In Patient 2, we recorded files from all possible cortical contact pairs (montage recordings, see Materials and Methods) during dyskinesia and observed that the dyskinesia-associated narrowband gamma power increase was spatially specific to one contact common to all bipolar pairs showing a strong gamma oscillation (contact 11, anterior commissure–posterior commissure coordinates: −38.02, −6.59, 64.21). This contact localized to the anterior part of precentral gyrus, extending over the precentral sulcus (Fig. 13).
STN stimulation entrains cortical gamma at half the stimulation frequency
STN stimulation is known to be able to suppress or exacerbate dyskinesia in PD, depending on exact stimulation parameters and contact locations (Zheng et al., 2010; Oyama et al., 2012). Both Patients 1 and 2 often experienced dyskinesia both on and off stimulation. When dyskinesia occurred in the on-stimulation state, the cortical gamma peak always occurred at half the stimulation frequency (Fig. 14A). This is unlikely to be due to stimulation artifact because it was not present when stimulation was delivered at the same settings, but in the absence of dyskinesia (Fig. 14B). Moreover, the distribution of phase angles between STN and motor cortex in the dyskinetic state, at half stimulation frequency, were similar for stimulation-on recordings compared to stimulation-off recordings. In contrast, a more narrow (presumed artifactual) phase angle distribution is observed at the actual frequency of stimulation and at folded subharmonics of stimulation frequency (Fig. 14D). These findings are consistent with reports of “partial entrainment” of neuronal discharge by STN DBS at therapeutic frequencies (Garcia et al., 2003; Hashimoto et al., 2003; Li et al., 2012; Agnesi et al., 2015). Our findings suggest that when this entrainment occurs in the frequency range associated with dyskinesia, DBS cannot suppress the dyskinetic state. As a proof of principle, we showed that the frequency at which the peak occurred could be moved by changing the stimulation frequency, in this example from 65 Hz with 130 Hz stimulation to 75 Hz with 150 Hz stimulation, (Fig. 14C), and peak shifts during stimulation wash in and wash out (Fig. 14E,F). These changes in the peak frequency of the gamma oscillation did not produce observable changes in the clinical phenomenology of the dyskinetic movements.
Discussion
We studied circuit mechanisms of dyskinesia using a novel, totally implanted chronic multisite brain-recording device in humans with PD. Motor cortex ECoG and STN LFPs were analyzed in two patients over 1 year, revealing a narrowband network oscillation between 60 and 90 Hz. This oscillation is closely associated with dyskinesia. When dyskinesia is present during DBS, this oscillation is entrained at a subharmonic of the stimulation frequency. These findings provide a mechanistic basis for a debilitating adverse effect of therapy in PD and suggest algorithms for feedback-controlled neurostimulation.
Narrowband gamma rhythms vs broadband gamma activity
The narrowband 60–90 Hz rhythm studied here should be distinguished from cortical “broadband gamma,” which is a wide-band phenomenon (typically 50–200 Hz) that tracks local activation and likely reflects underlying spiking activity rather than a narrowband oscillatory rhythm (Ray et al., 2008; Manning et al., 2009; Scheffer-Teixeira et al., 2013). Broadband gamma activity may play a role in the pathophysiology of the off-medication parkinsonian state, where akinesia is prominent and dyskinesia is absent. In this state, there is elevated resting broadband gamma in motor cortex (Crowell et al., 2012; Rowland et al., 2015) as well as excessive coupling of broadband gamma to the phase of the beta rhythm (de Hemptinne et al., 2013). Broadband gamma is also modulated by voluntary movement (Crone et al., 1998a; Miller et al., 2007). In contrast, the narrowband gamma oscillation elucidated in the present work is likely a manifestation of neuronal oscillations, is relatively unchanged by voluntary movements, and is associated with a hyperkinetic, rather than a bradykinetic, state. For illustration of this distinction see Figure 12. We use the terms “oscillation” or “rhythm” to make the distinction between narrowband phenomena versus arrhythmic broadband phenomena.
Functional role of gamma oscillations and relationship to hyperkinetic movements
Previous work has shown narrowband gamma oscillations to be a normal feature of cortical function, implicated in numerous cognitive processes, which may act to “bind” cortical regions together based on their phase relationships (Fries, 2009; Sohal, 2012). Alterations in cortical gamma rhythms have been implicated in diseases including autism and schizophrenia (Sohal, 2012). Altering the balance of excitatory and inhibitory activity in the cortex regulates these rhythms (Yizhar et al., 2011). Interactions of inhibitory fast-spiking interneurons play a critical role (Sohal, 2012; Salkoff et al., 2015), and these generators may be regulated by neuromodulators including acetylcholine (Teles-Grilo Ruivo and Mellor, 2013) and serotonin (Puig et al., 2010).
In subcortical regions, invasive human LFP recordings from temporarily externalized DBS leads have detected gamma rhythms in STN (Cassidy et al., 2002; Alonso-Frech et al., 2006; Trottenberg et al., 2006), globus pallidus (Weinberger et al., 2012), and thalamus in both PD patients and in patients with several nonparkinsonian conditions (Kempf et al., 2009). In PD patients, the subcortical gamma oscillation is present mainly in the on-medication state. However, the relationship to dyskinesia has been unclear. A combined STN LFP/magnetoencephalography study in PD detected a cortical gamma rhythm and STN–cortical coherence at the onset of voluntary movement, supporting a prokinetic role for the gamma oscillation (Litvak et al., 2012). Here, we are able to characterize a gamma oscillation that propagates through motor cortex and basal ganglia and is strongly associated with dyskinesia. Given the presence of gamma oscillations in the thalamus in several nondyskinetic conditions (Kempf et al., 2009), it is possible that dyskinesia arises when subcortical gamma rhythms are excessively propagated throughout the basal ganglia–thalamocortical loop, with a prominent representation in motor cortex. In support of this view, motor cortex recordings in a rodent model of parkinsonism showed a gamma oscillation during dyskinesia remarkably similar to that reported here (Halje et al., 2012). Together, these finding support the hypothesis that oscillations play a fundamental role in brain network dynamics and that alternations of these oscillations may manifest as disease (Voytek and Knight, 2015).
The mechanism by which gamma oscillations lead to dyskinesia is speculative. Oscillations can bias the probability of spike discharge, such that neuronal spiking tends to occur at a preferred oscillatory phase. Gamma oscillations tend to synchronize cortical neuronal pools so that common inputs to a cell arrive in close temporal succession, facilitating activation (Fries et al., 2001; Fries, 2009). In primary motor cortex, “fragments of movement” appear to be encoded by small groups of neurons, as demonstrated by the induction of complex movements by microstimulation (Graziano et al., 2002; Hatsopoulos et al., 2007). Thus, the coordination of these neuronal pools by an exaggerated gamma oscillation could release locally encoded fragments in rapid progression to produce choreiform activity.
Striatal mechanisms of levodopa-induced dyskinesia: possible links to gamma oscillations
The cellular origin of levodopa-induced dyskinesia remains controversial, but most current theories, developed in rodent models, emphasize changes in striatal microcircuitry induced by dopamine denervation followed by its unregulated restoration by medication (Fieblinger and Cenci, 2015). Pulsatile dopamine release may be mediated by serotonergic neurons (Carta et al., 2007). Striatal changes include morphological alternations in glutamatergic synapses onto striatal medium spiny neurons and concomitant changes in long-term potentiation and long-term depression, which strengthen corticostriatal connections (Fieblinger and Cenci, 2015). These changes might favor activation of striatal cells originating the prokinetic “direct” intrinsic basal ganglia pathway (Thiele et al., 2014).
A challenge in this field is relating the changes in striatal physiology to oscillatory phenomenon. We and others have proposed that striatal changes in the parkinsonian off state have the effect of reducing the basal ganglia “filter” of cortical activity, such that normal rhythms, including the motor beta rhythm, are excessively transmitted through the basal ganglia–thalamocortical loop, resulting in aberrant beta synchronization (Weinberger and Dostrovsky, 2011; de Hemptinne et al., 2013). A similar mechanism could underlie dyskinesia, again due to exaggerated “propagation” of an otherwise normal physiological rhythm through abnormally strengthened corticostriatal synapses.
Entrainment of gamma oscillations by stimulation
Clinically, dyskinesia is not only associated with levodopa in PD, but may be induced by DBS in both PD (Yelnik et al., 2000) and non-PD disorders (Mouton et al., 2006; Mallet et al., 2008; Ostrem et al., 2011). The entrainment of gamma rhythms by DBS (Fig. 14) offers a potential explanation for this. DBS is often delivered at 120–180 Hz, approximately twice the typical frequency of the dyskinesia-associated gamma oscillation. DBS entrains axonal activity to stimulus pulses, but this entrainment does not occur after every pulse (Li et al., 2012). Frequent failures of entrainment (in this case, every other stimulation pulse) could readily result in driving axonal or neuronal activity at 60–90 Hz. Our results suggest that irregular stimulation paradigms may be less prodyskinetic than the constant frequency stimulation that is used currently.
Control signals for closed-loop DBS
The finding of a brain rhythm reliably associated with dyskinesia has translational potential as a control signal in closed-loop DBS (Rosin et al., 2011; Little et al., 2013). Commercially available DBS systems are “open-loop” devices that do not respond to patients' symptom fluctuation. Stimulation-induced dyskinesia can be an important dose-limiting effect of DBS (Mouton et al., 2006; Mallet et al., 2008; Ostrem et al., 2011). To address this, a system could be designed that sensed the gamma rhythm and dynamically adjusted stimulation parameters to keep motor cortical spectral power in the relevant bandwidth below a specified level, mitigating hyperkinetic adverse effects. The cortical sensor provides a signal that is easily separable from normal cortical oscillatory activity and has less stimulation artifact than signals detected from the DBS lead. Control signals for closed-loop DBS related to akinesia, based on activity in the beta band, are under exploration (Little et al., 2013). However, beta activity is strongly affected by voluntary movement (Crone et al., 1998b; de Hemptinne et al., 2015), which may present a challenge for its use in closed-loop DBS.
Limitations
Determining the specificity of the gamma oscillation for dyskinesia depends on the accuracy of dyskinesia detection. In-clinic scoring was performed by a movement disorders neurologist, but the relative insensitivity of the scoring system and difficulty of continuously assessing dyskinesia during brain recording precluded precise examination of the timing of neural activity relative to changes in dyskinesia severity. Since dyskinesia is not always present during study visits, we chose to increase our number of recordings by allowing patients to trigger recordings at home. A cognitively intact PD patient with a history of dyskinesia is capable of scoring dyskinesia as present or absent, but we do not have independent verification of dyskinesia for home recordings. Additionally, because we used an implant without microrecording capability, our study focused on network rhythms recorded from macroelectrodes and did not address the relationship between dyskinesia and single unit discharge. However, it is likely that gamma oscillations exert a strong effect on neuronal discharge by entrainment to a preferred phase of the gamma rhythm (Trottenberg et al., 2006; Halje et al., 2012).
Conclusions
Using multisite chronic brain recordings in humans with PD, we characterize the association of a gamma oscillation in the basal ganglia–cortical motor network with dyskinesia. These findings illuminate the network dynamics underlying dyskinesia and suggest a strategy for feedback-controlled DBS.
Footnotes
This study was supported by the University of California President's Postdoctoral Fellowship (N.C.S.), the National Institute of Neurological Disorders and Stroke (R01 NS090913-01; P.A.S.), and the Michael J. Fox Foundation (C.D.H.). We thank Dr. Alexandra Nelson for her helpful comments and critical review of this manuscript and Andrew Miller for his help with the cortical reconstructions in Figure 1B.
University of California, San Francisco has filed a preliminary patent application based on results from this manuscript. N.C.S., C.D.H., J.L.O., and P.A.S. are coinventors on this patent. Devices and technical support were provided by Medtronic Inc. under a research agreement.
- Correspondence should be addressed to Dr. Nicole C. Swann, University of San Francisco, 513 Parnassus Health Science East 8, Room 823, San Francisco, CA 94143. Nicole.Swann{at}ucsf.edu