Cortical D1 and D2 dopamine receptor availability modulate methylphenidate-induced changes in brain activity and functional connectivity

Dopamine signaling plays a critical role in shaping brain functional network organization and behavior. Prominent theories suggest the relative expression of D1- to D2-like dopamine receptors shapes excitatory versus inhibitory signaling, with broad consequences for cognition. Yet it remains unknown how the balance between cortical D1R versus D2R signaling coordinates the activity and connectivity of functional networks in the human brain. To address this, we collected three PET scans and two fMRI scans in 36 healthy adults (13 female/23 male; average age 43 ± 12 years), including a baseline D1R PET scan and two sets of D2R PET scans and fMRI scans following administration of either 60 mg oral methylphenidate or placebo (two separate days, blinded, order counterbalanced). The drug challenge allowed us to assess how pharmacologically boosting dopamine levels alters network organization and behavior in association with D1R-D2R ratios across the brain. We found that the relative D1R-D2R ratio was significantly greater in high-level association cortices than in sensorimotor cortices. After stimulation with methylphenidate compared to placebo, brain activity (as indexed by the fractional amplitude of low frequency fluctuations) increased in association cortices and decreased in sensorimotor cortices. Further, within-network resting state functional connectivity strength decreased more in sensorimotor than association cortices following methylphenidate. Finally, in association but not sensorimotor cortices, the relative D1R-D2R ratio (but not the relative availability of D1R or D2R alone) was positively correlated with spatial working memory performance, and negatively correlated with age. Together, these data provide a framework for how dopamine-boosting drugs like methylphenidate alter brain function, whereby regions with relatively higher inhibitory D2R (i.e., sensorimotor cortices) tend to have greater decreases in brain activity and connectivity compared to regions with relatively higher excitatory D1R (i.e., association cortices). They also support the importance of a balanced interaction between D1R and D2R in association cortices for cognitive function and its degradation with aging.

T he neurotransmitter dopamine modulates the function of widespread brain networks, influencing a diverse range of behaviors including reward, motivation, motor, and cognitive function 1 . Dopamine exerts influence through two main types of receptors, the D1-like and the D2-like receptors (D1R/ D2R), which activate and inhibit adenylate cyclase, respectively, producing opposing effects on neuronal firing and fine tuning each others actions 2 . Yet it remains unclear how the balance between dopaminergic D1R and D2Rs signaling in the cortex impacts human brain function. The most widely used tool for studying human brain function relies on fMRI to measure bloodoxygenation level dependent (BOLD) signals as an index of brain activity and functional connectivity 3 . A better understanding of how dopamine signaling systematically influences the BOLD signal would benefit research on normal brain function as well as research on conditions with abnormal dopaminergic signaling, including substance use disorders, attention deficit/hyperactivity disorder, schizophrenia, and Parkinson's disease.
Early studies found that task-evoked dopamine release increased striatal BOLD signals 4,5 , but could not speak to how D1R/D2R signaling contributed to the results. Further, although midbrain dopamine neurons have direct projections to striatum, their signals cascade throughout the brain 6,7 based on topographically-organized cortico-striatal 'loops' 8 . These larger patterns of signaling play a fundamental role in shaping healthy behavior, especially cognition 9 , and have been largely understudied. Most studies in this space have used drug challenges to manipulate DA receptor signaling while observing changes in fMRI signals (i.e., pharmarcological fMRI), or correlated static PET measures of dopamine receptor availability with fMRI [10][11][12][13][14] . While these studies provided important insights, it is difficult to identify an overarching pattern of how D1R and D2R signaling in cortical and subcortical regions impacts brain function, because these studies tended to focus on specific seed regions that differed across studies. In addition, measuring receptor availability is only one part of the picture; it remains unknown how dopamine increases (e.g., from pharamacological manipulation) would affect brain activity and connectivity based on D1R and D2R availability and to their relative expression (i.e., the D1R/D2R ratio) across the brain. Some of the best evidence to date comes from nonhuman primate studies. Using the knowledge that D1Rs are lower-affinity than D2Rs and only stimulated at high dopamine levels, Mandeville et al. found that high versus low doses of amphetamine produced increases and decreases in the striatal fMRI signal, respectively. These results informed a model whereby a high relative occupancy of D1Rs to D2Rs should produce increases in the fMRI signal, and vice versa.
We sought to test the predictions of this model in humans and explain dopamine-induced changes in brain activity and connectivity across the neocortex. To do this, we collected baseline measures of D1R and D2R availability using PET and resting functional connectivity using fMRI; we also administered methylphenidate (60 mg oral) prior to a second D2R PET scan and fMRI sessions (Fig. 1). Based on the model of Mandeville et al., we hypothesized that networks with a higher relative D1R/ D2R ratio (and presumably a greater capacity for excitatory responses to dopaminergic stimulation) would show increases in brain activity to the methylphenidate challenge. Conversely, networks with lower D1R/D2R ratios would show decreases in activity to methylphenidate. We also expected that changes in within-network functional connectivity would parallel the changes in activity, since studies utilizing both glucose metabolism and fMRI have found strong correspondence between local activity and functional connectivity [15][16][17] . Finally, due to a broad literature on the topic, we tested how relative D1R/D2R ratios were associated with age and with cognitive performance focusing on spatial working memory. Based on consistent preclinical evidence that dopamine receptor signaling in prefrontal and parietal cortex is critical for spatial working memory 18 , we hypothesized that the D1R/D2R ratio in the association (as opposed to sensorimotor) cortices would be especially relevant for cognition and would decrease with aging.

