Abstract
The striatum comprises of multiple functional territories involved with multilevel control of behavior. Disinhibition of different functional territories leads to territory-specific hyperkinetic and hyperbehavioral symptoms. The ventromedial striatum, including the nucleus accumbens (NAc) core, is typically associated with limbic input but was historically linked to high-level motor control. In this study, performed in female Long–Evans rats, we show that the NAc core directly controls motor behavior on multiple timescales. On the macro-scale, following NAc disinhibition, the animals manifested prolonged hyperactivity, expressed as excessive normal behavior, whereas on the micro-scale multiple behavior transitions occurred, generating short movement segments. The underlying striatal network displayed population-based local field potential transient deflections (LFP spikes) whose rate determined the magnitude of the hyperactivity and whose timing corresponded to unitary behavioral transition events. Individual striatal neurons preserved normal baseline activity and network interactions following the disinhibition, maintaining the normal encoding of behavioral primitives and forming a sparse link between the LFP spikes and single neuron activity. Disinhibition of this classically limbic territory leads to profound motor changes resembling hyperactivity and attention deficit. These behavioral and neuronal results highlight the direct interplay on multiple timescales between different striatal territories during normal and pathological conditions.
SIGNIFICANCE STATEMENT The nucleus accumbens (NAc) is a key part of the striatal limbic territory. In the current study we show that this classically limbic area directly controls motor behavior on multiple timescales. Focal disinhibition of the NAc core in freely behaving rats led to macro-scale hyperactivity and micro-scale behavioral transitions, symptoms typically associated with attention deficit hyperactivity disorder. The behavioral changes were encoded by the striatal LFP signal and single-unit spiking activity in line with the neuronal changes observed during tic expression following disinhibition of the striatal motor territory. These results point to the need to extend the existing parallel functional pathway concept of basal ganglia function to include the study of limbic-motor cross-territory interactions in both health and disease.
Introduction
The striatum, the main input structure of the corticobasal ganglia pathway, receives converging glutamatergic inputs organized in multiple mostly-parallel pathways: sensorimotor afferents terminate primarily in the dorsolateral striatum, whereas limbic afferents target the ventromedial striatum and specifically the nucleus accumbens (NAc; McGeorge and Faull, 1989; Heilbronner et al., 2016). The NAc is further subdivided based on the anatomical division to the core and shell (Voorn et al., 1989; Heimer et al., 1991). Although the separation into territories is widely accepted, the specific functional role of individual territories and sub-territories is under debate. Recent studies of the NAc have focused on limbic functions such as reward, decision making, and addiction (Al-Hasani et al., 2015; Volkow and Morales, 2015; Floresco et al., 2018). Specifically, the core is associated primarily in goal directed behavior, whereas the shell is associated with motivational values (West and Carelli, 2016). However, earlier work emphasized the NAc role in bridging between the limbic and motor systems (Mogenson et al., 1980) and more specifically, control over overall activity in general, and locomotion in particular (Mogenson and Nielsen, 1984). This control over locomotion may be experimentally modulated by focal pharmacological manipulation of glutamate (Arnt, 1981; Swanson and Kalivas, 2000), dopamine (Costall and Naylor, 1975; Jones et al., 1981; Costall et al., 1984), and the extensive GABAergic network of the NAc (Morgenstern et al., 1984; Yoshida et al., 1991). Moreover, this modulation is affected by the interactions between the different pharmacological agents (Pycock and Horton, 1979; Dalia et al., 1998) leading to various hyperactive and hypoactive states. The NAc control over movement is primarily modulated by the NAc core as demonstrated by both lesion and pharmacological manipulations (Maldonado-Irizarry and Kelley, 1994, 1995).
These functional differences between striatal territories are apparent in the expression of a myriad of hyperkinetic and hyperbehavioral conditions involving the different territories. Tourette syndrome, a childhood onset neurological disorder characterized by the expression of motor and vocal tics, is primarily associated with abnormalities in the sensorimotor regions of the corticobasal ganglia pathway (Sowell et al., 2008; Tremblay et al., 2015). Attention deficit hyperactivity disorder (ADHD), a neurodevelopmental disorder characterized by attention difficulties, hyperactivity, and impulsivity, primarily involves limbic parts of the circuit (Giedd et al., 2006; Qiu et al., 2009). These two disorders show high comorbidity rates, in that >50% of all Tourette syndrome patients express comorbid ADHD (Schlaggar and Mink, 2003), hinting that the two disorders may share some common underlying pathophysiology. Disinhibition of different functional territories of the striatum by focal blocking of GABAA transmission in rodents and primates supports this hypothesis. Disinhibition leads to the expression of territory specific hyper-kinetic and hyper-behavioral symptoms; namely, motor tics following disinhibition of the motor territory (Tarsy et al., 1978; McCairn et al., 2009; Bronfeld et al., 2013) and hyperactivity following similar manipulation in the limbic territories (Morgenstern et al., 1984; Yoshida et al., 1991; Worbe et al., 2009).
In this study we use disinhibition to deconstruct the behavioral elements controlled by the NAc core on multiple timescales and identify the underlying neurophysiological mechanism encoding these elements. The relation between behavior and striatal activity has been explored during in the dorsolateral striatum during normal behavior (Klaus et al., 2017; Markowitz et al., 2018) and disinhibition (Bronfeld et al., 2011; Israelashvili and Bar-Gad, 2015; Oran and Bar-Gad, 2018) but is lacking for the ventral striatum during either states. This neuronal-behavioral relationship may serve as a basis for establishing the unique role of the NAc in the control of movement in normal and hyperkinetic conditions.
Materials and Methods
Animals
Fifteen adult female Long–Evans rats (weighing 288.5 ± 28.4 g, mean ± SD) were used in this study. The animals were kept in a temperature and humidity controlled room with a 12 h light/dark cycle and had ad libitum access to food and water. Experiments were performed during the light phase in a lighted setup. All procedures were approved and supervised by the Institutional Animal Care and Use Committee and were in accordance with the National Institute of Health Guide for the Care and Use of Laboratory Animals and the Bar-Ilan University Guidelines for the Use and Care of Laboratory Animals in Research. This study was approved by the National Committee for Experiments in Laboratory Animals at the Israel Ministry of Health.
Surgical procedures
The animals were implanted unilaterally with an injection cannula targeting the NAc core (AP: 1.5, ML: 2.2, DV: 7.0) and a recording electrode array (custom-made movable bundles of 16/32 Formvar isolated nichrome microwires, 25 μm diameter, A-M Systems) in the ventromedial striatum (AP: 0.5, ML: 2.2, DV: 5.0–7.0). In four animals, additional recording electrodes were placed in the dorsolateral striatum (AP: 0.5, ML: 4.0, DV: 3.5–5.5). All the surgical procedures were described in detail in previous studies by our group (Bronfeld et al., 2013; Yael et al., 2013; Israelashvili and Bar-Gad, 2015; Oran and Bar-Gad, 2018).
Experimental sessions
Following a recovery period lasting at least 7 d post-operation, the animals were connected to a wireless recording system and were placed in one of three behavioral arenas allowing free natural behavior: (1) open field (80 × 120 cm), (2) limited field (80 × 120 cm, containing six 20 × 20 cm cubes), and (3) linear field (8 × 180 cm). The animals' neuronal activity was continuously recorded and their behavior was monitored using a 50/60 fps rate video camera and 9-axis kinematic sensors located within the wireless acquisition system. Following a 30 min recording period during the naive state, the animals were removed from the behavioral arena and bicuculline methiodide was injected unilaterally into their NAc core. After the injection, the animals were returned to the arena, and their neuronal activity and behavior were recorded for an additional 30–45 min. The behavioral changes induced by the microinjection followed the temporal dynamics of the bicuculline spread and clearance from the tissue, in line with previous studies of disinhibition using bicuculline in the sensorimotor striatum (Bronfeld et al., 2013). At the end of the experiment and after return to normal behavior, the animals were returned to their home cage.
Microinjections
Unilateral micro-injections of the GABA-A antagonist bicuculline methiodide into the NAc core were performed with an acute injector (30 G stainless steel tube attached to a flexible Tygon microbore tube), which was manually inserted into a chronically implanted guiding cannula (25 G stainless steel tube) located 2 mm above the injection target. Bicuculline methiodide (Sigma-Aldrich) was dissolved in artificial CSF [concentrations (in mm): 145 NaCl, 15 HEPES, 2.5 KCl, 2 MgCl2, 1.2 CaCl2, PH 7.4 with NaOH] at a concentration of 1 μg/μl. A volume of 0.5 μl was injected at a rate of 0.5 μg/min using a 10 μl syringe (Hamilton) attached to a programmable syringe pump (NE-1000, New Era Pump Systems). One minute after the injection termination, the injector was removed and the guiding cannula was sealed by a cannula-dummy (28 G stainless steel wire). The injection procedure, volume and concentration were similar to those previously used during bicuculline injections to the sensorimotor striatum (Bronfeld et al., 2013).
Neuronal recordings and preprocessing
Electrophysiological signals were continuously recorded during the session using either a 16- or 32-channel wireless acquisition system (NeuroLog-16 and RatLog-64, Deuteron Technologies). The analog signal was bandpass filtered (1 Hz one-pole to 7000 or 10,000 Hz three-pole Butterworth filter), amplified (×200), and digitized (12 or 16 bits with a resolution of 3.3 or 0.2 μV, sampling rates: 29.297 or 32 kHz, respectively). Off-line, the recording files were concatenated, filtered (300–6000 Hz, 4-pole bidirectional Butterworth filter), and sorted for the extraction of single-unit spike times (Offline sorter, Plexon). This approach, in which the recording files are concatenated to create a single long file including the full experimental phase, allowed sorting the whole period using the same criteria and examination of the unit stability throughout the session. The suspected spikes were cut out of the streams in 2 ms windows and consequently projected over multiple parameters (such as the PCs and maximal and minimal values). Only units which formed a clear cluster, displayed a clean refractory period (<0.1% of the spikes occurring within short intervals) and which maintained the same waveform were included in the final dataset. The LFP signal was extracted using off-line low-pass filtration of the raw digitized signal (cutoff frequency: 300 Hz, 4-pole bidirectional Butterworth filter). To avoid filter-based distortion of waveforms (Yael and Bar-Gad, 2017) and timing information (Yael et al., 2018), all the off-line filters used in the study were bidirectional (zero-phase) filters. All subsequent analyses were performed using custom-written code (MATLAB, MathWorks).
Data analysis
Rat location detection.
Video recordings were off-line frame-by-frame analyzed to extract the 2-dimensional location of the animals. Detection used custom written MATLAB code based on a background subtraction algorithm. After subtraction of the background frame, individual frames were converted from RGB into grayscale coding. Using Otsu's method (Otsu, 1979), frames were further reduced from grayscale to a binary BW image. The binary image was then smoothed by eroding with a 2-dimensional disk shaped structuring object. Noise was subtracted by the removal of connected components <600 pixels. Finally, centroids were extracted and confirmed to be within a radius of 200 pixels from the centroid detected in the previous frame. The software and documentation are available via GIT (https://github.com/neurintlab/ratDetect).
Distance and velocity calculation.
After the extraction of the 2-dimensional rat location, the coordinates over time were convolved with a Gaussian window (SD = 2 frames). The Euclidean distance between successive coordinates was calculated, converted from pixels to metric units, and normalized by the frame rate to evaluate the corresponding velocity. The velocity was then smoothed using Gaussian window.
Auto and cross-correlations.
Spike train-based auto and cross-correlations were calculated for pairs of simultaneously recorded neurons, for a stable period of 10 min starting 10 min from the beginning of the session (pre-injection and postinjection states). The cross-correlation functions were only calculated for neurons recorded on different channels to avoid artificial correlations (Bar-Gad et al., 2001). The cross-correlation functions were calculated within a time window of ±4.5 s with a time resolution of 1 ms and were smoothed using a Gaussian window (SD = 40 ms). Z-scores were calculated based on a baseline period located ±2–4 s from zero time lag. The correlation was considered significant if the integral of a time window of ±10 ms normalized by the number of bins crossed the confidence levels (p < 0.01). Based on this significance test, the cross-correlation functions were then classified into one of three groups: positive, negative, or no correlation.
Raster plots, PLTHs, and joint PLTHs.
To evaluate the spiking activity around local field potential transient deflections (LFP spikes), peri-LFP-spike rasters were calculated by aligning 2000 ms windows extracted from the single-unit spike trains to the LFP spike onset. To calculate the peri-LFP-spike time histograms (PLTHs), the rasters were averaged and smoothed by a Gaussian window (SD 25 ms).
To estimate the temporal relationship between two spike trains around the LFP spike times we calculated the joint PLTHs (JPLTHs). Pairs of spike trains (bin size: 1 ms) based on simultaneously recorded neurons were used to create a matrix in which all the correlated firing combinations within ±1000 ms window for each event were summed. This matrix was corrected by subtraction of the shift-predicator matrix (calculated by multiplication of the corresponding PLTHs; Aertsen et al., 1989) and smoothed using 2D convolution with a Gaussian window (SD = 50 ms). To evaluate the time dependent coactivation of the single-unit pairs, we extracted the main diagonal of the corrected matrix corresponding to the coactivation of the neurons around the LFP within the window with zero time lags. Sparse matrices were manually removed from downstream calculations to avoid possible artifacts. The significance of the response was estimated by calculating the confidence lines (p < 0.001) based on baseline activity (±400–900 ms from the LFP spike onset). JPLTHs were considered to have a significant response if they crossed the significance lines in a time window of ±50 ms around the LFP spikes onset.
Movement termination detection.
Movement terminations were detected using a threshold-crossing algorithm combined with trough detection based on the animals' velocity (smoothed by a 100 ms SD Gaussian window). The threshold was set to the mean velocity − 0.2 × SD of the velocity. Events below threshold and local minima >1 SD of the velocity were identified when the inter-event intervals were limited to 200 ms. The detected movement termination train was used to calculate peri-LFP spike-triggered stop histograms.
Event triggered velocity.
The spike-triggered velocity was calculated for a period of 15 min starting 10 min after the initiation of pre- and postinjection states. The single unit/LFP spike train and its corresponding velocity (smoothed by a 100 ms SD Gaussian window) were synchronized. The event-triggered velocity within a time window of ±5 s was calculated by averaging the z-scores of the velocity around individual events. Confidence levels were calculated by a baseline period located in the first 1000 ms of the window. The relationship to behavior for significant event-triggered velocities was determined by the deviation of the velocity around the event onset (positive deviation, movement initiation; peak, peak velocity; negative, movement termination).
Statistical analysis.
The statistical hypotheses were assessed using either normal distribution tests or nonparametric tests in cases of non-normal distributions (such as intervals) depending on the specific dataset. The values are presented as either mean ± SD or mean ± SEM as stated in the text. The resulting p values are shown in raw format. In all cases, the tests were performed using their standard MATLAB implementation without any modifications.
Histology
At completion of experimental sessions, the animals were deeply anesthetized using a mixture of ketamine (100 mg/kg), xylazine (10 mg/kg), and morphinhydrochlorid (0.15 mg/kg). The microinjection location within the NAc core was marked using an injection of 0.5 μl of 1% toluidine blue and was used as the inclusion criterion. The rats were transcardially perfused with 0.9% saline and 10% paraformaldehyde and their brains were removed and stored in a fixation solution of 30% sucrose and 10% paraformaldehyde. The implanted structures were preserved and examined to confirm their relative location using 50 μm coronal sections.
Results
Behavioral effect: hyperactivity
The GABAA antagonist bicuculline was unilaterally injected into the NAc core of freely behaving rats (Fig. 1A). We conducted a total of 34 recording sessions in 15 rats, which included a naive (pre-injection) period followed by the microinjection, and the subsequent postinjection period. The movement of the rats was assessed with a combination of automatic detection using a single video stream and 3D motion sensors (accelerometers and gyroscopes). The injection led to a hyperactive state, which was expressed as increased normal behavior without abnormal motor symptoms (Movie 1). The overall ambulation of the rat increased significantly following the injection (Fig. 1B). The overall distance traveled increased significantly (n = 23, pre: 54.9 ± 5.9 m, post: 156.8 ± 13.7 m, mean ± SEM over a period of 25 min, p = 3.6 × 10−9, paired t test). This increase was highly dependent on the activity of the animal before the injection, such that highly active naive animals tended to increase their activity after the injection to a larger extent than less active ones (Pearson's correlation coefficient: 0.65, p = 0.0007; Fig. 1C). The bicuculline injection led to a prolonged exploration period in which the animals continuously ambulated in the arena. This contrasted sharply with the naive state, in which following a short period of exploration the level of ambulation decreased (Fig. 1D–F). Following the injection, the animals spent a higher fraction of time in active periods (Fig. 1G; n = 23, Pre: 0.48 ± 0.05, Post: 0.90 ± 0.02, mean ± SEM over a period of 25 min, p = 8.4 × 10−3, paired t test) and moved faster during these active periods (Fig. 1H; n = 23, Pre: 3.65 ± 0.15 cm/s, Post: 6.41 ± 0.46 cm/s, mean ± SEM over only the active periods, p = 6.6 × 10−7, paired t test) compared with the pre-injection state. The animals' velocity postinjection was positively correlated with their velocity before the injection (Pearson's correlation coefficient on logarithmic values: 0.59, p = 0.0032; Fig. 1H). This hyperactivity was expressed bilaterally, thus, the animals did not display a preference for the ipsilateral or contralateral sides of the injection hemisphere (turns/min: ipsilateral: 22.27 ± 0.99; contralateral: 23.10 ± 1.04, mean ± SEM, paired t test: p = 0.56; Fig. 1I). In a few sessions, the animals developed short tetanic seizures or unilateral symptoms such as circling ipsilateral to the injection side, and motor tics contralateral to the injection hemisphere, potentially due to diffusion of bicuculline along the cannula (McCairn et al., 2009). Sessions in which the animals expressed motor tics (n = 4) and sessions with no apparent behavioral effect (probably because of malfunction of the injection process, n = 7), were not included in the dataset.
Behavioral effects of NAc core disinhibition. A, The location of bicuculline injection within the NAc marked on a coronal section of the rat brain (right) and the corresponding atlas-based section (AP: 1.44; left). B, Rat tracks >25 min pre- (left) and post- (right) injection of a single session. C, The overall distance traveled over 25 min pre- and post-injection during all sessions (dashed red line, linear regression). D, Cumulative distance traveled pre- (blue) and post- (black) injection during a single session (dashed black line, point of separation). E, The rat track pre- (blue) and post- (black) injection separated to before (left) and following the (right) point of separation. F, Velocity over time during all the sessions, pre- (left) and post- (right) bicuculline injection normalized to the maximal velocity within each session. G, H, Relation between pre- and post-injection periods: (G) fraction of active state of the whole session and (H) median velocity. The solid black lines in C, G, and H mark the identity relationship. I, Number of ipsiversive and contraversive turns postinjection calculated for a period of 25 min for all sessions.
Effect of NAc bicuculline injection on behavior and neuronal activity. Behavior of the rat during the pre- (naive) and post- (hyperactive) injection states in the same session. A synchronized example of a single-channel LFP and spiking activity and the corresponding accelerometer and gyroscope signals, which are overlaid on the video stream.
Neuronal correlates: hyperactivity
The wideband filtered extracellular electrophysiological signal was continuously recorded from the striatum throughout the session. For the analysis of single-unit activity, the wideband signal was off-line high-pass filtered and single-unit spike times were extracted. Based on the action potential waveform, firing rate and pattern, the single units were classified as either spiny projection neurons [SPNs; also termed medium spiny neurons; Pre: n = 106, Post: n = 91], tonically active neurons (TANs; Pre: n = 31, Post: n = 19), or fast-spiking interneurons (FSIs; Pre: n = 13, Post: n = 17). The overall activity patterns of individual neurons and their interactions did not change following the injection despite the dramatic increase in the animals' behavior. The firing rates (Pre: SPN: 0.85 ± 0.09 TAN: 5.06 ± 0.52 FSI: 25.38 ± 3.03, Post: SPN: 1.15 ± 0.13 TAN: 6.07 ± 0.89 FSI: 23.17 ± 2.56 spikes/s, mean ± SEM; t test: SPN: p = 0.06, TAN: p = 0.3, FSI: p = 0.6; Fig. 2A), firing patterns, as assessed by the coefficient of variation (CV; Pre: SPN: 2.45 ± 0.09 TAN: 2.17 ± 0.13 FSI: 1.69 ± 0.23, Post: SPN: 2.55 ± 0.10 TAN: 2.10 ± 0.19 FSI: 1.55 ± 0.13, mean ± SEM, t test: SPN: p = 0.43, TAN: p = 0.78, FSI: p = 0.58; Fig. 2B) and autocorrelation functions (Fig. 2C) were unaltered following the bicuculline injection. In line with the preservation of single-unit activity properties, the inter-striatal correlations were not significantly modulated by the bicuculline injection (SPN-SPN: Pre: n = 315 pairs, Positive correlation: 93, negative correlation: 6; Post: n = 396 pairs, Positive correlation: 127, negative correlation: 4; SPN-TAN: Pre: n = 231 pairs, Positive correlation: 56, negative correlation: 13, Post: n = 175 pairs, Positive correlation: 29, negative correlation: 14; SPN-FSI: Pre: n = 84, Positive correlation: 18, negative correlation: 8, Post: n = 175 pairs, Positive correlation: 32, negative correlation: 14; χ2 values: SPN-SPN: 1.43, SPN-TAN: 4.02, SPN-FSI: 0.62, for all groups: p > 0.05; Fig. 2D).
Bicuculline-induced changes in single neuron activity and interactions. A, Mean firing rates, (B) CVs, and (C) autocorrelation functions of SPNs (blue), FSIs (red), and TANs (green), pre- (light colors) and post- (dark colors) bicuculline injection (error bars, 1 SEM). D, Cross-correlation functions of all the simultaneously recorded single-unit pairs pre- (left) and post- (right) injection. Significance lines of positive (bottom) and negative (top) correlations are marked by dashed black lines.
In a sharp contrast to the preserved properties of individual neurons, the injection led to changes in the LFP, which were expressed as large deflections (Fig. 3A). These deflections, termed LFP spikes, were characterized by a stereotypic shape which gradually developed over a few minutes following the injection, stabilized for a prolonged period, and finally decayed (Fig. 3B). The shape of the LFP spike was highly stereotypic during the stable period (Fig. 3C). During the stable periods, the LFP spike occurrences typically followed a Poisson distribution, which was reflected as an exponential distribution of the inter LFP spike intervals, with a refractory period (Fig. 3D). The rate of the LFP spikes varied across animals with a mean rate of 0.44 Hz (Fig. 3E). The rate of LFP spikes was highly reflective of the animals' state: the LFP spike rate and the behavioral effect strength (ratio of the distance traveled postinjection to the distance traveled pre-injection) were highly correlated (Pearson's correlation coefficient, r = 0.83, p = 0.001; Fig. 3F).
Bicuculline-induced changes in the LFP signal. A, Example of the LFP signal pre- (top) and post- (bottom) bicuculline injection. B, Shapes of all the LFP spikes recorded over a single session. C, Mean shape of the LFP spike during the stable period (solid black line) from the session shown in B. D, Inter-LFP-spike interval histogram taken from a single session. E, Medians of the inter-LFP-spikes interval histograms (black line, mean of medians). F, LFP spike rate as a function of behavioral effect strength (dashed line, linear regression).
The activity of individual neurons was phase-locked to the LFP spikes across all the subpopulations (Fig. 4A,B). Although all the recorded neurons displayed a modulation in their activity, the temporal structure of the modulation was diverse, including multiple phases of excitation and inhibition with a variety of time constants (Fig. 4B). In most cases, the neurons displayed a prolonged inhibition, which was preceded, in some cases, by a short early excitation. This typical response was also evident in the mean response across the different neuronal populations (Fig. 4B). Despite the common coactivation with the LFP spikes, the neuronal activity was not correlated in its responses. We calculated the JPLTHs of simultaneously recorded pairs of single units (n = 230 pairs with non-sparse matrices; Fig. 4C). Only a small fraction displayed correlated activity around the time of the LFP spike (15/230, 6.5%; Fig. 4D) indicating sparse locked neuronal activation around individual LFP spikes. These results demonstrate the translation of LFP information into a sparse finely-timed representation in the single-unit activity which is sufficient to drive the significant behavioral changes without changing in the striatal firing rates.
Single-unit activity around LFP spikes. A, Examples of peri-LFP-spike raster plots (top) and histograms (bottom) of individual SPN (top, blue), TAN (middle, green), and FSI (bottom, red) neurons. The inset within each histogram shows the mean spike shape using zero-phase (left) and nonlinear phase (right) filters. B, Matrices of single-unit peri-LFP-spike histograms (top) and their mean shape (bottom). C, Joint peri-LFP-spike histogram. Scale: Y, 2 SDs; X, 200 ms. D, Co-response types based on the JPLTHs calculated for all pairs.
Microscale behavior and neuronal activity
On the macro-scale, NAc core disinhibition caused an increase in normal behavior prolonging the active state. A closer examination of the dynamics of behavior revealed a complex temporal microstructure consisting of bouts of rapid movements separated by brief cessation of locomotion, which was apparent in both the video-based and the sensor-based movement analysis (Fig. 5A). The neuronal activity was correlated to this micro-structure in that individual LFP spikes appeared at the time of locomotion termination (Movie 2; Fig. 5B). The peri-LFP spike movement stop histograms (Fig. 5C) and the LFP spike-triggered velocity (Fig. 5D) reflected this correlation between the behavioral changes and the LFP signal. The correlation between the LFP spikes and movement stops was only partial, whereas not all of the LFP spikes were accompanied by movement termination (Fig. 5B). To test whether the LFP spikes that were related or unrelated to movement termination were similar, we compared the temporal structure of the two populations (Fig. 5E), and calculated their mean shape (F) and interval distribution (G,H). None of these measures was significantly different for the two conditions (Interval distribution, Wilcoxon signed rank test, p = 0.05; Mean shape: t test performed on the first principal component, p > 0.05 for all cases), suggesting that they represented events emerging from similar neuronal populations.
Micro-scale behavior and neuronal activity. A, Example of the video-based velocity and sensor-based acceleration pre- (blue) and post- (black) injection. B, Example of the animal's velocity (blue, top) and the corresponding LFP signal (black, bottom). Arrows mark the LFP time related (black arrow) and unrelated (gray arrow) to movement termination. C, Peri LFP-spike movement-termination histogram. D, LFP spike-triggered velocity (top) and the mean LFP triggered velocity (bottom). E, LFP spike shapes separated into movement-termination related (blue, bottom) and unrelated (black, top) and their (F) mean shape (related, blue; unrelated, black). G, Histograms of the intervals between movement-termination related (blue) and unrelated (black) LFP spikes and their (H) medians calculated for 11 sessions (red, mean of the sessions).
Behavior microstructure and its relationship to the LFP. The time series representing the LFP spike times is presented as a series of auditory beeps, which are overlaid on the video stream.
The behavioral micro-structure was accompanied by changes in the activity of individual neurons. Whereas the LFP spikes were highly associated with movement termination, individual neuron activity was correlated with additional components of movement, including movement initiation and the peak velocity (Fig. 6A,B) covering multiple parts of each of the behavioral bouts. Correlations of this type have only been reported in neurons from the sensorimotor territories directly encoding motor behavior (Schultz and Romo, 1988; Cui et al., 2013) and were expanded to describe discrete actions or behavioral syllables (Klaus et al., 2017; Markowitz et al., 2018). This relationship to different behavioral components was found in all the striatal subpopulations (SPN: n = 66 related: 36, 54%; TAN: n = 23, related: 13, 56%; FSI: n = 31, related: 18, 58%; fraction of movement related neurons: SPN, initiation: 14% peak velocity: 56%, termination: 22%, other: 8%; TAN: initiation: 15%, peak velocity: 46%, termination: 31%, other: 8%; FSI: initiation: 6%, peak velocity: 61%, termination: 22%, other: 11%; Fig. 6C). The fraction of neurons related to different parts of behavior did not differ pre- and post-bicuculline injection (Pre: n = 58 neurons, movement related: 35, 59% movement unrelated: 23, 41; Post: n = 64, movement related: 34, 53%, movement unrelated: 30, 47%; χ2 test = 0.65, p > 0.05; Fig. 6D). Thus, despite the major change in behavior following the injection, the encoding of the behavioral microstructure by individual neurons did not change.
Single-unit activity relation to micro-scale behavior. A, Spike trains of three ventromedial striatum neurons plotted below the locomotion velocity. B, Velocity aligned to the spike times (top) and the mean spike-triggered velocity (bottom) of the spike trains of the session shown in A showing that each neuron encodes a different component of movement (i, movement initiation; ii, peak velocity; iii, movement termination). C, Fraction of neurons associated with different components of movement, grouped by cell type (blue, SPN; green, TAN; red, FSI). D, Fraction of movement related and unrelated neurons pre- (left) and post- (right) injection.
Single-unit activity was correlated with both LFP spikes and behavior. We explored this three-way interaction by partitioning the LFP spikes based on their relationship to movement termination and calculating the peri-LFP spike histograms for each group. We examined the change in the spiking activity in each PLTH given the relationship between neuronal activity and behavior (Fig. 7A,B). A minority (12/39, 31%) of the neurons differed in their LFP spike locked activity between the two LFP groups, indicating a behavioral context dependency (Fig. 7A,C). The largest group from out of these neurons was associated with movement termination (Fig. 7D). Most of the neurons displaying LFP spike-locked activity were independent of movement termination (69%, 27/39; Fig. 7B,C). These neurons were associated with other components of movement, and only a small minority was related to movement termination (Fig. 7D), suggesting that individual neurons encode both normal movement and abnormal LFP spike-dependent changes independently of each other.
Three-way interaction between single-unit activity, LFP spikes and micro-scale behavior. A, B, Examples of behavioral context-dependent (A) and independent (B) LFP-spike locked activity of two single units illustrated by the peri-movement termination related (left) and unrelated (middle) LFP-spike histograms and the corresponding spike-triggered velocity (right). C, Fraction of behavioral context dependent and independent LFP-spike locked activity. D, Fraction of the neurons related to different components of movement in each group (green, behavioral context-dependent; gray, behavioral context independent LFP-spike locked activity).
Discussion
In this study we explored the effect of NAc core disinhibition on the rat behavior and the underlying neuronal activity. We used a wireless acquisition system, which enabled the animals to move unrestrained within a large behavioral arena while continuously recording both electrophysiological and kinematic signals. This approach enabled precise synchronization of the different constituents of the neurophysiological signal and the corresponding behavior on multiple temporal scales. On the macro-scale, disinhibition resulted in hyperactivity without any abnormal movements: the animals moved within the behavioral arena without the transition to rest typical of naive animals. The magnitude of hyperactivity was animal-dependent, because the activity level during the hyperactive macro-state was dependent on the activity baseline level during the naive state of the same animal. The macro-scale behavioral state transition was not associated with major changes in either the tonic activity of individual neurons or their interactions. A major change, however, was observed on the neuronal population level: the LFP displayed large deflections (LFP spikes) whose frequency was linearly related to the magnitude of the behavioral effect. On the micro-scale, the rats' activity did not consist of continuous locomotion but rather included transitions between movement and non-movement microstates (episodes). The transitions between the microstates tended to be encoded by the neuronal population as individual LFP spikes accompanied by finely-tuned sparse activity modulations of single neurons, encoding different microstates related to movement while adhering to the same distribution as seen in the neuronal population of the naive animals.
Although the pathophysiology underlying ADHD is still largely unknown, multiple converging evidence from human patients studies point to the prominent role of the corticobasal ganglia pathway (Giedd et al., 2006). Anatomical evidence indicates a decrease in the volume of the striatum and in particular the ventral striatum (Carmona et al., 2009; Qiu et al., 2009; Greven et al., 2015). Physiological evidence (fMRI) suggests differences in activation (Durston et al., 2003) and responsiveness (Scheres et al., 2007; Plichta and Scheres, 2014) of these brain areas. In line with these results, early studies in animals showed that experimental manipulation of the NAc core changes their overall activity. Local application of GABA antagonists (Wachtel and Andén, 1978; Morgenstern et al., 1984), as well as glutamate agonists (Arnt, 1981) induces hyperactivity whereas glutamate antagonists (Mogenson and Nielsen, 1984; Swanson and Kalivas, 2000) and GABA agonists (Jones et al., 1981; Wong et al., 1991) induce hypoactivity.
Hyperactivity, a key symptom in many ADHD patients (American Psychiatric Association, 2013), was expressed in this study as enhanced normal behavior including locomotion, exploration, and other behaviors which were indistinguishable from the repertoire of the naive animals. This persistent hyperactivity expressed for a prolonged period following the disinhibition differed from the short-term novelty-associated hyperactivity of spontaneous hypertensive rats (Knardahl and Sagvolden, 1979). The normal movements expressed following disinhibition differed from the abnormal movements expressed by animal models after dopaminergic modulation, which are associated with dyskinesia and stereotypy (Fog, 1970; Costall and Naylor, 1976). The hyperkinetic symptoms subsequent to disinhibition were expressed solely as hyperactivity and there were neither abnormal individual movements nor abnormal movement sequences. The ongoing behavior during the macro-scale hyperactivity state was assembled from multiple micro-scale transitions between movements separated by brief non-movement periods. Frequent shifts between activities are common in many ADHD patients and are related to difficulties in sustaining attention. Similar attention-deficit-like symptoms have been reported in primates using disinhibition in the globus pallidus externus (GPe; Grabli et al., 2004). Transitions between behavioral states also exist during these animals' normal behavior, as part of the normal repertoire (Markowitz et al., 2018). This supports the notion that behavioral transitions, which are attributed to the attentional state, may merely be another aspect of normal behavior, which is overexpressed during ADHD.
Disinhibition of different functional territories of the basal ganglia leads to a wide range of excessive movement (hyperkinetic) and behavioral (hyperbehavioral) symptoms (Tarsy et al., 1978; Morgenstern et al., 1984; Grabli et al., 2004; Worbe et al., 2009). Hyperactivity in the form of increased locomotion has been reported in studies using disinhibition of limbic territories in rats, including the NAc (Wachtel and Andén, 1978; Jones et al., 1981; Morgenstern et al., 1984) and the ventral pallidum (Mogenson and Nielsen, 1983; Austin and Kalivas, 1990), and in primates (Worbe et al., 2009). In addition, disinhibition of the NAc shell led to emotional dysregulation (Stratford and Kelley, 1997; Lopes et al., 2012), which may contribute indirectly to the observed changes in locomotion. Multiple different hyperkinetic states have been evoked using similar manipulations in the sensorimotor basal ganglia: disinhibition of the sensorimotor striatum leads to motor tics (Tarsy et al., 1978; McCairn et al., 2009; Bronfeld et al., 2013; Pogorelov et al., 2015; Klaus and Plenz, 2016), downstream disinhibition in the sensorimotor GPe leads to chorea (Grabli et al., 2004; Bronfeld et al., 2010) and disinhibition of the striatum in cats (Yamada et al., 1995), and the motor thalamus in primates (Guehl et al., 2000) leads to dystonia. The disinhibition of sensorimotor territories share common properties: they lead to the expression of different hyperkinetic symptoms in the form of abnormal movements, which are not part of the normal animal's repertoire, are specific to body part, and are expressed unilaterally contralateral to the injected hemisphere. Similar manipulations in the limbic part of the striatum lead to hyperactivity expressed as increased normal behavior, which is nonspecific to a body part and is expressed bilaterally. Disinhibition of the sensorimotor territory in the striatum was shown to lead to the appearance of LFP spikes and phase-locked individual neuronal activity which were correlated with the appearance of individual tics (Bronfeld et al., 2013; Israelashvili and Bar-Gad, 2015; Pogorelov et al., 2015; Klaus and Plenz, 2016). The timing of individual LFP spikes was correlated with unitary behavioral events, individual tics, resembling the correlation to individual movement terminations in the current study. Our results support the hypothesis that striatal disinhibition leads to the release of unwanted motor plans (Albin and Mink, 2006). The disinhibition hypothesis, which in the sensorimotor striatum results in the expression of motor tics, may be extended to the NAc in which similar loss of inhibition leads to behavioral switches and the increased performance of motor plans. The source of this abnormal disinhibition, i.e., SPN collaterals (Chuhma et al., 2011), FSI feedforward inhibition (Koos et al., 2004), or GPe feedback inhibition (Mallet et al., 2012), however, cannot be distinguished by this study. The similarity of network manipulation across the striatal territories and the resulting underlying neuronal activity fits well with the high comorbidity (∼50%) of tic expression and hyperactivity in Tourette syndrome patients (Schlaggar and Mink, 2003), hinting at similar fundamental neuronal changes underlying the two sets of symptoms. These different symptoms also support the functional distinction between territories in which the sensorimotor territory is associated with the form of individual movement primitives (Wiltschko et al., 2015; Klaus et al., 2017; Markowitz et al., 2018), whereas the limbic territory is associated with the number of movement primitives and switches between them.
The mechanism by which GABAergic disinhibition leads to macro- and micro-scale changes in behavior may be attributed to changes in the properties of either (1) the intrastriatal network, or (2) the extra-striatal inputs or their processing by the striatal neurons. The lack of changes in macro-scale activity or the interaction of the striatal neurons, together with similar neuronal encoding on the micro-scale do not support the likelihood of tonic alteration of the intrastriatal network as the basis for macro-scale state change. Studies manipulating the tonic inputs to the NAc have reported similar behavioral results using both glutamatergic (Arnt, 1981) and dopaminergic (Jones et al., 1981) agents. These results lend weight to the notion of a mechanism based on alteration of the input to the striatum or on the responsiveness to such input. This mechanism was addressed by an in vitro study investigating the NAc synaptic transmission using extracellular field recordings (Buckby and Lacey, 2001), which demonstrated that cortical input forms the principal drive of regenerative field changes, resembling the observed LFP spikes in our study following administration of GABAA antagonists. A recent study demonstrated, in vivo, that the timing of unitary behavioral events associated with LFP spikes (motor tics) was determined by the cortical input to the striatal network (Israelashvili and Bar-Gad, 2015), supporting the idea of an extra-striatal input-based mechanism. However, further studies are required to fully determine the source of the LFP spikes and their relation to the observed symptoms.
Disinhibition of the NAc core leads to neuronal changes consistent with those observed after disinhibition of the dorsal striatum but with different behavioral symptoms; i.e., macro-scale hyperactivity and micro-scale activity switches. These results emphasize the unique role of the NAc in controlling the overall level of activity and the switches between movement primitives, a function that is synergistic with the control of the individual movement primitives by the dorsolateral striatum. These complementary roles in normal behavior may also underlie the expression of the very different, yet highly comorbid symptoms of Tourette syndrome and ADHD.
Footnotes
This study was supported in part by an Israel Science Foundation Grant (297/18) and a BSF-NSF Collaborative Research in Computational Neuroscience Grant (2016744).
The authors declare no competing financial interests.
- Correspondence should be addressed to Izhar Bar-Gad at izhar.bar-gad{at}biu.ac.il