Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Directed Communication between Nucleus Accumbens and Neocortex in Humans Is Differentially Supported by Synchronization in the Theta and Alpha Band

  • Jörn M. Horschig ,

    jhorschig@gmail.com

    Affiliation Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands

  • Ruud Smolders,

    Affiliation Department of Psychiatry, Academic Medical Center, Amsterdam, The Netherlands

  • Mathilde Bonnefond,

    Affiliation Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands

  • Jan-Mathijs Schoffelen,

    Affiliations Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands, Max-Planck Institute for Psycholinguistics, Nijmegen, The Netherlands

  • Pepijn van den Munckhof,

    Affiliation Department of Neurosurgery, Academic Medical Center, Amsterdam, The Netherlands

  • P. Richard Schuurman,

    Affiliation Department of Neurosurgery, Academic Medical Center, Amsterdam, The Netherlands

  • Roshan Cools,

    Affiliations Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands, Department of Psychiatry, Radboud University Medical Center, Nijmegen, The Netherlands

  • Damiaan Denys ,

    ‡ These authors are senior authors on this work and also contributed equally to this work.

    Affiliations Department of Psychiatry, Academic Medical Center, Amsterdam, The Netherlands, The Netherlands Institute for Neuroscience, Royal Netherlands Academy of Arts and Science, Amsterdam, The Netherlands

  • Ole Jensen

    ‡ These authors are senior authors on this work and also contributed equally to this work.

    Affiliation Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands

Abstract

Here, we report evidence for oscillatory bi-directional interactions between the nucleus accumbens and the neocortex in humans. Six patients performed a demanding covert visual attention task while we simultaneously recorded brain activity from deep-brain electrodes implanted in the nucleus accumbens and the surface electroencephalogram (EEG). Both theta and alpha oscillations were strongly coherent with the frontal and parietal EEG during the task. Theta-band coherence increased during processing of the visual stimuli. Granger causality analysis revealed that the nucleus accumbens was communicating with the neocortex primarily in the theta-band, while the cortex was communicating the nucleus accumbens in the alpha-band. These data are consistent with a model, in which theta- and alpha-band oscillations serve dissociable roles: Prior to stimulus processing, the cortex might suppress ongoing processing in the nucleus accumbens by modulating alpha-band activity. Subsequently, upon stimulus presentation, theta oscillations might facilitate the active exchange of stimulus information from the nucleus accumbens to the cortex.

Highlights

  • Both theta (4–7 Hz) and alpha (9–14 Hz) facilitate corticostriatal interactions.
  • The nucleus accumbens communicates with cortical oscillations mainly in the theta-band.
  • Alpha oscillations represent control from the neocortex over nucleus accumbens.
  • Upon stimulus processing, the theta-band connectivity increases.

Introduction

There is growing evidence for a key contribution of communication between striatum and neocortex in human cognition, attention and behavior. However, little is known about the neuronal dynamics supporting such cortico-striato-cortical communication. Recordings from deep brain electrodes implanted in humans for the treatment of psychiatric disease provide a rare opportunity to investigate such subcortico-cortical dynamics. Deep-brain stimulation of the nucleus accumbens, a part of the ventral striatum, has been successful in treatment of refractory-resistant obsessive compulsion disorder [13], major depressive disorder [49] and also drug addiction [10,11]. Earlier studies reported on theta (4–7 Hz) and alpha (9–14 Hz) band oscillations in the human nucleus accumbens [1215], as well as top-down directed synchrony between nucleus accumbens and frontal electrodes in the low frequencies between 1 and 10 Hz [16]. However, none of these studies explored the differential roles of alpha and theta oscillations in communication between neocortex and the nucleus accumbens.

Anatomically, the nucleus accumbens receives input from different neocortical areas, mainly from temporal and prefrontal areas. Indirectly the nucleus accumbens projects back to the prefrontal cortex and parietal cortex via the globus pallidus, subthalamic nucleus and the thalamus [1719]. Computational models by Frank & O’Reilly [2022] propose that the striatum acts as a gatekeeper by deciding whether stimulus information should or should not be passed on to the prefrontal cortex. Moreover, the nucleus accumbens has recently been suggested to be involved in detection of visual information, and actively modulates the degree of neocortical frontoparietal connectivity [23]. In this study, we investigated whether the nucleus accumbens, as part of the ventral striatum, and the neocortex are employing theta and alpha band oscillations in anticipation and during processing of visual stimulation, and whether the nucleus accumbens contributes to frontoparietal connectivity in these two frequency bands.

Materials and Methods

Participants

Seven right-handed patients (one male diagnosed with chronic major depressive disorder, one male with cocaine and opiate addiction, and five female and one male patient with obsessive-compulsive disorder; 22–55 years of age).participated in the experiment. The experiment was approved by the local Medical Ethical Committee of the Academic Medical Center, University of Amsterdam. All patients provided written informed consent according to the Declaration of Helsinki and the local Medical Ethical Committee of the Academic Medical Center, University of Amsterdam prior to the experiment. All patients underwent surgery for implantation of deep brain electrodes in bilateral nucleus accumbens (NAc) between 2010 and 2012 (see [1,3] for more information about the exact procedure). The most ventral contact point was located in the core of the NAc, and the three other contact points were extending into the ventral part of the anterior limb of the internal capsule. Data from the male OCD patient was unsuited due to excessive movement during the experiment, and low signal-to-noise ratio particularly in intracranial electrodes.

Stimulus presentation and experimental paradigm

Equipment.

Stimulus presentation was performed using Presentation (Version 14.5; Neurobehavioural Systems, Inc.) and a laptop (HP 6730b) computer screen on a 15.4 inch display at a resolution of 1024 by 768 pixels (refresh rate of 60 Hz). The distance from the screen to the participants was kept around 60 cm.

Paradigm.

We adapted a covert attention switching paradigm also described in [24], see Fig 1. Squares were flashed on each side and subjects had to report the color of the attended square by a button press. When subjects detected a color change at the unattended side (signaling a switch trial), they had to report the color of the unattended square (but not the currently attended square) and switch attention to the unattended side in future trials.

thumbnail
Fig 1. The paradigm.

The attended side was initially indicated by a cue. Subjects had to fixate at the central cross and by button press indicate the color of the attended squares (left button for red and right button for blue). The 1200 ms prestimulus period was followed by the colored stimuli flashed for 33 ms. Subjects had to respond within 2500 ms. If there was a color change in the square of the unattended hemifield, attention had to switch to that direction (‘switch-trial). After the response there was a 1000 ms window for eye blinking. A Example of an explicit cue followed by a repeat trial. The subject had to covertly attend to the left and subsequently report the color of the stimuli by pressing the corresponding button (here: blue, right button). B Example of a switch trial. In the previous repeat trials, the subject had to attend to the left, because of the initially shown spatial cue. Upon stimulus presentation, the subject correctly switched attention and indicated so by reporting the color of the stimulus at the formerly unattended side (here: right, red color). If the subject responded according to the formerly attended side (here: left, blue), the switch trial would repeat up to four times. Repetitions of switch trials were removed from the analysis. If the subject did not switch after the fourth consecutive switch trial, another explicit spatial cue pointing to the formerly unattended side was presented (here: a rightward pointing arrow).At the beginning of each block, subjects were explicitly cued to which side to attend. From then on, the attended side was determined by stimuli properties alone. A central fixation point was presented during the entire experiment. Colored squares were flashed 1200ms after the beginning of each trial for about 33 ms (two frames = 2/60Hz). These stimuli were presented with six degrees eccentricity and two degrees lower than the fixation cross (measured from the fixation cross to the center of the stimuli). The squares were two degrees wide.