Results
We first mapped the relative availability of D1R, D2R, and the D1R/D2R ratio across the neocortex, and tested if these measures significantly differed across sensorimotor versus association cortices, using paired t-tests (map of the D1R/D2R ratio is shown in Fig. 2a). There were robust differences across canonical functional brain networks (Fig. 2b), such that relative D1R availability (t (35) = 46.34, p < 0.0001) and relative D2R availability (t (35) = 32.25, p < 0.0001) were much greater in association as compared to sensorimotor cortices. However, the regional differences were greater in D1R than D2R availability, which showed as a significantly higher relative D1R/D2R ratio in association as compared to sensorimotor cortices t (35) = 15.07, p < 0.0001; Fig. 2c). All relative PET measures met the assumption of a normal distribution (Kolmogorov-Smirnov distances < 0.12; all p's > 0.10). Inspection of individual networks revealed that D1R and D2R availability was highest in the cingulo-opercular network and lowest in the sensory hand network; meanwhile, the Study design and fMRI Pipeline. a Study design. On one day, participants underwent a [ 11 C]NNC-112 scan to assess D1R availability at baseline. Then, participants took the oral medication (either 60 mg methylphenidate, MP, or placebo, PLA) and underwent a [ 11 C]raclopride scan to assess D2R availability, and a resting fMRI scan to assess brain activity and connectivity. On a second day, the D2R and fMRI scans were repeated after the other oral drug was given (MP/PLA session order was counterbalanced). The second set of fMRI scans allowed us to estimate changes in brain activity and connectivity from MP-induced dopamine increases, relative to placebo. b fMRI pipeline. Multi-echo echo-planar imaging (EPI) scans were acquired and preprocessed (see Methods section). Then, signals were extracted from 10 cortical networks of the relative D1R-D2R ratio was highest in the fronto-parietal network and lowest in the sensory mouth network (Fig. 2d) 20 . There was moderate-to-strong correspondence (linear r = 0.484, quaratic r = 0.759; see Supplementary Fig. 1), and results from both radiotracers showed higher D2R availability in association as compared to sensorimotor cortices.
We also examined whether non-relative striatal D1R and D2R availability, D1R/D2R ratio, and methylphenidate-induced dopamine increases were associated with changes in brain activity and connectivity. To replicate prior work, we confirmed that methylphenidate significantly decreased D2R availability in ventral (t (34) = 4.835, p < 0.0001) and dorsal striatum (t (34) = 6.213, p < 0.0001), reflecting dopamine increases. At an uncorrected threshold, D1R in dorsal and ventral striatum were positively correlated with baseline (placebo session) brain activity in sensorimotor cortices (for each dorsal and ventral striatum, r 2 = 0.115, p = 0.043); however these results would not survive correction for multiple comparisons. None of the other baseline PET measures significantly correlated with baseline measures of activity/connectivity in association nor sensorimotor cortices, even at an uncorrected threshold. Further, methylphenidateinduced changes in D2R availability (dopamine increases) in ventral and dorsal striatum did not significantly correlate with methylphenidate-induced changes in brain activity/connectivity. Finally, we tested whether the relative D1R/D2R ratio was significantly associated with age and spatial working memory performance. As hypothesized, we showed a significant negative correlation with age in association (r 2 = 0.206, uncorrected p = 0.006; Fig. 4a) but not sensorimotor cortices (r 2 = 0.013, p = 0.508; Fig. 4b) such that there was an age-related decline in relative D1R/D2R ratio in the association cortices. Additionally we showed that the relative D1R/D2R ratio correlated positively with spatial working memory performance-that is, negatively correlated with task errors-in association (r 2 = 0.143, uncorrected p = 0.027; Fig. 4c) but not in sensorimotor cortices (r 2 = 0.030, p = 0.330; Fig. 4d) such that higher ratios were linked with better performance. Critically, neither the relative D1R nor  relative D2R availability on their own was significantly associated with age or SWM performance (all p's > 0.05), suggesting that the relative D1R/D2R ratio was more sensitive than either receptor alone.

Discussion
Using PET, fMRI, and a pharmacological challenge, we found evidence that methylphenidate-induced dopamine increases (and presumably also norepinephrine increases) triggered changes in brain function that differed between association and sensorimotor networks, and were associated with the relative availability of dopamine D1/D2Rs in these regions. Association networks, with a high relative concentration of excitatory D1Rs to inhibitory D2Rs, showed strong drug-induced increases in resting brain activity. Meanwhile, sensorimotor networks, with comparatively lower D1R/D2R ratios, showed decreases in resting brain activity and functional connectivity; these data largely confirmed predictions from an extant model 21 . Across individuals, a higher D1R/D2R ratio in association cortices was negatively correlated with age and positively correlated with spatial working memory performance. These data provide a possible framework for how medications like methylphenidate differentially affect resting functional networks, and highlight the importance of the relative balance of cortical D1R to D2R signaling in cognitive function and in brain aging.
Using data from a sample of healthy adults, we mapped the relative availability of D1R and D2R across the cortex, as previously done with PET measures of glucose metabolism using FDG 22,23 . Though [ 11 C]NNC-112 has long been used to map D1Rs in the neocortex, only recently have studies found that [ 11 C]raclopride can also be used to observe a reliable D2R signal extrastriatally 24,25 . By taking advantage of this information, we were able to identify regional patterns of the two receptor densities and, critically, the relative signaling of D1R to D2R ratios. We found that, while both D1Rs and D2Rs show a much higher concentration in association as compared to sensorimotor cortices, the difference was much greater for D1Rs. This resulted in a much higher D1R/D2R ratio in association cortices than in sesnorimotor cortices. To our knowledge this finding has not been described previously, but this may fit with human postmortem studies on D1R and D2R densities in different parts of the striatum, and what we know of these striatal subregions' anatomical projections. The head of caudate, which has strong reciprocal connections with association cortices including frontal and parietal cortices, has relatively higher D1R binding and lower D2R binding than the posterior putamen, which is highly interconnected with somatosensory cortices 8,26,27 . If the D1R/D2R ratio indeed follows these topographic cortico-striatal connections, we speculate that the high relative prevalence of D1Rs in association cortices has especial relevance for cognition: for instance, Alzheimer's disease, marked by aberrant brain function in predominantly association cortices (including the default mode network; refs. 28,29 ) is associated with a reduction in D1R but not in D2R availability, relative to controls 30 .
We hypothesized that, if D1R signaling is predominantly excitatory and D2R signaling inhibitory 2 , then dopamine increases would increase activity/connectivity in regions with higher D1R/D2R ratios, whereas regions with lower D1R/D2R ratios would decrease activity/connectivity. To test this, we administered a high dose of oral methylphenidate and found that Fig. 4 Associations of relative D1R/D2R ratio with age and spatial working memory performance. Correlations between relative D1R/D2R ratio and age (a, b) and cognitive performance (c, d). Trend lines reflect the regression line of best fit. WM working memory. **p < 0.01; *p < 0.05; ns not significant (p > 0.05).
our hypotheses were largely confirmed: resting fMRI activity (measured with fALFF) strongly increased after methylphenidate compared to placebo in association cortices (with higher D1R/ D2R ratios), but strongly decreased in sensorimotor cortices (with lower D1R/D2R ratios). These data bear out in humans the predictions of models for how D1-like and D2-like signaling affects neurovascular responses in non-human primates 21 . They also broadly align with studies showing that D1R binding in amygdala positively correlated with BOLD responses in the amygdala to fearful faces 31 , and with findings that D2R antagonists increase hemodynamic responses more in regions with higher D2R density 32 .
Within-network connectivity followed a similar, though not identical pattern: connectivity also significantly decreased in sensorimotor cortices but a much weaker decrease in association cortices. These findings largely confirmed our hypothesis and correspond well with prior FDG/fMRI studies which found a strong coupling between local activity and functional connectivity [15][16][17] . The methylphenidate-induced decreases in sensorimotor withinnetwork connectivity replicate prior work 33,34 and might be relevant to its therapeutic benefit: in adults with attention-deficit/ hyperactivity disorder (ADHD), symptoms of hyperactivity and restlessness are associated with higher functional connectivity in these circuits 35 . Of note, methylphenidate also boosts central norepinephrine in addition to dopamine, raising the possibility that these findings could be due to increases in norephinephrine. However, a study using atomoxetine, a selective norepinephrine transporter blocker, found a somewhat different pattern from that in the current study: in response to an oral 40 mg dose of atomoxetine, functional connectivity tended to decrease across the brain along a posterior-to-anterior gradient, decreasing most in regions with high connectivity at baseline 36 . The different findings should be interpreted with caution, since some of the methodology differs between studies; nevertheless, it appears that methylphenidate and atomoxetine produce fairly different effects on regional brain function.
Lastly, the relative D1R/D2R ratio in association cortices negatively correlated with age and positively correlated with spatial working memory performance across individuals. It is well-documented that various markers of dopaminergic function decline with age [37][38][39] . Although one human postmortem study suggested that striatal D1Rs and D2Rs appear to decline with age at a roughly equal rate, leaving the ratio unchanged 40 others found that D1Rs decline with age more rapidly than D2Rs 41 , and a recent meta-analysis of PET studies confirmed a more rapid decline in D1Rs (roughly 13-14% per decade) than D2Rs (roughly 8-9% per decade), both in frontal cortex and striatum 42 .
Here we replicate the finding of a declining D1R/D2R ratio with age and find significant effects in association (including prefrontal) cortices, but not in sensorimotor cortices. Age-related decline in dopamine receptors has been linked to reductions in prefrontal metabolism 43 and cognitive performance 44,45 . Further, prefrontal D1R and D2R signaling have each been specifically implicated in spatial working memory performance 18,46,47 . Here we highlight that the relative ratio of these two receptors may be a critical marker for understanding cognitive function throughout the lifespan. However, these data are cross-sectional, thus we cannot determine if age-related decline in the D1R/D2R ratio causes the lower spatial working memory performance.
There were several limitations to the current study. PET data were collected on two different scanners; however, we used methods to correct for any average differences in the PET images based on scanner, and more importantly, the primary analyses in this manuscript were conducted within-subjects. Other limitations include the radiotracers used. Specifically, for [ 11 C]NNC-112, though the signal is predominantly due to D1R, roughly 20% of the cortical signal may reflect binding to serotonin 5HT2a receptors 48 . More specific radiotracers are needed to obviate this confound. In the case of [ 11 C]raclopride, while it predominantly binds to D2-type receptors, its affinity is relatively low, thus specific to nonspecific signals from cortical regions are lower than other radiotracers such as [ 11 C]fallypride 19,49 . Still, consistent evidence suggests that a reliable and detectable specific signal exists 50 and our work follows a number of recent investigations that have successfully used [ 11 C]raclopride to study extrastriatal D2R [51][52][53] . Nevertheless, due to the somewhat lower specific signal of [ 11 C]raclopride, the individual difference findings should be interpreted with some caution and viewed as preliminary data in need of replication. Additionally, while many prior studies have argued that fALFF is a good proxy for brain activity, the exact significance of fALFF remains unknown. Studies using concurrent electrical recording of neuronal activity and fMRI in monkeys and rodents have found that low-frequency BOLD fluctuations are strongly associated with gamma-band local field potential activity and to a lesser extent multiple-unit neuronal spiking activity [54][55][56] . Nevertheless, the precise relationship between BOLD amplitude and neuronal activity remains under debate. Therefore, the interpretation of fALFF as reflecting 'activity' should be done cautiously.
Future studies could expand on this work by utilizing simultaneous PET-pharmacological fMRI 57 . There are notable advantages to this approach. Compared to sequential acquisition of PET and fMRI images, simultaneous measurement would likely yield a stronger relationship between pharmacologically-induced changes in PET and fMRI measures, since it would eliminate temporal differences in drug bioavailability. Further, it would eliminate temporal differences in an individual's mental/physical state, which is especially relevant for fMRI, since some non-trivial component of resting fMRI measures appears to be state-based (studies estimate that the intra-class coefficient for resting fMRI measures varies between 0.5 and 0.8 58,59 ). Finally, simultaneous acquisition allows one to leverage temporal patterns of PET and fMRI signals to model molecular adaptations to drugs such as receptor internalization, as recently performed in the non-human primate 60 . These efforts hold the promise of gaining richer insight into how dopamine receptor function shapes brain activity and connectivity.
In conclusion, we observed that relative availability of D1R and D2Rs and D1R/D2R ratio systematically varied across canonical sets of resting state networks. Methylphenidate produced a pattern of changes in brain activity and connectivity in a regional pattern aligning with these receptor densities, which was predicted by extant models of dopamine receptor signaling and fMRI responses 21 . Specifically, association cortices with higher relative D1R/D2R ratios, which presumably have greater potential for excitatory responses upon dopaminergic stimulation, tended to have greater increases in activity (and weaker decreases in connectivity) after a methylphenidate challenge, as compared to sensorimotor cortices with lower D1R/D2R ratios where methylphenidate decreased activity and connectivity. The D1R/ D2R ratio in association cortices also seems to have relevance for aging and cognitive function, which highlights the potential of studying these systems to understand neuropsychiatric disorders including substance use, attention-deficit/hyperactivity disorder and Parkinson's disease 61,62 .

Methods
Participants. Data from 36 healthy adults were included in the study (23 male, 13 female, age 22-64; for participant characteristics see Table 1). All participants provided written informed consent, and the Institutional Review Board committee of the National Institutes of Health approved the study. Participants were excluded if they had a history of substance abuse or dependence (other than nicotine) or a PET acquisition and drug administration. PET scans were used to measure D1R availability with [ 11 C]NNC-112 and to measure D2R availability with [ 11 C] Raclopride. For each individual, studies were conducted on one of two scanners: a high-resolution research tomography (HRRT) scanner (n = 17; 7 female; Siemens AG; Germany) or a Biograph PET/CT scanner (n = 19; 6 females; Siemens AG; Germany). The use of two different scanners was necessary due to scheduling limitations at our site. The methods for correcting differences between scanners are described in the PET analysis section below. All [ 11 C]NNC-112 scans were conducted at 10AM, in a baseline state, without any drug manipulation. [ 11 C]raclopride scans were conducted on two separate days: once 1 h after administration of an oral placebo pill (to assess baseline D2R availability) and once 1 h after administration of 60 mg oral methylphenidate (single blind; counterbalanced session order). Raclopride scans were conducted at the same time of day (1 PM) and in the same scanner for a given participant.
For [ 11 C]NNC-112, emission scans were started immediately after a maximum injection of 555 MBq. Twenty-one dynamic emission scans were obtained from time of injection up to 90 min after. For [ 11 C]raclopride, emission scans were started immediately after a maximum injection of 370 MBq. Twenty-two dynamic emission scans were obtained from time of injection up to 60 min after. Dynamic emission scan images were evaluated before analyses to ensure that motion artifacts or misplacements were not included.
PET analysis. PET images were coregistered to the high-resolution MRI T1 and T2 structural images. We used the minimal preprocessing pipelines of the Human Connectome Project for the spatial normalization to the stereotactic MNI space of the structural and PET scans 63 . Differences in geometry and PSF between cameras (PET/CT = 4 mm PSF; HRRT = 2.7 mm PSF) resulted in systematic voxelwise differences in signal intensity between PET/CT and HRRT images. To correct for these scanner-specific scaling effects and harmonize the data we used a voxelwise approach based on grand-mean scaling. We used an updated version of the ComBat Harmonization technique implemented in the ENIGMA study 64,65 . Originally proposed by Johnson et al. (2007) and implemented in the surrogate variable analysis (sva) package in R 66 , ComBat uses an Empirical Bayes framework to estimate the distribution scanner effects. It was shown to be superior to other harmonization methods for varieties of data types including DTI 67 and cortical thickness 68 . We conducted ComBat separately for each tracer for the PET measure of interest (i.e., receptor availability) to harmonize the data across scanners. For [ 11 C]raclopride measures, since we had placebo and methylphenidate treatments, we used drug, age, sex (male/female), and race (4 groups; white, black, Asian, and others) as covariates in the model. For [ 11 C]NNC measures, since there was no drug manipulation, we used only age, sex, and race as covariates.
D1R/D2R availability: striatum. Time-activity curves in the dorsal striatum (caudate and putamen), NAc, and cerebellum were used to obtain the distribution volume ratios using a Logan reference tissue model 70,71 . The accumbens-tocerebellum and the dorsal striatum-to-cerebellum distribution volume ratios correspond to BPnd+1, which was used to quantify D1R and D2R receptor availability and the ratio of the availability of D1R to D2R (D1R/D2R). We averaged the values for caudate and putamen to create one 'dorsal striatum' ROI, since caudate and putamen BPnd are highly correlated with one another (across participants, r ≈ 0.9). We also used the D2R availability estimates to compute 'dopamine increases' based on previous work: Relative D1R/D2R availability: neocortex. To assess regional differences in D1R/ D2R availability, we computed relative images by normalizing the D1R/D2R DISTRIBUTION VOLUME RATIO maps to the whole brain mean (FSL's MNI_T1_2mm_brain_mask image). D1R/D2R availability in the neocortex as measured by PET is dominated by a global signal with very little regional variability: across participants the correlation between receptor availability in different cortical networks is very high (e.g., ref. 13 ). Since our measures of brain function with fMRI (activity and connectivity) vary greatly across the networks, the relative measures of D1R/D2R availability enhance regional differences and are better suited for comparison with regional fMRI measures. Note that, although [ 11 C]Raclopride has been traditionally used only to measure striatal dopamine receptor availability due to lower D2R estimates in neocortex, recent work has found that it also produces reliable measures of exstrastriatal D2R 24,25 . Nevertheless, this remains under some debate and since [ 11 C]Raclopride does have lower sensitivity to detect extrastriatal D2R than some other radiotracers 49 , we took several precautions in the current study. First, we used spatial correlation to compare D2R availability in the neocortex between the group average Raclopride map in the current study and the group average [ 18 F]Fallypride map of 25 healthy adults from a publicly availably atlas 20 , since Fallypride has higher extrastriatal D2R sensitivity 19 , and demonstrated moderate-to-strong correspondence (see Supplementary Fig. 1). Second, we limited our analysis to large-scale sets of networks (sensorimotor versus association cortices; see 'Brain network parcellation scheme' section below) to reduce the possibility of insufficient signal-to-noise ratio when examining small individual regions with Raclopride. Finally, we only examined estimates of receptor availability and did not attempt to examine changes in receptor availability in neocortex to the methylphenidate challenge (i.e., dopamine increases), since drug-induced changes in Raclopride binding may not be reliably detected in some regions of neocortex 19 .
Brain network parcellation scheme: sensorimotor and association cortices. For resting fMRI network activity/connectivty and relative D1R/D2R maps, we followed procedures from prior work 22,72 and extracted data from broad sets of sensorimotor and association regions delineated using a popular resting fMRI parcellation scheme 73 (Fig. 1b). We used the consensus parcellation scheme (264 regions of interest) and created a 5 mm sphere at each region from which we extracted data. After computing each measure of interest (see processing sections below) we aggregated across networks by taking the average value across all sensorimotor and all association networks, following the precedent of prior studies 72 : Hand somato-motor, mouth somato-motor, visual, and auditory networks contributed to the sensorimotor system, and default mode, fronto-parietal, ventral attention, dorsal attention, cingulo-opercular, and salience networks contributed to the association system.
normalization to the stereotactic MNI space. 0.008-.09 Hz band-pass filtering was used to assess the low frequency fluctuations in the resting fMRI data. Signals from the white matter and CSF were regressed out of the data. Framewise displacements (FD) were computed from head translations and rotations using a 50 mm radius to convert angle rotations to displacements. Scrubbing was used to remove time points excessively contaminated with motion. Specifically, time points were excluded if the root mean square change in the BOLD signal (DVARS) from volume to volume met the criteria: DVARS > 0.5% and FD > 0.5 mm. For a measure of activity, we used the fractional amplitude of low-frequency fluctuations (fALFF), which maps spontaneous fluctuations in the 0.01-0.10 Hz frequency band. For a measure of connectivity, we used standard measures of within-network functional connectivity: Pearson correlation coefficients were calculated between pairs of ROI-averaged time courses after 0.01-0.10 Hz bandpass filtering and normalized to z-scores using the Fisher transformation.
Cogniton: spatial working memory performance. To assess cognitive performance, we used the spatial working memory task from the Cambridge Neuropsychological Test Automated Battery suite 74 . Participants completed this task in a baseline state, i.e. without having received a dose of methylphenidate or placebo. In this task, the participant sees an arrangement of colored squares on the screen. The computer hides a yellow token in one of the squares, and the participant must find it. Once the participant finds the token, the computer will hide it in another square, and it will never use the same square twice. There is a black bar on the side of the screen, where the participant puts each token once they have found it. The task is complete when the participant has found all of the tokens and filled the black bar. The primary outcome measure of this task is total number of errors, where fewer errors indicates better spatial working memory performance.
Statistics and reproducibility. Analyses were performed in R version 3.6.2 and in GraphPad Prism version 8.0.1.
To test for regional differences in the PET measures (relative D1R, relative D2R, and relative D1R/D2R ratio), we perfomed paired t-tests (association versus sensorimotor).
To test for methylphenidate-induced changes in the fMRI measures (brain activity/connectivity), we perfomed paired t-tests (placebo vs. methylphenidate), each for association and sensorimotor cortices, and Bonferroni-corrected for two comparisons. Then, to see if the pattern of methylphenidate-induced changes in brain activity/connectivity differed by cortical regions, we performed a two-way repeated measures ANOVA, with drug (placebo versus methylphenidate) and network (association versus sensorimotor) as factors, and examined the interaction effect.
As a control analysis, we also examined traditional measures of striatal receptor availability. We tested for methylphenidate-induced changes in the D2R availability (i.e., dopamine increases) and tested whether D1R, D2R, D1/D2 ratio, and dopamine increases correlated with baseline fMRI activity and connectivity, as well as methylphenidate-induced changes in activity and connectivity, both in sensorimotor and association cortices, using Pearson correlation.
Finally, we tested if the relative D1R/D2R ratio, each for association and sensorimotor cortices, was significantly associated with age and spatial working memory performance, using linear regression in Graphpad Prism. We hypothesized significant associations would be observed in the association cortices, based on a large literature showing that dopamine receptor signaling in association regions such as prefrontal/parietal cortices is critical for spatial working memory 18 . As an exploratory analysis, we repeated these tests in regression models that included sex and IQ as factors, using the lm function in R. Based on these tests, sex and IQ did not appear to play a major role in the results and these additional models are presented in Supplementary Note 1.
Due to poor image quality, fMRI data from one participant was removed. After this, fMRI analyses were performed twice: once with the entire sample (n = 35) and once after we removed participants with high levels of motion during resting fMRI (n = 5 removed due to >15% of timepoints scrubbed, reminder of particpants: n = 30). Since general findings did not change we report results from the full (n = 35) sample here.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
Summary data used to produce primary results are in a publicly available repository at: https://github.com/pmanza/Cortical-D1D2. Source data underlying figures are also available in Supplementary Data 1-2.
Code availability R scripts (R version 1.2.5019) used to produce primary results are in a publicly available repository 75 at: https://github.com/pmanza/Cortical-D1D2.