Single-dose effects of methylphenidate and atomoxetine on functional connectivity during an n-back task in boys with ADHD

Rationale Working memory deficits and associated neurofunctional abnormalities are frequently reported in attention-deficit/hyperactivity disorder (ADHD). Methylphenidate and atomoxetine improve working memory performance and increase activation of regions under-functioning in ADHD. Additionally, methylphenidate has been observed to modulate functional networks involved in working memory. No research, however, has examined the effects of atomoxetine or compared the two drugs. Objectives This study aimed to test methylphenidate and atomoxetine effects on functional connectivity during working memory in boys with ADHD. Methods We tested comparative effects of methylphenidate and atomoxetine on functional connectivity during the n-back task in 19 medication-naïve boys with ADHD (10–15 years old) relative to placebo and assessed potential normalisation effects of brain dysfunctions under placebo relative to 20 age-matched neurotypical boys. Patients were scanned in a randomised, double-blind, cross-over design under single doses of methylphenidate, atomoxetine, and placebo. Controls were scanned once, unmedicated. Results Patients under placebo showed abnormally increased connectivity between right superior parietal gyrus (rSPG) and left central operculum/insula. This hyperconnectivity was not observed when patients were under methylphenidate or atomoxetine. Furthermore, under methylphenidate, patients showed increased connectivity relative to controls between right middle frontal gyrus (rMFG) and cingulo-temporo-parietal and striato-thalamic regions, and between rSPG and cingulo-parietal areas. Interrogating these networks within patients revealed increased connectivity between both rMFG and rSPG and right supramarginal gyrus under methylphenidate relative to placebo. Nonetheless, no differences across drug conditions were observed within patients at whole brain level. No drug effects on performance were observed. Conclusions This study shows shared modulating effects of methylphenidate and atomoxetine on parieto-insular connectivity but exclusive effects of methylphenidate on connectivity increases in fronto-temporo-parietal and fronto-striato-thalamic networks in ADHD. Supplementary Information The online version contains supplementary material available at 10.1007/s00213-023-06422-7.