https://doi.org/10.1371/journal.pone.0138685.g001

Subjects had to report the color of the square on the attended side by pressing a button with their left (for red) or right hand (for blue). On the unattended side, the square was either grey (repeat trials) or colored in blue or red (switch trials). Subjects had to respond within 2500 ms. After responding, the fixation cross turned grey, indicating that the subject could blink or move the eyes in a 1000 ms period. Then the fixation cross turned white again indicating the start of the next trial. Subjects had to keep attention to one hemifield (repeat trial) and report the color of the square on that side until they detected a colored stimulus in the unattended hemifield (switch stimulus). The switch stimulus was detected if the color of the unattended target was correctly reported (detected switch trial). Subjects then had to keep attending the formerly unattended hemifield until a next switch trial was detected. If a subject failed to detect the switch stimulus (undetected switch trial), it was repeated with a random color (blue or red) up to four times. The probability of a switch trial was increasing with the number of trials since the last switch trial, but is not of further importance for this study (see [24] for more information).

Stimulus intensity.

In order to make the task sufficiently difficult, the intensity of the stimuli was varied across subjects and early trials in an adaptive staircase-like procedure on a 20-step scale (1:darkest; 20:brightest), starting at 10 for both the neutral and the colored stimuli. This was done in the first block (i.e. until the 15th detected switch trial). Trials from this block were not included in the later analyses. The intensities of the neutral and colored stimuli were modulated according to different criteria. Repeat trials should be sufficiently demanding while keeping response errors as low as possible. Therefore we reduced the brightness of the colored stimuli to a level in which the subject could perform the discrimination with no errors. The color intensity was reduced by one step after each detected switch trial and error-free performance. After a response error to repeat trials, the intensity was increased by one step again. The intensity of the neutral stimulus was adapted to manipulate the difficulty of switch trials. A large intensity difference between colored and neutral stimuli makes detection of switch trials easier (pop-out effect), whereas a similar intensity results in a harder task and less detected switch trials. We aimed at a correct response rate to switch trials between 25% and 75%. After a detected switch trial, the neutral stimulus was increased in intensity by one when less than 25% of all switch trials were detected. The intensity was decreased by one step when more than 75% of all switch trials were detected. The twenty levels of stimulus luminance were visually and mathematically matched according to the CIELAB specifications [25,26]. Note that the adaptation procedure was only applied in the first block until the 15th detected switch trial.

Instructions.

Prior to the experiment, participants received written and verbal task instructions. Subjects were instructed to prioritize accuracy rather than speed, but were informed that they should respond within 2500 ms. Subjects were instructed to detect the color at the cued side, but switch attention to the uncued side if the color at that side turned from grey to either blue or red. Subjects were informed that they would receive no response feedback. The instructions did not inform about other task properties. After the instructions, subjects had to complete a short tutorial on the computer, which explained the paradigm and introduced the stimuli. Afterwards the experimental task started. In our analysis, as already reported above, we discarded all trials up to the first experimental break.

Data acquisition

The ongoing brain activity was recorded using a 64-EEG-electrode recording system (Advanced Neuro Technology B.V.) at a sampling rate of 512Hz following the international 10–10 system. Eight channels were used to record the deep brain electrodes near the nucleus accumbens (four electrodes per lead), four to derive horizontal and vertical eye-movement, leaving 52 Ag/AgCl electrode distributed evenly on the scalp while leaving out areas covered by post-surgery bandages. Data was online common average referenced and offline re-referenced to the common electrodes across patients. For the final analysis, only electrodes that were common for all patients were included. These were the following 42 electrodes: AF7, AF8, C1, C2, C3, C4, C5, C6, CP1, CP2, CP3, CP4, CP5, CP6, CPz, Cz, F1, F7, F8, FC5, FC6, FCz, FT7, FT8, Fp1, Fp2, Fpz, Fz, O1, O2, Oz, P1, P2, P3, P4, P5, P6, PO3, PO4, PO7, POz and Pz. Purely for visualization purposes, missing electrodes were interpolated using a spherical spline interpolation [27].

Data preprocessing and analysis

Preprocessing.

The EEG data were analyzed using the Matlab-based FieldTrip toolbox, developed at the Donders Institute for Brain, Cognition and Behaviour [28]. We first visually inspected activity in all electrodes and rejected electrodes with a low signal-to-noise ratio. Then, we re-referenced all electrodes to a common average reference. Muscle and ocular artifacts were detected in a semi-automatic fashion, which included visual inspection and trial rejection based on variance and other measures as implemented in FieldTrip. We rejected all trials with an artifact in the interval -1000 ms to 500 ms relative to stimulus onset.

To estimate activity of the left and right nucleus accumbens, we computed the bipolar derivative of the most ventral and the second most ventral deep brain electrode (which was located in the core of the NAc) for the left and right lead, respectively. We will refer to these as left and right nucleus accumbens. For the main manuscript, only data from the left nucleus accumbens were used. Thus the term “nucleus accumbens” refers to the bipolar derivative of the electrodes in the core of the left nucleus accumbens and the most ventral electrode in the ventral part of the anterior limb of the internal capsule. Data from the right contact points showed generally more noise in most patients, and in one patient no data for the right nucleus accumbens was recorded. We therefore decided to conduct no further analysis on the data from the right nucleus accumbens.

Spectral analysis.

We computed Fourier coefficients from 0 Hz to 256 Hz for three time windows: the whole trial window from -1000 ms to 500 ms relative to stimulus onset, the prestimulus window from -500 ms to 0 ms, and the poststimulus window from 0 ms to 500ms. All data windows were zero-padded by 4500 ms to increase the spectral resolution and thereby smooth the Fourier estimates to avoid sudden transitions, especially around the 50 Hz noise-band. This resulted in a frequency resolution of 0.16 Hz for the whole trial window and 0.2 Hz for the pre- and poststimulus windows. For visualization, we show only the content from 2 to 18 Hz.

We computed the time-frequency representations (TFRs) of power from 0 to 256 Hz (0.5 Hz increments) for each trial from a -1000 ms to 500 ms interval around the stimulus onset. Spectral content was estimated using a 500 ms time window, which was multiplied with a Hanning window prior to applying a fast Fourier transform. For visualization, we only show the low frequencies up to 30Hz.

Coherence analysis.

We compute the squared coherence for each combination of the nucleus accumbens and extracranial EEG electrodes. Coherence Cxy(f) between two signals x and y is computed using the power spectra of both electrodes Sxx(f) and Syy(f) and their cross-spectrum Sxy(f) by the following equation [29]: (1)

Coherence is the frequency domain version of the linear correlation coefficient, and quantifies the consistency of the phase difference and amplitude correlations across observations. A coherence Cxy(f) of 0 indicates that the signals x and y have no consistent linear phase relationship, whereas a value of 1 indicate that the two signals are fully phase-coherent.

Non-parametric Granger causality.

Granger causality in the time-domain is well defined and has been well-established throughout the last century [30,31]. From the time-domain formulation of Granger causality, one can directly infer a formulation for non-parametric Granger causality in the frequency domain [3235]. We summarize the most important aspects of this derivation here and outline its relation to coherence as explained above.

A signal x is said to “Granger cause” signal y if the future of signal y can be better explained by incorporating knowledge about the past of signal x, above and beyond the prediction of the future of signal y based on its own past alone.

