Abstract
Deep brain stimulation (DBS) of the subthalamic nucleus is a symptomatic treatment of Parkinson’s disease but benefits only to a minority of patients due to stringent eligibility criteria. To investigate new targets for less invasive therapies, we aimed at elucidating key mechanisms supporting deep brain stimulation efficiency. Here, using in vivo electrophysiology, optogenetics, behavioral tasks and mathematical modeling, we found that subthalamic stimulation normalizes pathological hyperactivity of motor cortex pyramidal cells, while concurrently activating somatostatin and inhibiting parvalbumin interneurons. In vivo opto-activation of cortical somatostatin interneurons alleviates motor symptoms in a parkinsonian mouse model. A computational model highlights that a decrease in pyramidal neuron activity induced by DBS or by a stimulation of cortical somatostatin interneurons can restore information processing capabilities. Overall, these results demonstrate that activation of cortical somatostatin interneurons may constitute a less invasive alternative than subthalamic stimulation.
Similar content being viewed by others
Introduction
Parkinson’s disease results from the neurodegeneration of the nigro-striatal dopaminergic neurons. The main symptomatic treatment for Parkinson’s disease consists in substituting lacking dopamine with levodopa and/or dopaminergic agonists, but after a typical “honeymoon” period with dopaminergic therapy, patients inevitably develop motor complications1. At this stage, deep brain stimulation at high frequency of the subthalamic nucleus (DBS) constitutes to date the most efficient symptomatic treatment2,3. However, due to its surgical invasiveness and strict eligibility criteria, DBS benefits only to a minority of patients (~5–10%). Several hypotheses have been proposed to explain the beneficial effects of DBS4,5,6. Notably, a growing body of evidence points towards a cortical effect of DBS in both parkinsonian rodent models7,8,9,10,11 and patients12,13,14,15,16,17. Here, we reasoned that mimicking the cortical effects of DBS should reproduce its therapeutic benefits, thus paving the way for less invasive approaches.
For this purpose, we (i) determined DBS effects on cortical cell-type specific populations using a combination of in vivo electrophysiological and optogenetic approaches, (ii) reproduced these effects using optogenetics in freely-moving parkinsonian mice, and (iii) explored mathematically how DBS and DBS-guided optogenetics could restore cortical information processing capabilities. We showed that DBS normalized pathological hyperactivity of motor cortex pyramidal cells, while concurrently inhibiting parvalbumin (PV)- and activating somatostatin (SST)-expressing GABAergic interneurons. Furthermore, reproducing these effects by direct opto-activation of cortical SST interneurons alleviates motor symptoms in a parkinsonian mouse model. Lastly, our computational model shows that the dampening of the firing activity of pyramidal cells by DBS and DBS-guided optogenetics restores cortical information processing capabilities. Overall, these results establish that cortical SST interneurons constitute a promising target for a less invasive alternative to DBS.
Results
DBS decreases pathological hyperactivity of pyramidal cells
To understand how DBS affects cortical activity, we first aimed at depicting the electrophysiological signature of Parkinson’s disease at the neuronal level in the primary motor cortex (M1). We performed in vivo single unit juxtacellular recordings of pyramidal neurons in deep layers of M1 in a rat model of Parkinson’s disease (Fig. 1a and Supplementary Fig. 1a) by unilateral stereotaxic injection of 6-hydroxydopamine (6-OHDA) in the substantia nigra pars compacta (SNc). In anesthetized 6-OHDA-lesioned rats, we observed an increase in the spontaneous firing rate of M1 pyramidal cells (n = 41) compared to sham animals (n = 36) (p = 0.0014) (Fig. 1b), in line with previous observations10,18. We then explored the effect of DBS (STN stimulation parameters: 2–4 V, 60 μs at 130 Hz during 2 min) on this pathophysiological hyperactivity of M1 pyramidal neurons. In 6-OHDA-lesioned rats, we found that the increased firing activity was diminished by DBS (p = 0.0311, n = 20) back to physiological firing rates (Fig. 1b), with 68% of pyramidal cells inhibited by DBS (Fig. 1c). DBS also decreased M1 neuron firing rate in sham animals (p = 0.0266, n = 19), and the proportion of inhibited neurons, as well as the change in firing rate, were similar in sham and 6-OHDA-lesioned rats (p = 1 and p = 0.7894, respectively) (Supplementary Fig. 1b, c). Overall, DBS diminished the firing rate of pyramidal cells in sham rats and normalized their pathophysiological hyperactivity in parkinsonian rats.
GABAergic circuits mediate DBS inhibition of pyramidal cells
We next investigated the mechanistic underpinnings related to the decreased activity of pyramidal cells under DBS. For this purpose, we performed single-cell in vivo intracellular recordings of electrophysiologically identified M1 pyramidal cells in anesthetized rats (Fig. 2a). Spontaneous firing of pyramidal cells was determined before and during 150 s of DBS. We further confirmed that DBS decreased the spontaneous firing rate (p = 0.0002, n = 20) of M1 pyramidal neurons (Fig. 2b). This decrease was accompanied by a hyperpolarization of −3.7 ± 1.0 mV of their membrane potential (p = 0.0007, n = 20) (Fig. 2b), and a decrease in their membrane time constant (p = 0.0390, n = 18), input resistance (p = 0.0497, n = 19) and AP threshold (p = 0.0263, n = 20), without affecting their Ih current (p = 0.2336, n = 17). In a subset of cells (n = 3) that were antidromically activated by STN stimulation, the antidromically-evoked action potentials were rapidly shunted and a marked hyperpolarization was observed (Supplementary Fig. 1d). We further characterized DBS-induced changes in pyramidal cell excitability by applying successive depolarizing current steps before and during continuous DBS (Fig. 2c). DBS induced a decrease in depolarization-evoked activity, associated with an increased rheobase (p = 0.0312, n = 17), without affecting the gain (p = 0.1746, n = 17) of pyramidal cell f–I curve. Overall, the decrease of several properties of pyramidal cell excitability under DBS could participate in the decrease of their firing activity.
To determine the evoked conductances in M1 pyramidal cells, we next applied single STN stimulations (2–4 V, 60 μs) (Fig. 2d, e). We observed a post-stimulation hyperpolarization, sufficient to delay the evoked firing activity (p = 0.0057, n = 6) (Fig. 2d). Analysis of the voltage-dependency showed that the early phase of the evoked postsynaptic responses reversed at −71.7 ± 2.7 mV (n = 9) (Fig. 2e), which corresponds to the chloride reversal potential19. These results suggest that DBS recruits GABAergic circuits responsible for the hyperpolarization of pyramidal cells and the increase in rheobase, leading to the reduced activity observed in both in vivo juxtacellular and intracellular recordings of M1 pyramidal cells.
Somatostatin interneurons are activated by DBS
In M1, GABAergic inhibition is provided by local neuronal populations mainly composed of PV and SST interneurons20. To identify the cell-type specific populations recruited by DBS, we used genetically modified mice expressing channelrhodopsin (ChR2) in either PV (Pv::ChR2 mice) or SST (Sst::ChR2 mice) interneurons (Fig. 3, Supplementary Figs. 2 and 3). We ensured that ChR2 was expressed in the targeted populations with minimal non-specific expression (Supplementary Fig. 2).
In mice, stereotaxic targeting of the STN for DBS is challenging. To ensure the proper placement of the stimulation electrode, we first lowered a microelectrode to electrophysiologically localize the STN based on typical STN neuron firing (Supplementary Fig. 3a), and in a second step implanted the DBS electrode at the same coordinates (Supplementary Fig. 3b). We then performed in vivo juxtacellular recordings of M1 neurons in anesthetized Pv::ChR2 and Sst::ChR2 mice to monitor DBS-evoked responses in opto-identified neuronal subpopulations (Figs. 3a, b and Supplementary 3c, d). Namely, once the DBS electrode was inserted in STN, an optical fiber placed on top of M1 shone light (100 ms at 0.5 Hz) and opto-responsive neurons were detected by a recording microelectrode lowered within M1. PV and SST interneurons were distinguished from pyramidal neurons by their responses to light in Pv::ChR2 and Sst::ChR2 mice, and by post-hoc clustering of their waveform characteristics based on principal component analysis (Supplementary Fig. 3c, d). A subset of PV, SST, and pyramidal neurons were juxtacellularly labeled with neurobiotin for immunohistochemical and morphological identification and used as ground truths for the principal component analysis (Fig. 3c and Supplementary Fig. 3c–e). Most (75%) of the recorded neurons were located in M1 layer V (n = 72 neurons subjected to DBS; Supplementary Fig. 3e). In agreement with in vivo recordings in anaesthetized mice21, SST interneurons exhibited a lower spontaneous activity compared to PV cells (p = 0.0002, PV and SST median firing rates were 1.23 and 0.05 Hz, respectively; Fig. 3d). Consistent with our observations in rats, DBS (60 µs, 120 µA at 130 Hz) decreased the firing activity of pyramidal cells in the transgenic mouse lines (p = 0.0297, n = 28; Fig. 3d–f). We then investigated DBS effect on PV interneurons. Surprisingly, DBS decreased their activity (p = 0.0074, n = 28; Fig. 3d–f). We next examined SST interneurons and observed that DBS caused the opposite effect, i.e. increasing SST firing (p = 0.0040, n = 26; Fig. 3d–f). These effects were robust across cell populations as 64% of the SST cells were excited and 70% of the PV cells were inhibited (Fig. 3f). The activity of pyramidal, PV, and SST cells was stable under DBS, and reversible (with a shorter delay after DBS offset for PV cells than for SST and pyramidal neurons). Interestingly, the kinetics of DBS effects on pyramidal cells mirrored those of SST interneurons: we observed similar delays between DBS onset and its effects on pyramidal and SST cells, and also at DBS offset (Fig. 3e). Consistently, the modulation of pyramidal cells was strongly correlated with SST interneuron activity during and after DBS (r = −0.91, p = 1.0 × 10−9). Therefore, DBS efficiently drives activation of M1 SST interneurons, concurrently inhibiting M1 pyramidal cells.
Activity of PV and SST cells in sham and parkinsonian mice
Since DBS both normalizes pyramidal cell hyperactivity and differentially recruits PV and SST interneurons, we further explored the firing rate modulation of interneurons in parkinsonian conditions. Indeed, hyperactivity of pyramidal cells observed in 6-OHDA-lesioned animals could be caused by a decreased activity of GABAergic interneurons. We performed in vivo juxtacellular recordings of neurons opto-identified and clustered by principal component analysis as PV or SST cells in anesthetized sham (PV cells = 27 and SST cells = 26) and 6-OHDA-lesioned (PV cells = 14 and SST cells = 15) Pv::ChR2 and Sst::ChR2 mice (Fig. 4a, b). The spontaneous firing activity of PV cells was similar in sham and parkinsonian Pv::ChR2 mice (p = 0.8391 with n = 32 in sham and n = 14 in 6-OHDA-lesioned mice). Similarly, the spontaneous firing rate of SST cells was not different in sham and 6-OHDA-lesioned mice (p = 0.2112 with n = 32 in sham and n = 15 in 6-OHDA-lesioned mice). Pyramidal hyperactivity in parkinsonian models is not likely due to changes in electrophysiological activity of cortical GABAergic populations.
Inhibitory evoked-responses by PV or SST opto-activation
To investigate the synaptic weight of the inhibitory inputs onto pyramidal cells, we performed in vivo patch-clamp whole-cell recordings in anesthetized Pv::ChR2 and Sst::ChR2 mice (Fig. 4c–f). We recorded pyramidal cells upon opto-activation of PV or SST interneurons. First, we characterized the evoked-PSP following single opto-stimulation (3–20 ms). The opto-activation of PV and SST cells evoked PSPs of similar amplitudes in pyramidal cells, regardless of the membrane potential (p = 0.3397, n = 8 pyramidal cells in Pv::ChR2 vs. 6 in Sst::ChR2 held at −95/−80 mV; p = 0.8632, n = 7 vs. 7 held at −65/–50 mV and p = 0.8838, n = 7 vs. 9 held at −40/−30 mV) (Fig. 4c). We analyzed the voltage-dependency of opto-PSPs, which reversed at −71.5 mV in Pv::ChR2 (n = 6 pyramidal cells) and −75.2 mV in Sst::ChR2 (n = 4) mice (p = 0.0542) (Fig. 4c), close to the calculated chloride reversal, −71.9 mV. Single-pulse opto-activation of PV (n = 7) and SST (n = 6) cells induced PSPs whose area, amplitude, rise time and delay to peak increased with increasing opto-pulse duration (2-way repeated-measures ANOVA; Fig. 4d). There was no difference between PSPs evoked by PV or SST opto-activation, except for the delay to peak amplitude which was shorter upon PV than SST opto-activation (F1,11 = 66.39, p < 0.0001). We then investigated the effects of opto-activation of PV or SST cells at 67 Hz for 5 s (see “Methods” and Supplementary Fig. 4a–d for opto-stimulation frequency selection) (Fig. 4e–i). We observed that opto-activation of PV or SST cells induced a hyperpolarization of similar peak amplitude (p = 0.2270, Fig. 4f), with a shorter delay upon PV activation (p = 0.0124; Fig. 4g) that strongly inhibited the firing activity of pyramidal cells (p = 0.0392, n = 5 in Pv::ChR2 mice, p = 0.001, n = 8 in Sst::ChR2 mice) (Fig. 4e). However, in the later part of the pulse, opto-activation of SST cells induced a less pronounced membrane hyperpolarization than opto-activation of PV cells (p = 0.0371; Fig. 4h), together with a partial release of the spiking inhibition only in Sst::ChR2 mice (p = 0.0326; Fig. 4e). In addition, the variance of the membrane potential was larger under opto-activation of SST than PV cells (p = 0.0371; Fig. 4i), denoting a strong shunting of synaptic activity reaching the soma by PV opto-activation, whereas some synaptic activity subsisted (leading eventually to spikes) under SST opto-activation.
M1 SST cell opto-activation alleviates parkinsonian symptoms
Since DBS activates SST interneurons, we hypothesized that cortical SST interneurons could constitute a target to efficiently mimic the effects of DBS. Therefore, we tested whether the direct activation of GABAergic interneurons would improve motor symptoms in freely moving parkinsonian mice. We first compared the respective effects of SST and PV opto-activation at 67 Hz (see “Methods” and Supplementary Fig. 4a–d for opto-stimulation frequency selection). For this purpose, unilaterally 6-OHDA-lesioned mice (Sst::ChR2, Pv::ChR2 and wild-type mice) were ipsilaterally implanted with an optical fiber in M1 (Fig. 5a). We monitored the impact of opto-stimulation on the asymmetrical locomotor behavior induced by unilateral 6-OHDA-lesioning in three different tasks: the open field (Fig. 5b–g), cylinder test (Fig. 5h–m) and cross-maze (Fig. 5n–s).
Opto-activation of SST interneurons in 6-OHDA-lesioned Sst::ChR2 mice successfully decreased their asymmetrical locomotor behavior. Indeed, in the open field, opto-activation of SST cells in 6-OHDA-lesioned Sst::ChR2 mice decreased spontaneous ipsilateral rotations (n = 11 mice; rotations/min: p = 0.001; Fig. 5c and normalized rotations: p = 0.0002; Supplementary Fig. 4e). This decrease in asymmetrical behavior did not result from a decreased locomotor activity, which remained unaffected by opto-activation of SST cells (p = 0.5721; Fig. 5c). In the cylinder test, the opto-activation of SST cells decreased the asymmetry in front paw usage preference (p = 0.0156, n = 9; Fig. 5i). In the cross-maze, 6-OHDA-lesioned Sst::ChR2 mice (n = 11) exhibited a strong bias towards ipsilateral turns (Fig. 5o), as expected for hemi-parkinsonian rodents, while sham-lesioned Sst::ChR2 mice (n = 5) displayed no turn preference (Supplementary Fig. 4g). During SST opto-activation, a decrease in the ipsilateral bias (p = 0.0019) and an increase in the straight choice (p = 0.0022; Fig. 5o), were observed.
PV interneuron opto-activation only partially reproduced the improvement in motor symptoms induced by SST opto-activation. In 6-OHDA-lesioned Pv::ChR2 mice, opto-activation decreased spontaneous ipsilateral rotations (n = 12; rotations/min: p = 0.0123; Fig. 5d and normalized rotations: p = 0.0125; Supplementary Fig. 4e) without affecting locomotor activity (p = 0.0737; Fig. 5d). However, opto-activation of PV cells did not decrease the asymmetrical locomotor behavior in the cylinder test (p = 0.9139, n = 9; Fig. 5j) or in the cross-maze task (n = 12; p = 0.1876 for ipsilateral, p = 0.3419 for straight, and p = 0.0960 for contralateral turn preference; Fig. 5p).
Interestingly, opto-activation of cortical interneurons in Sst::ChR2 or Pv::ChR2 sham-mice did not induce an asymmetrical behavior contralateral to the opto-activation neither in the open field (Supplementary Fig. 4f) nor in the cross-maze task (Supplementary Fig. 4g). This indicates that the reduced asymmetry observed upon SST opto-activation in 6-OHDA-lesioned mice is due to an improvement of pathological symptoms rather than a generic effect of unilateral opto-activation of cortical interneurons.
We also ensured that opto-stimulation in wild-type (non-opsin expressing) 6-OHDA-lesioned mice, did not affect spontaneous ipsilateral rotations (n = 11; rotations/min: p = 0.5591; Fig. 5e and normalized rotations: p = 0.2540; Supplementary Fig. 4e) or locomotor activity (p = 0.6289; Fig. 5e) in the open field, asymmetry in front paw usage preference in the cylinder test (p = 0.3841, n = 9; Fig. 5k) and asymmetrical locomotor behavior in the cross-maze (n = 11; p = 0.8292 for ipsilateral, p = 0.9755 for straight, and p = 0.3409 for contralateral turn preference; Fig. 5q).
We next evaluated the effects of DBS and levodopa treatments in wild-type 6-OHDA-lesioned mice in open-field, cylinder and cross-maze tasks, for comparison with opto-activation of SST and PV cells. In the open-field, DBS did not reduce the spontaneous ipsilateral rotations (n = 12; rotations/min: p = 0.2082; Fig. 5f and normalized rotations: p = 0.2456; Supplementary Fig. 4e) but increased the locomotor activity (p = 0.0438; Fig. 5f), in line with a recent report22. DBS did not modify asymmetrical locomotor behavior for the cylinder test (n = 12, p = 0.9768; Fig. 5l) or in the cross-maze (n = 12; p = 0.2040 for ipsilateral, p = 0.7675 for straight, and p = 0.3388 for contralateral turn preference; Fig. 5r). Levodopa treatment (6 mg/kg) in 6-OHDA-lesioned mice efficiently reversed the spontaneous ipsilateral rotations into contralateral rotations (n = 10; rotations/min: p = 0.0004; Fig. 5g and normalized rotations: p = 2.0 × 10−5; Supplementary Fig. 4e) and increased the locomotor activity (p = 0.0135) in the open-field. The asymmetry in front paw usage preference was reversed to contralateral preference (n = 10, p = 0.0002; Fig. 5m) in the cylinder test, as well as the ipsilateral bias (n = 7, p = 0.0029) in favor of contralateral (p = 0.0029) turn preference in the cross-maze (Fig. 5s).
Therefore, specific opto-activation of cortical SST interneurons alleviates the asymmetrical behavior while electrical DBS increases locomotor activity in parkinsonian mice.
M1 L5 model recapitulates DBS effects on cortical firing
To better understand the beneficial effects of DBS and SST opto-activation, we theoretically tested the information processing capabilities of cortical networks in control and parkinsonian conditions under DBS and when activating PV or SST interneurons. To this purpose, we built a simplified spiking neural network model of L5 motor cortex, including pyramidal, PV, and SST cells modeled as adaptive exponential integrate-and-fire neurons (Fig. 6a). The choices of intrinsic parameters defining each population (Supplementary Table 1), network architecture and synaptic strengths (Fig. 6a) were guided by experimental23,24,25,26,27,28,29 and modeling30 data. In particular, we modeled PV and SST interneurons as fast-spiking and low-threshold spiking neurons, respectively. As for the network architecture, we set lower conductances for excitatory synapses compared to inhibitory ones31,32. Both PV and SST interneurons inhibited pyramidal cells, with similar connection probability but a stronger synaptic weight from PV to pyramidal cells25,26. Importantly, only PV cells received excitatory feedback24. The choice of asymmetric inhibitory strength between PV and SST neurons, as well as the respectively weak and strong self-inhibition for SST and PV cells are based on previous reports26,31. The hyperexcitability of pyramidal cells in the parkinsonian condition10 was considered by decreasing their firing threshold. Finally, DBS-mediated currents were modeled as instantaneous excitatory inputs to all cells, with no preconception of the pathway involved to reach each population (directly via antidromic or orthodromic connections from STN to cortex, or indirectly through basal ganglia loops), thus representing the sum of DBS network effects.
Consistent with our experimental findings, we observed a decrease in the firing rate of both pyramidal and PV cells while SST interneurons were activated at all tested frequencies and the impact was stronger with increasing stimulation frequency (p < 0.001; Fig. 6a). This phenomenon remained robust to changes in the parameters of 130 Hz DBS stimulation (pulse duration and amplitude; Supplementary Fig. 5a) and linearly scaled with the overall amount of current injected in the network (Supplementary Fig. 5b). The impact of DBS on network activity was almost immediate (15–20 ms lag) since DBS acts as an instantaneous current in our simplified model.
We next evaluated the relative contributions of PV and SST interneurons for driving changes in pyramidal cell activity under DBS. We found a larger decrease of pyramidal cell activity in response to DBS-induced stimulations of SST compared to PV interneurons (Fig. 6b and Supplementary Fig. 5c). In particular, no decrease in pyramidal cell activity was observed in the absence of DBS-induced current on SST interneurons (Fig. 6b and Supplementary Fig. 5c). We hypothesized that the asymmetrical efficacy of cortical interneurons in reducing pyramidal cell firing under DBS could be due to their asymmetrical connectivity with pyramidal cells. Indeed, in the motor cortex, L5 SST interneurons engage in feedforward inhibition, receiving no or little excitation from L5 pyramidal cells24, in contrast to PV interneurons. This hypothesis was validated when evaluating the impact of adding an excitatory feedback from pyramidal cells to SST neurons on pyramidal cell firing: once the strength of this connection became comparable to the excitatory drive from pyramidal cells to PV neurons, the modulation of pyramidal cell activity under DBS was no longer detectable (Fig. 6c).
Finally, we modeled the opto-activation of cortical interneurons by adding an external current to half of PV or SST neurons, constituted by a series of pulses repeated at the same frequency as those used experimentally. These simulated opto-activations of PV and SST interneurons at 67 and 130 Hz both led to a remarkable inhibition of pyramidal cell activity (Fig. 6d), yielding a very sparse residual activity of pyramidal cells, as experimentally observed (Supplementary Fig. 4c, d).
Cortical information processing under DBS or opto-stimulation
Based on our simplified network model, we investigated how pyramidal cell activity reduction, by DBS or SST opto-activation, could account for the improved motor function. The impact of the modifications of the network activity on information transmission in these conditions is not heuristically obvious: in Parkinson’s disease, the increased excitability of pyramidal cells could make the system more responsive to stimuli, but the higher spontaneous activity may interfere with stimulus-evoked activity. In contrast, the lower spontaneous activity induced by DBS and DBS-guided opto-activation may avoid signal degradation, but could bring the system to a less responsive state. We first tested whether pyramidal cells conserve the capacity to generate specific patterns of activity in response to stimuli despite sparse activity under DBS or opto-activation. We applied 28 inputs with different motifs (deterministic constant and ramping stimuli, and stochastic Ornstein–Uhlenbeck processes; Supplementary Table 2) and intensities to subsets of pyramidal cells (Fig. 7a, b), and systematically evaluated the information conveyed by the network. We found that for medium to high amplitude stimuli, pyramidal cell responses were less correlated with the stimulus in parkinsonian condition compared to control. Both DBS and SST opto-activation increased the correlation coefficients when compared to parkinsonian condition, with smoother response profiles to each stimulus (Fig. 7b, c and Supplementary Fig. 6). PV opto-activation led to weak correlation coefficients (Fig. 7b, c and Supplementary Fig. 6). Importantly, the differential impact of PV and SST opto-activation remained consistent over a large range of current intensities. The reduced efficacy of PV opto-activation was not only due to the higher silencing level of pyramidal cells but also to the induction of oscillatory-like responses (Supplementary Fig. 7a, b). Finally, low amplitude stimuli yielded weak correlations that were not significantly different across conditions and thus did not confer any advantage to the parkinsonian condition (Fig. 7c).
We then tested how accurately M1 activity patterns could be decoded by downstream areas, such as the striatum. We used four machine-learning decoding algorithms (nearest centroid classifier, multinomial logistic regression, linear discriminant analysis and support vector machines; Supplementary Table 3) to discriminate network responses to 20 different stimuli (Fig. 7a, d). Binned spiking patterns from all pyramidal cells were densified using a random matrix to mimic highly convergent cortico-striatal information transfer and the resulting matrices were used as inputs to the classifiers. The classifier accuracy at determining the stimulus identity was higher in control, DBS and SST opto-activation conditions relative to the parkinsonian condition; this result was consistent across the four classifiers (p < 0.001; Fig. 7d) and the efficacy of cortical interneurons opto-activation broke down for very high stimulation intensities, once pyramidal cells were too sparsely responsive to each stimulus (Supplementary Fig. 7c).
Thus, our modeling results indicate that reducing pyramidal cell hyperactivity by DBS or preferentially by SST opto-activation enables the cortical network to efficiently encode stimulus content and may facilitate the extraction of this information by downstream areas. These results are robust to variations in the stimulation parameters (Supplementary Fig. 7) and choices of network parameters. Importantly, the asymmetrical connectivity of SST interneurons with L5 pyramidal cells24 accounts for their higher efficiency at inhibiting pyramidal cells under DBS, when all three populations receive an external periodic pulse-shaped excitatory current; furthermore, this asymmetrical connectivity has important consequences in structuring the temporal dynamics of pyramidal cells under SST opto-activation as it does not impose an oscillatory-like pattern as PV opto-activation does in our model.
Discussion
In Parkinson’s disease, although DBS has been widely described as acting on distinct basal ganglia nuclei4,5,6, recent studies have also pointed towards an impact of DBS at the cortical level7,8,11,12,13,14,15,16,17. Here, we show that M1 pyramidal cells display a hyperactivity in parkinsonian rodents, which is counteracted by DBS. Our in vivo intracellular and juxtacellular recordings indicate that DBS recruits cortical GABAergic networks. Furthermore, opto-activation of cortical interneurons, and in particular SST interneurons, alleviates key motor symptoms in parkinsonian mice. Our theoretical modeling reveals that both DBS and SST interneuron opto-activation increase information processing capacity by counteracting the cortical hyperactivity, without interfering with the network temporal dynamics.
As reported here and in agreement with data collected from patients12,13,14,15,16,17, hyperactivity of pyramidal neurons may be a neuropathological hallmark in Parkinson’s disease. It likely results from abnormal basal ganglia dynamics as well as substantial dopaminergic denervation of the motor cortex possibly altering the excitation-inhibition balance17. It has been suggested that the therapeutic actions of levodopa may arise from dampening M1 hyperactivity, since it reduces M1 glutamate levels and increases the inhibitory tone in 6-OHDA-lesioned rats, and reduces M1 blood-oxygenation and glucose metabolism in parkinsonian patients13,17.
We report that DBS reduces pyramidal hyperactivity and exerts a dualistic action on cortical GABAergic interneurons, namely activation of SST and inhibition of PV interneurons. This agrees with a study showing that in a genetic parkinsonian rodent model, DBS reduced the abnormal M1 activity, hypothesized to result from the activation of cortical GABAergic interneurons33. The correlated time courses of the decreased hyperactivity of pyramidal cells and SST activation upon DBS combined with our modeling result suggest that the main beneficial effects of DBS on normalization of cortical activity and processing capabilities seem to be conveyed by the activation of SST cells. Several pathways could account for SST recruitment by DBS: antidromic activation of the cortico-subthalamic fibers, orthodromic activation of subthalamo-cortical fibers, and/or basal ganglia-thalamo-cortical loops4,11. The activation of SST interneurons by antidromic axonal reflex is unlikely to occur since cortico-STN pyramidal cells lack collaterals to SST interneurons24. The orthodromic STN-cortex pathway could recruit SST neurons, since it mainly targets superficial cortical layers34, where pyramidal cells directly connect to L5 SST interneurons35. Concurrently, cortical effects of DBS could occur via recruitment of other basal ganglia-thalamo-cortical loops. Other interneurons might be engaged during DBS, such as VIP cells36, since their inhibition would disinhibit SST cells. The concurrent recruitment of these anatomical pathways likely results in a new cortico-basal-ganglia wide state of equilibrium. This DBS-induced state stabilizes over time, which would explain the observed onset and offset delays in pyramidal and SST responses to DBS.
We mimicked DBS effects by opto-activating M1 SST interneurons and revealed this was sufficient to alleviate parkinsonian symptoms. Interestingly, both in experimental and modeling data, while opto-activation of PV interneurons was more efficient than SST to inhibit pyramidal cell activity, the beneficial outcome on motor symptoms and cortical processing was weaker compared to SST opto-activation. Recent studies have highlighted the key roles operated by SST interneurons in cortical information processing35,37,38,39. SST opto-activation in the somatosensory cortex ameliorates symptoms in a mouse model of neuropathic pain40, and improves affective state discrimination in the prefrontal cortex41. Under DBS, the need for SST interneuron activation to efficiently reduce cortical hyperactivity stems from their asymmetrical connectivity with pyramidal cells (Fig. 6b, c). Indeed, within L5 motor cortex, PV interneurons are involved in strong feedback loops with pyramidal neurons, while SST interneurons, in part recruited by L2/3 pyramidal neurons, provide disynaptic feedforward inhibition24. Furthermore, this asymmetrical connectivity has important consequences on the temporal structuring of the network activity patterns. Additional mechanisms could further reinforce this effect. The position of SST terminals onto pyramidal dendritic branches, by leading to a direct modulation of excitatory inputs, may explain why SST activation improves information processing compared to PV, which exerts a massive shunting (Fig. 4e, i) through somatic-targeted inhibition37,42,43. Overall, SST interneurons could thereby be more efficient than PV cells at exerting a sustainable inhibitory influence by filtering excitatory inputs onto pyramidal neurons (rather than simply silencing pyramidal cell output by shunting the inputs), thus resulting in an ameliorated performance of motor functions. A more detailed model, including a multi-compartmental description of pyramidal cells as well as short-term plasticity properties, would allow exploring how these differences in PV and SST physiology could account for the different behavioral outcomes of their opto-activation.
We compared the behavioral effects of SST opto-activation with typical treatments of parkinsonian symptoms: DBS and levodopa. In our unilateral model of parkinsonism, levodopa has strong effects, in particular on the asymmetrical symptoms, where the ipsilateral bias is not only reduced but reversed, thus creating a strong contralateral bias. This is in line with previous studies44 and thought to result from the supersensitivity of ERK signaling in the 6-OHDA-lesioned hemisphere45. Our DBS results mirror a recent study in 6-OHDA-lesioned mice22, reporting improved locomotion with no effect on rotational bias. Surprisingly, DBS in mice does not recapitulate the typical motor improvements observed in rats8,46,47, although results in rats also vary depending on the intensity and the duration of stimulation48,49,50. Several studies, including ours, choose an intensity below the dyskinetic/dystonia threshold, as used in clinical DBS. In mice, dyskinesia is reported in the absence of effect on rotational bias over a large range of parameters studied22. The small size of the mouse STN could make it difficult to stimulate efficiently enough to reduce asymmetry without deleterious impact on the pyramidal tract51. Importantly, opto-stimulating cortical interneurons avoids this issue, thus explaining why SST opto-activation was efficient on asymmetrical behavior. Yet unaffected hypolocomotion under opto-activation suggests that non-cortical effects of DBS (i.e. direct modulation of basal ganglia structures4,5,6, unlikely to be mimicked by cortical opto-activation) also contribute to the beneficial locomotor effects. Alternatively, the magnitude or spatial extent of DBS cortical effects (milder stimulation of SST neurons reaching potentially larger regions of frontal cortex than single fiber opto-stimulation) could contribute to thedifferent behavioral outcomes observed. Importantly, M1 SST opto-activation decreases asymmetrical behavior in 6-OHDA-lesioned mice but has no effect in sham animals. This argues for a therapeutic effect, rather than a net sum of two unrelated biases (pathological ipsilateral bias, and physiological contralateral bias from the inhibition of ipsilateral M1). Yet, the similar firing of interneurons in (anesthetized) sham and 6-OHDA animals suggests that this therapeutic effect does not stem from a pathological hypoactivity of interneurons in parkinsonian conditions, though, anesthesia effects on SST firing patterns could occlude differences between sham and lesioned animals. Similarly, our data pointing to SST recruitment by DBS was obtained in anaesthetized rodents. While we can reasonably hypothesize that in the absence of the dampening effect of anesthesia, DBS-activation of SST interneurons would be even higher, hence exerting an even stronger inhibition onto the hyperactive pyramidal cells, this remains to be confirmed.
The methodology used for testing the information processing capabilities of the network (Fig. 7) could be applied to experimental recordings: differences in the correlation between specific motor patterns and large-scale in vivo M1 single-unit activities as well in the responses reliability over multiple trials could be found by comparing sham, 6-OHDA lesioned mice, with or without DBS or SST opto-activation. Moreover, since the segregation of sensorimotor maps is degraded in the basal ganglia in the parkinsonian condition52,53, neuronal decoding algorithms trained on a subset of trials in which different types of movements are initiated could test our model prediction, i.e. that DBS or SST opto-activation would restore neuronal selectivity by decreasing cortical hyperactivity.
Overall, our study reveals the cortical SST interneurons as a promising therapeutic target for Parkinson’s disease. Therefore, increasing the inhibitory drive into the motor cortex could represent a useful strategy to improve motor symptoms. While there is an emerging therapeutic potential of optogenetics54,55,56,57,58,59,60, its medical application is still in its infancy, with the first human clinical trial underway to treat retinitis pigmentosa (clinicaltrials.gov, NCT03326336 and NCT02556736). More work is necessary to fully explore and optimize safe transfection and light delivery in humans58,60,61. Nevertheless, targeting cortical GABAergic networks with pharmacology or non-invasive brain stimulation61 such as transcranial stimulation62 could provide less invasive strategies than DBS, thus benefiting a larger population of parkinsonian patients.
Methods
Animals
All experiments were performed in accordance with the guidelines of the local Ethics Committee (CEEA-59) and EU directive (2010/63/EU). Adult male OFA rats (n = 31, 175–200 g) (Charles River, L’Arbresle, France) were used for electrophysiology. Adult mice of both sexes (n = 158, 2–12 months) were used for unitary extracellular recordings combined with opto-tagging, behavior, and immunostaining: wild type C57Bl6 (Charles River), and hybrid transgenic Pv::ChR2 and Sst::ChR2. The hybrid transgenic mice were heterozygous for both genes, obtained by mating a cre-driver transgenic line ensuring specific expression of the cre recombinase in PV or SST neurons (PV-cre: Jackson Laboratory stock #008069 or SST-cre: Jackson Laboratory stock #013044, Charles River) with a transgenic reporter line containing channelrhodopsin-2(H134R) (ChR2: Jackson Laboratory stock #012569, Charles River). Animals were housed in an approved animal facility under standard 12-hour light/dark cycles, with food and water available ad libitum, and nesting materials provided.
Rodent model of Parkinson’s disease: 6-OHDA lesions
Rats were anesthetized with pentobarbital (30 mg/kg, i.p.), and adjusted to a surgical plane with ketamine (27.5 mg/kg, i.m.; Imalgène, Mérial, Lyon, France). Mice were anesthetized with an i.p. mix of ketamine (67 mg/kg) and xylazine (13 mg/kg, Rompun, Bayer, Puteaux, France). All animals received a bolus of desipramine (Tocris, Bio-Techne, Lille, France) dissolved in saline (rats: 25 mg/kg, i.p.; mice: 20 mg/kg, i.p.) 30 min before the injection of 6-OHDA (Sigma, Steinheim, Germany) to prevent neurotoxin-induced damage of noradrenergic neurons. For rats, 7 µL of vehicle (saline containing 0.01% w/v ascorbic acid, Sigma) or 6-OHDA (2.5 mg/ml in vehicle) was injected at 16 µL/h in the substantia nigra pars compacta (from the interaural line: AP 3.7 mm, ML 2.1 mm, 7.55 mm depth from the cortical surface). For mice, 0.5 µL of vehicle (saline containing 0.02% w/v ascorbic acid) or 0.5 µL of 6-OHDA (6 mg/ml in vehicle) was injected at 3 µL/h in the right medial forebrain bundle (from bregma: AP 1.2 mm, ML 1.2 mm, DV 4.9 mm). Mice who failed to display spontaneous ipsilateral rotations (i.e., were seen to explore indifferently towards ipsi- and contralateral sides when observed in their homecage and on the scale during daily weighting, without amphetamine challenge) and weight loss in post-lesion recovery were excluded from the study.
In vivo intracellular electrophysiological recordings
After 3 weeks of recovery, we performed intracellular recordings in 6-OHDA-lesioned and sham rats maintained in narcotized and sedated states with fentanyl (4 µg/kg, i.p.; Janssen-Cilag, Issy-Les-Moulineaux, France), immobilized with gallamine triethiodide (40 mg, i.m., every 2 h; Specia, Paris, France) and artificially ventilated (UMV-03, UNO, Zevenaar, Netherlands). Craniotomies were drilled above M1 (from the interaural line: AP 12.5 mm; ML 3.8 mm) and STN (AP 5.2 mm; ML 2.5 mm). Body temperature was maintained at 36.5 °C with a homeothermic blanket. Intracellular recordings were performed using glass micropipettes filled with 2 M potassium acetate (40–70 MΩ) and the active bridge mode of an Axoclamp 2B amplifier (Molecular Devices, Union City, CA). Data were sampled at 25 kHz via a CED 1401 interface using the Spike2 data acquisition program (Cambridge Electronic Design, Cambridge, UK). The STN ipsilateral to the recorded pyramidal cells was stimulated with a bipolar electrode (SNE-100; Rhodes Medical Instruments, Woodlands Hill, CA) at 8.1 mm depth. Electrical stimulation consisted in either single-shock STN stimulations or 130 Hz DBS (60 µs, 2–4 V), using a stimulus isolator (DS2A, Digitimer, WelWyn Garden City, UK) driven by a pulse stimulator (Pulsemaster A300, WPI, Hitchin, UK).
In vivo juxtacellular electrophysiological recordings
Rats
Unilaterally 6-OHDA-lesioned or sham-lesioned rats were anesthetized with chloral hydrate (400 mg/kg, i.p., supplemented by continuous injection delivered at 60 mg kg−1 h−1). Craniotomies: see intracellular recordings. Single-unit activity was recorded extracellularly in deep layers of M1 (depth range of recorded neurons: 1422–2622 µm) using glass micropipettes (10–15 MΩ) filled with NaCl (1 M) and identified as pyramidal based on their waveform features (with waveform width from depolarization to trough >1 ms; Supplementary Fig. 1a). Single neuron action potentials were recorded using the active bridge mode of a Axoclamp-2A amplifier (Molecular Devices), amplified with a DAM50 (WPI). Spikes were sampled at 10 kHz via a CED 1401 interface using Spike2 (Cambridge Electronic Design). The STN ipsilateral to the recorded pyramidal cells was stimulated as described in intracellular recordings. Neuronal spontaneous activities were recorded for at least 120 s before, during, and after DBS (60 µs, 2–4 V at 130 Hz).
Mice
Mice were anesthetized using urethane (1.5 g/kg i.p.). A craniotomy was made over the STN (from bregma: AP −1.8 mm, ML 1.5 mm). Recording electrodes (15-25MΩ) contained NaCl (0.5 M) and neurobiotin (1.5%). Electrode signals were amplified using an Axoclamp-2B amplifier (Axon Instruments), digitized at 25 kHz and stored in Spike2. STN position (depth of ∼4.3 mm) was based on characteristic discharges of STN neurons in urethane anesthesia63. Bursts in STN neurons were detected using the Poisson surprise method. To confirm STN position, some STN neurons were electroporated and filled with neurobiotin. Once the STN was located, the recording electrode was removed and a bipolar electrode (CBDSF75, FHC, Bowdoin, ME) was targeted into the STN using the same coordinates. For opto-identification of neuronal populations, a craniotomy was made over M1 ipsilateral to the STN (from bregma: +1.5 mm, ML 1.5 mm). An optical fiber connected to the light source (Plexbright 465 nm, Plexon, Dallas, TX) was placed on the cortical surface. Glass electrodes were lowered through the cortex while flashing light pulses (100 ms, 0.5 Hz, 10 mW). After opto-identification, spontaneous activity was defined with 2 min of stable recording before launching DBS (60µs, 120 µA at 130 Hz during 2 min) using pulse isolator and stimulator (IsoFlex + Master-8, A.M.P.I, Jerusalem, Israel). Finally, neurons underwent a series of 1 min opto-activations (3 ms) at 13, 67, and 130 Hz, with at least 1 min of recovery between each series. Neurons were juxtacellularly labeled, and electrolytic lesions of the STN were performed (30 µA DC for 30 s) for posthoc confirmation. Animals with lesions outside the STN were excluded. Because Pv::ChR2 and Sst::ChR2 mouse strains display a small proportion of off-target expression of ChR2 (Supplementary Fig. 2), we used cluster analysis to combine opto-activation and waveform characteristics for a better identification (Supplementary Fig. 3). Principal component analysis (Matlab) was applied to waveform characteristics (time and amplitude of the trough and of the second peak relative to the first peak, on the filtered waveform) and opto-activation success rate. Neurons were clustered using hierarchical and k-means clustering algorithms (Matlab), testing all combinations of the first 2, to all principal components, with 2–7 target clusters. The resulting cluster silhouettes were ranked and the best score was selected (mean silhouette = 0.73, with 2 principal components, 2 clusters, k-means algorithm), resulting in two clearly separated clusters. Morphological and immunohistochemical identity of labeled neurons (n = 18) allowed us to identify the two clusters as pyramidal cells and interneurons. Two cells where the morphological and clustered identity did not match were excluded. Opto-response characteristics: the success rate was the percentage of light pulses eliciting at least a spike (occurring during the pulse); the opto-response latency was the time elapsed between the onset of the light pulse and the first spike within the light pulse (mean, SD and CV of the latency were calculated). Neuronal responses to DBS were classified as non-modulated if the variation of their spontaneous activity was <5% compared to baseline activity.
In vivo patch-clamp recordings
Mice were anaesthetized with urethane as described for juxtacellular recordings. A 0.5 × 0.5 mm craniotomy was performed to expose M1 (from bregma: AP +0.5 mm, ML 1.5 mm). Borosilicate glass pipettes (3.5–5 MΩ) contained (mM): K-gluconate (125), KCl (8), HEPES (10), phosphocreatine (10), ATP-Mg (4), GTP-Na (0.3) and EGTA (0.3) (adjusted to pH 7.35 with KOH). Signals were amplified using an EPC10-2 amplifier (HEKA Elektronik, Lambrecht, Germany) and sampled at 10 kHz, with Patchmaster v2x32 (HEKA Elektronik). ChR2 was activated using a 475 nm-laser diode light source (FLS-475, DIPSI, Cancale, France) with single light pulses (3–20 ms, 10 mW) at 0.33 Hz, or with 3 ms-pulse trains at 67 Hz during 5 s. The optic fiber was placed on top of the pia with an angle of 15°. Pyramidal cells were recorded at 790–1350 µm from the pia.
Behavior
Optic fiber or DBS electrode implantation
For opto-stimulated mice, an optic fiber cannula (200 µm core, 0.39 NA, CFML 12U-20, Thorlabs, Maisons-Laffitte, France) was implanted in M1 (from bregma: AP +1.2 mm, ML 1.5 mm, 0.2–0.4 mm depth from the surface). For DBS-stimulated mice: see juxtacellular recordings. After closing the craniotomy with silicone (Kwik-Cast, WPI), the implant was cemented to the skull and protected by a coppermesh circlet or a plastic cylinder. Mice recovered for 10 days before behavioral testing.
Optogenetic and DBS stimulations
Blue light (475 nm) was delivered from a laser diode light source (FLS-475, DIPSI) through a rotary joint (FRJ_1x2i_FC-2FC, Doric Lenses, Québec, Canada) connected to the mouse optic fiber cannula, controlled by a Master-9 (A.M.P.I.) delivering 3 ms light pulses (10 mW max at the tip) at 67 Hz. We chose a stimulation frequency (67 Hz) that was close to the opsin tau-off (~15–20 ms64). Indeed, higher frequencies such as 130 Hz (the same frequency as electrical DBS) would lose efficiency quickly, while lower frequencies such as 13 Hz, able to entrain all the neurons at each pulse (Fig. S4a–d), might impose pathological-like synchronous activation65. DBS stimulation: electrical current (60 µs, 120 µA, 130 Hz). The light/electrical stimulation and the synchronized video were recorded with KJE-1001 (Amplipex, Szeged, Hungary).
Levodopa treatment
In the first session, mice were injected with saline (6 µL/g, i.p.), and went through behavioral tests after 30 min in that order: open field, cross-maze, cylinder test, with 5-min rest in their homecage between each test. In the second session, mice were injected (6 µL/g, i.p.) with levodopa (6 mg/kg) with benserazide (3.75 mg/kg) and underwent the same tests in the same order. All the behavioral testing occurred within 30-100 min after the saline or levodopa injection.
Open field
Mice were placed in a 40 × 25 × 20 cm arena for 30 s habituation period. For opto- and DBS-stimulation sessions consisted of alternating opto- or DBS-stimulation periods off and on for 1 min, for 10 min. Saline and levodopa sessions lasted for 5 min. Head position was tracked using DeepLabCut algorithm66, then custom Matlab scripts.
Cylinder test
Mice were placed in a glass cylinder (12 cm diameter). For opto- and DBS-stimulated mice, sessions consisted in 3 epochs (before, during and after stimulation) of 10 rearing episodes each. Saline and levodopa sessions consisted in 10 rearing episodes each. The asymmetry of paw contacts on the cylinder surface was measured as % of contacts using paw ipsilateral to the lesion − % of contacts using the contralateral paw.
Cross-maze
Mice were left free to explore the cross-maze (arms: 37 × 7 cm) until they performed at least 30 turns (left, right or straight) within 30 min. Right turns: turns towards the ipsilateral side of the dopaminergic lesion. For opto- and DBS-stimulation, sequential independent sessions were performed: off, on and then off during the whole sessions. For saline and levodopa, 2 sessions (≥1-day interval) were recorded. The percentage of turn in each direction was calculated on the first 30 turns.
Histology
Assessment of 6-OHDA lesion
Coronal slices (50–80 µm thickness) of striatum and SNc were processed for Tyrosine hydroxylase(TH)-immunohistochemistry, after endogenous peroxidases inactivation, using a rabbit anti-TH (1/1000, 36 h at 4 °C, AB152, Merck Millipore), and biotinylated goat antibody against rabbit IgG (1/200, BA 1000, Vector Laboratories, CliniSciences, Nanterre, France), followed by revelation using an ABC kit (PK-6100, Vector Laboratories).
Efficacy and specificity of transgenic mouse lines
Coronal slices (30–50 µm) from Pv::ChR2 or Sst::ChR2 mice were processed for immunostaining of PV or SST, and YFP (ChR2 expression reporter). For SST/YFP immunostaining, an extra step of antigen retrieval was performed in citrate buffer for 15 min at 75 °C. Primary antibodies: rabbit anti-PV (1/2500, PV25, Swant, Marly, Switzerland), rat anti-SST (1/100, MAB354, Merck Millipore), chicken anti-GFP (1/1000, AB13970, Abcam, Paris, France), incubated 24 h (PV/YFP) or 72 h (SST/YFP) at 4 °C; and secondary antibodies: donkey anti-rabbit coupled to Alexa647 (1/500, 711-605-152, Jackson ImmunoResearch), donkey anti-rat coupled to Alexa647 (1/500, 712-605-153, Jackson ImmunoResearch), and goat anti-chicken coupled to Alexa488 (1/500, 103-545-155, Jackson ImmunoResearch) incubated 1–2 h at RT.
Juxtacellular labeling
Floating coronal sections (60 µm) were incubated (2 h at RT) in Alexa 488-conjugated streptavidin (1/250, S11223, Invitrogen, Thermofisher) in PBS containing 0.2% Triton X-100. Sections were processed for immunostaining of SST and PV as described above.
Images were taken with a stereozoom fluorescence microscope (Axiozoom, Zeiss) or confocal microscope (SP5, Leica), and analyzed using ImageJ software.
Statistics
Unless otherwise stated, normal data are displayed as mean ± SEM, and differences between groups were assessed using two-tailed unpaired or paired t-test. Non-normal data are presented as boxplots (Matlab built-in boxplot function defaults) where the center line is the median, the box represents the 25% and 75% quartiles (Q1 and Q3), and the whiskers extends to the last data point within 1.5* the interquartile range outside of the box range (bottom whisker: Q1–1.5*(Q3–Q1), and top whisker: Q3 + 1.5*(Q3–Q1)) (data outside the whisker range are not shown, but are included in the statistical analysis), and differences between groups were assessed using two-tailed Mann–Whitney U test (unpaired data) and Wilcoxon’s signed rank test (paired data). Normality of each dataset was tested using D’Agostino and Pearson’s Omnibus K² test. 2-way ANOVAs were performed using Prism 5.0 (GraphPad, San Diego, CA, USA).
Computational model: spiking network model
We built a simplified spiking model of layer 5 of the motor cortex: the network consisted of 800 pyramidal, 120 PV and 80 SST cells, consistent with the 1:5 ratio between cortical excitatory and inhibitory cells and the larger proportion of PV vs. SST neurons in L525. Neuron of index i was modeled as an adaptive exponential integrate-and-fire neuron67: its activity depends upon the dynamics of a fast voltage variable \(v^{i}\) and a slow adaptation variable \(w^{i}\), according to the Eqs. (1) and (2):
The value of the voltage \(v^{i}\) is determined by intrinsic dynamics and input integration mechanisms. The intrinsic dynamics consists of a linear relaxation towards the resting state and an exponential spiking current active once the spike threshold Vthres is reached. The voltage also integrates excitatory/inhibitory synaptic inputs from the network as well as external stimulations: external constant currents \(I_{{\mathrm{ext}}}\), specific stimuli \(I_{{\mathrm{stimulus}}}\), DBS-induced currents \(I_{{\mathrm{DBS}}}\), as well as noisy currents modeled as independent Gaussian white noise \({(\xi^{i})}\). These stochastic inputs encompass the variety of sources of fluctuations of the voltage68, whose amplitude is determined by the parameter \({\sigma}\). The slower adaptation variable \(w^{i}\) serves as a negative feedback, whose dynamics is linear in \(v^{i}\) and controlled by the parameter \(a\) determining the strength of the adaptation-voltage coupling.
A spike is triggered when the voltage approaches the threshold Vthres. Upon firing, the variable \(v^{i}\) is reset to a fixed value, equal to the leak reversal potential, Vreset = Eleak = −60mV67, whereas \(w^{i}\) is increased by a fixed amount b, corresponding to the spike-triggered adaptation.
Choices of intrinsic parameters (Supplementary Table 1) were guided by experimental data. First, pyramidal cells had a larger leak conductance gleak and capacitance C compared to interneurons. The spike threshold of pyramidal cells was adjusted in parallel with the amount of constant external current Iext received such that, in absence of stimulus, their firing rate remained around 1 Hz23. Secondly, PV cells were modeled as fast-spiking interneurons, with a sharp spiking onset and no adaptation25 whereas the profile of SST interneurons was characterized by a lower spike threshold and the existence of spike-frequency adaptation26,67,69. Adaptation parameters were chosen following Brette and Gerstner67 and Naud et al.69. The parkinsonian model incorporated the hyperexcitability of pyramidal neurons through a decreased spike threshold10.
Synaptic connections were randomly distributed, with a probability p and fixed synaptic weight w (Fig. 6a). We opted for a simple conductance-based description of synaptic currents: the excitatory and inhibitory conductances gexc and ginh display a discrete jump following spikes (arriving at times texc or tinh), after a delay of 1 ms, and decay exponentially with time constants τexc and τinh and reversal potentials Eexc and Einh, respectively50 (Eqs. (3) and (4)).
To mimic the somatic impact of DBS on cortical neurons, we added an external current IDBS to the equation of the voltage variable in every cell of the network. More precisely, considering the periodic nature of DBS-induced somatic currents, IDBS corresponds to a series of square pulses (2 ms and 120 pA, at 130 Hz). Considering the different putative activation pathways, IDBS impact neurons after specific delays, chosen heterogeneous: 0 ms for half of pyramidal and PV cells (fast antidromic pathways), 2 ms for the other half of pyramidal and SST cells (orthodromic loops). The absence of such delays does not compromise any of the modeling results. For optogenetic stimulations, a series of square pulses (3 ms, 600 pA) was applied to half of PV or SST neurons. This fraction was chosen considering the fact that approximately 70% of targeted neurons express the channelrhodopsin, among which about 70% respond to the light pulse. These optogenetic-like pulses were repeated at a variable frequency, which unless stated otherwise was set to the experimental value of 67 Hz. Additional tests with different current amplitudes, ranging from 200 to 800 pA were also simulated (Supplementary Fig. 7).
Simulations of the network activity were done using a custom code developed in Matlab R2016 (The Mathworks).
Computational model: analysis
Firing rates
Average firing rates were computed for each population from 1 s simulations of network activity. When mimicking PV or SST opto-activation, the average firing rates of PV or SST interneurons were calculated only based on the neurons that were directly activated.
Heatmaps of pyramidal cell firing rate modulation were obtained by varying the intensity (from 0 to 120 pA) of IDBS (with default parameters: square pulses of 2 ms, repeated at 130 Hz) received by two populations at the same time, while keeping the default intensity (120 pA) for the third population. Pyramidal cell firing rate modulation corresponds to the average firing rate of pyramidal cells under DBS, divided by the average firing rate of pyramidal cells in the parkinsonian condition.
A linear regression of pyramidal cell firing rates as a function of the amount of IDBS (corresponding to the total amount of current in pA injected over one second, equal to: stimulation intensity x pulse duration x stimulation frequency) was performed.
Pearson correlation coefficients
To test the capacity of the network to discriminate and respond specifically to various inputs, an additional current was injected to a subset of pyramidal cells (200 randomly chosen cells). We explored the responses of the network to three types of stimuli, presented during 500 ms: constant input, linear ramping input (decreasing or increasing amplitude with time) or noisy inputs (Ornstein–Uhlenbeck processes xt with different mean μ, variance σ and time constant τ, according to \(dx_{t} = {\textstyle{{\left( {\mu - x_{t}} \right)} \over \tau }}dt + \sigma \,dW_{t}\) with x0 = 0). These either deterministic or stochastic inputs (Supplementary Table 2) were chosen such as to mimic some activity patterns observed in pyramidal neurons in the motor cortex70.
For each condition, we used two different methods for quantifying the capacities of extracting information from the network activity patterns. We first measured for each trial the Pearson correlation coefficient ρ between the moving spike count of all pyramidal cells across time for ramping and stochastic stimuli. The moving spike count was calculated based on the sum of the number of spikes across all pyramidal cells, with a moving time interval of 10 ms. The first and last 10 ms of stimulus presentation were discarded to avoid boundary effects. For each stimulus, an average correlation coefficient was obtained by averaging the Pearson correlation index across 100 independent trials (with different intrinsic noises, but identical connectivity matrices, to explore specifically the variability of the responses of a given network to the same input).
Classification of stimuli
We also estimated the efficiency with which downstream neurons might discriminate the network responses to various stimuli, beyond the sole knowledge of the mean firing rate. We used four supervised-learning algorithms (Supplementary Table 3) to classify the responses of pyramidal cells to 20 stimuli. Our approach consisted of the following steps (Fig. 7a): our network responses are a time-binned matrix M in which each row corresponding to a given pyramidal cell contains the number of spikes emitted for every 10 ms. The first 200 rows corresponded to the responses of the pyramidal cells directly activated by the stimulus. In order to reproduce the highly convergent cortical motor inputs received by striatal neurons, we densified the responses by contracting these 800 × 50 M matrices into a 100 × 50 S matrix, defined as: S = W.M, where the weight matrix W is a random matrix, identical for all stimuli presentations, with each element generated from the uniform distribution on the interval [0, 1] (Fig. 7a). These S matrices were then used as inputs for the supervised-learning algorithms (nearest centroid classifier, multinomial logistic regression, linear discriminant analysis and support vector machines71,72). For each condition, the dataset consisted of 100 repetitions for each of the 20 stimuli (with the same network parameters and convergence matrix W, but independent realizations of the intrinsic noise ξ). The classifiers were trained to discriminate the population response given the stimulus on 80% of the data sample, using stratified k-fold cross validation (k = 5). This procedure was repeated using 5 independent seeds. The training and testing of the classifiers were run using scikit-learn and keras packages in Python 3.5 (Python Software Foundation, www.python.org).
Statistics
Normally distributed data are displayed as mean ± SD, and differences between conditions were assessed using one-way ANOVA for each feature and Tukey–Kramer post-hoc test. Since some distributions of correlation coefficients were not normal, all data is summarized using the median and presented in details as boxplots (using the same convention as described above; Supplementary Fig. 6) and differences between conditions were assessed using the Kruskal–Wallis test (with Tukey–Kramer post-hoc test) on every stimulus. The statistics for the Pearson correlation coefficients and classifier accuracies were performed considering the full model (Supplementary Fig. 7).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
The data sets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Code availability
Custom Matlab codes used for the computational model are publicly available on GitHub using the following link: https://github.com/cpiette95/Information_processing_DBS_Cortex.
Change history
30 July 2020
The original version of this Article was updated shortly after publication following an error that resulted in the ORCID ID of Laurent Venance being omitted.
References
Obeso, J. et al. Past, present, and future of Parkinson’s disease: a special essay on the 200th Anniversary of the Shaking Palsy. Mov. Disord. 32, 1264–1310 (2017).
Limousin, P. et al. Effect on parkinsonian signs and symptoms of bilateral subthalamic nucleus stimulation. The Lancet 345, 91–95 (1995).
Ashkan, K., Rogers, R., Bergman, H. & Ughratdar, I. Insights into the mechanisms of deep brain stimulation. Nat. Rev. Neurol. 13, 548–554 (2017).
Deniau, J. M., Degos, B., Bosch, C. & Maurice, N. Deep brain stimulation mechanisms: beyond the concept of local functional inhibition. Eur. J. Neurosci. 32, 1080–1091 (2010).
Wichmann, T., Bergman, H. & DeLong, M. Basal ganglia, movement disorders and deep brain stimulation: advances made through non-human primate research. J. Neural Transm. 125, 419–430 (2017).
Aum, D. & Tierney, T. Deep brain stimulation foundations and future trends. Front. Biosci. 23, 162–182 (2018).
Li, S., Arbuthnott, G., Jutras, M., Goldberg, J. & Jaeger, D. Resonant antidromic cortical circuit activation as a consequence of high-frequency subthalamic deep-brain stimulation. J. Neurophysiol. 98, 3525–3537 (2007).
Gradinaru, V., Mogri, M., Thompson, K., Henderson, J. & Deisseroth, K. Optical deconstruction of parkinsonian neural circuitry. Science 324, 354–359 (2009).
Li, Q. et al. Therapeutic deep brain stimulation in parkinsonian rats directly influences motor cortex. Neuron 76, 1030–1041 (2012).
Degos, B., Deniau, J. M., Chavez, M. & Maurice, N. Subthalamic nucleus high-frequency stimulation restores altered electrophysiological properties of cortical neurons in parkinsonian rat. PLoS ONE 8, e83608 (2013).
Anderson, R., Farokhniaee, A., Gunalan, K., Howell, B. & McIntyre, C. Action potential initiation, propagation, and cortical invasion in the hyperdirect pathway during subthalamic deep brain stimulation. Brain Stimul. 11, 1140–1150 (2018).
Cunic, D. et al. Effects of subthalamic nucleus stimulation on motor cortex excitability in Parkinson’s disease. Neurology 58, 1665–1672 (2002).
Payoux, P. Subthalamic nucleus stimulation reduces abnormal motor cortical overactivity in Parkinson disease. Arch. Neurol. 61, 8 (2004).
Haslinger, B., Kalteis, K., Boecker, H., Alesch, F. & Ceballos-Baumann, A. Frequency-correlated decreases of motor cortex activity associated with subthalamic nucleus stimulation in Parkinson’s disease. NeuroImage 28, 598–606 (2005).
Fraix, V., Pollak, P., Vercueil, L., Benabid, A. L. & Mauguière, F. Effects of subthalamic nucleus stimulation on motor cortex excitability in Parkinson’s disease. Clin. Neurophysiol. 119, 2513–2518 (2008).
Kim, S. et al. Effects of subthalamic nucleus stimulation on motor cortex plasticity in Parkinson disease. Neurology 85, 425–432 (2015).
Lindenbach, D. & Bishop, C. Critical involvement of the motor cortex in the pathophysiology and treatment of Parkinson’s disease. Neurosci. Biobehav. Rev. 37, 2737–2750 (2013).
Pelled, G., Bergman, H. & Goelman, G. Bilateral overactivation of the sensorimotor cortex in the unilateral rodent model of Parkinson’s disease - a functional magnetic resonance imaging study. Eur. J. Neurosci. 15, 389–394 (2002).
Dégenètais, E., Thierry, A. M., Glowinski, J. & Gioanni, Y. Synaptic influence of hippocampus on pyramidal cells of the rat prefrontal cortex: an in vivo intracellular recording study. Cereb. Cortex. 13, 782–792 (2003).
Tremblay, R., Lee, S. & Rudy, B. GABAergic interneurons in the neocortex: from cellular properties to circuits. Neuron 91, 260–292 (2016).
Ma, W. et al. Visual representations by cortical somatostatin inhibitory neurons—selective but with weak and delayed responses. J. Neurosci. 30, 14371–14379 (2010).
Schor, J. S. & Nelson, A. B. Multiple stimulation parameters influence efficacy of deep brain stimulation in parkinsonian mice. J. Clin. Invest 130, 3833–3838 (2019).
Gentet, L., Avermann, M., Matyas, F., Staiger, J. & Petersen, C. Membrane potential dynamics of GABAergic neurons in the barrel cortex of behaving mice. Neuron 65, 422–435 (2010).
Apicella, A., Wickersham, I., Seung, H. & Shepherd, G. Laminarly orthogonal excitation of fast-spiking and low-threshold-spiking interneurons in mouse motor cortex. J. Neurosci. 32, 7021–7033 (2012).
Tanaka, Y., Fujiyama, F., Furuta, T., Yanagawa, Y. & Kaneko, T. Local connections of layer 5 GABAergic interneurons to corticospinal neurons. Front. Neural Circuits 5, 12 (2011).
Pfeffer, C., Xue, M., He, M., Huang, Z. & Scanziani, M. Inhibition of inhibition in visual cortex: the logic of connections between molecularly distinct interneurons. Nat. Neurosci. 16, 1068–1076 (2013).
Fino, E. & Yuste, R. Dense inhibitory connectivity in neocortex. Neuron 69, 1188–1203 (2011).
Kiritani, T., Wickersham, I., Seung, H. & Shepherd, G. Hierarchical connectivity and connection-specific dynamics in the corticospinal-corticostriatal microcircuit in mouse motor cortex. J. Neurosci. 32, 4992–5001 (2012).
Weiler, N., Wood, L., Yu, J., Solla, S. & Shepherd, G. Top-down laminar organization of the excitatory network in motor cortex. Nat. Neurosci. 11, 360–366 (2008).
Chadderdon, G. et al. Motor cortex microcircuit simulation based on brain activity mapping. Neural Comput. 26, 1239–1262 (2014).
Cavallari, S., Panzeri, S. & Mazzoni, A. Comparison of the dynamics of neural interactions between current-based and conductance-based integrate-and-fire recurrent networks. Front. Neural Circuits 8, 12 (2014).
Destexhe, A. Self-sustained asynchronous irregular states and Up–Down states in thalamic, cortical and thalamocortical networks of nonlinear integrate-and-fire neurons. J. Comput. Neurosci. 27, 493–506 (2009).
Carron, R. et al. Early hypersynchrony in juvenile PINK1−/− motor cortex is rescued by antidromic stimulation. Front. Sys. Neurosci 8, 95 (2014).
Degos, B., Deniau, J. M., Le Cam, J., Mailly, P. & Maurice, N. Evidence for a direct subthalamo-cortical loop circuit in the rat. Eur. J. Neurosci 27, 2599–2610 (2008).
Muñoz, W., Tremblay, R., Levenstein, D. & Rudy, B. Layer-specific modulation of neocortical dendritic inhibition during active wakefulness. Science 355, 954–959 (2017).
Fu, Y., Tucciarone, J. M., Espinosa, J. S., Sheng, N., Darcy, D. P., Nicoll, R. A., Huang, Z. J. & Stryker, M. P. A cortical circuit for gain control by behavioral state. Cell 156, 1139–1152 (2014).
Urban-Ciecko, J. & Barth, A. Somatostatin-expressing neurons in cortical networks. Nat. Rev. Neurosci. 17, 401–409 (2016).
Yavorska, I. & Wehr, M. Somatostatin-expressing inhibitory interneurons in cortical circuits. Front. Neural Circuits 10, 76 (2016).
Scheyltjens, I. & Arckens, L. The current status of somatostatin-interneurons in inhibitory control of brain function and plasticity. Neural Plast. 2016, 1–20 (2016).
Cichon, J., Blanck, T., Gan, W. & Yang, G. Activation of cortical somatostatin interneurons prevents the development of neuropathic pain. Nat. Neurosci. 20, 1122–1132 (2017).
Scheggia, D. et al. Somatostatin interneurons in the prefrontal cortex control affective state discrimination in mice. Nat. Neurosci. 23, 47–60 (2020).
Ramaswamy, S., Colangelo, C. & Muller, E. B. Distinct activity profiles of somatostatin-expressing interneurons in the neocortex. Front. Cell Neurosci. 11, 273 (2017).
Cardin, J. A. Inhibitory interneurons regulate temporal precision and correlations in cortical circuits. Trends Neurosci. 41, 689–700 (2018).
Lundblad, M., Picconi, B., Lindgren, H. & Cenci, M. A. A model of L-DOPA-induced dyskinesia in 6-hydroxydopamine lesioned mice: relation to motor and cellular parameters of nigrostriatal function. Neurobiol. Dis. 16, 110–123 (2004).
Santini, E. et al. L-DOPA activates ERK signaling and phosphorylates histone H3 in the striatonigral medium spiny neurons of hemiparkinsonian mice. J. Neurochem. 108, 621–633 (2009).
Darbaky, Y., Forni, C., Almaric, M. & Baunez, C. High frequency stimulation of the subthalamic nucleus has beneficial antiparkinsonian effects on motor functions in rats, but less efficiency in a choice reaction task. Eur. J. Neurosci. 18, 951–956 (2003).
Maesawa, S. et al. Long-term stimulation of the subthalamic nucleus in hemiparkinsonian rats: neuroprotection of dopaminergic neurons. J. Neurosurg. 100, 679–697 (2004).
Shi, L. H., Woodward, D. J., Luo, F., Anstrom, J., Schallert, T. & Chang, J. Y. High frequency stimulation of the subthalamic nucleus reverses limb-use asymmetry in rats with unilateral 6-hydroxydopamine lesions. Brain Res. 2013, 98–106 (2004).
Gubellini, P., Eusebio, A., Oueslati, A., Melon, C., Kerkerian-Le Goff, L. & Salin, P. Chronic high-frequency stimulation of the subthalamic nucleus and L-DOPA treatment in experimental parkinsonism: effects on motor behaviour and striatal glutamate transmission. Eur. J. Neurosci 24, 1802–1914 (2006).
Fang, X., Sugiyama, K., Akamine, S., Sun, W. & Namba, H. The different performance among motor tasks during the increasing current intensity of deep brain stimulation of the subthalamic nucleus in rats with different degrees of the unilateral striatal lesion. Neurosci. Lett. 480, 64–68 (2010).
Tommasi, G. et al. Pyramidal tract side effects induced by deep brain stimulation of the subthalamic nucleus. J. Neurol. Neurosurg. Psychiatry. 79, 813–819 (2008).
Rothblat, D. S. & Schneider, J. S. Alterations in pallidal neuronal responses to peripheral sensory and striatal stimulation in symptomatic and recovered parkinsonian cats. Brain Res. 705, 1–14 (1995).
Cho, J., Duke, D., Manzino, L., Sonsalla, P. K. & West, M. O. Dopamine depletion causes fragmented clustering of neurons in the sensorimotor striatum: evidence of lasting reorganization of corticostriatal input. J. Comp. Neurol. 452, 24–37 (2002).
Paz, J. et al. Closed-loop optogenetic control of thalamus as a tool for interrupting seizures after cortical injury. Nat. Neurosci. 16, 64–70 (2013).
Chen, B. et al. Rescuing cocaine-induced prefrontal cortex hypoactivity prevents compulsive cocaine seeking. Nature 496, 359–362 (2013).
Walsh, J. et al. Optogenetic manipulation of ventral tegmental area (VTA) neurons that project to the nucleus accumbens (NAc) and medial prefrontal cortex (mPFC). Protocol Exchange https://doi.org/10.1038/protex.2012.049 (2012).
Creed, M. C., Pascoli, V. & Lüscher, C. Refining deep brain stimulation to emulate optogenetic treatment of synaptic pathology. Science 347, 659–664 (2015).
Adamantidis, A. et al. Optogenetics: 10 years after ChR2 in neurons-views from the community. Nat. Neurosci. 18, 1202–1212 (2015).
Rajasethupathy, P., Ferenczi, E. & Deisseroth, K. Targeting neural circuits. Cell 165, 524–534 (2016).
Viana Magno, L. A. et al. Optogenetic stimulation of the M2 cortex reverts motor dysfunction in a mouse model of Parkinson’s disease. J. Neurosci. 39, 3234–3248 (2019).
Lüscher, C. & Pollak, P. Optogenetically inspired deep brain stimulation: linking basic with clinical research. Swiss Med. Wkly. 146, w14278 (2016).
Bouthour, W., Krack, P. & Lüscher, C. A deeply superficial brain stimulation. Mov. Disord 32, 1326–1326 (2017).
Magill, P., Bolam, J. & Bevan, M. Dopamine regulates the impact of the cerebral cortex on the subthalamic nucleus–globus pallidus network. Neuroscience 106, 313–330 (2001).
Lin, J. Y., Lin, M. Z., Steinbach, P. & Tsien, R. Y. Characterization of engineered channelrhodopsin variants with improved properties and kinetics. Biophys. J. 96, 1803–1814 (2009).
Singh, A. Oscillatory activity in the cortico-basal ganglia-thalamic neural circuits in Parkinson’s disease. Eur. J. Neurosci. 48, 2869–2878 (2018).
Mathis, A. et al. DeepLabCut: markerless pose estimation of user-defined body parts with deep learning. Nat. Neurosci. 21, 1281–1289 (2018).
Brette, R. & Gerstner, W. Adaptive exponential integrate-and-fire model as an effective description of neuronal activity. J. Neurophysiol. 94, 3637–3642 (2005).
Faisal, A., Selen, L. & Wolpert, D. Noise in the nervous system. Nat. Rev. Neurosci. 9, 292–303 (2008).
Naud, R., Marcille, N., Clopath, C. & Gerstner, W. Firing patterns in the adaptive exponential integrate-and-fire model. Biol. Cybern. 99, 335–347 (2008).
Costa, R. M., Cohen, D. & Nicolelis, M. A. Differential corticostriatal plasticity during fast and slow motor skill learning in mice. Curr. Biol. 13, 1124–1134 (2004).
Bishop, C. Pattern Recognition and Machine Learning (Springer, New York, 2006).
Hastie, T., Friedman, J. & Tibshirani, R. The Elements of Statistical Learning (Springer, New York, 2017).
Acknowledgements
We thank the L.V. lab members, Christian Giaume and Pierre Magistretti for helpful suggestions and critical comments. We thank Sylvie Perez, Alessandra Romei, Farah Hadj-Idris, Angèle Roudeau and Charlotte Branco for technical assistance for 6-OHDA lesioned-rodents and immunohistochemistry. We thank France Maloumian for her help in the schematic illustrations in Figs. 3 and 5. This work was supported by grants from Fondation de France (Maladie de Parkinson) CRCPEN, Fondation du Collège de France, France Parkinson, Fondation Patrick Brou de Laurière, INSERM, Collège de France and CNRS. C.P. is supported by Ecole Normale Supérieure, G.G by Fondation Patrick Brou de Laurière, and S.V., W.D. and M.V. by the Collège de France.
Author information
Authors and Affiliations
Contributions
Conceptualization: L.V., B.D., and J.T.; S.V., B.D., A.A.A., and W.D. performed in vivo extracellular electrophysiological experiments; S.V. performed in vivo juxtacellular electrophysiological, opto-identification and immunohistochemistry experiments; B.D. performed in vivo intracellular electrophysiological experiments; W.D. performed in vivo patch-clamp recordings; M.V. and GG performed in vivo optogenetic, behavioral and immunohistochemistry experiments; M.V., B.D., S.V., C.P., W.D., G.G., and L.V. performed analysis; B.D., S.V., G.G., and M.V. performed 6-OHDA lesions; J.T. and C.P. carried out the conception and the design of the mathematical model; C.P. performed the acquisition and analysis of data from the mathematical model; L.V., B.D., J.T., C.P., M.V., and S.V. wrote the paper and all authors have edited and corrected the paper; Funding acquisition: B.D. and L.V.; L.V. supervised the whole study.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Un Kang, Alex Kwan and Giuseppe Sciamanna for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Valverde, S., Vandecasteele, M., Piette, C. et al. Deep brain stimulation-guided optogenetic rescue of parkinsonian symptoms. Nat Commun 11, 2388 (2020). https://doi.org/10.1038/s41467-020-16046-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-16046-6
This article is cited by
-
Advancement in modulation of brain extracellular space and unlocking its potential for intervention of neurological diseases
Med-X (2024)
-
Neuronal and synaptic adaptations underlying the benefits of deep brain stimulation for Parkinson's disease
Translational Neurodegeneration (2023)
-
High frequency DBS-like optogenetic stimulation of nucleus accumbens dopamine D2 receptor-containing neurons attenuates cocaine reinstatement in male rats
Neuropsychopharmacology (2023)
-
Thalamus drives vocal onsets in the zebra finch courtship song
Nature (2023)
-
Cerebellar stimulation prevents Levodopa-induced dyskinesia in mice and normalizes activity in a motor network
Nature Communications (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.