Abstract
Functional connectivity (FC) analyses of resting-state fMRI data allow for the mapping of large-scale functional networks, and provide a novel means of examining the impact of dopaminergic challenge. Here, using a double-blind, placebo-controlled design, we examined the effect of l-dopa, a dopamine precursor, on striatal resting-state FC in 19 healthy young adults. We examined the FC of 6 striatal regions of interest (ROIs) previously shown to elicit networks known to be associated with motivational, cognitive and motor subdivisions of the caudate and putamen (Di Martino et al., 2008). In addition to replicating the previously demonstrated patterns of striatal FC, we observed robust effects of l-dopa. Specifically, l-dopa increased FC in motor pathways connecting the putamen ROIs with the cerebellum and brainstem. Although l-dopa also increased FC between the inferior ventral striatum and ventrolateral prefrontal cortex, it disrupted ventral striatal and dorsal caudate FC with the default mode network. These alterations in FC are consistent with studies that have demonstrated dopaminergic modulation of cognitive and motor striatal networks in healthy participants. Recent studies have demonstrated altered resting state FC in several conditions believed to be characterized by abnormal dopaminergic neurotransmission. Our findings suggest that the application of similar experimental pharmacological manipulations in such populations may further our understanding of the role of dopaminergic neurotransmission in those conditions.
Introduction
Functional connectivity (FC) analyses of resting-state fMRI data (Biswal et al., 1995) have recently provided novel insights into the brain's functional architecture. Region of interest (ROI)-based analyses elicit detailed maps of FC, which recapitulate functional distinctions within regions such as the anterior cingulate (Margulies et al., 2007), striatum (Di Martino et al., 2008) and hippocampus (Kahn et al., 2008). Furthermore, disrupted FC has been observed in several clinical conditions, including ADHD, schizophrenia and Alzheimer's disease, as well as in normal aging (Greicius et al., 2004; Andrews-Hanna et al., 2007; Zhou et al., 2007; Castellanos et al., 2008). Resting state FC approaches thus provide an ideal complement to, but not a replacement for, task-based approaches for mapping brain circuits in healthy and pathological populations.
Dopaminergic modulation of neuronal communication is crucial to normal functioning of brain circuits subserving thought and behavior (Montague et al., 2004; Robbins, 2005). Abnormal dopaminergic neurotransmission, particularly in striatal circuits, is associated with the alterations in thought, mood and behavior observed in many pathological conditions, as well as in normal aging (Bäckman et al., 2006; Graybiel, 2008). Resting state FC approaches offer an ideal means for interrogating large-scale striatal networks. Furthermore, the disruption of FC in populations characterized by dopamine (DA) abnormalities prompts the hypothesis that DA may play a role in the regulation of FC in striatal circuits.
We evaluated this hypothesis, through an examination of the effects of DA challenge on striatal FC. Previous studies suggest that DA manipulations can modulate both task-based (Honey et al., 2003; Nagano-Saito et al., 2008) and resting state FC (Li et al., 2000; Achard and Bullmore, 2007). Recently, we mapped the FC of six ROIs placed throughout caudate and putamen (Di Martino et al., 2008). We observed distinct patterns of FC for each seed, which paralleled motivational, cognitive and motor subdivisions predicted by models of striatal circuitry (Alexander et al., 1986; Haber, 2003), and which were consistent with a meta-analysis of corticostriatal coactivation across 126 task-based neuroimaging studies (Postuma and Dagher, 2006). Here, we applied the same systematic seeding approach to examine the effects of the DA metabolic precursor, l-dihydroxyphenylalanine (l-dopa) on resting state FC in 19 young adults, using a double-blind, placebo-controlled design. To our knowledge, this is the first study to use such a design to investigate the effects of a DA manipulation on seed-based resting-state FC.
We predicted that the previously observed patterns of striatal FC observed would be replicated in the placebo condition, but altered by l-dopa. Studies in PD patients suggest l-dopa improves or normalizes function within motor networks (Haslinger et al., 2001; Buhmann et al., 2003). We hypothesized that l-dopa would increase FC between striatal motor regions (putamen) and subcortical and cortical regions associated with motor function. Given studies demonstrating that l-dopa modulates task-evoked activation within frontostriatal circuits (Mattay et al., 2002; Rowe et al., 2008), and can improve performance on tasks that rely on those circuits (Cools, 2006), we hypothesized that l-dopa would also alter FC between caudate and prefrontal cortex.
Materials and Methods
Participants
Twenty healthy adults (8 females; mean age = 27.3 years), free from a history of neurological or psychiatric illness and psychotropic medications, contraindications to MRI or Madopar (l-dopa), participated in this study. Volunteers completed a medical questionnaire and interview with a recruiting nurse. One participant showed abnormally high l-dopa serum levels in the placebo condition (16,400 pmol/L), and was excluded from the study (19 remaining participants: mean age = 26.2 years, 7 females). Participants were reimbursed Au$30 to cover their travel costs. All procedures were reviewed and approved by the institutional review board of Queensland University.
Design
Participants were assigned to receive either placebo or l-dopa in session 1 and the other treatment in session 2 at least 2 weeks later, in a randomized double-blind, placebo controlled crossover design. l-dopa was administered in the form of a Madopar capsule containing 100 mg of levodopa (l-dopa) and 25 mg of benserazide (a peripheral decarboxylase blocker) by a research nurse. The dose was identical to that used previously in psychopharmacological studies reporting behavioral effects (e.g., Kischka et al., 1996; Angwin et al., 2004; Copland et al., 2009). Participants were scanned early in the morning and were advised not to eat breakfast since food may affect absorption rates. Placebo was administered in identical capsules. Functional imaging took place ∼45 min after ingestion of the capsule, consistent with peak plasma l-dopa levels (Olanow et al., 2000). Blood pressure and heart rate were measured before capsule ingestion and after the imaging session. Blood was taken immediately before the imaging session and plasma DA levels measured via gas chromatography/mass spectrometry. Participants also completed a standardized self-report measure of subjective feelings (e.g., alertness, lethargy, tension) before the imaging session (Bond and Lader, 1974).
Functional imaging
Functional imaging data were acquired using a 4T Bruker MedSpec system equipped with a transverse electromagnetic head coil (Vaughan et al., 2002). We collected one “resting state” scan, comprising 200 contiguous echo planar imaging (EPI) whole-brain functional volumes (TR = 2100 ms; TE = 30 ms; flip angle = 90°, 36 slices, matrix = 64 × 64; acquisition voxel size = 3.6 × 3.6 × 3.6 mm). During this scan participants were asked to relax with their eyes open. Before the resting-state scan, a point-spread function (PSF) mapping sequence was acquired to correct geometric distortions in the EPI volumes (Zaitsev et al., 2003). A high-resolution T1-weighted anatomical image was also acquired using a magnetization prepared gradient echo sequence (TR = 2500 ms; TE = 3.83 ms; TI = 1500 ms; flip angle = 8°; 256 slices, acquisition voxel size = 0.9 × 0.9 × 0.9 mm).
Image preprocessing and individual analyses
AFNI (Cox, 1996) was used to perform the initial preprocessing steps of slice timing, correction for interleaved acquisition (using Fourier-space time-series phase-shifting) and motion correction (by aligning each volume to a “base” image using Fourier interpolation). The first five volumes were discarded to allow for magnetization equilibration. All other data processing was performed using FSL (www.fmrib.ox.ac.uk). Further image preprocessing comprised spatial smoothing (using a Gaussian kernel of FWHM 5 mm), mean-based intensity normalization of all volumes by the same factor [each subject's entire four-dimensional (4-D) dataset was scaled by its global mean], and temporal bandpass filtering (high-pass temporal filtering: Gaussian-weighted least-squares straight line fitting, with σ = 100.0 s; Gaussian low-pass temporal filtering HWHM 2.8 s). Registration of high-resolution structural images to the MNI152 template (Montreal Neurological Institute) with 1 mm3 resolution was performed using the FSL linear registration tool FLIRT (Jenkinson and Smith, 2001; Jenkinson et al., 2002). Transformation to MNI152 standard space was then further refined using FNIRT nonlinear registration (Andersson et al., 2007a,b). Linear registration of each participant's functional timeseries to the space of the high-resolution structural image, was also performed using FLIRT.
Movement parameters
Movement in each of the cardinal directions (X, Y and Z), and rotational movement around 3 axes (pitch, yaw and roll) were calculated for each participant. Data were also visually inspected for movement-related artifacts.
Functional connectivity
Nuisance signal regression.
To control for the effects of physiological processes (such as fluctuations related to cardiac and respiratory cycles), and motion, we removed signal associated with several nuisance covariates. Specifically, we regressed each subject's 4-D volume on nine predictors that modeled nuisance signals from white matter (WM), CSF, the global signal, and six motion parameters. Correction for time series autocorrelation (prewhitening) was performed. The six motion covariates were generated as part of the output of the AFNI motion correction program 3dvolreg. The global signal regressor was generated by averaging across all voxels within the brain. Finally, to generate the WM and CSF nuisance covariates, we first segmented each individual's high-resolution structural image, using the FAST segmentation program provided by FSL. The resulting segmented WM and CSF images were then thresholded to ensure 80% tissue type probability. These thresholded masks were then applied to each individual's time series, and a mean time series was calculated by averaging across all voxels within the mask.
This nuisance signal regression step produced prewhitened, 4-D residuals for each participant. Voxelwise scaling was performed on the 4-D residuals, by dividing each voxel's time series by its SD. Performing this step ensures that the ensuing FC estimates represent partial correlations rather than regression parameter estimates (e.g., Vincent et al., 2006), and removes potential between-condition differences in the magnitude of BOLD fluctuations (Sorg et al., 2007). As a final preprocessing step, each participant's time series was spatially normalized by applying the previously computed transformation to MNI152 standard space, with 1 mm3 resolution. Specifically, the linear function describing the transformation from functional to high-resolution structural space was applied to each participant's scaled residuals volume using FLIRT. Subsequent to this initial linear coregistration step, the nonlinear warping transformations derived for that participant's high-resolution structural image were applied to the functional volume, to produce a scaled residuals volume in MNI152 space, with 1 × 1 × 1 mm resolution.
ROI selection and seed generation.
We used the same 12 striatal seeds (six in each hemisphere) described by Di Martino et al. (2008). Each seed ROI was approximately spherical (volume = 257 × 1 mm3 voxels, radius ∼4 mm). These were located (in MNI152 space) in the (1) inferior ventral striatum (VSi) (±9 9 −8, corresponding to the nucleus accumbens), (2) superior ventral striatum (VSs) (±10 15 0, corresponding to the ventral caudate), (3) dorsal caudate (DC) (±13 15 9), (4) dorsal caudal putamen (dcP) (±28 1 3), (5) dorsal rostral putamen (drP), (±25 8 6), and (6) the ventral rostral putamen (vrP) (±20 12 −3). As outlined in detail by Di Martino et al. (2008), these seed locations are consistent with both anatomical and functional subdivisions of the striatum (Heimer and Alheid, 1991; Drevets et al., 1999; Postuma and Dagher, 2006). The 12 seed ROIs were applied to each participant's prewhitened, 4-D residuals, and a mean time series was calculated for each seed by averaging across all voxels within the seed.
To further confirm the anatomical validity of our findings, and to control for potential interindividual anatomic variability, we conducted additional probability-weighted seed-based analyses of the striatum. We created right- and left-hemisphere masks for the nucleus accumbens, caudate, and putamen, as defined by the Harvard–Oxford Structural Atlas, a validated probabilistic atlas included with FSL, which parcellates each hemisphere into anatomically distinct cortical and subcortical regions (Kennedy et al., 1998; Makris et al., 1999). In this atlas, every voxel is assigned a value that corresponds to its probability of belonging to a given anatomical region (e.g., caudate). To control for potential interindividual anatomic variability for each striatal region, we used these probability values to weight each voxel's time series within that region (Stark et al., 2008). Mean time series were then extracted for each striatal region by averaging across all probability-weighted time series within that region, for each participant.
Seed-based functional connectivity analyses
We performed 12 separate multiple regression analyses in which we regressed each subject's 4-D residuals volume on each of the 12 seed time series, using the FSL program FEAT. These analyses produced subject-level maps of all voxels that were positively and negatively correlated with the seed time series. The same procedure was followed for the six probability-weighted ROIs.
Group-level analyses
Analyses of group-level FC for the l-dopa and placebo conditions, and pairwise l-dopa vs placebo comparisons were performed using a repeated-measures mixed-effects ordinary least-squares model as implemented in FSL. This group-level analysis produced thresholded Z-score maps (“networks”) of positive and negative connectivity for each striatal ROI and for each condition. Direct voxelwise condition comparisons produced thresholded Z-score maps of those voxels that showed significant l-dopa-related changes in FC with each ROI. For all analyses, corrections for multiple comparisons were performed at the cluster level using Gaussian random field theory (min Z >2.3; cluster significance: p < 0.05, corrected).
Results
Blood plasma DA levels were higher on l-dopa (mean 7145 pmol/L) compared with placebo (mean 302 pmol/L). Self-rated subjective mood did not differ significantly on l-dopa versus placebo, nor did blood pressure or heart rate (p > 0.05).
Movement
The root mean square (rms) of the movement parameters did not exceed 1 mm or 1° in any of the cardinal directions or rotational axes. There were no significant differences between conditions in the mean, maximum or rms of the movement parameters.
Striatal functional connectivity: placebo
Relative to our previous study of striatal FC (Di Martino et al., 2008), the present data were acquired from a different set of participants, in a different magnet of different field strength. Nonetheless, the patterns of functional connectivity observed for the placebo condition were remarkably consistent with those previously observed by Di Martino et al. (Figs. 1⇓⇓–4; supplemental Table S1, available at www.jneurosci.org as supplemental material). The three putamen seeds exhibited positive correlations with secondary motor cortical and subcortical regions, as well as prefrontal associative regions. Within the putamen, the network associated with the ventral seed (vrP) comprised more extensive regions of PFC than did the dorsal seeds (dcP and drP), which demonstrated correlations with superior premotor areas, the supramarginal and angular gyri that were not observed for the vrP seed (Figs. 1, 2).
Left hemisphere putamen functional connectivity. Thresholded Z-score maps (min Z >2.3; cluster significance: p < 0.05, corrected) of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity. Images are displayed according to radiological convention (left is right).
Right hemisphere putamen functional connectivity. Thresholded Z-score maps (min Z >2.3; cluster significance: p < 0.05, corrected) of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity. Images are displayed according to radiological convention (left is right).
In the caudate, we replicated the ventral-dorsal distinction previously observed by Di Martino et al. (2008). Whereas the more dorsal caudate seeds (DC and VSs) were associated with more dorsal areas of lateral and medial PFC, the VSi seed was primarily correlated with limbic regions (Figs. 3, 4). In supplemental Results A (available at www.jneurosci.org as supplemental material), we address in greater detail differences between the present set of analyses and those of Di Martino et al. (2008). Furthermore, because we did not have specific predictions concerning the effect of l-dopa on the negative FC of the striatum, we do not explicitly discuss the patterns of negative functional connectivity in the main text. However, the patterns of negative FC observed were consistent with those previously observed by Di Martino et al. (2008), and are shown in supplemental Figures S1–S4 (see also supplemental Results A and Fig. S5, available at www.jneurosci.org as supplemental material). Note that no significant effects of l-dopa on negative functional connectivity were observed, for any of the ROIs examined.
Left hemisphere caudate functional connectivity. Thresholded Z-score maps of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity. Images are displayed according to radiological convention (left is right).
Right hemisphere caudate functional connectivity. Thresholded Z-score maps of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity. Images are displayed according to radiological convention (left is right).
l-Dopa versus placebo: putamen
We performed direct comparisons of FC in the l-dopa and placebo conditions using a repeated-measures mixed-effects model. These analyses produced thresholded Z-score maps of those voxels that showed significant l-dopa-related changes in FC with each ROI (min Z >2.3; cluster significance: p < 0.05, corrected). Although the broad patterns of connectivity were highly similar for both the placebo and l-dopa conditions, we observed significantly stronger connectivity between all six putamen seeds and the anterior cerebellum and cerebellar vermis, in the l-dopa condition (Figs. 1⇑⇑–4; supplemental Fig. S6, available at www.jneurosci.org as supplemental material; Table 1) (cerebellar nomenclature derived from Schmahmann et al., 2000). Additionally, the right hemisphere ventral putamen seed, vrP exhibited increased connectivity with the dorsal midbrain and ventral brainstem (proximal to the pons), and all three left hemisphere putamen seeds exhibited increased connectivity with this ventral brainstem region. No regions of l-dopa-related decreases in functional connectivity were observed, nor were any significant differences in negative connectivity for any of the putamen seeds observed.
l-Dopa versus placebo: caudate
The group-level map of DC connectivity in the l-dopa condition (Figs. 3, 4) showed a striking reduction in connectivity between right DC and regions of the default mode network (see also supplemental Fig. S7, available at www.jneurosci.org as supplemental material). This was supported by direct comparisons, which revealed stronger correlations between the right DC and the posterior cingulate, ventromedial and dorsomedial PFC for placebo, relative to l-dopa. Significantly greater connectivity was also observed between both right and left VSs and posterior cingulate cortex (PCC), and between right and left VSi and PCC and medial occipital cortex (Figs. 3, 4; Table 2), for placebo, relative to l-dopa.
In the group-level maps of FC for the caudate ROIs, increased connectivity with lateral PFC in the l-dopa condition, relative to placebo is apparent (Figs. 3, 4, which show the group-level connectivity for the l-dopa and placebo conditions separately, see the top two rows for each seed). In direct comparisons, these differences only reached statistical significance for the left VSi seed: we observed increased FC between left VSi and right inferior prefrontal cortex, as well as between left VSi and the bilateral anterior and right posterior cerebellum, ventral brainstem and globus pallidus (Fig. 3, Table 2).
Because the default mode network is hypothesized to have a negative or even competitive relationship with regions related to attentional control, such as PFC (e.g., Fox et al., 2005; Fransson, 2005; Kelly et al., 2008), we further examined the changes in functional connectivity between left VSi and right vlPFC and between each of the caudate seeds and PCC (a core region of the default mode network) as a result of l-dopa administration. Specifically, we computed the mean difference in FC between the l-dopa and placebo conditions in the right vlPFC region that exhibited significant (i.e., min Z >2.3; cluster significance: p < 0.05, corrected) l-dopa-related increase in FC with left VSi (Fig. 5). We also computed the mean difference in FC between the l-dopa and placebo conditions in each of the PCC regions (clusters) that exhibited significant l-dopa-related decreases in FC with left VSi, right DC, right VSs and left VSs. We then calculated the correlation, across participants, between the l-dopa-related increase in FC and each of the l-dopa-related decreases. For each pairing, we observed an inverse correlation between the l-dopa-related changes in connectivity (Fig. 5): the greater the l-dopa-related increase in connectivity between left VSi and vlPFC, the greater the l-dopa-related decrease in connectivity for the same participant, between left VSi and PCC (r = −0.47, p < 0.05), between right DC and PCC (r = −0.63, p < 0.005), between right VSs and PCC (r = −0.47, p < 0.05). The same relationship just escaped significance for left VSs (r = −0.45, p = 0.051). This inverse relationship was not present for any of the other clusters which demonstrated l-dopa-related increases in connectivity with left VSi (i.e., cerebellum, brainstem and globus pallidus).
A–D, l-Dopa-related increase in VSi-vlPFC connectivity is inversely related to the l-dopa-related decrease in connectivity with PCC. We observed a significant inverse correlation between the l-dopa-related changes in connectivity: the greater the l-dopa-related increase in connectivity between VSs and PFC, the greater the l-dopa-related decrease in connectivity between the VSs and PCC for the same participant. Across participants, there were significant inverse correlations between the l-dopa-related changes in connectivity: the greater the l-dopa-related increase in FC between VSi and vlPFC, the greater the l-dopa-related decrease in FC for the same participant between left VSi and PCC (A), between right DC and PCC (B), between right VSs and PCC (C), and between left VSs and PCC (D).
Probability-weighted ROIs
Figures 6 and 7 show the patterns of FC observed for the six probability-weighted ROIs, corresponding to the nucleus accumbens, caudate and putamen. The patterns of FC associated with these larger ROIs were highly consistent with those observed for the smaller striatal seeds, albeit with less fine-grained detail, as would be expected for larger ROIs. The effects of l-dopa were also consistent with those observed for the smaller seeds, although more circumscribed (see Table 3). Both right and left nucleus accumbens ROIs exhibited reduced connectivity with posterior cingulate and medial occipital cortices as a result of l-dopa administration. Left nucleus accumbens also exhibited increased FC with (primarily posterior) cerebellum and brainstem, but not with right inferior prefrontal cortex, as was observed for left VSi. The right caudate ROI exhibited reduced FC with PCC, as was observed for the right DC seed. Finally, the putamen ROIs exhibited increased connectivity with regions of left cerebellum, and left putamen also exhibited increased connectivity with the ventral brainstem (left only). No significant differences in negative functional connectivity were observed.
Left hemisphere probability-weighted functional connectivity. Thresholded Z-score maps of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity, for each of the three probability-weighted ROIs. Images are displayed according to radiological convention (left is right).
Right hemisphere probability-weighted functional connectivity. Thresholded Z-score maps of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity, for each of the three probability-weighted ROIs. Images are displayed according to radiological convention (left is right).
Iterative connectivity analyses
To further explore the striking reductions in connectivity between the caudate and the default mode network in the l-dopa condition, we performed an additional FC analysis in which we examined the connectivity of the PCC region that demonstrated reduced connectivity with the right hemisphere probability-weighted caudate ROI. We first identified the region within PCC which exhibited the greatest connectivity with the right hemisphere probability-weighted caudate ROI in the placebo, relative to the l-dopa condition (Fig. 7). This peak coordinate (x = 0, y = −42, z = 38) then formed the center of a spherical seed ROI (volume = 257 × 1 mm3 voxels, radius ∼4 mm). Maps of significant positive FC for this PCC seed, and of significant condition-related differences in FC were then generated in the same manner as for each of the primary striatal seed ROIs.
The results of this iterative FC analysis support our primary findings of reduced connectivity between the caudate and regions of the default mode network. We observed a drastic reduction in the connectivity of the default mode network in the l-dopa condition, relative to placebo, particularly for correlations between the caudate and the PCC (Fig. 8). Specifically, l-dopa significantly reduced connectivity between the PCC seed and proximal PCC regions, rostral ACC/medial PFC, ventral striatum and dorsal caudate.
l-Dopa reduces functional connectivity within the default mode network. The iterative connectivity analysis, which examined the FC of the PCC region that demonstrated reduced connectivity with the right hemisphere probability-weighted caudate ROI (seed location shown on left of figure). The figure displays the thresholded Z-score maps of positive functional connectivity for each condition, and significant l-dopa-related increases and decreases in functional connectivity for this PCC seed ROI. A significant reduction in the connectivity of the default mode network in the l-dopa condition, relative to placebo, was observed. Images are displayed according to radiological convention (left is right).
Discussion
The present work is among the first resting state fMRI studies to examine the effects of a placebo-controlled dopaminergic manipulation on FC within large-scale networks (e.g., Achard and Bullmore, 2007). In the placebo condition, we replicated the findings of Di Martino et al. (2008), observing patterns of FC consistent with hypothesized motor (dorsal putamen), cognitive (ventral putamen, dorsal caudate, superior ventral striatum) and affective (inferior ventral striatum) subdivisions of the striatum. In line with our predictions, l-dopa increased FC in motor networks joining putamen and cerebellum and brainstem. In contrast, l-dopa reduced connectivity between all three caudate seeds and regions of the “default mode” network, particularly PCC. Furthermore, l-dopa-related decreases in FC between caudate and PCC were inversely correlated with an l-dopa-related increase in FC with left vlPFC (Fig. 5). Highly similar patterns of FC and l-dopa related changes in FC were observed when we used a set of anatomically defined, probability-weighted ROIs (Figs. 6, 7). Further exploratory examination of the posterior cingulate FC revealed a striking reduction of connectivity within the default mode network as a result of l-dopa administration (Fig. 8). Not only did l-dopa reduce striatal involvement within the default mode network, but it also reduced connectivity within PCC and between PCC and medial PFC.
l-Dopa and motor and cognitive functional connectivity
Acute administration of a low dose (100 mg) of l-dopa increased connectivity between putamen and brainstem and cerebellum, relative to placebo. This observation is consistent with neuroimaging studies demonstrating that l-dopa and other dopaminergic manipulations improve motor performance and normalize task-evoked activity in primary and secondary cortical motor areas, brainstem and cerebellum in both healthy and PD populations (Haslinger et al., 2001; Buhmann et al., 2003; Hershey et al., 2004; Asanuma et al., 2006; Palmer et al., 2009). Although earlier loop models of striatal motor networks (Alexander et al., 1986; Haber, 2003) did not consider subcortical interactions between striatum and cerebellum, recent studies have provided evidence of cerebellar input to the striatum, via a cerebello-thalamo-striatal pathway (Ichinohe et al., 2000; Hoshi et al., 2005). Moreover, interactions between striatal and cerebellar networks are increasingly implicated in movement disorders (Yu et al., 2007; Neychev et al., 2008) as well as in psychiatric conditions such as ADHD (Krain and Castellanos, 2006; Brennan and Arnsten, 2008).
l-Dopa also increased FC between the VSi, corresponding to the nucleus accumbens, and left vlPFC, relative to placebo. This observation mirrors the findings of a recent study (Nagano-Saito et al., 2008), in which DA depletion was associated with decreased connectivity between right dlPFC and dorsal caudate, and between both right and left vlPFC and the ventral striatum during a set-switching task. DA depletion also eliminated a positive relationship between frontostriatal FC and task performance, observed in the control (undepleted) condition. Conversely, studies have demonstrated improved performance on cognitive tasks subserved by frontostriatal circuits as a result of increased DA levels (Cools, 2006). Together, these studies suggest that FC and neuronal function within frontostriatal circuits is subject to modulation by dopaminergic signals from the striatum.
l-Dopa and the default mode network
Previous studies have demonstrated that activity (Shulman et al., 1997; McKiernan et al., 2003) and connectivity (Fransson, 2006) within the default mode network are reduced during task performance. The reduction in activity or connectivity typically occurs in association with increased activation in task-relevant regions, such as dorsal ACC and lateral PFC (Fransson, 2006; Tomasi et al., 2006). Failure to suppress default mode activity during task performance, coupled with decreased task-evoked activation in prefrontal areas has been related to impaired behavioral performance (slower RTs, increased errors) (Polli et al., 2005; Weissman et al., 2006). Dopaminergic manipulations may have similar effects. Nagano-Saito et al. (2008) demonstrated that under conditions of DA depletion, task-related suppression of default mode network activity was reduced, as was FC between PFC and striatum. They also observed a trend toward slower task performance in the DA depleted condition, relative to the undepleted condition. In our study, we augmented, rather than depleted DA, and, as would be predicted, our findings mirrored those of Nagano-Saito et al. Not only did l-dopa reduce FC within the default mode network, the reduction of FC between PCC and three of the caudate seeds (right dorsal caudate, right superior ventral striatum and left inferior ventral striatum) was significantly inversely related to the extent to which l-dopa increased FC between left nucleus accumbens (VSi) and right vlPFC.
How might DA bring about this bidirectional effect? Bamford and colleagues (2004a,b) propose that DA alters the signal-to-noise ratio within corticostriatal loops by reinforcing active corticostriatal synaptic connections and inhibiting inactive ones. Similar hypotheses have been proposed regarding the role of DA throughout the brain (Rolls et al., 1984; Kischka et al., 1996; Cepeda and Levine, 1998; Nicola et al., 2004), and regarding the role of DA in combination with norepinephrine (Arnsten, 2007). Speculatively, reduced connectivity within the default mode network, coupled with increased striatal-PFC or striatal-cerebellar connectivity might reflect such increased signal-to-noise. Striatal DA circuits may provide a mechanism for the active suppression of the default mode network under conditions that require increased processing of external stimuli (signal) relative to internal, self-directed processing (noise).
Limitations
Although we observed l-dopa-related increases in FC between dorsal putamen, a region hypothesized to serve primarily motor functions, and regions of the brainstem and cerebellum, we did not observe increased FC between putamen and primary or secondary motor cortical regions. One possible explanation is that we did not observe the same primary motor connectivity for the dorsal putamen seeds as did Di Martino et al. (2008; see supplemental Results A, available at www.jneurosci.org as supplemental material, for further details). A second possibility is that, as a subtherapeutic dose of l-dopa was administered, the effects of l-dopa on striatal motor networks were subtle. Larger effects extending to primary and secondary motor cortices may be observed at different doses, which should be assessed in a future study.
We demonstrated robust effects of l-dopa on striatal FC that we have hypothesized were related to the modulatory effects of DA on striatal motor and cognitive function. We and others have observed significant alterations in performance on behavioral tasks with the same dosage of l-dopa (Kischka et al., 1996; Angwin et al., 2004; Copland et al., 2009). However, there are interindividual differences in the effects of DA challenge on a given cognitive or motor function that reflect interindividual differences in baseline DA levels. Furthermore, inverse-U-shaped relationships exist both between DA levels and BOLD/PET activation levels, and between DA levels and behavioral performance (e.g., Cools, 2006; Rowe et al., 2008; Cools et al., 2009). Because we had no external measures of the impact of l-dopa on motor or cognitive performance, we cannot know whether the observed changes in FC were associated with improvement or decrements in performance. In fact, l-dopa-related reductions of FC between PCC and other default mode network components (e.g., medial PFC) may even have a negative impact on task performance, as suggested by a recent study relating working memory performance to default mode connectivity (Hampson et al., 2006). Future studies should address the potential effects of smaller and larger doses of l-dopa on both FC and performance on motor and cognitive tasks subserved by striatal circuits, as well as how individual differences impact the dose–response relationship.
Controversy surrounding the nature of negative relationships between networks (anticorrelations) in FC analyses persists (Skudlarski et al., 2008; Murphy et al., 2009). However, recent work examining the test-retest reliability of anticorrelations (Shehzad et al., 2009) and the emergence of anticorrelations in a computational model of cortical connectivity (Ghosh et al., 2008) encourages confidence in their fidelity. Although we did not explicitly discuss the patterns of negative functional connectivity associated with each our ROIs, they were consistent with those previously observed by Di Martino et al. (2008), and are shown in supplemental Figures S1–S4, available at www.jneurosci.org as supplemental material. Interestingly, we observed no significant effects of l-dopa on negative connectivity, suggesting that relationships within the functionally connected networks we observed were more sensitive to l-dopa than relationships between those networks and their negatively correlated counterparts. Further studies examining the impact of pharmacological manipulations on the negative relationships among networks are required to test this hypothesis. Such studies may shed light on the neurophysiological bases of negative relationships among networks.
Implications and future directions
Our observation of an l-dopa-related reduction of FC within the default mode network may help advance our understanding of previous studies that have demonstrated significant differences in default mode FC in conditions such as ADHD (Castellanos et al., 2008), and in normal aging (Andrews-Hanna et al., 2007; Damoiseaux et al., 2008). Dysregulation of catecholaminergic neurotransmission is thought to play a role in these conditions (Laruelle et al., 2003; Bäckman et al., 2006; Brennan and Arnsten, 2008; Tripp and Wickens, 2008; but see Gonon, 2009). Here, we have provided initial support for the hypothesis that DA contributes to the regulation of FC in the brain, bolstering the suggestion that the abnormalities in default mode connectivity in these conditions may be related to DA dysregulation. One way in which future studies might further examine this hypothesis is to examine whether a DA challenge normalizes FC of the default mode network in a population such as older adults.
Summary
We observed robust effects of acute low-dose l-dopa administration on striatal networks implicated in motor and cognitive function, in a sample of healthy young adults. Our results suggest that, although large-scale patterns of resting state FC are spatially and temporally stable, alterations of network coherence may result from changes in neurochemical state. Connectivity within the default mode network was particularly sensitive to DA manipulation. These findings are consistent with a role for DA in the modulation of FC observed during task performance, as well as in the altered FC observed in pathological conditions.
Footnotes
-
This work was supported by an Australian Research Council (ARC) Discovery grant (DP0452264) awarded to D.A.C. and K.M., and by a National Institute on Drug Abuse grant (1R03DA024775-01) awarded to C.K.
-
The authors declare no competing interests.
- Correspondence should be addressed to either of the following: Dr. Katie McMahon, Centre for Magnetic Resonance, University of Queensland, Brisbane, QLD 4072, Australia, katie.mcmahon{at}cmr.uq.edu.au; or Dr. Michael P. Milham, New York University Child Study Center, 215 Lexington Avenue, 14th Floor, New York, NY 10016, michael.milham{at}nyumc.org