Total interdependence fx,y(f) between of two signals x and y is computed by their power spectra Sxx(f) and Syy(f) and their cross spectra Sxy(f) as (2) where ln denotes the natural logarithm, and S(f) = Sxx(f)Syy(f) − Sxy)(f)Syx(f). From Eqs (1) and (2), one can infer that (3)

Furthermore, total interdependence can be decomposed into three different causality terms: (4)

Where fx→y denotes the causal influence from signal x to signal y, fy→x denotes the causal influence from signal y to causal x and fx·y denotes the instantaneous causality, i.e. contributions by non-linear interactions between x and any contributions exogenous to signal x and y or by common input from a third signal.

From the derivations of [3233] it follows that (5) (6) where Σ2 and Γ2 denote the variance of the noise-term in the bivariate model of signal x and y, respectively (i.e. can be obtained from their noise covariance matrix Σ), denotes the spectral transfer function of signal x after normalization, and * the complex conjugate. The instantaneous causality is thus the difference between the total interdependence of x and y and the sum of the two causality terms fx→y and fy→x. Note that the spectral representations Sxx and Syy can also be fully estimated by the spectral transfer function Hxx and Hyy and the noise-covariance matrix Σ. For a more thorough derivation see [3233].

The normalized spectral transfer function of signal x can be derived using a matrix factorization method on the spectral representation of signal x [36]. Note that the spectral representations Sxx and Syy can also be fully estimated by the normalized spectral transfer function xx and yy and the noise-covariance matrix Σ.

Statistical analysis.

Statistical significance of the coherence of the neural data was assessed using a non-parametric cluster-based permutation test [37,38]. In the cluster-based permutation test, notational significant clusters in space are detected using a parametric test-statistic. We z-transformed the coherence values to correct for the variable number of degrees of freedom (trials numbers) across patients. In accordance to [38], if there are more than 20 degrees of freedom and if the squared population coherence is between 0.4 and 0.95, then tanh−1 (Cxy(f)) is approximately normal distributed with mean tanh−1 (Cxy(f)) − 1/(df − 2) and variance 1/(df − 2). Thus, a standard z-transformation can therefore be performed to obtain Zxy as follows: (7)

Note that although above mentioned requirements are not met in our data, this is no reason for concern as we applied the here described nonparametric test. This transformation was only done to account for the different amount of degrees of freedom across patients.

We compared the observed z-transformed coherence Zxy with a surrogate distribution of maximal z-transformed coherence cluster sizes. One instance of the surrogate distribution was computed as follows: For each patient, we randomly shuffled trials in the intracranial data, while remaining the order of trials of the EEG data. Then, we averaged over time and frequency of interest (here -1s to 0.5s, and 4 Hz to 14 Hz). After we computed the z-transformed coherence Zxy for each patient of the shuffled data, we computed grand-average z-transformed coherence across patients. If the summed cluster size exceeded 3.92, which is twice the threshold at alpha = 2.5%, we defined the cluster as a notational cluster. This was conducted 500 times to obtain a surrogate distribution. We then compared the size of notational clusters in the original observation with the size of notational clusters found in the surrogate distribution. We assessed significance in the original observation by a 5%-alpha level, i.e. 95% of the surrogate data had to contain smaller cluster sizes than those in the original observation for the test to be statistically significant. A cluster-based permutation test controls for multiple comparisons as only one comparison is made (the cluster size of the actual observation versus the surrogate distribution of cluster sizes).

Further statistical tests were performed by analysis of variance (ANOVA) and post-hoc paired t-tests. As parametric tests require normally distributed data, we log-transformed data prior to performing ANOVA or paired t-tests if data were not normally distributed as assessed by a Kolgomorov-Smirnov test. Effect sizes were computed by assessing η2 = SSeffect/SStotal and subsequently transformed to Cohen’s f (see e.g. [39]) as (8)

While η2 give an estimate of the explained variance, Cohen’s f is a standardized measure of effect sizes and widely used throughout different fields (psychometric testing, psychological tests, neuroscientific investigations, etc.)

In addition, we performed a permutation test on Granger directionality estimates. For each of the six subjects, we obtain one Granger estimate for the cortex-to-nucleus accumbens direction and one Granger estimate for the nucleus accumbens-to-cortex direction. For each subject, we subtracted these two directionality estimates to obtain the difference. In line with permutation testing, we permuted the directionality estimates for each subject, which effectively was achieved by swapping the sign of the difference in directionality. With two possible directions in six subjects, there are in total 26 = 64 permutations possible, of which one is the original observation. We computed the grandaverage difference in directionality for the original observation and for all other 63 permutations. Then, we assessed how many permutations resulted in a more extreme grandaverage difference in directionality with a binomial test. This number divided by the total number of permutations results in a p-value. We assessed statistical significance at an alpha level of 0.05.

Results

We investigated data from six patients (four with OCD, one with major depression disorder and one with drug addiction) with four-contact deep brain stimulation electrodes implanted bilaterally in the ventral anterior internal capsule, with the deepest of the four electrode contact points situated in the nucleus accumbens. Prior to deep brain stimulation, the deep brain electrodes were externalized, and we simultaneously recorded data from the intracranial electrodes and the scalp EEG. All the reported results are from the left nucleus accumbens. Data from the right nucleus accumbens were more noisy and incomplete in some patients. Patients performed a demanding attention task without any response feedback or reward, which did not include symptomatic related stimuli. Patients had to continuously attend one hemifield and report the color of shortly presented boxes, until a stimulus change at the unattended side occurred. The onset of the stimuli was predictable: after a 1200 ms prestimulus period, two boxes were shortly presented for 33 ms. Patients had to discriminate the color of the box on the attended side, while the unattended box was grey. In some trials, the box at the unattended side was colored as well, and the subjects had to switch attended side and report the color of the box presented on that side. From then on the patients had to attend that side. In this study we specifically investigated how theta-band and alpha-band oscillations supported the functional connectivity between the nucleus accumbens and the cortex, both in anticipation of and in response to the stimulus.

Behavioral results

The total duration of the experiment was 30 minutes. On average, the participants completed 446±(SD) 59 trials, of which 224±34 trials were on attention to the left and 223±26 on attention to the right. 120±29 trials were switch trials, and 326±39 trials were repeat trials. Participants responded correctly in 96±2% of all repeat trials, with no significant difference between attention left and attention right repeat trials (t(5) = -1.7, p = 0.16). Participants detected switch trials correctly in 74%±18% (switch rate) of all switch trials. Notably, participant 4 (the depressed patient) had a switch rate of only 48%, while all other participants had a switch rate between 76% and 94%. There was no significant difference between switch trials to the left or the right hemifield within participants (t(5) = 0.66, p = 0.54). These rates are comparable with the one we reported in [24], where healthy participants performed the same task in a more demanding version (further eccentricity of the stimuli) and correctly responded to 91±4.6% of all repeat trials and 64.2±19% of all switch trials.

Participants responded to repeat trials in 758±185 ms, with no significant difference between attention left and right hemisphere (t(5) = 0.54, p = 0.6). In switch trials, participants responded in 1119±278 ms with no significant difference between switch trials to the left and right hemifield (t(5) = 0.19, p = 0.86) and no significant difference between correct and incorrect responses to switch trials (t(5) = 0.09, p = 0.93).

Spatial topography of coherence

