Abstract
Dopamine (DA) levels in the striatum are increased by many therapeutic drugs, such as methylphenidate (MPH), which also alters behavioral and cognitive functions thought to be controlled by the PFC dose-dependently. We linked DA changes and functional connectivity (FC) using simultaneous [18F]fallypride PET and resting-state fMRI in awake male rhesus monkeys after oral administration of various doses of MPH. We found a negative correlation between [18F]fallypride nondisplaceable binding potential (BPND) and MPH dose in the head of the caudate (hCd), demonstrating increased extracellular DA resulting from MPH administration. The decreased BPND was negatively correlated with FC between the hCd and the PFC. Subsequent voxelwise analyses revealed negative correlations with FC between the hCd and the dorsolateral PFC, hippocampus, and precuneus. These results, showing that MPH-induced changes in DA levels in the hCd predict resting-state FC, shed light on a mechanism by which changes in striatal DA could influence function in the PFC.
SIGNIFICANCE STATEMENT Dopamine transmission is thought to play an essential role in shaping large scale-neural networks that underlie cognitive functions. It is the target of therapeutic drugs, such as methylphenidate (Ritalin), which blocks the dopamine transporter, thereby increasing extracellular dopamine levels. Methylphenidate is used extensively to treat attention deficit hyperactivity disorder, even though its effects on cognitive functions and their underlying neural mechanisms are not well understood. To date, little is known about the link between changes in dopamine levels and changes in functional brain organization. Using simultaneous PET/MR imaging, we show that methylphenidate-induced changes in endogenous dopamine levels in the head of the caudate predict changes in resting-state functional connectivity between this structure and the prefrontal cortex, precuneus, and hippocampus.
Introduction
Methylphenidate (MPH; e.g., Ritalin) is one of the most common psychostimulants prescribed to treat children with attention deficit hyperactivity disorder (ADHD) (Arnsten, 2006). It binds to the dopamine (DA) transporter (DAT) and norepinephrine (NE) transporter and blocks neurotransmitter reuptake (Krause, 2008; Hannestad et al., 2010), thereby increasing their extracellular levels (Volkow et al., 2002), which alters behavioral/cognitive functions thought to be under the control of the PFC dose-dependently (Sprague and Sleator, 1977; Arnsten, 2011; Rajala et al., 2012). The pattern of effects is described by an inverted-U dose–response function (Sprague and Sleator, 1977; Gamo et al., 2010; Rajala et al., 2012). Prior studies at the systems-wide level have examined the effects of MPH on the networks that connect different brain regions, but questions remain about the relationship between changes to the organization of the brain's functional networks and changes to extracellular levels of DA brought about by the administration of oral doses of the drug (Zimmer, 2009).
Changes in extracellular concentration of DA can be measured with PET exploiting the sensitivity of select PET radioligands to endogenous neurotransmitter competition for receptor binding. For example, the binding of the D2/D3 DA receptor antagonist [18F]fallypride (Mukherjee et al., 1997) can be measurably reduced following treatment with d-amphetamine or cocaine (Schlaepfer et al., 1997) as a result of increased synaptic DA competition for the binding sites. Similar changes in PET ligand binding have been reported with therapeutic doses of MPH (Volkow et al., 2002; del Campo et al., 2013), which binds to DAT, yielding significant increases in extracellular DA in the striatum, where DAT availability is highest (Ouchi et al., 1999). We chose [18F]fallypride because of its high affinity for the DA D2/D3 receptor, which allows for the measurement of changes in endogenous DA levels in areas where receptor density is lower, such as the cerebral cortex (Christian et al., 2000). Microdialysis measurements in monkeys using doses of MPH equivalent to those used to treat humans (Doerge et al., 2000) yielded similar results. However, increased levels of DA were observed in the PFC only at much higher, cognitively impairing doses (Kodama et al., 2017). Thus, the question remains as to how increased levels of DA in the striatum resulting from therapeutic doses of MPH could lead to changes in PFC function.
Changes in the functional organization of the brain measured with fMRI in humans in response to MPH have shown alterations in different brain networks, particularly the default mode and frontostriatal networks. Changes in functional connectivity (FC) in frontostriatal circuitry following administration of MPH have been mixed, however, with some studies reporting increases in frontostriatal connectivity (Nagano-Saito et al., 2008; Kelly et al., 2009; Rubia et al., 2009), and others reporting decreases (Farr et al., 2014; Mueller et al., 2014; Hong et al., 2015), which may be explained by the fact that most studies used a single dose of the drug and different experimental conditions. Task activation data also show increased activation in frontal regions and basal ganglia following MPH administration (Peterson et al., 2009; Tomasi et al., 2011).
In this study, we bridged the gap between neurochemistry and functional organization by simultaneously measuring changes in extracellular DA using PET imaging of the D2/D3-specific neuroligand [18F]fallypride, and FC with resting-state fMRI (rsfMRI) at various doses of MPH in awake rhesus monkeys (Fig. 1). Our results identify, for the first time, an MPH dose-dependent link between [18F]fallypride nondisplaceable binding potential (BPND) in the head of the caudate (hCd) and FC between this area, used as a seed, and prefrontal, hippocampal, and motor regions, which are involved in attention, memory, and motor functions affected by MPH.
Schematic representation of the study. MPH blocks DAT, which causes extracellular levels of DA to increase. Red “?,” Goal of the study: to link changes in extracellular levels of DA to changes in FC (ΔFC). Extracellular levels of DA were inferred by measuring changes in BPND of [18F]fallypride brought about by different doses of MPH in specific brain regions. Changes in FC were measured with rsfMRI. The link between BPND and MPH was established using a GLM.
Materials and Methods
Animals and procedures.
Three male rhesus monkeys (Macaca mulatta;10–12 years of age, 8–15 kg in weight), purchased from the Wisconsin National Primate Research Center, were used. The animals were individually housed in the same room and received identical husbandry and veterinary care; no other monkeys were tested. All experimental procedures were approved by the University of Wisconsin Institutional Animal and Care and Use Committee and were in accord with the National Institutes of Health Guide for the care and use of laboratory animals.
Behavioral training.
All training and behavioral procedures were based on positive reinforcement, which consisted of pieces of fruit and/or liquid. The animals were first trained to enter a custom-built primate chair, then trained to accept being placed in the sphinx position using a custom-feature of the chair that allowed for gradual forward tilting (15 degree steps). When the full sphinx position was achieved, the animals inside the chair were placed on the bed of a mock scanner for acclimation and exposure to the scanner's noise. Subsequently, they underwent a surgical procedure to attach an MR-compatible head post. Briefly, under sterile conditions, the top of the head of the monkeys was prepared for surgery. An incision was made rostrocaudally, at the center the head, the soft tissues retracted, the bone exposed, and 14 ceramic screws (Thomas Recording, Type SI06) were inserted in the skull, and a custom head post was attached to the screws with Justi Cold-cure ortho resin (American Tooth Industries). The post was placed caudally on the skull, centered approximately at interaural level, to minimize interference with MR and PET signals from the frontal part of the brain. The post was used with a coupler to restrain the subject's head using a feature of the MR-compatible primate chair. After the surgery, a period of 8 weeks was allowed to pass before restraining the subjects' heads to ensure that the implant area had healed properly, and that bone had grown around the ceramic screws. During this period, the animals continued to be placed in the primate chair and were trained on tasks that did not require restraining the head. Emphasis was placed on training the animals to move as little as possible.
Behavioral setup for scanning experiments.
In addition, the animals were also trained to enter a custom tabletop restraining device, and to accept the insertion of a removable 20 GA catheter, 0.25 or 2 inch length (Surflash IV, Terumo) into the cephalic vein of the leg. Once inserted, the catheter was secured in place with tape and vet wrap. The monkey was then guided to the primate chair for transport to the scanning facility. This was an essential step for the administration of the PET tracer, which was also accomplished with positive reinforcement.
Upon arriving to the scanner facility, the monkey was placed in the sphinx position and onto the bed in front of the scanner, where its head was restrained with a custom head-post holder, and a 3 inch receive-only MR loop coil was placed above the subject's head for imaging. A water spout was placed in front of the subject's lips to provide rewards during the preparation. In addition, the catheter was exposed and connected to a line filled with heparinized saline, and secured in place. The animal was inserted into the bore of the scanner tail first. A video-based eye tracker (Eyelink 1000Plus, SR Research, RRID:SCR_009602) was placed 150 cm away from the head, immediately in front of an MR-compatible computer screen used to present a dot for the monkey to fixate during data acquisition.
Experimental timeline.
After administering the randomly selected dose of MPH on the scanner bed, the setup of the subject and equipment was completed (5–10 min). A series of structural scans and a fieldmap for distortion correction were acquired (∼20 min) before starting the 10 min rsfMRI scans. Before each rsfMRI scan, the signals from the eye tracker were calibrated; during each 10 min run, the subject was presented with a central fixation dot using Experiment Builder (SR Research), fixation of which was not enforced. No rewards were administered during the resting-state scans. Approximately 45–60 min after oral administration of MPH, the tracer was injected through the catheter extension tube, which was positioned to the rear of the scanner during setup. A series of 10–12 10 min rsfMRI runs were collected for 3 h after tracer injection, simultaneous with the 3 h PET scan (Fig. 2).
Drug delivery and dosing.
MPH hydrochloride (Sigma-Aldrich) was administered orally, dissolved in 0.5 ml grape-cranberry juice after the subject was fully set up on the scanner bed, but before it was moved into the bore of the magnet. The doses used in this study (1.5, 3, and 6 mg/kg) were used previously (Rajala et al., 2012) and were chosen based on the work of Doerge et al. (2000) who determined that an acute oral dose of 3 mg/kg MPH in the monkey resulted in plasma levels similar to those used therapeutically in humans treated for ADHD. Drug dosing was calculated in milligrams per kilogram based on the animals weight the day of the experimental session. Vehicle (juice only) was administered in a similar manner as a control. Each subject underwent four scanning sessions, one with each of the three doses plus vehicle.
PET imaging.
[18F]Fallypride was synthesized as described previously (Mukherjee et al., 1995); 3 mCi was injected in an intravenous bolus, and data were acquired for 180 min (Fig. 2). Images were reconstructed by ordered subset expectation maximization to 1 min frames with corrections for attenuation and scatter using 511 keV attenuation coefficient maps based on a combination of MR images of the subjects and a coregistered CT image of the housing, MR coil, and other apparatus within the PET FOV.
PET image alignment.
PET images were aligned to a rhesus brain template as follows (Jenkinson et al., 2002; Moirano et al., 2018). To better visualize whole brain, sum images 0–15 min after injection of tracer were created. These 0–15 min images were aligned by 6 degrees of freedom (shifts and rotations) to the template, and the resulting transformation matrices were applied to the dynamic, 0–180 min, images. The dynamic images were motion corrected with respect to the image frame at 15 min using the FMRIB Software Library's (FSL) mcflirt. The 0–15 min sum images were again created and aligned within subject as follows. For each subject, images were first aligned by 6 degrees of freedom to the template, averaged, and aligned to the average. The average was again aligned to the template by 9 degrees of freedom (shifts, rotations, and zooms). The resulting transformation matrices were concatenated and applied to the dynamic images. This yielded motion-corrected images aligned within subject to the template space.
Pharmacokinetic modeling.
Time activity curves were calculated for five ROIs based on the Paxinos rhesus atlas: hCd, putamen, PFC, substantia nigra/ventral tegmental area (SN/VTA), and a cerebellar reference region (Paxinos et al., 2009; Moirano et al., 2018). Binding potentials with respect to nondisplaceable uptake (BPND) were estimated using the simplified reference tissue model (Lammertsma and Hume, 1996; Innis et al., 2007; Vernaleken et al., 2011). The binding potentials' relationship to MPH dose was examined by fitting the following mixed-effects GLM as follows:
where β0 and β1 are the intercept and slope fixed effects for the MPH dose of the kth (1 ≤ k ≤ 4) observation, aj and bj are the random effects specific to subject j (1 ≤ j ≤ 3), BPjk is the subject-specific binding potential response, and ejk is the random error term (Gelman and Hill, 2007). Mixed-effects models were fit using the MATLAB (RRID:SCR_001622) Statistics and Machine Learning Toolbox (fitlmematrix).
fMRI analysis.
All fMRI data analyses were performed using the AFNI software package (Cox, 1996) (RRID:SCR_005927) unless otherwise indicated. Specific programs used are provided in parentheses. T1-weighted structural MRI scans were aligned to a template brain (Calabrese et al., 2015) using a 12-parameter affine transformation (AFNI's @auto_tlrc script). Functional MRI echo-planar time series were first corrected for B0-field distortions using an acquired B0-field map and FSL's fugue. Data were then corrected for motion using rigid-body registration (3dvolreg) and corrected for slice-timing differences (3dTshift). Time points with a volume-to-volume motion (sum squared difference of the 6 realignment parameters) exceeding 0.2 mm, as well as the subsequent 10 TRs, were censored and not used in the subsequent analyses. On average, 7.6% of volumes were censored (median = 5.1%). The corrected EPI data were then aligned to the T1-weighted structural MRI using a 12-parameter affine transformation (align_epi_anat.py). This transformation was combined with the T1 to template transformation and applied to the EPI data in original space to warp the EPI data into template space with only one resampling step. The average signals in the white matter and CSF, as well as their temporal derivatives, were then regressed from the data (3dTproject). Finally, the data were spatially smoothed by 1 mm within a mask of the brain.
First, we examined the FC between the hCd, the region where PET [18F]fallypride BPND was most strongly and significantly associated with MPH dose, and the entire PFC gray matter (on average). The hCd was defined using a digitized version of the Paxinos atlas (Paxinos et al., 2009; Moirano et al., 2018). The template brain was segmented into gray matter, white matter, and CSF using FSL's fast. A PFC mask was created by taking all gray matter voxels anterior to the genu of the corpus callosum and dilating this mask by 1 voxel (2 mm) into the brain to account for the lower-resolution EPI and applied spatial smoothing. The goal of this analysis was to address the hypothesized effects of increased extracellular DA levels resulting from MPH on the frontal lobe, in general, in the control of cognitive functions. The following GLM was fit:
where β0 and β1 are the intercept and slope fixed effects for the fallypride binding potential BP of the kth observation, aj and bj are the random effects specific to Subject j, FCjk is the subject specific FC response, and ejk is the random error term. In this case, (1 ≤ k ≤ Kj), where Kj = 4 MPH doses × number of runs at each dose (KGO = 50, KGR = 48, and KMI = 49).
Voxelwise FC maps were generated by computing the average signal intensity time course in the caudate head and then computing the temporal correlation with the signal intensity time course in all other brain voxels (AFNI's 3dDeconvolve). Areas that showed significant association between [18F]fallypride BP and FC were then determined using a linear mixed-effects analysis (3dLME) across all subjects and imaging runs, with [18F]fallypride BP as a covariate. The 3 subjects were treated as random effects in the mixed-effects model. Voxelwise FC maps were corrected for multiple comparisons using a cluster-based approach and an individual-voxel p value threshold of 0.001. This cluster-based approach estimates the spatial autocorrelation function (as the sum of a Gaussian and exponential) of the following form:
with a = 0.82, b = 2.07, and c = 8.41 (using AFNI's 3dFWHMx), and then runs a Monte Carlo simulation (3dClustSim) to determine the likelihood of getting clusters of a certain size (Cox et al., 2017).
Results
Data were collected from 3 adult male rhesus monkeys (Subjects GO, GR, and MI) trained to sit quietly in the bore of the PET/MR scanner for up to 5 h to be imaged while at rest. Scans were performed following oral administration of MPH (Fig. 2) at four dose levels (0, 1.5, 3.0, and 6.0 mg/kg).
Experiment timeline. Upon arrival to the scanner, the first 30 min were used to set up the eye tracking and reward delivery equipment and position the subject on the scanner bed. Once the subject was ready, MPH was administered orally, just before moving the subject inside the bore of the magnet. The camera and display monitor were set up, and then MRI scanning was started. Structural MRI attenuation correction and fieldmap correction scans were collected. Forty-five to 60 min after MPH administration, MRI acquisition was stopped and the PET scan started 1 min before a bolus of [18F]fallypride was given via a catheter line from behind the scanner. Subsequently, a series of 10–12 10 min rsfMRI scans were acquired simultaneously with the PET scan for a total scan time of 180 min.
PET
Our first goal was to obtain an index of changes in endogenous DA concentration in striatal and cortical regions of the brain. BPND of the highly selective D2/D3 antagonist [18F]fallypride was measured by dynamic PET imaging. Elevated levels of D2/D3 specific binding were readily observed in the striatum (Fig. 3). Estimates of BPND were obtained from the time course of radiotracer concentrations in anatomically identified ROIs using the simplified reference tissue model (Fig. 4) (Lammertsma and Hume, 1996; Vernaleken et al., 2011) and the cerebellum as the reference region of negligible specific binding. This yielded a BPND for each of the 3 subjects at the four dose levels of MPH, given orally ∼1 h before the administration of [18F]fallypride (Fig. 2; Table 1). BPND was found to correlate negatively with MPH dose in the hCd, where BPND decreased by 1.37% (p = 0.005) per milligrams per kilogram of MPH, suggesting increased extracellular DA competition (to reduce [18F]fallypride binding) resulting from the administration of MPH (Fig. 5).
PET imaging of D2/D3 receptor availability. [18F]Fallypride binding map of Subject GR at baseline exhibiting high concentrations of available D2/D3 receptors in striatum. Green represents hCd. Sagittal, coronal, and axial slices shown at left striatum. Binding is shown on a nonlinear scale to visualize uptake in regions outside of striatum.
Pharmacokinetic modeling of [18F]fallypride binding. Time-activity curves for ROIs in Subject GR at baseline. +, Individual dynamic measurements. Solid lines indicate simplified reference tissue model fits to data. y axis indicates standardized uptake value (SUV), which was computed as radioactivity concentration/(injected radioactivity/body weight). Cb, Cerebellar reference region.
[18F]Fallypride binding (BPND) in ROIs
[18F]fallypride binding potentials from model fits to head of caudate time-activity curves, plotted as a function of administered MPH dose. This shows a significant dose dependence, which is expected to result from competition with increasing concentrations of extracellular DA (slope = −0.21; fractional change −1.37% per milligrams per kilogram, p = 0.0045). A linear mixed-effects model describes the relationship between BPND and MPH dose, where the intercept and slope are fixed effects and subjects were treated as random effects to model individual differences (intercepts and slopes). The conditional predicted responses and 95% point-wise CIs are shown for 3 individual subjects (solid color). Dotted line indicates the marginal fixed effect.
fMRI
Because MPH affects cognitive functions thought to be under the control of the PFC, our next goal was to determine whether the significant decrease in [18F]fallypride BPND observed in the hCd, as measured by PET (Fig. 5; Table 1), was correlated with the FC between the hCd and the PFC. Thus, we first examined the influence of MPH on the FC between the hCd and PFC gray matter as a whole. The data were fit with a mixed-effects model (Eq. 2). The analysis showed that the hCd-PFC FC was significantly (negatively) correlated with [18F]fallypride BPND in the hCd (p = 0.019, Fig. 6) for the multiple MPH doses. That is, lower [18F]fallypride BPND, indicating increased competition from endogenous DA, was associated with higher hCd-PFC FC.
a, FC was computed between two ROIs: the head of the caudate (hCd; left, green) and the entire PFC gray matter (right, green ROI). b, FC between hCd and PFC versus the [18F]fallypride binding potential. Each point is the FC computed from one 10 min rsfMRI run. A linear mixed-effects model describes the relationship between FC and BPND, where the intercept and slope are fixed effects and subjects were treated as random effects to model individual differences (intercepts and slopes). The conditional predicted responses and 95% point-wise CIs are shown for 3 individual subjects (solid color). Dotted line indicates the marginal fixed effect. hCd-PFC FC is negatively related to [18F]fallypride BPND (slope = −0.023, p = 0.019).
Subsequently, we focused on the association between [18F]fallypride BPND and hCd FC on a voxelwise basis across the entire brain, with corrections for multiple comparisons. This analysis revealed that hCd [18F]fallypride BPND was significantly (negatively) associated with the FC between the hCd and three clusters in the brain (Fig. 7): the dorsolateral PFC (dlPFC) (p = 0.0014), the hippocampus (p = 0.0003), and the precuneus area (p = 0.0032).
a, ROI (green) in the hCd that was used as a seed for computing voxelwise FC. b, Colored regions represent areas where hCd FC is significantly (negatively) correlated with [18F]fallypride binding potential (BPND) in hCd. There is a significant negative correlation between [18F]fallypride BPND and the connectivity between hCd and dlPFC (p = 0.0014), hippocampus (p = 0.0003), and precuneus (p = 0.0032). All p values are corrected for multiple comparisons.
Last, the results of the FC analyses for the three individual regions identified with the voxelwise approach show that the effects are not driven by outliers from a single subject (Fig. 8). The linear regression lines in these plots were computed separately for each set; that is, they are not the output of the voxelwise mixed-effects model and are considered descriptive, not inferential. The statistical significance of these effects is as indicated above by the voxelwise p values, corrected for multiple comparisons.
Plot showing individual runs of FC versus [18F]fallypride binding potential (BP) in hCd for the three significant clusters from Figure 7, computed using hCd as a seed region and (a) dlPFC (slope: −0.12023; intercept: 2.025), (b) hippocampus (slope: −0.094; intercept: 1.430), and (c) precuneus (slope: −0.161; intercept: 2.374). Each point represents FC computed from one 10 min rsfMRI run. Colors of the points represent each of the 3 subjects. The linear regression lines in these plots were computed separately for each set and are considered descriptive, not inferential.
Discussion
Central to the quest for understanding the mechanisms underlying the effects of therapeutically relevant doses of MPH on the function of the nervous system has been to determine how changes in cognitive functions thought to be under the control of the PFC result from increases in dopaminergic function in the striatum (Volkow et al., 2002; Kodama et al., 2017), where DAT expression is highest (Ouchi et al., 1999). This study sheds much needed light on this issue and points to a potential, biologically viable mechanism.
The data show that increasing doses of MPH led to decreasing [18F]fallypride BPND, indicative of increasing levels of extracellular DA competing for D2/D3 DA receptor binding sites, particularly in the hCd nucleus. These decreases in BPND in the hCd are negatively related to the FC of this structure to the entire frontal lobe (Fig. 6). Subsequent voxelwise analyses revealed that decreases in BPND in the hCd are negatively related to the dlPFC, the hippocampus, and an area in the precuneus (Fig. 5). Thus, increasing levels of extracellular DA in the hCd are associated with increased connectivity of the hCd to these areas.
The dlPFC is a key component of frontostriatal circuitry; thus, increased hCd-dlPFC connectivity may be a mechanism by which MPH, acting primarily on the striatum, affects prefrontal cortical activity thought to underlie improvements in cognitive performance resulting from the administration of oral doses of MPH (Grace et al., 2007; Krugel et al., 2009). More specifically, the dlPFC is a key component of the attention network, involved in sustained attention. Increased hCd-dlPFC connectivity may therefore be a mechanism for the observed improvement in attention resulting from MPH.
[18F]Fallypride BPND was also significantly correlated with the connectivity between hCd and the most rostral region of the precuneus, an area thought to be involved in sensory motor function because it is functionally connected to the superior parietal cortex, paracentral lobule, and motor cortex (Margulies et al., 2009). Thus, increased connectivity of the hCd to this area may be a mechanism by which MPH reduces hyperactivity.
Last, the connectivity between the caudate and hippocampus has previously been shown to be associated with DA function (Farr et al., 2014). Our findings are consistent with this earlier study that also showed increased caudate-hippocampus FC following the administration of a single dose of MPH. Due to the central role of the hippocampus in memory, it is possible that this connection is a mechanism by which MPH influences memory function.
At the neuronal circuit level, these increases in FC could result from improvements in signal-to-noise ratio produced by increased levels of DA in the hCd, as demonstrated by the effects of iontophoretic application of DA in the immediate vicinity of single units recorded in the striatum of behaving monkeys (Rolls et al., 1984) and rats (Kiyatkin and Rebec, 1996). High-frequency broadband activity may serve as a proxy for local multiunit activity (Nir et al., 2008), which also covaries with spontaneous fluctuations in BOLD signals (Mantini et al., 2007; Shmuel and Leopold, 2008). DA is thought to reduce background firing rates more than that of the signal, such that information flow in the corticostriatal pathway is maximized (Bamford et al., 2004; Nagano-Saito et al., 2008; Bubl et al., 2015). To the degree that MPH further augments an increase in striatal signal to noise, the temporal correlation with functionally connected areas would also be expected to increase, consistent with our observations. In addition, it must be noted that MPH may also affect cortical function by other mechanisms, such as norepinephrine receptors (Arnsten, 2007), which was not studied here.
Several factors were essential for obtaining the present results. First, to directly relate the PET and fMRI measurements, the subjects needed to be imaged in the same physiological state. This was accomplished by using a simultaneous PET/MR system. This approach also allowed us to estimate [18F]fallypride BPND, a more biologically relevant method to measure changes in extracellular levels of DA than MPH dose alone given that the behavioral effects of the drug can vary considerably between subjects and between days (Rajala et al., 2012, 2015). In addition, FC, despite its robustness at delineating brain networks, fluctuates considerably at time scales of minutes to days, likely reflecting the state of the brain at each moment in time (Allen et al., 2014). Therefore, having simultaneous measurements of changes in extracellular DA and functional brain organization is key to linking these two measures. This approach will make it possible for future work to directly investigate the role of changes in DA levels on cognitive task performance and the underlying changes in brain activation.
Second, we chose an animal model over humans to accommodate the administration of various doses of MPH and multiple PET sessions. The awake rhesus monkey (M. mulatta) was selected because of its similarity to humans in the gross organization of the frontal lobes (Wise, 2008), cortical distribution of dopaminergic terminals (Berger et al., 1991), and functional resting-state networks (Hutchison et al., 2011). The awake state ensured that the results would not be compromised by potential spurious effects of anesthetics on the nervous system (Populin, 2005) and cardiovascular function (Brevard et al., 2003). Third, changes in DA levels brought about by a single dose of MPH may be different from day to day in the same subject, as evidenced by observed behavioral variability with the same dose level (Rajala et al., 2012, 2015); thus, the use of an animal model allowed for repeated measures in a within-subjects design. In addition, the administration of several doses was required to capture the pattern of changes occurring with increased dose that could explain the inverted-U profile of MPH dose–response functions seen with behavioral measures (Gamo et al., 2010; Rajala et al., 2012). Finally, the use of research animals allowed for more stringent control of outside factors that could influence these measures.
It is important to note that our measurements of baseline [18F]fallypride binding are in agreement with those reported previously (Christian et al., 2009), and the fractional change in [18F]fallypride BPND versus MPH dose is consistent with that observed in humans (del Campo et al., 2013) after accounting for species differences in MPH pharmacokinetics (Wargin et al., 1983; Doerge et al., 2000; Gill et al., 2012; Soto et al., 2012).
Thus, the present results demonstrate, for the first time, significant MPH-induced increases in extracellular DA levels using PET imaging of D2/D3 receptor availability (Mach et al., 1997; Gill et al., 2012; Allen et al., 2014) (Fig. 5) in the hCd of the rhesus monkey and, using simultaneous fMRI imaging, directly link these increases in DA to changes in FC between the hCd and prefrontal, hippocampal, and motor regions.
Footnotes
This work was supported by the University of Wisconsin–Madison Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation, and the Defense Advanced Research Projects Agency Biological Technologies Office Targeted Neuroplasticity Training Program under the auspices of Doug Weber and Tristan McClure-Begley through the Space and Naval Warfare Systems Command Systems Center Pacific Grant N66001-17-2-4010. We thank the staff of the Wisconsin National Primate Research Center and the imaging technologists of the University of Wisconsin Radiology Department. Research reported in this publication was supported in part by the Office Of The Director, National Institutes of Health under Award Number P51OD011106 to the Wisconsin National Primate Research Center, University of Wisconsin-Madison. This research was conducted in part at a facility constructed with support from Research Facilities Improvement Program grant numbers RR15459–01 and RR020141–01. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
The authors declare no competing financial interests.
- Correspondence should be addressed to Luis C. Populin at lpopulin{at}wisc.edu