Abstract
Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) have been used to study the neural correlates of reward anticipation, but the interrelation of EEG and fMRI measures remains unknown. The goal of the present study was to investigate this relationship in response to a well established reward anticipation paradigm using simultaneous EEG-fMRI recording in healthy human subjects. Analysis of causal interactions between the thalamus (THAL), ventral-striatum (VS), and supplementary motor area (SMA), using both mediator analysis and dynamic causal modeling, revealed that (1) THAL fMRI blood oxygenation level-dependent (BOLD) activity is mediating intermodal correlations between the EEG contingent negative variation (CNV) signal and the fMRI BOLD signal in SMA and VS, (2) the underlying causal connectivity network consists of top-down regulation from SMA to VS and SMA to THAL along with an excitatory information flow through a THAL→VS→SMA route during reward anticipation, and (3) the EEG CNV signal is best predicted by a combination of THAL fMRI BOLD response and strength of top-down regulation from SMA to VS and SMA to THAL. Collectively, these findings represent a likely neurobiological mechanism mapping a primarily subcortical process, i.e., reward anticipation, onto a cortical signature.
Introduction
Functional magnetic resonance imaging (fMRI) studies have identified a neural circuit underlying reward anticipation in humans and nonhuman species. This circuit includes the ventral striatum (VS) as a pivotal structure that is innervated by mesostriatal dopamine and linked to the anterior cingulate cortex (ACC), supplementary motor area (SMA), insula, and thalamus (THAL; Knutson et al., 2001a,b; Kirsch et al., 2003; Knutson and Cooper, 2005; Liu et al., 2011). Anticipatory VS activation has been shown to be a reliable fMRI assay within and across subjects (Plichta et al., 2012), sensitive to genetic variants (Kirsch et al., 2006; Forbes et al., 2009; Hahn et al., 2011) and linked to impulsivity and (an)hedonic states in clinical and nonclinical populations (Kirsch et al., 2006; Scheres et al., 2007; Hahn et al., 2009; Plichta et al., 2009; Forbes et al., 2010; Juckel et al., 2012).
Despite considerable fMRI research on the neural correlates of reward anticipation, there are only relatively few studies focusing on the time-resolved neurophysiological signatures of this process. Previous studies have identified several event-related potentials (ERPs) derived from the electroencephalogram (EEG) relevant to reward processing (Broyd et al., 2012a) among those the contingent negative variation (CNV), a slow frontocentrally negative electrical brain activity, which is directly related to expectation and preparation of motivated responses (Walter et al., 1964).
The CNV amplitude has been shown to increase with reward in animal and human studies (Boyd et al., 1979; Pierson et al., 1987; Kotani et al., 2001; Broyd et al., 2012b; Capa et al., 2013) and source analysis (Gómez et al., 2003) suggested important generators of the CNV in the ACC, SMA, and posterior regions including the inferior parietal cortex (IPC). This has been validated in a simultaneous EEG-fMRI study (Nagai et al., 2004).
In the context of reward processing, ERP and fMRI reward anticipation markers have been studied only separately. Therefore, the interrelation between them is still unknown. The primary goal of the present study was to investigate the relationship between EEG and fMRI signals in response to a well established fMRI reward anticipation paradigm with simultaneous EEG-fMRI data recording.
By using conventional general linear model (GLM) analysis and analysis of causal interactions between the THAL, VS, and SMA based on mediator analysis and dynamic causal modeling (DCM) we show that (1) THAL fMRI BOLD activity is mediating intermodal correlations between the EEG CNV signal and the fMRI BOLD signal in SMA and VS, (2) the underlying causal connectivity network consists of top-down regulation from SMA to VS and SMA to THAL along with an excitatory information flow through a THAL→VS→SMA route during reward anticipation, and (3) the EEG CNV signal is best predicted by a combination of THAL fMRI BOLD response and strength of top-down regulation from SMA to VS and SMA to THAL.
Materials and Methods
Subjects.
Eighteen healthy right-handed subjects participated in the experiment. For technical reasons, data from four subjects had to be discarded: one subject had an insufficient number of artifact-free EEG trials; two subjects showed significant outlier mean CNV values; and one subject produced significant outlier fMRI results (for outlier analyses, see below). The final sample consisted of 14 subjects (8 males and 6 females) with a mean age of 24 years (SD = 4.24).
All participants were informed about the nature of the experiment before providing written informed consent. The combined EEG-fMRI investigation of healthy participants and the entire experimental procedure was in accordance with the Declaration of Helsinki and was approved by the local ethics committee of the Medical Faculty Mannheim of Heidelberg University.
Paradigm.
We used a modified version of a monetary reward anticipation task (Kirsch et al., 2003) that has been shown to reliably and robustly activate the VS (Plichta et al., 2012). In this paradigm, the subject must respond sufficiently fast to a light-flash (duration = 100 ms) on the visual display screen. The upcoming flash was signaled by a stimulus S1 that informed the subject about the consequences of their response to the flash stimulus S2. Stimulus onset asynchrony between S1 and S2 was variable (3–5 s) with a mean of 4 s. After the response of the subject to the flash and an additional time lag of 200 ms, the outcome was reported on the screen for 1.5 s. After an intertrial interval of 2 s plus an reaction time (RT)-dependent value that ensured equal trial durations [i.e., maximum of RT window (1 s) minus actual RT], the next trial started.
The two conditions were: (1) monetary reward anticipation (S1: smiley) where the subject could win 0.50 Euros if the response was sufficiently fast, or (2) control condition (S1: scrambled smiley) where only a verbal feedback was given. During the outcome phase the participants saw the outcome of the latest trial (“+0.50€”; “+0.00€”, or “Fast response!”; “Too slow”) as well as their cumulative sum total at that point in time. In this event-related paradigm, each of the above conditions was presented 50 times in a pseudo-randomized order. Additionally, we used so-called boost trials to maximize reward processing through uncertainty and to ensure task engagement (Riba et al., 2008), i.e., after successful performance the subject won unexpectedly 2 Euros instead of 0.50 Euros. Boost trials occurred after eight wins on average. The RT window for a sufficiently fast response was adaptively tailored to the individual RTs to collect a comparable number of winning trials across subjects (∼60%). After the final trial presentation, a black screen was presented for 12 s. The total run length was 14 min 52 s.
Data acquisition.
The EEG was recorded inside the scanner using a 5 kHz sampling rate, 32 mV input range and 0.1–250 Hz bandpass filters. The EEG data were transmitted from two MR-compatible BrainAmp amplifiers (Brain Products) via optic fibers outside the scanner room. Sixty equidistantly spaced silver/silver chloride (Ag/AgCl) scalp electrodes recorded the EEG signal using MR-compatible caps (Easycap) with twisted and fixed electrode cables. EEG montage was based on a selection of 10–20 system positions. To account for different head circumference, three different cap sizes were available. F1 served as recording reference, F2 was the ground electrode. An additional electrooculogram (EOG) electrode was placed below the outer canthus of each eye and two further electrodes (one right to the sternum and one on the left chest close to the heart) to record the electrocardiogram (ECG). Electrode impedances were kept <20 kΩ, except ECG, and EOG electrodes (<30 kΩ) as well as Reference and Ground (<10 kΩ). The EEG was monitored while scanning using online correction software (RecView, Brain Products).
EEG analysis.
EEG data were corrected for the MRI gradient (Allen et al., 2000) and cardioballistic artifacts (Allen et al., 1998) using standard template subtraction procedures as implemented in the Brain Vision Analyzer software 2.0 (Brain Products). EEG data were digitally low-pass filtered (70 Hz) and down sampled to 500 Hz. After exclusion of physical artifacts via raw data inspection, infomax independent component analysis (ICA) was used to remove ocular artifacts (blinks, movements), remaining cardioballistic artifacts (Debener et al., 2008) and residual artifacts related to gradient modulation (Brain Vision Analyzer software 2.0, Brain Products). EEG data were re-referenced to an average reference, baseline corrected by subtracting the average over a 500 ms prestimulus interval and low-pass filtered with a cutoff of 30 Hz. Segmentation into ERP epochs of 4.5 s began 500 ms before stimulus onset and ended 4000 ms poststimulus onset. ERP averages for both conditions (monetary reward anticipation and control) were calculated within each and across all participants. The CNV was captured at electrode Central zero (Cz) within two different time windows (CNV1: from 1 to 2 s and CNV2: from 2 to 3 s poststimulus presentation). Dividing the CNV into successive 1 s time window preceding the (earliest) target onset was based on previous work (van Leeuwen et al., 1998; Nagai et al., 2004) because these CNV phases may reflect partly different neural processes (Rohrbaugh et al., 1976). The CNVs were exported as mean amplitudes for each time window.
fMRI image acquisition.
All MRI sequences were performed on a 3-tesla whole-body scanner (Magnetom Trio, Siemens Medical Solutions). Functional data were acquired using T2*-weighted single shot gradient echo planar imaging (EPI) sequences with the following scanning-parameters: TR/TE = 2210/28 ms; flip angle = 90°; 36 axial slices (slice thickness = 3 mm + 1 mm gap) ascending, FOV = 220 mm × 220 mm, matrix = 64 × 64. Slices were oriented ∼20° steeper than AC-PC-plane. The first four EPI volumes were discarded to allow for T1 equilibration. After recording of the functional images, a high-resolution T1-weighted 3D MRI sequence was conducted (192 sagittal slices, slice thickness 1 mm, FOV = 256 mm × 256 mm, matrix = 256 × 256). These parameters provided whole brain coverage, including cerebellum, scalp, eyes, and nose, to avoid wrap-around artifacts. The paradigm was presented to the subject via video goggles (Resonance Technologies) controlled by the software Presentation (Neurobehavioral Systems). Behavioral data were recorded with response buttons (Current Designs).
fMRI data analyses.
The fMRI data were analyzed using statistical parametric mapping (SPM8; Wellcome Trust Centre for Neuroimaging, London, UK). Preprocessing of the fMRI data included slice-time correction, motion correction, spatial normalization into Montreal Neurological Institute [MNI] space, resampling to 2 × 2 × 2 mm3, and spatial smoothing with an 8 mm full-width at half-maximum Gaussian kernel. Spatial normalization was performed by calculating linear (12-parameter affine) and nonlinear transformations of the mean EPI image from each time series to the SPM EPI template in MNI space, and then applying these same transformation parameters to the time series.
Statistical analyses comprised first level temporal modeling within a GLM framework to generate a 3D map corresponding to estimated regressor response amplitudes. Regressors of interest (monetary reward anticipation, control), both modeled as events with a duration of zero seconds, were convolved with the default SPM hemodynamic response function (HRF) computed as a two-parameter gamma function. Eight regressors of no-interest were included: two regressors coding the feedback and the responses (convolved with the HRF) and six motion parameters (not convolved with the HRF). For all three tasks, a high-pass filter with a cutoff frequency of 1/128 Hz was used to attenuate low-frequency components. All analyses were corrected for serially correlated errors by fitting a first-order autoregressive process (AR[1]) to the error term.
For the combined analysis of EEG and fMRI data we used two different approaches. First, we calculated the mean CNV parameters per subject (average amplitude over predefined CNV1 and CNV2 windows; see above) and used them as covariates in the second-level analyses, revealing the correlation between the ERP difference amplitudes and fMRI brain activation during reward anticipation (as contrasted to the control condition). Second, we applied a trial-to-trial analysis: here we used the single-trial ERP amplitudes of each participant as a parametric modulator of the condition-dependent BOLD responses at the single-subject level. Contrast images from these first-level models were then introduced into a one-sample t test at the second level to test for group effects.
ROI definitions.
For the combined analyses of EEG and fMRI data, we defined the following regions-of-interest (ROIs). First, we defined the ventral striatum as a fusion of the “caudate head” mask taken from the WFU-PickAtlas (human-atlas TD Brodmann areas+) and the “accumbens” mask from the Harvard–Oxford Subcortical Structural Atlas (implemented in FSLView 3.1.8; http://www.cma.mgh.harvard.edu/fsl_atlas.html; probability threshold was set to 50%). The left and right VS were treated as separate ROIs. Second, we used an anatomical definition of the thalamus taken from the WFU-PickAtlas (human-atlas aal). Finally, for the SMA we used an empirical mask based on binarized second-level activation maps calculated from a published independent subject sample (Plichta et al., 2012) that was scanned twice (N = 25) using the original reward paradigm and the identical contrast (monetary reward anticipation > control). This was done because the atlas definition of the SMA is much larger than typically observed activation foci, hence probably heterogeneous in function and likely to result in low ROI mean response values due to a mixing of strongly and weakly responding voxels.
Mediator analysis.
A mediator analysis was conducted to assess whether the neural regions that show significant correlation with the ERP explain similar variance or whether certain regions subsume the impact of others. Specifically, THAL is assumed to be involved relatively early in the reward anticipation process, whereas the SMA, as a motor execution relay, might represent a relatively late component in information flow. We formally tested mediation by means of the Sobel test (Baron and Kenny, 1986) accompanied by a bootstrapping method with N = 5000 bootstrap samples (Preacher and Hayes, 2004) using the PROCESS procedure (Hayes, 2012) for SPSS version 20 (IBM).
Dynamic causal modeling.
To infer the causal structure of the network between the ROIs, we additionally applied bilinear DCM of the network between THAL, VS, and SMA. DCM regards the brain as a nonlinear dynamic system of interconnected neural nodes, and models experimental conditions as perturbations of the system's dynamics (Friston et al., 2003; Stephan et al., 2007, 2010). After fitting a set of a priori defined models (“model space”) to the data, models are compared and the winning model with the highest log evidence relative to the others is assumed to describe the network structure with the highest accuracy and the lowest complexity (Stephan et al., 2009, 2010).
We extracted the first eigenvariate from 8 mm spheres around the individual peak voxels of the GLM contrast “monetary reward anticipation > control” within the functional masks of THAL, VS, and SMA on the left hemisphere and constructed a model space covering all possible causal models. In all models, the monetary reward anticipation and the control condition were fixed as driving inputs into the THAL, resulting in a model space of 32 models. In addition, we assumed that both conditions activate the same network but might lead to differences in connectivity strength. Thus, the monetary reward anticipation condition was allowed to modulate the strength of all existing connections in each model to account for changes in connectivity between conditions.
After estimation, all models were entered in a group analysis with random effects Bayesian model selection (RFX BMS) and the parameters of the winning model were extracted from the individual subjects (Stephan et al., 2009).
Due to the lateralization of the CNV data which is consistent with Sobotka et al. (1992) and the right-handed task nature we focused our DCM analyses on the left hemisphere.
Data quality check.
EEG and fMRI data were checked for significant outliers separately. Two strategies were used: (1) Boxplot analysis with a criterion for data exclusion of 1.5 times the interquartile range, and (2) extreme studentized deviate test with a criterion of α = 0.05. Both approaches were performed in repeated cycles of (1) outlier detection, (2) outlier exclusion, and (3) reanalysis until no more outliers were apparent. Both approaches identified the identical set of three subjects (two EEG dataset and one fMRI dataset outlier).
Results
Behavioral effects
Participants reacted significantly faster after the presentation of a monetary stimulus (monetary reward anticipation: 167.35 s; control: 196.04 s; t = −4.99; df = 13; p < 0.001) and reached the reaction time window significantly more often during monetary trials compared with the control condition (monetary reward anticipation: 30.21; control condition: 21.07 (absolute counts); t = 4.52; df = 13; p = 0.001). All participants gained money with a mean payoff of 22.78 € (SD, 2.54; range, 18–26.50 €).
fMRI main effects
Figure 1A shows the group-level results for the contrast monetary reward anticipation > control. The expected reward anticipation network including VS, THAL, and SMA was clearly detectable (Table 1 shows fMRI statistics).
A, Significant group-level fMRI results for the contrast “monetary reward anticipation>control”. The typical BOLD activation pattern of reward anticipation was detected including VS, thalamus, SMA, and other structures. B, Grand average ERPs showing CNV waves (and standard error) at electrode position Cz developing after cue presentation signaling either monetary reward anticipation or control. Significant differences between the conditions occur both in CNV1 (time window 1–2 s) and CNV2 (2–3 s). Scalp topographies for monetary reward anticipation and control condition (time window 1–3 s) are shown in C. Electrode position Cz is marked by an asterisk.
fMRI statistics; main effects for the contrast (monetary reward anticipation > control)
EEG main effects
Figure 1B,C shows the grand-average waveforms at the Cz electrode site of the CNV for monetary reward anticipation and control conditions along with the EEG channel scheme and scalp topography of the CNV. Highly significant condition effects (monetary reward anticipation > control) were found for both CNV phases (Table 2 shows EEG statistics).
EEG statistics; main effects
Joint EEG-fMRI analysis
Analyzing the correlation of the CNV measures with the BOLD response amplitudes across subjects revealed a network comprising the cerebellum, THAL, mid-brain structures, insula, SMA, sensorimotor areas, and parietal cortices. Hypothesis-driven analyses of the correlation between the CNV and BOLD responses in the predefined ROIs revealed significant correlations in VS, THAL, and SMA (Fig. 2A–C).
Significant correlation of the CNV amplitude (difference score: monetary reward anticipation > control) and BOLD fMRI signal (monetary reward anticipation > control) within the predefined ROIs (A, SMA; B, VS; C, THAL). The scatterplots show ROI average BOLD signal of all voxels within VS, SMA, and THA ROI. Note that results for CNV1 are shown and that the effects for CNV2 are consistently weaker (data not shown).
Trial-to-trial analyses (see Materials and Methods) revealed a similar correlation network including the identical structures as listed above. However, the effects were relatively small compared with the second-level correlations and did not survive correction for multiple testing.
Intermodal correlation effects were consistently stronger for CNV1 compared with CNV2. Therefore, further analyses focused on CNV1.
Mediation analyses showed that the relationship between both ROIs (VS and SMA) and the CNV1 measure was fully mediated by thalamic activation (Fig. 3; Table 3). Importantly, the converse was not true, i.e., the relation between thalamic activation and CNV1 did not significantly change after the inclusion of VS or SMA as mediator variables.
Mediation analyses. Relation between the EEG CNV and the fMRI response (A) in the VS and (B) in the SMA is fully mediated by the THAL, i.e., the significance of bivariate relations (gray numbers) fully vanishes (red numbers) when THAL thalamus is considered as a mediator. Note that the converse is not true, i.e., the relation between the thalamic fMRI response and the EEG CNV is not significantly reduced when either the VS (C) or the SMA (D) is included as a mediator. Numbers are correlation coefficients; *p < 0.05; **p < 0.01.
Mediation analysis
DCM analyses
RFX BMS revealed one clear winner model describing the underlying neuronal activity best: the fully connected model had a model exceedance probability of 0.71 in a model space of 32 models where all probabilities sum up to 1 (Fig. 4A). This supports the hypothesis that all three nodes in the network are reciprocally connected and are influencing each other.
Results of the DCM analyses. The exceedance probabilities resulting from the BMS for the 32 models are shown in A. The fully connected model is the clear winner (p = 0.71). B, The endogenous connectivity, i.e., connectivity regardless of task demands. C, Modulation of connectivity due to monetary reward anticipation.
Extracting the parameters of the winning model revealed a strong intrinsic excitatory influence of the THAL on both the VS and SMA. In addition, VS showed an inhibitory influence on the SMA (Fig. 4B). To further elucidate whether a physiologically plausible surrogate of the revealed full mediation by THAL exists, the modulation of the network by the monetary reward anticipation condition was analyzed. The monetary reward anticipation condition led to a decrease in the drive from THAL to SMA and to an increase in its drive to the VS, i.e., the influence of the VS on the SMA became excitatory during the monetary reward anticipation condition (Fig. 4C).
DCM-CNV relationship
To explore whether DCM parameters explain additional variance of the CNV, stepwise multiple regression (inclusion criterion p < 0.10) was applied with the BOLD contrast estimates of THAL, VS, and SMA together with the six DCM parameters for the modulated network as predictors for the CNV. This procedure identified a model consisting of THAL activation and two DCM parameters: SMA→VS and SMA→THAL connectivity modulation (model statistics: F(3,10) = 11.6; p < 0.001; adjusted R2=.71). This model which combines BOLD activation and DCM connectivity parameters explained significantly more variance than THAL BOLD effects alone. These results suggest that beyond a strong thalamic activation impact, the scalp-recorded CNV signal contains a significant amount of SMA-mediated top-down regulatory processes (Fig. 5).
Partial regression plots showing the relationship of the modulated DCM parameters from SMA→VS (left) and SMA→THAL (right) with the CNV amplitude (difference score: monetary reward anticipation > control) beyond pure THAL BOLD activation effects.
Discussion
In the present simultaneous EEG-fMRI study on reward anticipation, we investigated intermodal correlations between two widely known reward anticipation signals, namely the fMRI BOLD response in the VS, and the CNV observed in EEG. Using a modified version of a highly reliable fMRI reward task (Kirsch et al., 2003; Plichta et al., 2012) we showed that (1) the simultaneously recorded CNV component is detectable and sensitive to the experimental conditions, (2) THAL fMRI BOLD activity is mediating intermodal correlations between the CNV1 signal and the fMRI BOLD response in SMA and VS, (3) the underlying causal connectivity network consists of top-down regulation from SMA to VS and SMA to THAL along with an excitatory information flow through a THAL→VS→SMA route during reward anticipation, and (4) the EEG CNV1 signal is best predicted by a combination of THAL fMRI BOLD response and strength of top-down regulation from SMA to VS and SMA to THAL.
The correlation of anticipatory BOLD signals and connectivity parameters with the scalp-recorded CNV1 component suggests that these two signals share similar neural mechanisms. This finding adds to the weight of evidence implied by many previous EEG studies (for review, see Broyd et al., 2012a) but direct empirical evidence has so far been lacking. However, it must be noted that such correlations do not imply that the measured EEG signals stem directly or exclusively from these subcortical regions (Cohen et al., 2009, 2011). Alternatively, top-down regulatory processes from the frontal cortex to the striatum (Cohen et al., 2012) as captured by the scalp-recorded CNV1 might include information about VS BOLD activity.
The revealed network presented in Figure 4 is in line with the understanding that thalamic nuclei supply striatal neurons with information of behaviorally significant events (Matsumoto et al., 2001). Behaviorally significant events transmitted along the thalamostriatal projections, in coordination with the motivational value of the dopamine inputs, are therefore a precursor for proper action selection through the basal ganglia–thalamocortical-striatal circuitry (Smith et al., 2004). Interestingly, direct THAL connectivity to motor areas including the SMA has been shown in the primate brain (Fang et al., 2006) and in the present study were found to be the most dominant excitatory link in the intrinsic connectivity model (Fig. 4B). An integrative view for the modulated network pattern would be that THAL is acting on VS via thalamostriatal connections (Groenewegen and Berendse, 1994; Kimura et al., 2004; Smith et al., 2004, 2009) signaling increased arousal and vigilance (Kinomura et al., 1996) associated with the anticipated monetary reward versus the control condition. The VS influence on the SMA, e.g., to enhance and/or suppress action related signals, does not necessarily imply a monosynaptic connection. Alternatively, the VS can act on frontal cortex via direct projections to the nucleus basalis in the basal forebrain which is the main source of cholinergic fibers (Haber, 2011). VS, in turn, is influenced by top-down-directed modulation from the medial prefrontal cortex as shown by Cohen et al. (2012). Termination of this signal flow in the SMA finally reactivates the thalamus via reciprocal corticothalamic projections (Haber and Calzavara, 2009).
From a practical perspective, the intermodal correlation we report here potentially provides a surface EEG, and therefore likely cortical, surrogate for neural processes during reward anticipation including modulation of VS activity. The human dopaminergic reward system can be directly stimulated by drugs of abuse such as cocaine (Breiter et al., 1997) and modulated by dopamine releasers like amphetamine (Knutson et al., 2004) or neuroleptic drugs (Abler et al., 2007). Therefore, anticipatory VS responsivity has become an attractive target for pharmaco-fMRI studies, i.e., studies wherein pharmacological compounds are administered to modulate the brain's response to task-based fMRI paradigms (Kassubek et al., 2011; Gilman et al., 2012; van Hell et al., 2012; van Wouwe et al., 2012). Given that EEG recordings are less demanding than fMRI, both technically and financially, identification of an EEG surrogate for reward anticipation would provide a useful alternative methodology to probe this brain system. To date, pharmaco-EEG studies have focused on feedback-related brain activity (Santesso et al., 2009) or CNV without reward-modulation (Kratz et al., 2012). Our results supply future pharmaco-EEG studies with an EEG surrogate more specific to reward anticipatory effects including modulation of VS activity. Again, the thalamic transmission in combination with top-down regulation from SMA to VS and THAL identified in our pathway analysis would be a likely neurobiological mechanism mapping a primarily subcortical process onto a cortical signature.
Although the CNV amplitude was significantly correlated with the fMRI signals measured in the VS and the THAL, it is highly unlikely that the measured EEG signal is generated in these subcortical regions. For example, neurons within the VS do not have a geometrically parallel organization like cortical pyramidal cells. Therefore, it is implausible that potentials from the striatum propagate to the scalp (Cohen et al., 2011). On the other hand, even brainstem structures generate EEG signals detectable through averaging (Jewett et al., 1970) and at least one recent positive report about scalp-measured subcortical EEG exists that needs further replication (Foti et al., 2011). However, even when the CNV signal, significantly correlated with fMRI VS signal, has a different origin the correlation itself is a predictor of VS activity and our DCM analysis support the assumption that the thalamic and VS signal modulate at least one important generator of the cortical CNV in the SMA. Correlation of CNV amplitudes with experimentally induced VS changes would be a critical prerequisite that could be tested by future studies.
Limitations
Using only two anticipation conditions and immediate feedback to ensure enough trials for EEG averaging in the applied paradigm still yielded typical reward anticipation effects unlike plain attentional or motor activation patterns. Studies with additional control conditions have identified contributions of valence and salience (Diekhof et al., 2012), and excluded confounds from unspecific activation (Kirsch et al., 2003), or arousal and motor preparation (Damen and Brunia, 1987; Kotani et al., 2001; Dillon et al., 2009). Specificity of our results by means of potential differences in monetary win versus loss anticipation cannot be tested. However, RT matched CNV data (Table 2) showed that pure motor preparation effects cannot explain the revealed differences between the experimental conditions (Sobotka et al., 1992).
Because of the subjects' motor responses and related artifacts we could not analyze target- or feedback-related EEG components (Broyd et al., 2012a). To achieve this, further modification of the reward paradigm would be necessary, i.e., the introduction of a larger delay period between response and feedback.
Joint EEG-fMRI trial-to-trial analysis might offer an interesting possibility of additional analysis (Bland et al., 2011) because of its' sensitivity for dynamic effects potentially revealing the involvement of additional brain regions (Scheibe et al., 2010). Recent studies have begun to focus on correlations of moment-to-moment fluctuations and behavior (Fox et al., 2007; Mars et al., 2008). However, because the trial-to-trial fluctuation effects in our study were rather small in size, further modifications of the applied reward paradigm would again seem to be necessary to make it more suitable for EEG recordings and this type of analysis. For example, single-trial analysis suffers from the superposition of task-relevant signals by task-unrelated “background” brain activity. Therefore, the resulting signal-to-noise ratio (SNR) of the observed single trial is expected to be rather low. Methodological studies should investigate the suitability of the reward task for single-trial analysis in combination using sophisticated statistical approaches (Debener et al., 2006). For the present study, our primary aim was to increase the trial number from 10 (Plichta et al., 2012) to 50 to improve the SNR in the EEG responses but other characteristics, such as the time interval between S1 and S2, could be reconsidered as well (Birbaumer et al., 1990).
Finally, the impact of different fMRI preprocessing strategies (e.g., smoothing kernel) on localization of reward-related brain activity (Sacchet and Knutson, 2012) was not tested here but should be considered in future studies.
Conclusions
In a simultaneous EEG-fMRI study on reward anticipation, we have demonstrated intermodal correlations between two widely known fMRI and EEG signals sensitive to reward anticipation, namely BOLD response within core structures of a reward anticipation network and the EEG CNV signal. We showed that (1) THAL fMRI BOLD activity is mediating intermodal correlations between the EEG CNV1 signal and the fMRI BOLD response in SMA and VS, (2) the underlying causal connectivity network consists of top-down regulation from SMA to VS and SMA to THAL along with an excitatory information flow through a THAL→VS→SMA route during reward anticipation, and (3) the EEG CNV signal is best predicted by a combination of THAL fMRI BOLD response and strength of top-down regulation from SMA to VS and SMA to THAL. The intermodal correlations enable predictions of the responsiveness of the reward anticipation network based on a scalp surface recorded EEG signal.
Footnotes
This work was supported by the German Ministry of Education and Research (BMBF, 01GQ1003B), the Deutsche Forschungsgemeinschaft (SFB 636), and NEWMEDS. The work leading to these results has received funding from the Innovative Medicines Initiative Joint Undertaking under Grant agreement no. 115008. We thank Dagmar Gass, Heike Schmidt, Vera Zamoscik, Christine Niemeyer, and Stefanie Rechtsteiner for assistance with data collection.
- Correspondence should be addressed to Dr Michael M. Plichta, Central Institute of Mental Health, Department of Psychiatry and Psychotherapy, Medical Faculty Mannheim/Heidelberg University, J5, 68159 Mannheim, Germany. Michael.Plichta{at}zi-mannheim.de