First, we investigated frontostriatal connectivity by computing coherence between the nucleus accumbens and the scalp EEG electrodes for the interval from -1000 ms to 500 ms relative to stimulus onset [29]. For all connectivity analyses, we pooled data of all trial types. Fig 2A shows two clearly distinct peaks, one in the theta-frequency range from 4 to 7 Hz and the other one in the alpha-frequency range from 9 to 14 Hz. The data were z-transformed to normalize over subjects prior to averaging. To assess whether the observed coherence was statistically significant, we shuffled the data over trials and recalculated the coherence and used the maximum coherence value in this surrogate distribution as a threshold. We identified coherence values between 4 Hz and 14 Hz to be above this threshold calculated from the shuffled data. To statistically assess spatial specificity of this frequency-band, we performed a cluster-based permutation with surrogate data to identify in what scalp EEG electrodes the coherence with nucleus accumbens in the frequency range from 4 to 14 Hz was strongest [38]. The permutation test controls for multiple comparisons over electrodes and frequency bins. There was a highly significant difference between the observed and the surrogate data (p<0.01), with one cluster over the frontal cortex and one cluster over parietal cortex, which we subsequently defined as regions of interest (ROIs). Figs 2B and 1C show topographic representations of the coherence in the theta- and the alpha-band. Coherence in the theta-band was in general higher than in the alpha-band, and coherence in the frontal ROI was generally stronger than the parietal ROI (see also Fig 2D).

thumbnail
Fig 2. Coherence between the nucleus accumbens (NAc) electrodes and the scalp EEG.

A Z-transformed coherence spectrum between nucleus accumbens and frontal scalp electrode Fz. The horizontal dotted line (Z = 2.8) indicates the threshold of the coherence calculated from data shuffled over trials (see Methods). Coherence values between 4 Hz and 14 Hz are above this threshold. Electrodes Fz and Oz are circled. B Coherence in the theta-band showed a strong frontal topography and a weaker parieto-occipital topography. C Coherence in the alpha-band was strongest around frontal electrodes. The *’s indicate scalp electrodes that were significantly coherent between 4–14 Hz with nucleus accumbens when controlling for multiple comparisons over electrodes and frequency bins using a cluster randomization approach. A strong frontal and a weaker parietal cluster emerged. D Example of coherence in the parietal cluster. The Z-transformed coherence spectrum between nucleus accumbens and occipitoparietal scalp electrode Oz is depicted.

https://doi.org/10.1371/journal.pone.0138685.g002

Spectrotemporal dynamics of the oscillatory coupling

Subsequently, we explored the spectrotemporal dynamics of the coherence in response to the stimulus. To this end we conducted a time-resolved coherence analysis from -1000 ms to 500 ms relative to stimulus onset using a sliding window technique with Hanning-tapered windows of 500 ms length). In Fig 3A we show the coherence between electrodes in the frontal cluster and the nucleus accumbens electrodes (marked in Figs 2C and 1D). In the prestimulus interval, we observed high, sustained coherence in the theta- and alpha-bands. During stimulus processing, around 300 ms after stimulus onset, alpha-band coherence appeared weaker whereas coherence in the theta-band became stronger. This speaks to the functional relevance and different functional roles of theta- and alpha-band functional connectivity during stimulus processing.

thumbnail
Fig 3. Time-frequency dynamics of cortico-striatal interactions between nucleus accumbens (NAc) and frontal cortex (FC).

A Time-frequency coherence spectrum. In the prestimulus period (-1000 ms to 0 ms), both alpha and theta band coherence were strong. After stimulus onset (0 ms), alpha band coherence decreased while theta band coherence increased. B Granger causality from frontal electrodes to nucleus accumbens was highest in the alpha range (9–14 Hz) and diminished upon stimulus processing (after 0 ms). C Granger causality from nucleus accumbens to frontal electrodes was strongest in the theta range and increased upon stimulus processing (after 0 ms).

https://doi.org/10.1371/journal.pone.0138685.g003

In Fig 3A, alpha-band coherence appeared to become weaker whereas coherence in the theta-band became stronger around 300 ms after stimulus onset. When comparing the prestimulus interval (-500ms to 0ms) with the poststimulus interval (0ms to 500ms), these changes were, however, not significant, neither the theta-band increase (around Fz, t(5) = 1.31, p = 0.25), nor the alpha-band decrease (around Fz, t(5) = -1.22, p = 0.28). When posthoc pre-selecting the 300ms time point, however, we found that theta-band activity was significantly higher than during stimulus onset (comparing coherence from 4–7Hz at 0ms vs. 300ms; t(5) = 2.91, p = 0.03). No effect at this time point for the alpha-band (9–14Hz) was found (t(5) = -1.96, p = 0.11).

To further investigate the functional connectivity, we conducted a time-resolved Granger causality analysis with the same parameters that we used to compute the temporal coherence spectrum [3335]. Fig 3B shows the spectrotemporal dynamics of the Granger causality measure applied to frontal EEG and intracranial nucleus accumbens electrodes. In order to assess whether the resulting frequency-bands of the Granger causality analysis in the two directions were different from each other, we assessed the peak frequency of the Granger values in each subject in the range between 4 Hz and 14 Hz. We compared peak frequency from NAc to cortex (6.16 Hz +/- 0.6 Hz) to the peak frequency from cortex to NAc (8.9Hz +/- 1.7Hz) and found them to be significantly different from each other (t(5) = 3.42, p = 0.02). Consistent with contemporary literature and the coherence results, the lower of these two frequency peaks is in the theta-band between 4 Hz and 7 Hz, and the higher of the two peak frequencies is around the alpha-band, defined between 9 Hz and 14 Hz.

In Fig 3B, we further observe that in anticipation of the stimulus, sustained activity in the NAc predicted alpha band-activity of the cortex (from here on referred to as ‘granger-causal influence’ or simply ‘connectivity’). Upon stimulus processing, the alpha-band connectivity became weaker. In the reverse direction (Fig 2C), we found that theta-band activity from the NAc predicted activity around the frontal cortex in anticipation of the stimulus, which further increased during processing of the stimulus. Results for the parietal ROI are comparable, therefore not shown.

One might argue that modulation of coherence and Granger causality measures can be explained by concurrent changes in power. Therefore, we compared power in the theta and alpha bands between the prestimulus (-500 ms to 0 ms relative to stimulus onset) and poststimulus intervals (0 ms to 500 ms relative to stimulus onset). There was no statistically significant changes in power between these intervals in the theta-band, neither in the nucleus accumbens (t(5) = -1.12, p = 0.31), nor the frontal cluster (t(5) = 1.21, p = 0.28) and a near significant theta power decrease in the parietal cluster (t(5) = 2.27, p = 0.07). We conclude that the increase in coherence at the theta-band in response to the stimulus cannot be explained by a concurrent theta power increase in the cortex nor the nucleus accumbens. In the nucleus accumbens, there was no statistically significant change in alpha-power in response to the stimulus (t(5) = 1.00, p = 0.36), but there was a significant alpha-power decrease in the frontal ROI (t(5) = 3.92, p = 0.01) and in the parietal ROI t(5) = 3.1, p = 0.03). Thus, we cannot exclude that the decrease in alpha coherence in response to the stimulus is affected by the strong decrease in cortical alpha power.

To access the robustness of the modulations in Granger causality we conducted a repeated-measures ANOVA with four factors: region of interest (frontal and parietal EEG electrodes), time (the -500–0 ms prestimulus window and the 0–500 ms poststimulus window), frequency (the 4–7 Hz theta band and the 9–14 Hz alpha band) and directionality (nucleus accumbens to neocortex and neocortex to nucleus accumbens). Granger estimates were log-transformed to approach a normal distribution (Kolgomorov-Smirnov test; non-log transformed data: p<0.001; log-transformed data: p>0.5). In the following we report only significant, or near-significant results from the ANOVA (p<0.1); for all results see Table 1. The ANOVA showed a main effect of frequency (F(1,5) = 7.68, p<0.00393, Cohen’s f = 0.39), caused by Granger causality in the theta-band generally being stronger than in the alpha-band. There was a trend towards a main effect of time (F(1, 5) = 4.78, p = 0.08, Cohen’s f = 0.17).