Functional magnetic resonance imaging (fMRI) studies of working memory in youth with ADHD relative to neurotypical individuals show consistent increases in fronto-parietal coupling (Massat et al. 2012;Bédard et al. 2014;Wu et al. 2017), along with connectivity increases in networks comprising posterior areas such as cuneus, precuneus, and occipital regions (Wu et al. 2017), between occipital and cortico-striato-cerebellar regions, and between cerebellum and brainstem (Massat et al. 2012).Similarly, adults with ADHD relative to neurotypical adults show connectivity increases in various networks that include inferior and superior frontal regions, dorsal cingulate, and cuneus (Wolf et al. 2009).Additionally, abnormal reductions in working memoryrelated connectivity have been observed within executive control (Wu et al. 2017) and fronto-cingulate networks in youth (Bédard et al. 2014), and within fronto-cinguloparieto-cerebellar networks in adults with ADHD (Wolf et al. 2009).This variability in connectivity differences may be further exacerbated by the known functional brain maturation processes taking place in adolescence and the variability of ages of participants recruited for research (Pfeifer and Allen 2021).
The catecholamine transporter blocker methylphenidate (Wilens 2008;Faraone 2018) and the non-stimulant, noradrenaline transporter blocker atomoxetine (Bymaster et al. 2002;Gallezot et al. 2011) are the most commonly used pharmacological treatments for ADHD showing good clinical efficacy (Cortese et al. 2018).However, their specific neurofunctional mechanisms of action are not fully understood.While both methylphenidate and atomoxetine have shown to improve working memory performance of individuals with ADHD (Gau and Shang 2010;Yang et al. 2012;Shang and Gau 2012;Ni et al. 2013;Coghill et al. 2014;Pievsky and McGrath 2018b;Rubio Morell and Expósito 2019), distinct neurofunctional effects of stimulants and non-stimulants have been observed (Cubillo et al. 2014a).A meta-regression of stimulant effects on brain activation in ADHD during n-back working memory paradigms found that treatment was associated with greater activation in middle and superior frontal regions (McCarthy et al. 2014).A direct comparison of single doses of methylphenidate and atomoxetine on brain activation during the n-back task in the same cohort of children with ADHD found both drug-specific and shared effects (Cubillo et al. 2014a).Both drugs increased fronto-insular, temporal, and striatal activation in patients compared to unmedicated controls.However, atomoxetine selectively normalised abnormalities observed under placebo in the right dorsolateral prefrontal cortex (DLPFC) and upregulated this region relative to methylphenidate, while methylphenidate selectively enhanced activity of the left inferior frontal cortex (IFC)/DLPFC relative to controls during the 2-back condition (Cubillo et al. 2014a).
Also, methylphenidate has shown to normalise functional connectivity differences between boys with ADHD and neurotypical controls in executive control, fronto-parietal, and auditory networks, and to increase the connectivity within the executive control network compared to placebo (Wu et al. 2017).Studies further showed that, in youth with ADHD, stimulants decrease fronto-striatal connectivity compared to no medication (Sheridan et al. 2010) and increase the connectivity within the fronto-parietal network, as well as this network's functional connections with other cortical regions, including anterior cingulate cortex (ACC), ventrolateral prefrontal cortex, and precuneus compared to placebo (Wong and Stevens 2012).Nonetheless, no studies explored the effects of atomoxetine on working memoryrelated connectivity or compared the network effects of the two drugs in ADHD.Given that adolescence is a period of progressive neurofunctional maturation and specialisation of brain networks (Rubia 2013;Pfeifer and Allen 2021), it is crucial to understand the effects medications may have during this developmental stage.
Consequently, the aim of this pseudo-randomised, double-blind, placebo-controlled, cross-over fMRI study was to test the comparative effects of single doses of methylphenidate and atomoxetine relative to placebo on functional networks involved in working memory in medication-naïve boys with ADHD, and to test for potential normalisation effects of the two drugs on brain connectivity differences during placebo relative to unmedicated age-matched neurotypical controls.This study extends a previously published analysis of shared and distinct effects of the two drugs on brain activation during working memory in the same cohort (Cubillo et al. 2014a) to the network level.The advantage of single-dose comparisons is studying drug effects and probing the role of catecholamines without the confounds of symptomatic improvement, side effects, or brain adaptation following long-term dosing (Konrad et al. 2007;Nakao et al. 2011;Fusar-Poli et al. 2012;Frodl and Skokauskas 2012;McCarthy et al. 2014;Lukito et al. 2020).
In line with previous reports of fronto-parietal hyperconnectivity during working memory (Massat et al. 2012;Bédard et al. 2014;Wu et al. 2017), we expected greater frontal and parietal connectivity during working memory compared to the control task condition in individuals with ADHD under placebo relative to controls.Considering evidence of normalisation/upregulation effects of methylphenidate and atomoxetine on working memory-related brain activation (Cubillo et al. 2014a;McCarthy et al. 2014) and of normalising/modulating effects of methylphenidate on working memory-related functional connectivity (Sheridan et al. 2010;Wong and Stevens 2012;Wu et al. 2017), we expected that both medications would minimise the connectivity differences between patients and controls.

Methods
The analysis of this data focusing on brain activation has been published previously (Cubillo et al. 2014a).The analysis described here focuses on working memory-related functional connectivity.

