Abstract
Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective therapy for the motor symptoms of Parkinson's disease (PD). However, the neural elements mediating symptom relief are unclear. A previous study concluded that direct optogenetic activation of STN neurons was neither necessary nor sufficient for relief of parkinsonian symptoms. However, the kinetics of the channelrhodopsin-2 (ChR2) used for cell-specific activation are too slow to follow the high rates required for effective DBS, and thus the contribution of activation of STN neurons to the therapeutic effects of DBS remains unclear. We quantified the behavioral and neuronal effects of optogenetic STN DBS in female rats following unilateral 6-hydroxydopamine (6-OHDA) lesion using an ultrafast opsin (Chronos). Optogenetic STN DBS at 130 pulses per second (pps) reduced pathologic circling and ameliorated deficits in forelimb stepping similarly to electrical DBS, while optogenetic STN DBS with ChR2 did not produce behavioral effects. As with electrical DBS, optogenetic STN DBS exhibited a strong dependence on stimulation rate; high rates produced symptom relief while low rates were ineffective. High-rate optogenetic DBS generated both increases and decreases in firing rates of single neurons in STN, globus pallidus externa (GPe), and substantia nigra pars reticular (SNr), and disrupted β band oscillatory activity in STN and SNr. High-rate optogenetic STN DBS can indeed ameliorate parkinsonian motor symptoms through reduction of abnormal oscillatory activity in the STN-associated neural circuit, and these results highlight that the kinetic properties of opsins have a strong influence on the effects of optogenetic stimulation.
SIGNIFICANCE STATEMENT Whether STN local cells contribute to the therapeutic effects of subthalamic nucleus (STN) deep brain stimulation (DBS) in Parkinson's disease (PD) remains unclear. We re-examined the role of STN local cells in mediating the symptom-relieving effects of STN DBS using cell type-specific optogenetic stimulation with a much faster opsin, Chronos. Direct optogenetic stimulation of STN neurons was effective in treating the symptoms of parkinsonism in the 6-hydroxydopamine (6-OHDA) lesion rat. These results highlight that the kinetic properties of opsins can have a strong influence on the effects of optogenetic activation/inhibition and must be considered when employing optogenetic to study high-rate neural stimulation.
Introduction
Deep brain stimulation (DBS) is a surgical therapy for the motor symptoms of advanced Parkinson's disease (PD). The subthalamic nucleus (STN) is a common DBS target, and STN DBS produces remarkable improvements in resting tremor, rigidity, bradykinesia, and akinesia (Krack et al., 1998; Deep-Brain Stimulation for Parkinson's Disease Study Group et al., 2001). Despite well-established clinical benefits, how STN DBS exerts its therapeutic effects remains unclear, and this limits treatment optimization, as well as potential applications of DBS in other neurologic and psychiatric disorders.
Brain tissue is anatomically heterogeneous and fiber-dense, and DBS could modulate the activity of multiple neuronal elements in target sites and adjacent areas (Alexander et al., 1990; Albaugh and Shih, 2014). Given this complexity, understanding the neural elements that contribute to therapeutically effective DBS will enable improvements to this treatment. Application of cell-type-specific optogenetic DBS led to the conclusion that the therapeutic effects of STN DBS do not result from stimulation of neurons in the STN, but rather result from antidromic activation of the (hyperdirect) cortical input to STN (Gradinaru et al., 2009). However, the finding that optogenetic activation of STN neurons did not ameliorate parkinsonian symptoms, and the subsequent conclusion that antidromic activation of the hyperdirect pathway is the therapeutic target of STN DBS, are unclear because the slow response kinetics of channelrhodopsin-2 (ChR2) are incapable of maintaining high-rate firing, a critical parameter for effective STN DBS (Boyden et al., 2005; Lin et al., 2009; Yizhar et al., 2011; Jun and Cardin, 2020). During high-rate stimulation, ChR2 expressing neurons respond with irregular patterns of activity with firing rates well below the rates of optical stimulation (Boyden et al., 2005; Lin et al., 2009; Hight et al., 2015). Importantly, random patterns of electrical STN DBS are not effective at treating parkinsonian symptoms in either the 6-hydroxydopamine (6-OHDA) lesioned rat (McConnell et al., 2016) or in humans with PD (Dorval et al., 2010). Thus, it remains unclear which neural elements and circuits are essential for the therapeutic benefits of STN DBS.
Here, we used optogenetic STN DBS with an ultrafast opsin, Chronos, to re-assess the effects of selective STN DBS on behavior and neural activity. Chronos exhibits much faster kinetics than ChR2 (Klapoetke et al., 2014; Jun and Cardin, 2020) and is able to follow high-rate [>100 pulses per second (pps)] stimulation with much greater fidelity (Hight et al., 2015). Chronos was packaged into an adeno-associated virus serotype 5 (AAV5) with a calcium/calmodulin-dependent protein kinase II (CaMKII) promoter to generate local cell-specific expression in STN in the 6-OHDA lesioned rat model of PD. We quantified the behavioral effects of optogenetic DBS of the STN on parkinsonian motor symptoms using circling behavior and forelimb adjusting steps. We conducted single unit recordings in STN, substantia nigra pars reticular (SNr), and globus pallidus externa (GPe) during different rates of optogenetic STN stimulation and quantified the effects of STN DBS on neural activity. Cell-specific high-rate STN optogenetic DBS alleviated pathologic circling and improved abnormal stepping, behavioral effects that were not achieved with similarly targeted optogenetic STN DBS with ChR2, and disrupted abnormal oscillatory activity in the STN and downstream neural circuits.
Materials and Methods
All animal care and experimental procedures were approved by the Institutional Animal Care and Use Committee (IACUC) at Duke University.
Animals and surgery
Nineteen female Sprague Dawley rats (250–300 g) were used for behavioral studies and electrophysiological recordings. Five rats were injected with Chronos and implanted with optical fibers in the STN; three rats were injected with ChR2 and implanted with fibers; five rats were implanted with metal microwire electrodes in the STN for electrical stimulation. Six rats were injected with Chronos and used for acute optogenetic stimulation and recording of single units. All rats received unilateral infusion of 6-OHDA into the medial forebrain bundle (MFB) to produce hemi-parkinsonism. Surgery was performed under general anesthesia with sevoflurane (induction at 7% in 2 l/min O2; maintained at 3.5% with 1.5 l/min O2). Dexamethasone (5.0 mg/kg) was injected subcutaneously to reduce potential for brain swelling during surgery. The body temperature was maintained at ∼37°C with a water heating blanket. Craniotomies were performed over the unilateral STN and MFB according to a stereotactic atlas of the rat brain (STN: AP −3.6 mm from bregma, ML 2.6 mm; MFB: AP −2.0 mm, ML 2.0 mm) as previously described (McConnell et al., 2012).
Viral injection and unilateral dopamine depletion
For optogenetic experiments, the viral vectors rAAV5-CaMKII-Chronos-GFP (1.8 × 1013 vg/ml) and rAAV5-CaMKII-hCHR2 (H134R)-EYFP (4.6 × 1012 vg/ml) were stereotyped and packaged by the University of North Carolina at Chapel Hill Vector Core. The virus (0.6 µl) was injected into the STN (DV −7.0 mm) at a rate of 0.1 µl/min through a 1-µl 32-Gauge Hamilton syringe, waiting 5 min after injecting every 0.2 µl, and the needle was left in the brain for another 10 min after the full injection. Metal microelectrode recordings were used during the surgery to localize the STN before the injection. Rats were then lesioned to cause unilateral degeneration of dopaminergic neurons in the substantia nigra pars compacta (SNc) and produce hemi-parkinsonism. The 6-OHDA (6 µl, 2.5 mg/ml in 0.2% ascorbic acid dissolved with saline, Sigma-Aldrich) was prepared immediately before use and infused into MFB (DV −7.5 mm) through 10-µl Hamilton syringe at a rate of 1 µl/min, waiting 5 min after every 2 µl, and the needle was left in the brain for another 10 min after the full injection. Thirty minutes before lesioning, the rats were administered 50 mg/kg pargyline (intraperitoneally, Sigma-Aldrich) and 5 mg/kg desipramine (intraperitoneally, Sigma-Aldrich) to inhibit monoamine oxidase and protect noradrenergic neurons.
Fiber implantation
For rats with optogenetic stimulation, an optical fiber (200-μm core diameter, 1.25-mm OD ceramic zirconia ferrule, Precision Fiber Products) was lowered into the STN (DV ∼−6.8 mm from bregma) and cemented in place with dental acrylic. Rats were then single housed for at least five weeks to allow time for viral expression and recovery from lesion and surgery before the behavioral paradigms began. For electrophysiological recordings with optical stimulation, the surgical openings were sutured after 6-OHDA and virus injection. The rats were allowed to recover for at least five weeks before the acute recordings.
Electrode implantation
Rats for electrical STN DBS were implanted with stimulating electrodes in STN, recording electrodes in SNr and a cannula in MFB for 6-OHDA injection. A stimulating electrode array (2 × 2, platinum-iridium, 75-µm diameter, 0.3-mm interelectrode spacing, and 10-kΩ impedance; Microprobes) was implanted in the STN. A 4 × 4 recording multichannel electrode array (platinum-iridium, 75-µm diameter, 0.25-mm interelectrode spacing, ∼500-kΩ impedance; Microprobes) was implanted in the SNr (from bregma in mm, AP −5.8, ML 2.3, DV −7.2). A 23-Gauge stainless steel cannula was placed in the MFB (DV −7.5 mm). The implanted arrays and cannula were secured using dental acrylic attached to stainless steel bone screws anchored to the skull. After five weeks of recovery, rats were lesioned with 6-OHDA under 3.5% sevoflurane anesthesia. The rats were pretreated with 50 mg/kg pargyline and 5 mg/kg desipramine injection (intraperitoneal) 30 min before the lesion. Then 6-OHDA (10 µl, 2.5 mg/ml) was infused through the cannula at a rate of 2 µl/min using a 10-µl Hamilton syringe. Rats were left for recovery for at least one week before any additional measurements.
Behavioral tests
Stimulation parameters
For optogenetic experiments, after at least five weeks for viral expression and nigral neuron degeneration, rats were connected to a 473-nm DPSS laser (Shanghai Laser) via fiber optic cables and placed inside a testing chamber. Stimulation was pulsed at a fixed rate (5, 20, 75, 100, and 130 pps) with a pulse width of 1 ms for Chronos-injected rats and 5 ms for ChR2-injected rats. We used 5-ms pulses for ChR2 stimulation to be consistent with and thereby enable direct comparison with previous studies of optogenetic STN DBS with ChR2 (Gradinaru et al., 2009). Since Chronos has much faster off-kinetics and is able to follow reliably high rate, we used a shorter pulse width of 1 ms. We confirmed that optogenetic STN DBS in Chronos-injected animals (n = 3) using 5-ms pulses produced effects on circling and linear speed that were comparable to those with 1-ms pulses. The fiber output was 10 mW measured with a PM20A photodetector (Thorlabs). For electrical DBS, stimulation with charge balanced biphasic pulses (pulse width of 90 μs/phase) was conducted using an isolated voltage to current converter (A-M System) controlled by MATLAB software. The appropriate amplitude of stimulation for each rat was determined before the start of experiments based on sustained motor responses including increased contralateral turning, decreased ipsilateral turning, increased activity, and increased rearing, and a lack of motor side effects including involuntary muscle contractions of the limbs and neck during 130-pps stimulation.
Circling test
To induce robust and protracted circling, unilaterally lesioned rats received a single injection of methamphetamine (1.875 mg/kg in 0.9% saline) 30 min before placement in a cylinder (So et al., 2012). The cylinder was placed in a dark sound attenuating chamber equipped with an infrared camera to record rat behavior. A rotating optical (Doric) or electrical commutator (PlasticsOne) was used to prevent cable twisting. Behavioral effects of DBS were examined using fixed 130 pps and randomized blocks of 5 stimulation rates: 5, 20, 75, 100, and130 pps with presentation order randomized within each block. For fixed 130-pps optical stimulation, pulses were delivered in trains of 6 trials (each trial was composed of 10 s on and 20 s off) for a total of 3 min. Each behavioral recording session started with 3 min control period, followed by a 3-min stimulation period and ended with 3 min light off. The entire session was 9 min long. For rate-randomized stimulation, pulse trains were delivered for 10 s at each rate followed by a 20-s off interval. Each session started with at least 2-min control period, followed by the stimulation period, and ended with 2-min light off. For each rat, six sessions were recorded for both 130 pps and the randomized rates, and the repeated sessions were performed at least 2 d apart. For rats with electrical DBS, four blocks of five different rates of DBS were delivered, with 60-s stimulation epochs spaced 120 s apart. Circling behavior was video recorded and tracked using behavioral analysis software (Clever Systems). The position of the nose and the base of tail over time were computed. Angular velocity and distance traveled per minute (linear speed) were calculated offline from the tracking data using MATLAB. Normalized angular velocity or linear speed for each trial was calculated by dividing the angular velocity or linear speed during the stimulation-on period by the average angular velocities or linear speeds during the prestimulation-off and poststimulation-off periods immediately before and after the stimulation-on period. Normalized angular velocities or linear speed for each stimulation rate were then calculated by averaging across all randomized blocks.
Adjusting steps test
Hemi-parkinsonian rats exhibit specific deficits in contralateral limb use, and deficits in forelimb adjusting steps is a validated measure of parkinsonian akinesia in rats (Schallert et al., 1978; Glajch et al., 2012). Each rat was held with their hind limbs elevated and moved backward at a steady rate traversing a 1-m glass corridor (Runway, CleverSys) over ∼3 − 4 s. The movement was videotaped and analyzed offline by counting the number of steps made with the contralateral and ipsilateral forelimbs. Continuous optical or electrical DBS at 20- or 130-pps DBS was delivered in addition to no stimulation control trials. Two to three trials were recorded in each session, and this was repeated across four to six sessions. The behavioral effects of DBS were quantified by the ratio of the number of steps with the contralateral forelimb to the number of steps with the ipsilateral forelimb.
Electrophysiological recording
After at least five weeks of recovery following viral injection in STN, electrophysiological recordings in STN, GPe, and SNr were conducted in ipsilaterally lesioned rats under urethane anesthesia (1.2 g/kg, subcutaneous). Following a midline incision of the scalp, soft tissue was retracted, and a craniotomy was drilled over STN, GPe, and SNr. The dura was removed at each craniotomy location, and the brain was covered with normal saline. A bone screw was implanted at the back of the head and served as electrical ground. Optical stimulation and unit recording in STN were conducted using a custom made optrode with an optical fiber (200 µm) attached to single tungsten microelectrode (0.5 MΩ, 75-µm diameter) with the tip of the fiber 0.3–0.4 mm above the electrode tip and angled so the emitted light encompassed the electrode tip (Sparta et al., 2011). Sixteen channel electrode arrays were used for single unit recordings from GPe or SNr during optical stimulation in the STN. Each electrode/array was slowly advanced using a micromanipulator. Separate recording depths for each penetration were spaced 300–400 μm apart. Neuronal activity was recorded using a multichannel acquisition processor system (Plexon Inc). The data were sampled at 40 ksamples/s. Stimulation at 130 pps or five randomized rates (5, 20, 75, 100, and 130 pps) occurred in a train of 10 trials (each trial consisted of 10-s stimulation on and 20-s stimulation off). Each session began with baseline recording for at least 60 s in which no stimulation was presented and ended with baseline recording of at least 60 s.
Electrophysiological data analysis
Electrophysiological data were analyzed by a combination of custom-written MATLAB scripts. Single units were extract and sorted with the Plexon Offline Sorter using principal component analysis (Plexon Inc). Spikes were detected by applying a bandpass filter (150 Hz to 9 kHz) to raw traces with an amplitude threshold of −4*SD and 800-µs duration snippets of spike waveforms were extracted. Only single-unit activity with a clear separation from noise was used for the analysis, and units with refractory period shorter than 1 ms were removed. Neural activity within a window of −10 − 20 s aligned to stimulus onset was averaged in 200-ms bins and averaged across trials to construct peristimulus time histograms (PSTHs). The PSTHs were then Z-scored by subtracting the mean baseline firing rate (measured from a baseline window from −10 to 0 s before trial initiation) and dividing by the standard deviation. To classify “excited” versus “inhibited” units, neural activity during 10-s stimulation on was compared with a baseline window from 10 s before the stimulus onset (two-tailed t test, p < 0.05). Peak latency was computed as time elapsed between the stimulus onset and the time point within 100 ms when the response reached its first peak. Spectral analysis of the spike times was used to quantify temporal structure in the data. Multitaper analysis was used to construct spectral estimators (Chronux version 2.12; Mitra and Pesaran, 1999) on 10-s windows using 5 Slepian data tapers. Spectra are presented on a logarithmic scale. Spectral power in the β band (12–30 Hz), low β band (12–20 Hz), and high β band (20–30 Hz) during the stimulation on period were calculated and compared with the average of the spectral power during the prestimulation off and poststimulation off periods immediately before and after the stimulation on period (two-tailed t test, p < 0.05). Power spectrograms were constructed within a window of −10–20 s aligned to stimulus onset using a 1-s sliding window and 0.1-s steps.
Statistical analysis
Statistical significance between conditions was determined using a one-way or two-way repeated measures ANOVA. The normality of data were judged by frequency distributions (histogram) and boxplots (Ghasemi and Zahediasl, 2012). Post hoc paired comparisons were used if the corresponding main effect or interaction was significant at p < 0.05. To control for multiple comparisons, pairwise post hoc tests were performed with the Tukey's HSD test with a p < 0.05 significance cutoff. For single comparisons, two-sample Mann–Whitney test or paired t test was performed. All statistical tests were performed using MATLAB. All results are presented as mean ± SE.
Histology
Following completion of experiments, rats were deeply anesthetized with urethane (1.8 g/kg, i.p.) and perfused transcardially with 0.1 M PBS followed by 4% paraformaldehyde in 0.1 M PBS. Brains were postfixed in 4% paraformaldehyde overnight at 4°C and then transferred to 30% sucrose (4°C). The brains were cut into 40- or 50-µm sections in the coronal plane (except one rat with electrical DBS was horizontally sectioned) using a cryostat (CM3050S, Leica Microsystems) and processed for staining. Tyrosine hydroxylase (TH) immunohistochemistry was used to determine the extent of degeneration of dopaminergic neurons in the SNc and electrode/fiber locations. Briefly, brain sections were first rinsed and blocked for 1 h in 10% goat serum, and then incubated in anti-TH antibody (AB152; 1:1000, Vector Laboratories) overnight at 4°C in solution with 10% goat serum and 0.25% Triton X-100 with PBS. The sections were incubated with biotinylated goat anti-rabbit secondary antibody (BA-1000, 1:250, Vector Laboratories) with 10% goat serum and 0.25% Triton X-100 in PBS for 1 h at room temperature. After rinsing, the sections were incubated in a VECTASTAIN Elite ABC kit (Vector Laboratories) solution for 1 h and then visualized using DAB solution. To visualize GFP positive cells, the neighboring sections were directly mounted with DAPI-FluoroMount-G (SouthernBiotech) after rinsing. To determine the cell-specific expression of Chronos-GFP, brain sections were immunostained for CaMKII or glutamate. Slices were incubated with either primary anti-CaMKII antibody (1:300, Millipore, catalog #05-532) or primary rabbit anti-glutamate antibody (1:500, Millipore, AB133) with 10% goat serum and 0.25% Triton X-100 overnight at 4°C. Secondary antibodies (Alexa Fluor 594 goat anti-mouse IgGI) were used to visualize CaMKII or glutamate (1:500, Life Technologies). All sections were mounted with DAPI-FluoroMount-G (SouthernBiotech). Images were acquired with Nikon Eclipse TE2000-U, Zeiss AxioPlan 2 or Zeiss 780 inverted confocal microscope. The locations of electrodes, fibers, and injections were compared with a rat brain atlas (Paxinos and Watson, 2005).
Results
Histologic verification
We unilaterally injected 6-OHDA into MFB and produced unilateral loss of nigral dopaminergic cells indicated by dramatic reductions in TH immunoreactivity in the SNc of the lesioned hemisphere (Fig. 1A) in all rats used for behavioral and electrophysiological studies (n = 19). To assess the behavioral effects of optogenetic STN DBS, we expressed Chronos (n = 5) or ChR2 (n = 3) and implanted optical fibers in the STN. To compare the behavioral effects with traditional electrical DBS, we also implanted stimulating electrodes in the STN in five rats. The locations of implanted optical fibers (n = 8) and stimulating electrodes (n = 5, at least one electrode of the bipolar pair used for stimulation) were confirmed to reside within the STN (Fig. 1B,C). We targeted Chronos-GFP driven by the CaMKII promotor to neurons in STN. The virus injection labeled neurons in the STN and their projection terminals in the GPe and SNr (Fig. 2A). We did not observe Chronos-GFP labeled neurons in the motor cortex (M1) that projects to STN, excluding retrograde transmission to M1. Chronos-GFP was co-expressed with CaMKII (Fig. 2B) and glutamate (Fig. 2C) demonstrating successful and specific targeting of excitatory neurons in the STN.
Optogenetic STN DBS improved pathologic circling
We quantified the effects of optogenetic and electrical DBS on methamphetamine-induced circling behavior. Optogenetic STN DBS at 130 pps in Chronos-injected animals reduced pathologic circling behavior as indicated by decreases in angular velocity compared with the pre-light- and post-light-stimulation periods (Fig. 3A,B, left; F(2,545) = 67.1, p < 0.001, one-way repeated measures ANOVA). This was consistent with the results of electrical DBS, which was highly effective in reducing circling (Fig. 3B, middle; F(2,59) = 66.4, p < 0.001, one-way ANOVA). However, optogenetic STN DBS at 130 pps in ChR2-injected animals did not influence pathologic circling (Fig. 3B, right; F(2,305) = 1.43, p = 0.240) consistent with previous results (Gradinaru et al., 2009). The amelioration of pathologic circling during electrical DBS was significantly greater than during optogenetic DBS using Chronos (Fig. 3C; p < 0.001, Mann–Whitney test). Further, neither electrical stimulation nor optogenetic stimulation in Chronos-injected or ChR2-injected rats produced changes in linear speed (Fig. 3D,E; Chronos: F(2,545) = 0.07, p = 0.935; electrical: F(2,59) = 0.04, p = 0.960; ChR2: F(2,305) = 1.06, p = 0.346), and there were no significant differences between the three conditions (p > 0.27, Mann–Whitney test). Importantly, these data indicated that the reductions in angular velocity during DBS did not result from reductions in overall movement or freezing during stimulation.
The behavioral effects of electrical STN DBS are strongly dependent on stimulation pulse rate (McConnell et al., 2012), and we applied five different rates of optogenetic STN DBS. Chronos optical DBS and electrical DBS decreased ipsilateral circling with higher stimulation rates reflected as a decrease in angular velocity compared with the off conditions [Fig. 4A,B; two-way ANOVA with significant main effects of stimulation rate (F(4,1265) = 40.91, p < 0.001) and stimulation type (F(2,1265) = 154.48, p < 0.001) and a significant interaction between these factors (F(8,1265) = 11.91, p < 0.001)]. Post hoc paired comparisons revealed that optical STN DBS in Chronos-injected rats at high rates (75, 100, 130 pps; all p < 0.001) relieved ipsilateral turning, while low rate (5 and 20 pps) DBS was not effective. This was consistent with the effects of electrical DBS (75, 100, 130 pps; all p < 0.01), but not optogenetic DBS in ChR2 expressing animals (all p > 0.05; Fig. 4B). Meanwhile, the stimulation rate had no significant main effect on linear speed (Fig. 4C; two-way ANOVA, F(4,1265) = 0.36, p = 0.839). There was no difference between the three stimulation types (stimulation type, F(2,1265) = 2.98, p = 0.051), and there was no significant interaction between stimulation rate and stimulation type (F(8,1265) = 0.97, p = 0.457).
Optogenetic STN DBS improved forelimb stepping
Optogenetic DBS at 130 pps in Chronos-injected rats corrected the bias for using the unimpaired forepaw, as indicated by the increases in the ratio of contralateral to ipsilateral steps, whereas optogenetic DBS at 20 pps did not improve impaired forelimb use (Fig. 5, left; F(2,251) = 115.5, p < 0.001, one-way ANOVA). Similarly, electrical DBS at 130 pps increased use of the contralateral limb and corrected the stepping bias, while 20 pps electrical DBS had no effect (Fig. 5, middle; F(2,178) = 11.1, p < 0.001). In contrast, optogenetic DBS in ChR2 expressing animals did not have any effect on the use of the impaired forepaw at either 130 or 20 pps (Fig. 5, right; F(2,124) = 0.64, p = 0.527).
Changes in neural activity in STN, GPe, and SNr during 130-pps optogenetic STN DBS
Single unit neural activity was recorded from STN (n = 34 units) in six rats, from GPe (n = 38) in three rats, and from SNr (n = 61) in three rats, and recording electrode locations were confirmed by postmortem histology (Figs. 6A, 7A, 8A).
Neuron firing rates
Optogenetic DBS at 130 pps increased or decreased neural activity in the STN (Fig. 6). Approximately 53% of STN single units exhibited an increase in firing rate during light stimulation (excited, n = 18/34). These neurons rapidly increased their firing rate at stimulus onset, and firing then stabilized at a higher rate for the duration of the stimulation. Meanwhile, 32% of STN neurons exhibited a decrease in firing rate (inhibited, n = 11/34). Inhibited neurons decreased their firing rates, and spiking stabilized at lower rate for the duration of the stimulation (Fig. 6C,D). The remaining 5% of neurons (n = 5) did not exhibit significant changes in firing rate during light delivery. Similarly, neurons in GPe and SNr exhibited both increases and decreases in firing rate during 130-pps optogenetic STN DBS. The majority of neurons in GPe exhibited increased firing (53%, n = 20/38), and only four neurons (11%) exhibited decreased firing in response to light stimulation (Fig. 7). Conversely, a majority of neurons (60%, n = 37/61) recorded in the SNr decreased their firing rates during STN stimulation, and only 20% of SNr neurons (n = 12/61) increased their firing rates in response to light stimulation (Fig. 8).
We quantified response latencies of units in STN, GPe, and SNr measured from stimulus onset to the response peak for excited units and to the response trough for inhibited neurons (Fig. 9A–C). Excited units in the STN responded rapidly to light stimulation with latencies of 12.5 ms (median; Fig. 9A), while inhibited units decreased their firing rates at longer latencies of 20 ms (median; Fig. 9A; p = 0.036, Mann–Whitney test). This difference in latencies suggested that inhibitory responses in the STN may result from recruitment of inhibitory neurons downstream of light-activated excitatory neurons. The latencies of excited GPe neurons exhibited a bimodal distribution with a median latency of 37.5 ms, which was significantly longer than STN excited neurons (p < 0.001), while the four inhibited GPe neurons had latencies of 10 ms (median). Both excited and inhibited neurons in the SNr exhibited longer latencies (excited, median = 25 ms, p = 0.009; inhibited, median = 20 ms, p = 0.015) than STN excited neurons. The longer latencies of inhibited neurons in the SNr suggested that the suppression of neural responses in the SNr may result from network mechanisms involving local inhibitory recruitment and external inhibition from another input, such as GPe.
To understand further the suppression of neural activity in the STN, SNr and GPe, we conducted a correlation analysis to examine the influence of baseline firing rate on response during optogenetic DBS (Fig. 9D–F). There was a significant negative correlation between the reduction of firing rate and baseline firing rate in inhibited STN neurons (r = −0.90, p = 0.0002, Pearson's linear correlation coefficient), but no significant correlation was found in excited STN neurons (r = 0.02, p = 0.935). Similarly, the reductions in firing rate of inhibited neurons in SNr and GPe were also inversely correlated with their baseline firing rates (SNr, r = −0.60, p = 0.0001; GPe, r = −0.99, p = 0.008). The increase in firing rate of excited neurons in SNr displayed significant dependence on their baseline firing rates (r = 0.61, p = 0.036), but this was not observed for excited GPe neurons (r = −0.20, p = 0.401).
Oscillatory activity
Changes in neural oscillatory activity in basal ganglia circuits are associated with the parkinsonian state and altered by electrical DBS. We conducted spectral analysis of single unit firing times to quantify the effects of optogenetic DBS on oscillatory activity. Before stimulation, the power spectra of STN neurons (26/34, 76%) showed a pronounced peak in the β band (12 − 30 Hz). This oscillatory activity was entirely eliminated during optical stimulation at 130 pps but returned to prestimulation levels after stimulation was turned off (Fig. 10A,D). 22/26 (85%) STN neurons exhibited decreased β band oscillatory activity during 130-pps optogenetic DBS (t(25) = 2.10, p = 0.044, paired t test). In GPe, 25/38 recorded neurons (63%) exhibited increased oscillatory activity in the β band in the parkinsonian state. In contrast to STN, GPe β oscillatory power was either unchanged or shifted to a slightly higher frequency during STN optogenetic DBS (Fig. 10B,E; t(24) = −0.35, p = 0.728), and only six GPe neurons (6/25, 24%) showed deceased β band activity during stimulation. The majority of SNr neurons (52/61, 85%) displayed pronounced β band oscillatory activity before STN DBS, and this oscillatory activity was clearly decreased during optogenetic STN DBS and recovered after cessation of DBS (Fig. 10C,F); 38/52 (73%) SNr neurons exhibited a reduced oscillatory peak in the β band during 130 pps optical DBS (t(51) = 6.59, p < 0.001, paired t test).
β Oscillatory activity can be further subdivided into low-frequency (12–20 Hz) and high-frequency (20–30 Hz) components, and differential changes in the low and high β bands may have differential effects on parkinsonian symptoms. There was reduction of high β power in STN neurons during optogenetic STN DBS at 130 pps (t(25) = 2.19, p = 0.039), but no changes were detected in low β band power (t(25) = 1.09, p = 0.289). Both low and high β activity in SNr were significantly suppressed during optical STN DBS (low β, t(51) = 5.25, p < 0.001; high β, t(51) = 6.07, p < 0.001), while no significant reductions were found in GPe for either β band (low β, t(24) = 0.29, p = 0.775; high β, t(24) = 0.13, p = 0.896).
To quantify the relationship between changes in neural firing rate and changes in oscillatory activity during STN DBS, we conducted two-way ANOVA of neuronal response type based on changes in firing rates (excited, inhibited, and non-responsive) and changes in oscillatory power (prestimulation vs stimulation). There were no significant interactions between firing rate response and changes in oscillatory power in the STN (F(2,57) = 197.8, p = 0.332), SNr (F(2,121) = 106.8, p = 0.102), or GPe (F(2,75) = 45.1, p = 0.542) during optogenetic STN DBS.
Effects of optogenetic STN DBS rate on STN and SNr neuronal activity
We recorded STN (n = 36) and SNr (n = 59) neurons during different rates of optogenetic STN DBS. 12/36 STN neurons were recorded during all five rates (5, 20, 75, 100, and 130 pps), while 24 STN neurons were recorded during DBS at four rates (20, 75, 100, and 130 pps). In SNr, 43 neurons were recorded during DBS with all five rates and 16 neurons were recorded during DBS with four rates.
Neural spiking
The changes in firing rates in STN and SNr, both increases and decreases in rate, were accentuated at higher rates of optogenetic STN DBS (Fig. 11). However, population data analysis showed that there were no differences in the mean firing rates between baseline and stimulation in STN across stimulation rates (two-way ANOVA with stimulus rate and stimulation on-off as factors, rate F(4,253) = 0.23, p = 0.920, stimulation F(1,253) = 1.06, p = 0.304 and no significant interaction between these factors F(4,253) = 0.12, p = 0.977). While there was a significant reduction in mean firing rate of SNr neurons during stimulation (two-way ANOVA, significant effect with stimulation F(1,419) = 6.06, p = 0.014), there was no significant effect of stimulation rate (F(4,419) = 0.07, p = 0.992), and no significant interaction (F(4,419) = 0.22, p = 0.927). However, since DBS produced both increases and decreases in firing rate in different neurons (Figs. 7, 8), there were significant effects of stimulation rate on the absolute changes of firing rate in responsive neurons in both STN (Fig. 11C; n = 29, one-way ANOVA with stimulation rate as factor, F(4,126) = 2.97, p = 0.02) and SNr (Fig. 11D; n = 45, one-way ANOVA, F(4,209) = 9.52, p < 0.001).
Oscillatory activity
We computed the spectral power of single unit firing to quantify the effects of optogenetic STN DBS rate on oscillatory activity (Fig. 12A,B). Before stimulation, the power spectra of 19 STN neurons (19/36, 53%) and 37 SNR neurons (37/59, 63%) showed a pronounced peak in the β frequency band (12 − 30 Hz). Most of neurons in STN [5 pps, 8/12 (67%); 20 pps, 14/19 (74%); 75 pps, 15/19 (79%); 100 pps, 15/19 (79%); 130 pps14/19 (74%) and SNr (5 pps, 16/25 (64%); 20 pps, 25/37 (68%); 75 pps, 25/37 (68%); 100 pps, 25/37 (68%); 130 pps 24/37 (65%)] exhibited a significant reduction of β band power during stimulation at different rates of optogenetic DBS (paired t test, all p < 0.05). Further, higher-rate optical stimulation produced stronger reductions in neural β band oscillatory activity in both STN (one-way ANOVA, F(4,90) = 4.99, p = 0.001) and SNr (F(4,180) = 8.26, p < 0.001) in comparison to lower rates, as indicated by the ratio of β band power during stimulation to prestimulation (Fig. 12C,D).
Discussion
We quantified effects on motor behavior and neuronal activity of optogenetic STN DBS in the unilateral 6-OHDA lesioned rat using an ultrafast opsin, Chronos, that can follow the high rates required for effective DBS (Hight et al., 2015). A CaMKII promoter generated robust expression of Chronos-GFP in glutamatergic STN neurons and their projections in GPe and SNr. High-rate optogenetic STN DBS, with the kinetically-appropriate Chronos construct, alleviated ipsilateral circling bias and increased use of the impaired forelimb. As with electrical DBS, behavioral effects of optogenetic STN DBS were strongly dependent on stimulation rate. Further, we confirmed prior results that high-rate optogenetic STN DBS with kinetically-slower ChR2 was not effective in treating parkinsonism.
Optogenetic STN DBS resulted in both increases and decreases in firing rates of neurons within STN, SNr and GPe, and the absolute changes in rates were accentuated at higher rates of stimulation. β Band oscillatory activity was observed in most neurons in STN, GPe, and SNr, and high-rate optogenetic STN DBS suppressed β band oscillatory activity within the STN and SNr. This suppression was DBS rate dependent, and high rates produced stronger suppression. Changes in β power were not correlated with changes in firing rates, implying that the efficacy of DBS is due to suppression of pathologic oscillatory activity, rather than generation of specific changes in firing rate.
Multiple neural elements may be affected by DBS and mediate therapeutic effects, while others may mediate side effects. Determining which circuits contribute to effective DBS requires manipulating selectively the activity of specific neural elements. Optogenetics provides cell-type specificity and was used previously to determine the neural elements contributing to symptom relief by STN DBS in parkinsonian rodents (Gradinaru et al., 2009). That study concluded that activation of STN local cells was neither necessary nor sufficient to relieve parkinsonian motor dysfunction. However, the effects of DBS are strongly dependent on the stimulation rate, and the slow kinetics of ChR2 precluded following high-rate stimulation (>100 pps) necessary for effectiveness, as high-rate laser pulses did not result in faithful high-rate neural activity (Klapoetke et al., 2014; Hight et al., 2015; Jun and Cardin, 2020). Therefore, it was critical to determine the role of STN local cells in the symptom-relieving effects of STN DBS using faster opsins. Optogenetic activation of Chronos-expressing STN neurons at 130 pps reduced pathologic circling and improved forelimb stepping, consistent with effects of electrical DBS. Chronos exhibits higher overall driven firing rates and greater firing rate fidelity than ChR2, as well as greater sensitivity to illumination (Yizhar et al., 2011; Klapoetke et al., 2014; Hight et al., 2015; Jun and Cardin, 2020). Thus, it is unlikely that the behavioral differences between ChR2 and Chronos resulted from differences in the intensity of light required to elicit neuron activation, but rather that they resulted from the differences in the kinetic properties of the opsins. These results highlight the importance of considering the kinetic limits of opsins when probing the effects and mechanisms of DBS (Gradinaru et al., 2009; Lin et al., 2011; Sanders and Jaeger, 2016; Servello et al., 2016; Ordaz et al., 2017).
The relief of pathologic circling during electrical DBS was greater than during DBS using Chronos, which may suggest that STN local cell-specific activation may not produce complete alleviation of motor symptoms. Antidromic cortical input may be engaged by electrical DBS and enhance therapeutic efficacy, consistent with the hypothesis of DBS antidromic activation in M1 (Gradinaru et al., 2009; Li et al., 2012; Sanders and Jaeger, 2016). In addition, optogenetic responses depend on opsin expression, and increasing the opsin expression level could enhance optogenetic excitability (Meng et al., 2019) and improve behavioral outcomes. Thus, the difference between optogenetic and electrical effects on circling may be due to insufficient opsin expression and/or limited light illumination area (Yizhar et al., 2011). Differences in baseline circling across animals may reflect differences in the extent of dopaminergic lesion and/or differences in sensitivity to methamphetamine (So et al., 2012).
Changes in neural firing rates and patterns in the basal ganglia are associated with parkinsonism and are altered by DBS (Raz et al., 2000; Goldberg et al., 2002, 2004; Bar-Gad et al., 2004; Hammond et al., 2007; Erez et al., 2009; Chiken and Nambu, 2013; Little and Brown, 2014). Previous studies reported different effects of electrical DBS on neural firing rate. Some studies showed that STN DBS resulted in an increase in GPe and GPi rates (Hashimoto et al., 2003), while others reported a reduction of GPi and SNr activity (Benazzouz et al., 1995, 2000). Consistent with our results with optogenetic DBS, several studies demonstrated that STN DBS produced both increases and decreases in rates in the GPe, SNr, and striatum, and that there were no significant changes in the population mean firing rates (Shi et al., 2006; Bosch et al., 2011; McConnell et al., 2012). Like electrical stimulation, optogenetic stimulation can generate complex responses at the stimulation site, despite selective activation of specific cell types (Han et al., 2009; Galvan et al., 2012), and optogenetic STN DBS both increased and decreased firing rates of STN neurons consistent with previous findings in optogenetic studies of striatal and cortical neurons (Han et al., 2009; Galvan et al., 2012). The inhibition of STN neurons occurred immediately after light stimulation without prior activation. Further, latencies of inhibited neurons were longer than those of excited neurons, and the reduced firing rates of inhibited neurons were significantly correlated with their baseline firing rates. These observations support that the reduction of firing rates in neurons during light stimulation resulted from synaptic activation of downstream inhibitory neurons that suppressed the recorded excitatory neurons.
Optogenetic STN stimulation generated both activation and inhibition in SNr and GPe neurons, with most SNr neurons being inhibited and most GPe neurons being activated. The latencies of both inhibited and excited SNr neurons were longer than those of STN excited neurons, consistent with effects in SNr being mediated transsynaptically. While excitation of SNr neurons was most likely caused by glutamatergic input from the STN, suppression of SNr neurons could result from network mechanisms involving both local and external inhibitory inputs. Similar to inhibited neurons in STN, inhibition of SNr neurons occurred promptly on light stimulation without prior activation, and reductions in firing rate were dependent on baseline firing rates, suggesting the involvement of local inhibitory networks. Further, activation of more neurons in GPe may generate additional GABAergic inhibitory input to SNr (Windels et al., 2005), and this hypothesis supports a role of GPe-SNr connections in mediating the effects of STN DBS (Windels et al., 2000, 2003).
Changes in neural firing patterns play an important role in the pathophysiology of PD (Hamani and Lozano, 2003; Wichmann and DeLong, 2003). Dopamine depletion induces abnormal neural activity in the basal ganglia, and β oscillatory activity is implicated in the pathophysiology of PD (Hammond et al., 2007). β Band activity can be subdivided into low (13–20 Hz) and high (20–30 Hz) frequencies (Priori et al., 2004), with low β band activity more strongly correlated with rigidity and bradykinesia, and high β power elevated in patients with freezing of gait (Toledo et al., 2014). Further, β activity appears to shift to lower frequencies with increased dopamine depletion (Connolly et al., 2015), and low β activity is more sensitive than high β to dopaminergic treatment (Priori et al., 2004; López-Azcárate et al., 2010). Collectively, these results suggest that differential changes in the low and high β bands may have differential effects on parkinsonian symptoms.
Previous studies with electrical DBS indicated that therapeutic effects result from overriding pathologic neural activity and disrupting abnormal neural firing patterns (Hashimoto et al., 2003; Grill et al., 2004; Dorval et al., 2008, 2010; Rosin et al., 2011; McCairn and Turner, 2015; McConnell et al., 2016; So et al., 2017). Consistent with these findings, optogenetic STN DBS suppressed β band activity in STN and SNr, and effects were dependent on the stimulation rate. High-rate optical DBS strongly reduced β oscillatory activity, while low rates had smaller effects or were ineffective. We quantified the effects of STN optogenetic DBS on both low and high β oscillatory activity. High β activity of STN neurons was reduced during optogenetic STN DBS, while low β activity was unchanged. Meanwhile, both low and high β activity of SNr neurons was dramatically reduced, while there were no changes in the β oscillatory activity in GPe. Interestingly, the oscillatory activity of some GPe and STN neurons shifted from low β to high β during STN optogenetic DBS, and this may contribute to symptom relief. Further, the changes in oscillatory activity in STN and SNr neurons were independent of changes in firing rate, supporting the hypothesis that the therapeutic effect of DBS resulted from suppressing pathologic oscillatory activity. Cell-specific optogenetic activation of STN local neurons was thus sufficient to alter neural firing patterns within basal ganglia circuits and ameliorate motor symptoms in parkinsonian rats.
Footnotes
This work was partially supported by the National Institute of Neurological Disorders and Stroke under Award Numbers National Institutes of Health R01 NS091236 and R37 NS040894 (to W.M.G.) and R03 NS108097 (to C.Y.).
The authors declare no competing financial interests.
- Correspondence should be addressed to Warren M. Grill at warren.grill{at}duke.edu