In addition, we found a significant time by frequency interaction (F(1, 5) = 20.92, p = 0.006, Cohen’s f = 0.17, see Fig 4A). A posthoc analysis revealed that there was a significant theta-band Granger causality increase in the poststimulus interval compared to the prestimulus interval (t(5) = 3.52, p = 0.017). In contrast, Granger causality in the alpha band showed no significant change in time (t(5) = -0.08, p = 0.94). Furthermore, theta-causality in the post-stimulus interval was stronger than alpha-causality (t(5) = 4.27, p = 0.008). This suggests that stimulus processing is associated with increased connectivity in the theta-band between nucleus accumbens and neocortex.

thumbnail
Fig 4. Significant interaction effects of the factorial ANOVA on the Granger estimates.

A The time by frequency interaction (F(1, 5) = 20.92, p<0.01, Cohen’sf = 0.17) indicates that there was a significant theta band increase poststimulus compared to prestimulus Granger causality (t(5) = 3.52, p<0.05). No such change was present in the alpha-band. In the poststimulus window, Granger causality in the theta band was significantly stronger than in the alpha band (t(5) = 4.27, p<0.01). B The ANOVA showed a significant frequency by directionality interaction (F(1, 5) = 9.68, p<0.05, Cohen’s f = 0.17). In post-hoc tests, we found that the granger-causal influence from nucleus accumbens on cortex was stronger in the theta-band than in the alpha-band (t(5) = 4.84, p<0.01). In contrast, there was no significant difference in Granger causality in the alpha-band between cortex-to-nucleus accumbens and nucleus accumbens-to-cortex direction (t(5) = 1.75, p = 0.14). However, for five out of six subjects Granger estimates in the alpha-band were stronger for the cortex-to-nucleus accumbens direction than the other way around (see Fig 5). Data are represented as mean ± SEM. * = p<0.05, ** = p<0.01. = p<0.01.

https://doi.org/10.1371/journal.pone.0138685.g004

Crucially, we found a significant interaction between frequency and directionality (F(1, 5) = 9.68, p = 0.00265, Cohen’s f = 0.17; see Fig 4B). Posthoc t-tests revealed that the difference in theta-band Granger causality between nucleus accumbens to cortex and cortex to nucleus accumbens was larger than the same difference in the alpha-band. Granger causality from nucleus accumbens to cortex was significantly stronger in the theta-band than in the alpha-band (t(5) = 4.84, p = 0.005), indicating that theta-oscillations represent the primary granger-causal influence from nucleus accumbens to cortex.

In addition, we found that alpha-band activity from the cortex to NAc provided more predictive information than in the reverse direction from nucleus accumbens to the cortex in five out of six patients, which indicates that cortical alpha-band oscillations are granger-causing alpha-oscillations in the nucleus accumbens. Using conventional posthoc t-tests, this difference was, however, not significant (t(5) = 1.75, p = 0.14). The lack of statistical significance might be explained by the low number of subjects and by the low frequency resolution (2Hz) due to the 0.5s time windows, which results in spectral leakage between the theta- into the alpha-band. To further investigate the directional-specificity of the alpha-band, we analyzed the whole trial time-window from -1s to 0.5s using a permutation test. We performed a permutation test of Granger causality difference in the alpha-band across participants. For each participant, we computed the average Granger causality in the alpha-band in the whole time-window (-1 s to 0.5 s) between the nucleus accumbens and the cortex (i.e. the electrodes in two ROIs) and in the reverse direction (see Fig 5). Then, we subtracted the values for the two directions from one another. We created a permutation distribution by swapping the Granger estimates of the directionalities within subjects and performed the formerly mentioned procedure for all 64 (= 26) possible permutation. We found that the average directionality difference was more extreme for only two out of the 64 permutations than in the actual observation (one-sided test; p<0.05). The same procedure for the theta band did not reveal any significant directionality preference (one-sided test; p>0.05). We therefore conclude that connectivity in the alpha-band reflects primarily communication from the cortex to the nucleus accumbens.

thumbnail
Fig 5. Granger causality values for the theta- and the alpha-band in the two directions calculated for the full time interval (-1 to 0.5 s).

The solid line indicates the main diagonal to help illustrating in how many participants the Granger results were stronger in the one compared to the other direction A In half of the subjects, theta-band oscillations were stronger from cortex to nucleus accumbens (NAc) than the other way around. B Alpha-band oscillations were stronger from the cortex to the nucleus accumbens than from nucleus accumbens to the cortex for five out of six subjects (permutation test across subjects, p<0.05). This indicates that cortical alpha-band activity granger-cause activity in the nucleus accumbens.

https://doi.org/10.1371/journal.pone.0138685.g005

The nucleus accumbens mediates long-range cortical communication

Recently it has been argued that subcortical regions, e.g. the pulvinar, have an impact on large scale cortical communication [40,41]. Recent findings by van Schouwenburg, den Ouden & Cools [23,42] showed that the striatum has a modulatory influence on frontoposterior connectivity. Therefore, we investigated whether frontoposterior coherence in the theta and alpha-band is influenced by nucleus accumbens activity. We did this by considering the fronto-parietal coherence in relation to the contribution from the nucleus accumbens. Partial coherence helps explaining how the connectivity pattern between two sources is influenced by a third source [43,44]. In line with the findings by van Schouwenburg and colleagues, we hypothesized that if the nucleus accumbens is coordinating frontoposterior connectivity, then partializing out the contributions of the nucleus accumbens will reduce the fronto-parietal coherence. While the frontal and occipital cortex coherence was relatively strong in all frequency bands due to volume conduction, we still observed two peaks around the theta (4–7 Hz) and the alpha-band (9–14Hz; Fig 5, black line). When partializing out the contributions of the nucleus accumbens (Fig 6, red line), we found that coherence in the alpha-band decreased significantly (t(5) = 3.00, p = 0.03), whereas this was not the case for theta band coherence (t(5) = 1.84, p = 0.12). Consistently, five out of six subjects showed a decrease in alpha-band coherence, whereas four out of six showed a decrease in theta-band coherence. This finding corroborates the idea that the striatum facilitates long-range cortical communication between the frontal and posterior cortex.

thumbnail
Fig 6. Partial coherence results.

Ordinary coherence between the frontal and parietal cortex showed clear peaks in the theta and alpha frequency band. When partializing out the coherence from the nucleus accumbens, the frontoposterior coherence in the 9–14 Hz alpha-band (t(5) = 3.00, p = 0.03) decreased significantly.

https://doi.org/10.1371/journal.pone.0138685.g006

Discussion

We investigated how connectivity in corticostriatal networks is implemented by neuronal oscillations in the theta- (4–7 Hz) and alpha-band (9–14 Hz). For the first time, we showed bi-directional temporospectral dynamics between the nucleus accumbens and the neocortex in humans. First, we showed that the nucleus accumbens controlled the cortex primarily in the theta-band, and only a weak influence from the nucleus accumbens to the cortex in the alpha-band was present. Furthermore, we showed that this connectivity was already present in anticipation of visual stimulation, and increased during stimulus processing. We therefore suggest that theta-band oscillations serve the active exchange of stimulus information. Second, cortical to nucleus accumbens connectivity is implemented by alpha-band oscillations. Alpha-band connectivity in the neocortex predicted activity in the nucleus accumbens, and was, just like theta activity, present in anticipation of a visual stimulus, but did not increase upon stimulus processing. In line with current ideas on the functional role of cortical alpha-oscillations this might be an inhibitory con, serving to inhibit ongoing processing in the nucleus accumbens. We conclude that theta- and alpha-band oscillations serve different functional roles in cortico-striatal communication.