Participants
Data from 19 medication-naïve boys (10-15 years old) with a diagnosis of ADHD combined presentation and 20 neurotypical control boys (10-15 years old) were included in this study.Patients were recruited from South London Child and Adolescent Mental Health Services.ADHD diagnosis was confirmed by an experienced child psychiatrist using the standardised Maudsley Diagnostic Interview (DSM-IV-TR criteria; Goldberg and Murray 2002).Boys with ADHD scored above the clinical threshold for ADHD symptoms on the strengths and difficulties questionnaire for parents (SDQ; Goodman and Scott 1999) and the Conners parent rating scale (CPRS-R; Conners 2008), and below the clinical threshold on the social communication questionnaire (SCQ) to exclude participants with high autism traits (Rutter et al. 2003).Where SCQ scores were inconclusive or not provided, the child's clinician was consulted to rule out autism comorbidity.Neurotypical boys were recruited through advertisement in the same South London area.They scored below the clinical cut-off on the SDQ, SCQ, and CPRS-R (Table 1).All participants were right-handed.
Exclusion criteria for all participants were MRI-related contraindications, mean framewise displacement during scanning >1mm, IQ <70 on the Wechsler abbreviated scale of intelligence (WASI; Wechsler 1999), history of substance abuse, neurological deficits, presence of psychiatric disorders (except for conduct disorder/oppositional defiant disorder in the ADHD group, N = 2), learning disability, reading, speech, or language disorder.Data for additional 11 participants were available (N ADHD = 1, N control = 10) but were not included in the current study due to excessive motion (N ADHD = 1), poor functional data normalisation (N control = 1), and control participant's age exceeding the age range of patients (N control = 9).
Welch's t-tests showed no significant group difference for age but a significant difference for IQ with ADHD participants scoring lower than controls, which is typical in this population (Rommel et al. 2015).Moreover, ADHD participants scored significantly higher on the SDQ and SCQ (Table 1).Missing data points (N SCQ = 9, N SDQ = 2) were mean imputed.
Participants were reimbursed £50 for each session.Parental/child written informed consent/assent were obtained.Ethical approval was granted by the Joint South London and Maudsley/Institute of Psychiatry, Psychology & Neuroscience Research Ethics Committee (07/H0807/84).2008).Both drugs and placebo were administered in identical capsules 1.5 h before scanning to allow maximum absorption (Chan et al. 1983;Witcher et al. 2003) and to preserve blinding.Similar dosages and time lapses between drug administration and scans have shown to be sufficient to observe changes in brain activation and performance in ADHD (Cubillo et al. 2014a(Cubillo et al. , 2014b;;Kowalczyk et al. 2019;Rubia et al. 2011aRubia et al. , 2011b;;Smith et al. 2013).Controls were scanned once, unmedicated.

N-back paradigm
The 6-min block design working memory task consisted of four conditions (Ginestet and Simmons 2011;Cubillo et al. 2014a).This parametric task with increasing working memory load was chosen to avoid floor or ceiling effects in task performance (Pongpipat et al. 2021).During 1-back, 2-back, and 3-back conditions, participants were presented with series of letters (1s duration, 2s inter-trial interval) and responded with their right thumb using a button box whenever the letter shown was the same as the letter presented one, two, or three before it, respectively (e.g.2-back: B/J/A/J).This requires both storage and continuous updating of stimuli held in working memory (Rac-Lubashevsky and Kessler 2016).In the baseline vigilance 0-back condition, participants responded to each 'X' that appeared on the screen (Fig. 1).The task consisted of 180 trials presented in 12 blocks.Before each block, written instructions (2s) indicated which condition would be shown.Each block comprised one condition lasting 30s and consisting of 14 stimuli presentations: 3 targets and 11 non-targets.Each condition was presented three times.Performance data were recorded during scanning.
Participants practised the task once before scanning.

MRI acquisition
Data were acquired on a GE Signa HDx 3T system (General Electric, Milwaukee, Wisconsin) with an 8-channel head coil at the Centre for Neuroimaging Sciences, King's College London.Functional T2*-weighted data comprised 39 interleaved slices acquired with a gradient EPI sequence bottom-to-top over 186 volumes depicting BOLD contrast covering the whole brain (TE = 30ms, TR = 2s, flip angle = 75°, in-plane voxel size = 3.75mm, slice thickness = 3.5mm, slice gap = 0.5mm).Wholehead structural image was acquired at the beginning of each scanning session using an MPRAGE protocol using parameters based on ADNI (TE = 2.85ms, TR = 6.99s, flip angle = 8°, in-plane voxel size = 1.02mm, 166 slices, slice thickness = 1.2mm, slice gap = 1.2mm).

Data preprocessing
Data were processed using SPM12 (v7487; https:// fil.ion.ucl.ac.uk/ spm; Penny et al. 2007) and MATLAB (v9.5.0;MATLAB 2018).The origin of all images was reset to the anterior commissure.Each participant's functional data were slice-time corrected, realigned, and co-registered with the participant's anatomical image, and smoothed at 6mm FWHM Gaussian kernel.Anatomical images were segmented into grey matter, white matter, cerebrospinal fluid, and non-brain components using unified segmentation.The resultant unsampled grey and white matter images were used to create a study-specific template using DARTEL (Ashburner 2007).The generated deformation flow fields were used in functional data normalisation to the MNI template.Visual quality checks were performed at all stages.

Seed region selection
Three seed regions were chosen for functional connectivity analysis based on their consistent responses to working memory paradigms in paediatric and adult meta-analyses (Owen et al. 2005 (rDLPFC;43, 31, 30;Owen et al. 2005).The seeds were created as 5mm radius spheres using MarsBar (Brett et al. 2002).

Generalised psychophysiological interaction
Functional connectivity during the n-back task was examined using generalised psychophysiological interaction (gPPI) method (McLaren et al. 2012).Separate gPPI models were created for each seed region.For each model, the first eigenvariate of the seed's time-series was extracted for each participant and deconvolved to obtain an estimate of the physiological activity within that region.Subsequently, five PPI regressors were created by multiplying the estimate of the seed's physiological activity by vectors representing each of the task conditions (instruction, 0-back, 1-back, 2-back, and 3-back).PPI regressors were convolved with a canonical haemodynamic response function and entered into a GLM together with a regressor corresponding to the seed's timeseries, 24 motion parameters (Friston et al. 1996), scrubbing regressors for volumes where framewise displacement (FD) > 0.5mm (Power et al. 2012(Power et al. , 2014)), and two nuisance regressors modelling mean white matter and cerebrospinal fluid signal.The 3-back condition was not included in contrasts, as participants in all groups were observed to make errors on >30% of trials.The contrast of interest explored the interaction of the seed with high working memory load condition compared to the baseline vigilance condition (2-back > 0-back).This contrast from first-level models was taken forward to group-level random-effects analysis.
To check the gPPI modelling quality, three one-sample t-tests of the interaction between the high working memory load contrast (2-back > 0-back) and each of the seeds were conducted in the control group (Supplementary Materials).
For hypothesis-testing, independent samples t-tests with age in months as a covariate were used to compare participants with ADHD under each drug condition with the neurotypical control group.Paired samples t-tests explored drug effects within the ADHD cohort.Cluster-based FWE correction (p < 0.05) with a cluster-forming threshold of p uncorrected < 0.001 was applied to account for multiple comparisons.

Conjunction analyses for the assessment of normalisation and upregulation effects
The goal of treatments in psychiatry is to minimise dysfunction observed in patients, potentially through normalisation of function.An absence of a significant difference, however, cannot be considered as evidence of normalisation (Ranganathan et al. 2015).A conjunction analysis between a comparison of neurotypical controls relative to patients under placebo and patients under a drug relative to patients under placebo would indicate regions where shared activation or connectivity in patients under a drug and neurotypical controls suggests a drug-related normalisation.Consequently, to test for potential normalisation effects of drugs on ADHDrelated neurofunctional abnormalities and any drug-related compensatory changes, we used confirmatory conjunction analyses.

Sensitivity analysis
Sensitivity analyses were conducted on performance and fMRI data excluding participants with comorbid oppositional defiant and/or conduct disorder (N = 2).

Scan order effects
Repeated measures ANOVAs were used to test for scan order effects in patients using accuracy, MRT, and SDRT.

Task performance
Descriptive statistics for performance are reported in Table 2.

Within-patients comparisons
Repeated measures ANOVAs revealed a significant effect of working memory load on accuracy, MRT, and SDRT.Increasing working memory load was associated with lower accuracy, along with slower and more variable responses across drug conditions.No effects of drug condition or interaction effects between drug condition and working memory load were observed (Table 2).

Case-control comparisons
Across all participants, mixed-measures ANOVAs showed that increasing working memory load was associated with reduced accuracy, longer MRT, and larger SDRT.Significant interaction effects between working memory load and group were observed for MRT when controls were compared to patients under placebo, due to patients responding slower during 0-back and 2-back, and faster during 1-back and 3-back; however, post hoc tests revealed these effects were statistically non-significant (0-back, p = 0.568; 1-back, p = 0.917; 2-back, p = 0.506; 3-back, p = 0.106).No other statistically significant effects were observed (Table 2).Covarying for age did not affect the results.

Within-patients comparisons
Paired samples t-tests showed no significant differences for any seeds in within-patients comparisons.

Case-control comparisons
Under placebo, independent samples t-test showed that participants with ADHD exhibited greater connectivity between rSPG and left-hemispheric insula/central operculum relative to controls.This effect was due to more negative beta values during 2-back compared to 0-back in controls, but no such difference in participants with ADHD under placebo (Fig. 2, Table 3).No differences between controls and participants with ADHD under placebo were observed for connectivity with rMFG or rDLPFC.
Under methylphenidate, participants with ADHD showed greater connectivity than controls between rMFG and (1) bilateral: insula, ACC, Heschl's gyrus, putamen, supramarginal gyri (SMG), parietal operculum; (2) left: postcentral gyrus and planum temporale/operculum; (3) right: thalamus, paracingulate gyrus, and temporal occipital fusiform cortex.Additionally, greater connectivity in participants with ADHD under methylphenidate between rSPG and bilateral ACC, as well as right parietal operculum/SMG was found compared to controls.Mean beta values indicated that neurotypical controls exhibited lower beta values during 2-back than 0-back.In contrast, participants with ADHD under methylphenidate showed the opposite effect, with higher beta values during 2-back compared to 0-back (Fig. 2, Table 3).No differences, however, were observed for rDLPFC connectivity between participants with ADHD under methylphenidate and controls.
Comparisons of participants with ADHD under atomoxetine and controls revealed no significant differences for connectivity with any seed.
To test for potential associations between connectivity and performance, 2-back accuracy (main index of performance in this task) was entered into group-level models showing case-control differences as a covariate of interest.Small volume correction was applied to restrict the test to regions showing significant differences between individuals with ADHD and controls.No statistically significant associations between accuracy and functional connectivity were observed.

Conjunction analyses for the assessment of normalisation and upregulation effects
To assess whether the lack of significant case-control differences in connectivity between rSPG and left-hemispheric insula/central operculum (observed under placebo) under methylphenidate and atomoxetine reflected drug normalisation effects, we conducted confirmatory conjunction analyses.Given that the within-patient comparisons of methylphenidate/atomoxetine and placebo did not reveal any statistically significant differences across the whole brain, here we restrict our search only to voxels within the left insula/central operculum that showed abnormally high connectivity with rSPG under placebo.In the absence of a separate ADHD group, we acknowledge the potential of bias of this approach.Two conjunction analyses were performed for the comparisons of (i) controls > individuals with ADHD under placebo and individuals with ADHD under methylphenidate > individuals with ADHD under placebo and (ii) controls > individuals with ADHD under placebo and individuals with ADHD under atomoxetine > individuals with ADHD under placebo.Small volume correction using a mask defined based on the insula/central operculum cluster of case-control differences under placebo was applied to both analyses.The conjunction analyses revealed no common areas of connectivity within this cluster in individuals with ADHD under either drug and controls, thereby failing to support full normalisation of rSPG-insula/operculum hyperconnectivity.
Additional conjunction analyses were conducted to test whether regions showing heightened connectivity in ADHD under methylphenidate relative to controls during 2-back > 0-back were also upregulated with methylphenidate relative to placebo within patients.Two conjunction analyses were conducted, one for connectivity with rMFG and another one for the connectivity with rSPG.These analyses included individuals with ADHD under methylphenidate > individuals with ADHD under placebo and individuals with ADHD under methylphenidate > controls.Small volume correction using masks of regions showing case-control differences when ADHD individuals were under methylphenidate was applied.Conjunction analyses revealed common upregulation with methylphenidate in ADHD relative to placebo and relative to neurotypical controls in connectivity between rMFG and SMG (cluster size = 12, peak-level Z = 3.35, p FWE = 0.026, MNI coordinates: 51, −36, 34) and rSPG and SMG (cluster size = 16, peak-level Z = 3.65, p FWE = 0.026, MNI coordinates: 52, −36, 36).Beta values revealed that while the change between vigilance and high working memory load was similar under both placebo and methylphenidate, under methylphenidate those regions had more negative values across task conditions.

Sensitivity analyses
Sensitivity analyses revealed no impact of oppositional defiant and/or conduct disorders on results of performance and fMRI analyses.

Discussion
This study investigated the effects of single doses of methylphenidate and atomoxetine relative to placebo on working memory-related functional connectivity in boys with ADHD and compared to neurotypical controls.Relative to controls, individuals with ADHD under placebo exhibited increased connectivity between rSPG and left central operculum/insula.This cluster was not observed when individuals with ADHD under methylphenidate or atomoxetine were compared to controls.Under methylphenidate, youth with ADHD compared to neurotypical controls showed increased connectivity between rMFG and bilateral insula, ACC, putamen, right thalamus, and other bilateral parieto-temporal regions, as well as between rSPG and predominantly righthemispheric ACC and parietal areas.No differences, however, were observed within patients across drug conditions or between ADHD participants under atomoxetine and controls.Participants with ADHD under each drug condition showed comparable performance to controls and no effects of drug condition on performance were observed within patients.
The hyperconnectivity in youth with ADHD relative to controls between rSPG and left central operculum/insula was no longer observed when individuals with ADHD were under either methylphenidate or atomoxetine.Combined with the lack of difference when comparing drug to placebo directly, the conjunction analysis supports a conclusion that neither medication was associated with Fig. 2 Functional connectivity during 2-back > 0-back and graphs showing mean beta values within each statistically significant cluster during 0-back and 2-back in neurotypical controls and participants with ADHD. A. Regions showing significantly higher functional connectivity with the rSPG in participants with ADHD under placebo compared to neurotypical controls.B. Regions of greater functional connectivity with rSPG (red) and rMFG (yellow) in the ADHD group under methylphenidate compared to controls.Axial slices are marked with the z coordinate.rMFG, right middle frontal gyrus; MPH, methylphenidate; rSPG, right superior parietal gyrus ◂ complete functional normalisation.Furthermore, under methylphenidate, additional widespread connectivity increases were observed relative to controls between rMFG and fronto-striato-temporo-parietal regions, and between rSPG and cingulo-parietal areas.Additionally, conjunction analyses revealed that methylphenidate-related upregulation of rMFG-SMG and rSPG-SMG connectivity was observed relative to both placebo and controls.Nevertheless, the lack of differences between methylphenidate and placebo sessions across the whole brain within patients complicates the interpretation of this finding.Although our study did not observe as widespread changes in working memory-related connectivity between placebo and methylphenidate sessions as past research (Wong and Stevens 2012;Wu et al. 2017), we report for the first time methylphenidate-induced increases of working memoryrelated connectivity in ADHD relative to neurotypical controls.These changes were due to controls showing a negative shift in coupling between rMFG/rSPG and functionally connected regions when moving from less to more demanding conditions, while participants with ADHD under methylphenidate showing an opposite positive shift in network engagement between less and more demanding task conditions.Additionally, SMG regions showing shared upregulation of connectivity under methylphenidate relative to controls and placebo revealed similar loaddependent changes under placebo and under methylphenidate; however, methylphenidate was associated with more negative connectivity across task conditions.The regions showing hyperconnectivity under methylphenidate relative to controls overlap with fronto-parietal networks comprising ACC, middle frontal, postcentral, supramarginal, and lingual gyri shown to be upregulated with stimulants relative to placebo in previous working memory studies in regularly treated youth with ADHD (Wong and Stevens 2012).Meta-regression analysis across fMRI studies of working memory has also reported that stimulant use was associated with increased middle frontal activation (McCarthy et al. 2014).Similarly, greater activation of insula, putamen, and ACC was associated with stimulant treatment in a meta-analysis of various cognitive paradigms (Rubia et al. 2014).Crucially, our results complement the previously published activation analysis of this dataset (Cubillo et al. 2014a), by demonstrating that methylphenidate is associated not only with increases in middle frontal, middle/superior temporal, striatal, and thalamic activation relative to controls but also with increases in the interconnectivity across these regions.
Similar to methylphenidate, atomoxetine led to a partial reduction of rSPG-operculum/insula hyperconnectivity observed in ADHD under placebo relative to controls, while not showing full normalisation or within-patients effects relative to placebo.In contrast to methylphenidate, which enhanced connectivity in several fronto-parieto-temporal and fronto-striato-thalamic regions, atomoxetine did not induce any changes relative to controls.This may be due to the single-dose design of this study.Although single doses of atomoxetine have been shown to modulate brain activation during working memory in this cohort (Cubillo et al. 2014a) and during other cognitive functions (Smith et al. 2013;Cubillo et al. 2014b;Kowalczyk et al. 2019), longer term administration is typically needed for full clinical benefits (Montoya et al. 2009), while methylphenidate shows clinical effects faster (Greenhill et al. 2001).Consequently, investigations of chronic atomoxetine treatment may be more suitable in better understanding atomoxetine's effects on neural networks in ADHD.
Our observation of abnormally increased rSPG-operculum/insula connectivity in ADHD under placebo compared to controls is consistent with past reports of hyperconnectivity in fronto-parietal networks including SPG in youth with ADHD (Massat et al. 2012;Bédard et al. 2014;Wu et al. 2017).Fronto-parietal networks are crucial to supporting maintenance and storage of information in working memory (Ekman et al. 2016) and increases in connectivity within them are associated with task engagement (Pongpipat et al. 2021) and improved task performance in neurotypical populations (Shen et al. 2015).Our finding of hyperconnectivity between rSPG and central operculum/insula in ADHD under placebo in the absence of working memory deficits may reflect a compensatory mechanism for the maintenance of working memory function, particularly considering consistent reports of widespread reductions of fronto-striatal and temporo-parietal activation in ADHD during working memory (Silk et al. 2005;Vance et al. 2007;Kobel et al. 2009;Cortese et al. 2012;Cubillo et al. 2014a;McCarthy et al. 2014;Chantiluke et al. 2015), including in the same cohort (Cubillo et al. 2014a), and during other executive functions (Cortese et al. 2012;Hart et al. 2013;Norman et al. 2016;Rubia 2018;Lukito et al. 2020).The differences were due to controls showing a negative shift in rSPG connectivity when switching between baseline vigilance and high working memory load, while patients showed similar network engagement in both conditions.The decreases of connectivity with increasing task load in neurotypical individuals were unexpected given the frequent reports of loaddependent increases of activation and connectivity in taskrelevant networks (O'Hare et al. 2008;van den Bosch et al. 2014;Vogan et al. 2016;Le et al. 2020) and given an association between higher task load and greater activation in task-relevant regions in our control sample (Supplementary Materials; Cubillo et al. 2014a).This may be due to ongoing maturation processes occurring in adolescence.In fact, greater load-dependent changes are present in adults than children (Vogan et al. 2016) and increasing age has been associated with decreased working memory-related engagement of superior frontal, postcentral, inferior parietal, and cingulate regions in adolescence (Andre et al. 2015).Consequently, the reductions of connectivity with increasing load in our control group may reflect the progressive specialisation of higher order regions through weakened connections of diffuse networks.Conversely, individuals with ADHD under placebo did not show load-dependent changes, potentially indicating less efficient adjustment to changing task demands in ADHD.This could have been tested by investigating the 3-back condition; however, we opted against it due to the large error rate confounding this comparison and increasing variance within blocks, thus reducing sensitivity to group differences.Furthermore, greater parietal connectivity in ADHD relative to neurotypical peers may relate to functional maturation delay of working memory networks (Andre et al. 2015) given evidence that younger children exhibit less load-dependent modulation of task-relevant networks (van den Bosch et al. 2014).This could suggest that connectivity differences between age-matched neurotypical youth and youth with ADHD are similar to a comparison between older and younger children.This would be in line with resting state connectivity studies showing maturation lag in the connections between fronto-parietal and attention networks and the default mode network in ADHD (Sripada et al. 2014).Finally, some of these findings may also be attributed to sex differences, given that our sample comprised only boys and males with ADHD tend to exhibit more pronounced differences from their neurotypical peers than females with ADHD (Valera et al. 2010;Dupont et al. 2022).
Lack of drug effects, case-control differences, or main effects of working memory load within controls on rDLPFC connectivity were unexpected.However, while this region was chosen due to its strong involvement in verbal n-back tasks, this seed was derived from a meta-analysis in adults (Owen et al. 2005).The functional maturation processes occurring in childhood and adolescence involve progressive shifting from diffuse to localised and specialised networks (Casey et al. 2005;Rubia 2013;Andre et al. 2015;Stevens 2016).Given that DLPFC is among the last regions to mature (Casey et al. 2005;Rubia 2013), lack of drug or working memory load effects on functional connectivity with rDLPFC might reflect the young age of our cohort.Although rDLPFC showed task-relevant activation during high working memory load relative to vigilance in our control group, it may be that the functional connections between rDLPFC and other task-relevant regions have not yet fully matured, and thus did not show the expected task-dependent modulation or case-control differences and drug effects.
Lack of methylphenidate and atomoxetine effects on task performance may be due to the absence of working memory deficits observed in our ADHD cohort.Although there is meta-analytic evidence for working memory performance deficits in ADHD, they seem larger in visuospatial than verbal domain (Martinussen et al. 2005;Willcutt et al. 2005;Kasper et al. 2012).Selection of a verbal n-back task was motivated by its wide use in neuroimaging (Owen et al. 2005; Yaple and Arsalidou 2018), facilitating comparisons with previous studies.Furthermore, neurofunctional differences between individuals with ADHD and controls, despite similar performance, have been seen in previous studies using n-back (Massat et al. 2012;Cubillo et al. 2014a;Wu et al. 2017) and other verbal working memory tasks (Hale et al. 2007;Wolf et al. 2009).While the interpretation of functional implications of drug-related modulation of brain networks in these cases is not straightforward, the lack of performance differences ensures that the observed neurofunctional changes are task-specific rather than a reflection of error monitoring.Nonetheless, future research should consider using visuospatial working memory tasks in investigations of medication effects given evidence that stimulants lead to greater improvements in the visuospatial domain (Bédard and Tannock 2008;Bédard et al. 2014).Furthermore, we cannot rule out that lack of performance differences was due to low sample size, given that this study was powered for fMRI analyses (Thirion et al. 2007).
This study is strengthened by a double-blind, placebocontrolled design and inclusion of medication-naïve individuals with ADHD.Considering previous reports of associations between long-term stimulant treatment and structural (Nakao et al. 2011;Frodl and Skokauskas 2012;Lukito et al. 2020), neurochemical (Fusar-Poli et al. 2012), and neurofunctional changes in ADHD (Konrad et al. 2007;Hart et al. 2013;McCarthy et al. 2014;Norman et al. 2016;Lukito et al. 2020), investigating drug effects without the confounds of previous pharmacotherapy is critical.
This study is not without limitations.Our participants were unmedicated at the time of the study randomisation potentially because many of them were newly diagnosed with ADHD and have not yet begun medication treatment.While recruitment of a medication-naïve cohort protected against the confounding impact of chronic pharmacological treatment, it may have biased our sample towards those with less severe symptoms (Hong et al. 2014).Nonetheless, this is unlikely given that the scores of our ADHD cohort on the SDQ questionnaire were above the clinical threshold for ADHD and similar to those typically reported in the ADHD literature (Becker et al. 2006;Hall et al. 2019).While the single-dose design avoided the influence of symptomatic improvement, side effects, or neural changes associated with long-term treatment, it likely favoured detection of methylphenidate effects.Methylphenidate offers immediate clinical benefits (Greenhill et al. 2001), while atomoxetine requires longer term administration (Montoya et al. 2009).Furthermore, although the within-patient design minimised the influence of between-subject variability on drug effects, it meant that patients were only truly medication-naïve on their first visit.There is, however, no evidence suggesting that a single dose of methylphenidate or atomoxetine can lead to lasting neural changes; furthermore, to minimise any carry-over effects, study visits were scheduled one week apart, corresponding to more than five half-lives of each medication (Dhariwal and Jackson 2003).The within-patient design also meant that patients completed the task three times, while controls only once, for ethical and financial reasons.Although the n-back task is associated with improved performance after repeated practice (Pergher et al. 2018;Chen et al. 2019), the lack of scan order effects in this study indicates it is unlikely the current findings are confounded by practice effects.Furthermore, conjunction analyses could have benefited either from two independent cohorts, to first define case-control differences and then separately test drug effects, or from a baseline scan for patients.Such data, however, was not available, and thus the placebo condition appears in both contrasts.Therefore, this analysis is not completely without bias.Consequently, the reliability of these findings should be investigated in a well-powered replication study.The greater prevalence of ADHD in boys meant that this study recruited a fully male sample.Consequently, given the sexual dimorphism of ADHD and evidence of distinct neurofunctional profiles in males and females with ADHD, the generalisability of the current findings is limited (Valera et al. 2010;Dupont et al. 2022).Finally, to maximise cohort homogeneity, we only recruited right-handed individuals with combined ADHD, meaning these findings may not generalise to other ADHD populations, such as adults, lefthanded individuals, or those with inattentive or hyperactiveimpulsive presentations.
Overall, this is the first study of comparative effects of methylphenidate and atomoxetine on working memoryrelated functional connectivity in medication-naïve boys with ADHD.We showed that while both drugs shared incomplete normalising effect of abnormally enhanced connectivity in parieto-insular regions, only methylphenidate, but not atomoxetine, led to widespread increases in connectivity within task-relevant networks in ADHD relative to controls.This evidence extends past research showing upregulating effects of methylphenidate on functional connectivity during working memory in ADHD (Wong and Stevens 2012;Wu et al. 2017), suggesting they are exclusive to methylphenidate.Furthermore, in the context of past research exploring the comparative effects of single-dose methylphenidate and atomoxetine (Cubillo et al. 2014a), this study suggests that while both drugs show shared upregulation of working memory-relevant activation and a downmodulatory effect on functional connectivity abnormalities present under placebo, indicative of network reorganisation, only methylphenidate leads to widespread connectivity increases in fronto-temporo-parietal and fronto-striatothalamic regions.

Fig. 1
Fig. 1 A schematic representation of the n-back block design task. A. A visual representation of the task block structure.B. Examples of n-back conditions of different working memory load (0-, 1-, 2-, and 3-back).ITI, inter-trial interval

Table 2
Performance measures and analyses for the n-back task for neurotypical controls and patients with ADHD under each drug condition

Table 3
Significant clusters of differences between neurotypical controls and participants with ADHD under each drug condition and their peak connectivity during the 2-back > 0-back contrast ADHD, Attention-Deficit/Hyperactivity Disorder; FWE, Family-Wise Error; L, Left; rMFG, Right Middle Frontal Gyrus; MNI, Montreal Neurological Institute; R, Right; rSPG, Right Superior Parietal Gyrus