Oscillatory coupling between the nucleus accumbens and the frontal cortex has been previously found in humans during reward anticipation and reward processing. Cohen and colleagues showed coherence between nucleus accumbens and frontal cortex in the theta-band upon reward processing [13] and during reward anticipation [16]. Earlier, Cohen et al. found preliminary evidence for bidirectional communication between nucleus accumbens and frontal cortex during reward anticipation and reward processing, but did not investigate frequency-specificity [45]. Here we extend these findings by showing that connectivity between the nucleus accumbens and the cortex is implemented by oscillatory activity. This connectivity is present without expectation of reward or feedback, and is present both in anticipation of visual stimulation and during stimulus processing. Our findings are not likely to be explained by pickup of activity originating in the nucleus accumbens. In order to be measurable in the EEG an active neural population needs to be in an open field configuration [46]. The nucleus accumbens mostly consists of medium spiny neurons (e.g. [14,47]), which is a closed field configuration and it thus does not produce electric currents that can be measured as potential differences at the scalp level.

The strong influence from the nucleus accumbens to the neocortex might be surprising, given the polysynaptic connections from the striatum to the neocortex via the globus pallidus, subthalamic nucleus and the thalamus [1719,48,49] and the hippocampus [50,51]. In contrast, the striatum receives monosynaptic input from the prefrontal areas [17]. Therefore one might expect a stronger influence from the neocortex to the striatum than the other way around. We, however, found that the control from nucleus accumbens over the cortex was about equally strong as the control from the cortex over the nucleus accumbens. In rats, McCracken and Grace have found that high-frequency stimulation of nucleus accumbens resulted in stronger low-frequency activity (0.5–4 Hz) and beta-activity (13–30 Hz) in orbitofrontal cortex [52,53], which also speaks to a strong, causal influence from nucleus accumbens to neocortex. Moreover, it was recently shown in monkeys that the striatum (here specifically the caudate nucleus) exerts a stronger influence in the beta-band (12–30 Hz) on the frontal cortex than the frontal cortex does on the striatum [54]. Likewise in humans, Smolders and colleagues have shown that deep brain stimulation of electrodes dorsally adjacent to the nucleus accumbens resulted in higher phase stability in the theta range over frontal scalp electrodes in OCD patients compared to when stimulation was turned off [55]. Our findings match these observations and also suggest that the polysynaptic connections from nucleus accumbens to the neocortex allow for a strong, causal influence.

We found a significant frequency-by-directionality interaction. In subsequent posthoc tests we found that this interaction indicated, that connectivity between NAc and Cortex was stronger in the theta than the alpha band, and also that there was a stronger difference between theta and alpha band Granger causality for the NAc-to-cortex than for the cortex-to-NAc direction. This, however, does not mean that theta oscillations exclusively subserve bottom-up processsing (from NAc to cortex). We tested this question by a permutation test and tested whether theta-band Granger causality is stronger in the NAc-to-cortex than in the cortex-to-NAc direction. Theta connectivity did not differ between directions. Our results thus suggest that the NAc is primarily, but not uniquely, employing theta oscillations to communicate with the cortex, probably reflecting an exchange of stimulus information (cf. [13,14,16]) as well as stimulus evoked activity.

While the theta-band connectivity might reflect a bi-directional exchange of stimulus information, we propose that alpha-band connectivity represents inhibitory control originating in the frontal cortex to actively suppress striatal processing. Functionally, this would be well in line with the recent view on cortical alpha-band oscillations (cf. [5661]). Further substantiating this idea, it has been shown that electrical stimulation of the frontal cortex exerts strong inhibition onto striatal processing via hippocampal and thalamic pathways [62]. In order to verify these hypotheses on the differential functional roles of theta- and alpha-band oscillations, future studies could test how alpha- and theta-band connectivity are modulated in response to stimuli with different behavioral relevance. In our design, behaviorally irrelevant stimuli were simultaneously presented with relevant stimuli. Presentation of behavioral irrelevant and relevant stimuli in isolation will allow investigating the functional roles of theta- and alpha-band connectivity in more detail. We hypothesize that upon detection of a behavioral relevant stimulus, the striatum allowing for gating of the stimulus to prefrontal cortex, reflected by increased cortitical to nucleus accumbens theta-band phase synchrony. In contrast, presentation of behavioral irrelevant stimuli should elicit an increase in alpha-band connectivity between the cortex and the nucleus accumbens in order to prevent information about the stimulus to be exchanged.

It has recently been found that the nucleus accumbens is actively involved in gating visual information between neocortical regions. This was demonstrated in an fMRI study showing that the nucleus accumbens actively modulates frontoparietal connectivity during a visual attention task comparable to the one employed here [23,42]. A similar striatal modulation of corticocortical connectivity has been found in response to prediction errors [63]. Similarly we also found that the nucleus accumbens modulates frontoparietal connectivity as assessed by alpha band coherence. This is well in line with computational models, which predict that the striatum is involved in stimulus selection [2022]. Concluding, our findings show the nucleus accumbens, the frontal and the parietal cortex constitute a functional entity during visual attention tasks, which are employing different frequency-bands for directional communication. In accordance with recent literature, this network might be involved in stimulus selection and processing [21,23,42].

The deep brain electrodes were implanted for clinical reasons to treat OCD, depression and addiction symptoms. Although we analyzed the most ventrally implanted electrodes, extending into the nucleus accumbens (see Methods), clinical effect in OCD patients are achieved by stimulation more dorsally, in the ventral part of the anterior limb of the internal capsule (vALIC), and not directly in the nucleus accumbens [3]. Deep brain stimulation of the vALIC modulates white matter tracts connecting basal ganglia to prefrontal and orbitofrontal cortical areas, which are likely to influence the interaction of cortex and the nucleus accumbens.

We investigated a mixed patient population of OCD, depression and drug addiction, which raises the question how our findings translate to the healthy brain. For example, OCD patients have been found to exhibit hyper-connectivity between the nucleus accumbens and prefrontal cortex (e.g. [64]). However, our results on the mixed group of participants suggest that our effects are not restricted to OCD patients, but extend to patients suffering from other pathologies. Also, we found no pathologic-specific differences between the patient groups, thus the found mechanism might be universal and not pathological. In addition, the task that the patients performed was kept free from pathological symbols or images to avoid symptom provoking stimulations. Moreover behavioral performance did not differ between patients and healthy subjects, who participated in an earlier study with a modified, slightly harder version of the task [24]. We therefore think that it is likely that similar corticostriatal connectivity patterns are present in healthy individuals as well.

In conclusion, our results shed new light on how functional connectivity in cortico-striatal networks is implemented by neuronal oscillations. Theta- and alpha-band oscillations reflected connectivity in different directions, and served different functional purposes. In line with recent investigations [23,42], we propose that the striatum actively gates behavioral relevant information to the neocortex during visual attention, here mediated by theta-oscillations, while behaviorally irrelevant information becomes suppressed as reflected by strong alpha-band connectivity. In future studies, the behavioral relevance of the spectrotemporal dynamics of cortical to nucleus accumbens connectivity could investigate this further by presenting distracters as well as target stimuli, which allows to specifically study the dynamics associated with gating of behaviorally relevant information.

Author Contributions

Conceived and designed the experiments: JH MB RC DD OJ. Performed the experiments: RS. Analyzed the data: JH MB JS RC OJ. Wrote the paper: JH RS MB JS PvdM PRS RC DD OJ. Performed the neurosurgery: PvdM PRS.

References

  1. 1. Denys D, Mantione M, Figee M, van den Munckhof P, Koerselman F, Westenberg H, et al. Deep brain stimulation of the nucleus accumbens for treatment-refractory obsessive-compulsive disorder. Arch Gen Psychiatry. 2010;67: 1061–1068. pmid:20921122
  2. 2. Huff W, Lenartz D, Schormann M, Lee S-H, Kuhn J, Koulousakis A, et al. Unilateral deep brain stimulation of the nucleus accumbens in patients with treatment-resistant obsessive-compulsive disorder: Outcomes after one year. Clin Neurol Neurosurg. 2010;112: 137–143. pmid:20006424
  3. 3. van den Munckhof P, Bosch DA, Mantione MHM, Figee M, Denys DAJP, Schuurman PR. Active stimulation site of nucleus accumbens deep brain stimulation in obsessive-compulsive disorder is localized in the ventral internal capsule. Acta Neurochir Suppl. 2013;117: 53–59. pmid:23652657
  4. 4. Bewernick BH, Kayser S, Sturm V, Schlaepfer TE. Long-term effects of nucleus accumbens deep brain stimulation in treatment-resistant depression: evidence for sustained efficacy. Neuropsychopharmacol Off Publ Am Coll Neuropsychopharmacol. 2012;37: 1975–1985.
  5. 5. Blomstedt P, Sjöberg RL, Hansson M, Bodlund O, Hariz MI. Deep brain stimulation in the treatment of depression. Acta Psychiatr Scand. 2011;123: 4–11. pmid:20961293
  6. 6. Grubert C, Hurlemann R, Bewernick BH, Kayser S, Hadrysiewicz B, Axmacher N, et al. Neuropsychological safety of nucleus accumbens deep brain stimulation for major depression: effects of 12-month stimulation. World J Biol Psychiatry Off J World Fed Soc Biol Psychiatry. 2011;12: 516–527.
  7. 7. Nauczyciel C, Robic S, Dondaine T, Verin M, Robert G, Drapier D, et al. The nucleus accumbens: a target for deep brain stimulation in resistant major depressive disorder. J Mol Psychiatry. 2013;1: 17. pmid:26019865
  8. 8. Schlaepfer TE, Bewernick BH. Deep brain stimulation for major depression. Handb Clin Neurol. 2013;116: 235–243. pmid:24112897
  9. 9. Giacobbe P, Mayberg HS, Lozano AM. Treatment resistant depression as a failure of brain homeostatic mechanisms: implications for deep brain stimulation. Exp Neurol. 2009;219: 44–52. pmid:19426730
  10. 10. Luigjes J, van den Brink W, Feenstra M, van den Munckhof P, Schuurman PR, Schippers R, et al. Deep brain stimulation in addiction: a review of potential brain targets. Mol Psychiatry. 2012;17: 572–583. pmid:21931318
  11. 11. Pierce RC, Vassoler FM. Deep brain stimulation for the treatment of addiction: basic and clinical studies and potential mechanisms of action. Psychopharmacology (Berl). 2013;229: 487–491.
  12. 12. Cohen MX, Axmacher N, Lenartz D, Elger CE, Sturm V, Schlaepfer TE. Good Vibrations: Cross-frequency Coupling in the Human Nucleus Accumbens during Reward Processing. J Cogn Neurosci. 2008;21: 875–889.
  13. 13. Cohen MX, Axmacher N, Lenartz D, Elger CE, Sturm V, Schlaepfer TE. Nuclei Accumbens Phase Synchrony Predicts Decision-Making Reversals Following Negative Feedback. J Neurosci. 2009;29: 7591–7598. pmid:19515927
  14. 14. Axmacher N, Cohen MX, Fell J, Haupt S, Dümpelmann M, Elger CE, et al. Intracranial EEG correlates of expectancy and memory formation in the human hippocampus and nucleus accumbens. Neuron. 2010;65: 541–549. pmid:20188658
  15. 15. Lega BC, Kahana MJ, Jaggi J, Baltuch GH, Zaghloul K. Neuronal and oscillatory activity during reward processing in the human ventral striatum. Neuroreport. 2011;22: 795–800. pmid:21975313
  16. 16. Cohen MX, Bour L, Mantione M, Figee M, Vink M, Tijssen MAJ, et al. Top-down-directed synchrony from medial frontal cortex to nucleus accumbens during reward anticipation. Hum Brain Mapp. 2012;33: 246–252. pmid:21547982
  17. 17. Alexander GE, DeLong MR, Strick PL. Parallel Organization of Functionally Segregated Circuits Linking Basal Ganglia and Cortex. Annu Rev Neurosci. 1986;9: 357–381. pmid:3085570
  18. 18. Haber SN. The primate basal ganglia: parallel and integrative networks. J Chem Neuroanat. 2003;26: 317–330. pmid:14729134
  19. 19. Middleton FA, Strick PL. Basal Ganglia Output and Cognition: Evidence from Anatomical, Behavioral, and Clinical Studies. Brain Cogn. 2000;42: 183–200. pmid:10744919
  20. 20. Frank MJ, Loughry B, O’Reilly RC. Interactions between frontal cortex and basal ganglia in working memory: a computational model. Cogn Affect Behav Neurosci. 2001;1: 137–160. pmid:12467110
  21. 21. Hazy TE, Frank MJ, O’Reilly RC. Towards an executive without a homunculus: computational models of the prefrontal cortex/basal ganglia system. Philos Trans R Soc B Biol Sci. 2007;362: 1601–1613.
  22. 22. O’Reilly RC, Frank MJ. Making working memory work: a computational model of learning in the prefrontal cortex and basal ganglia. Neural Comput. 2006;18: 283–328. pmid:16378516
  23. 23. van Schouwenburg MR, Ouden HEM den, Cools R. The human basal ganglia modulate frontal-posterior connectivity during attention shifting. J Neurosci Off J Soc Neurosci. 2010;30: 9910–9918.
  24. 24. Horschig JM, Jensen O, van Schouwenburg MR, Cools R, Bonnefond M. Alpha activity reflects individual abilities to adapt to the environment. NeuroImage. 2014;89: 235–243. pmid:24361665
  25. 25. Rubner Y, Tomasi C, Guibas LJ. A metric for distributions with applications to image databases. Sixth International Conference on Computer Vision, 1998. 1998. pp. 59–66.
  26. 26. Ruzon M., Tomasi C. Color edge detection with the compass operator. Computer Vision and Pattern Recognition, 1999 IEEE Computer Society Conference on. 1999. p. -166 Vol. 2.
  27. 27. Perrin F, Pernier J, Bertrand O, Echallier JF. Spherical splines for scalp potential and current density mapping. Electroencephalogr Clin Neurophysiol. 1989;72: 184–187. pmid:2464490
  28. 28. Oostenveld R, Fries P, Maris E, Schoffelen J-M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci. 2011;2011: 156869. pmid:21253357
  29. 29. Rosenberg JR, Amjad AM, Breeze P, Brillinger DR, Halliday DM. The Fourier approach to the identification of functional coupling between neuronal spike trains. Prog Biophys Mol Biol. 1989;53: 1–31. pmid:2682781
  30. 30. Granger CWJ. Investigating Causal Relations by Econometric Models and Cross-spectral Methods. Econometrica. 1969;37: 424–438.
  31. 31. Wiener N. The theory of prediction. Modern mathematics for the engineer. New York: McGraw-Hill; 1956.
  32. 32. Geweke JF. Measures of Conditional Linear Dependence and Feedback Between Time Series. J Am Stat Assoc. 1984;79: 907–915.
  33. 33. Ding M, Chen Y, Bressler SL. Granger Causality: Basic Theory and Application to Neuroscience. In: Schelter B, Winterhalder tthias, Timmer J, editors. Handbook of Time Series Analysis. Wiley-VCH Verlag GmbH & Co. KGaA; 2006. pp. 437–460. Available: http://onlinelibrary.wiley.com/doi/10.1002/9783527609970.ch17/summary
  34. 34. Dhamala M, Rangarajan G, Ding M. Analyzing information flow in brain networks with nonparametric Granger causality. NeuroImage. 2008;41: 354–362. pmid:18394927
  35. 35. Dhamala M, Rangarajan G, Ding M. Estimating Granger causality from fourier and wavelet transforms of time series data. Phys Rev Lett. 2008;100: 018701. pmid:18232831
  36. 36. Wilson GT. The Factorization of Matricial Spectral Densities. SIAM J Appl Math. 1972;23: 420–426.
  37. 37. Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods. 2007;164: 177–190. pmid:17517438
  38. 38. Maris E, Schoffelen J-M, Fries P. Nonparametric statistical testing of coherence differences. J Neurosci Methods. 2007;163: 161–175. pmid:17395267
  39. 39. Steiger JH. Beyond the F test: Effect size confidence intervals and tests of close fit in the analysis of variance and contrast analysis. Psychol Methods. 2004;9: 164–182. pmid:15137887
  40. 40. Saalmann YB, Pinsk MA, Wang L, Li X, Kastner S. The pulvinar regulates information transmission between cortical areas based on attention demands. Science. 2012;337: 753–756. pmid:22879517
  41. 41. Ketz NA, Jensen O, O’Reilly RC. Thalamic pathways underlying prefrontal cortex–medial temporal lobe oscillatory interactions. Trends Neurosci. 0.
  42. 42. van Schouwenburg MR, Ouden HEM den, Cools R. Selective Attentional Enhancement and Inhibition of Fronto-Posterior Connectivity by the Basal Ganglia During Attention Switching. Cereb Cortex N Y N 1991. 2013;
  43. 43. Kocsis B, Bragin A, Buzsáki G. Interdependence of Multiple Theta Generators in the Hippocampus: a Partial Coherence Analysis. J Neurosci. 1999;19: 6200–6212. pmid:10407056
  44. 44. Lopes da Silva FH, Vos JE, Mooibroek J, Van Rotterdam A. Relative contributions of intracortical and thalamo-cortical processes in the generation of alpha rhythms, revealed by partial coherence analysis. Electroencephalogr Clin Neurophysiol. 1980;50: 449–456. pmid:6160987
  45. 45. Cohen MX, Axmacher N, Lenartz D, Elger CE, Sturm V, Schlaepfer TE. Neuroelectric Signatures of Reward Learning and Decision-Making in the Human Nucleus Accumbens. Neuropsychopharmacology. 2008;34: 1649–1658. pmid:19092783
  46. 46. Nunez PL. Toward a quantitative description of large-scale neocortical dynamic function and EEG. Behav Brain Sci. 2000;23: 371–398; discussion 399–437. pmid:11301576
  47. 47. Tepper JM, Bolam JP. Functional diversity and specificity of neostriatal interneurons. Curr Opin Neurobiol. 2004;14: 685–692. pmid:15582369
  48. 48. Haber SN, Kunishio K, Mizobuchi M, Lynd-Balta E. The orbital and medial prefrontal circuit through the primate basal ganglia. J Neurosci Off J Soc Neurosci. 1995;15: 4851–4867.
  49. 49. Haber SN, Calzavara R. The cortico-basal ganglia integrative network: The role of the thalamus. Brain Res Bull. 2009;78: 69–74. pmid:18950692
  50. 50. Tort ABL, Kramer MA, Thorn C, Gibson DJ, Kubota Y, Graybiel AM, et al. Dynamic cross-frequency couplings of local field potential oscillations in rat striatum and hippocampus during performance of a T-maze task. Proc Natl Acad Sci. 2008;105: 20517–20522. pmid:19074268
  51. 51. Battaglia FP, Benchenane K, Sirota A, Pennartz CMA, Wiener SI. The hippocampus: hub of brain network communication for memory. Trends Cogn Sci. 2011;15: 310–318. pmid:21696996
  52. 52. McCracken CB, Grace AA. Nucleus Accumbens Deep Brain Stimulation Produces Region-Specific Alterations in Local Field Potential Oscillations and Evoked Responses In Vivo. J Neurosci. 2009;29: 5354–5363. pmid:19386932
  53. 53. McCracken CB, Grace AA. High-Frequency Deep Brain Stimulation of the Nucleus Accumbens Region Suppresses Neuronal Activity and Selectively Modulates Afferent Drive in Rat Orbitofrontal Cortex In Vivo. J Neurosci. 2007;27: 12601–12610. pmid:18003839
  54. 54. Antzoulatos EG, Miller EK. Increases in Functional Connectivity between Prefrontal Cortex and Striatum during Category Learning. Neuron. 2014;83: 216–225. pmid:24930701
  55. 55. Smolders R, Mazaheri A, van Wingen G, Figee M, de Koning PP, Denys D. Deep Brain Stimulation Targeted at the Nucleus Accumbens Decreases the Potential for Pathologic Network Communication. Biol Psychiatry. 2013;74: e27–e28. pmid:23623188
  56. 56. Foxe JJ, Snyder AC. The Role of Alpha-Band Brain Oscillations as a Sensory Suppression Mechanism during Selective Attention. Front Psychol. 2011;2: 154. pmid:21779269
  57. 57. Jensen O, Bonnefond M, VanRullen R. An oscillatory mechanism for prioritizing salient unattended stimuli. Trends Cogn Sci. 2012;16: 200–206. pmid:22436764
  58. 58. Jensen O, Mazaheri A. Shaping functional architecture by oscillatory alpha activity: gating by inhibition. Front Hum Neurosci. 2010;4: 186. pmid:21119777
  59. 59. Klimesch W. EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis. Brain Res Brain Res Rev. 1999;29: 169–195. pmid:10209231
  60. 60. Klimesch W, Sauseng P, Hanslmayr S. EEG alpha oscillations: the inhibition-timing hypothesis. Brain Res Rev. 2007;53: 63–88. pmid:16887192
  61. 61. Klimesch W. α-band oscillations, attention, and controlled access to stored information. Trends Cogn Sci. 2012;16: 606–617. pmid:23141428
  62. 62. Calhoon GG, O’Donnell P. Closing the Gate in the Limbic Striatum: Prefrontal Suppression of Hippocampal and Thalamic Inputs. Neuron. 2013;78: 181–190. pmid:23583113
  63. 63. Ouden HEM den, Daunizeau J, Roiser J, Friston KJ, Stephan KE. Striatal Prediction Error Modulates Cortical Coupling. J Neurosci. 2010;30: 3210–3219. pmid:20203180
  64. 64. Figee M, Luigjes J, Smolders R, Valencia-Alfonso C-E, van Wingen G, de Kwaasteniet B, et al. Deep brain stimulation restores frontostriatal network activity in obsessive-compulsive disorder. Nat Neurosci. 2013;16: 386–387. pmid:23434914