The neural correlates of visuo-spatial working memory in children with autism spectrum disorder: effects of cognitive load

Background Research on the neural bases of cognitive deficits in autism spectrum disorder (ASD) has shown that working memory (WM) difficulties are associated with abnormalities in the prefrontal cortex. However, cognitive load impacts these findings, and no studies have examined the relation between WM load and neural underpinnings in children with ASD. Thus, the current study determined the effects of cognitive load on WM, using a visuo-spatial WM capacity task in children with and without ASD with functional magnetic resonance imaging (fMRI). Methods We used fMRI and a 1-back colour matching task (CMT) task with four levels of difficulty to compare the cortical activation patterns associated with WM in children (7–13 years old) with high functioning autism (N = 19) and matched controls (N = 17) across cognitive load. Results Performance on CMT was comparable between groups, with the exception of one difficulty level. Using linear trend analyses, the control group showed increasing activation as a function of difficulty level in frontal and parietal lobes, particularly between the highest difficulty levels, and decreasing activation as a function of difficulty level in the posterior cingulate and medial frontal gyri. In contrast, children with ASD showed increasing activation only in posterior brain regions and decreasing activation in the posterior cingulate and medial frontal gyri, as a function of difficulty level. Significant differences were found in the precuneus, dorsolateral prefrontal cortex and medial premotor cortex, where control children showed greater positive linear relations between cortical activity and task difficulty level, particularly at the highest difficulty levels, but children with ASD did not show these trends. Conclusions Children with ASD showed differences in activation in the frontal and parietal lobes—both critical substrates for visuo-spatial WM. Our data suggest that children with ASD rely mainly on posterior brain regions associated with visual and lower level processing, whereas controls showed activity in frontal lobes related to the classic WM network. Findings will help guide future work by localizing areas of vulnerability to developmental disturbances.

Previous neuroimaging, electrophysiology and neurochemical studies in ASD have identified atypical white and grey matter volumes [15][16][17][18], functional connectivity [10,19], cortical sulcal and gyral anatomy [20], brain lateralization [21], neural perfusion [22] and serotonin synthesis capacity [23] compared to those of typically developing (TD) individuals, with the frontal cortex implicated in a number of these differences [24]. Although there is substantial evidence for developmental anatomical abnormalities of the frontal lobes in ASD, associations with cognitive performance are lacking [25]. It remains to be determined if executive dysfunction in ASD is related to functional as well as anatomical frontal lobe aberrations. Further, protracted maturation of the frontal lobes makes the functions they support, such as WM, susceptible to developmental disturbances [26,27] but amenable to therapeutic interventions. WM allows for the temporary storage, rehearsal and maintenance of information. It is important for learning, social ability [28], academic achievement [29] and many complex cognitive operations [30,31]. Impairments in WM have been reported in ASD, yet our understanding of its development and neural correlates is still very limited in this population. We investigated the neural systems underlying visuo-spatial WM capacity using functional magnetic resonance imaging (fMRI) to determine the role of cognitive load on possible functional differences between children with and without ASD. Cognitive load is a multidimensional construct referring to the processing resources that performing a particular task imposes on one's cognitive system [32]-in the present study, the level of WM task difficulty was used as a measure of cognitive load.
Extant research on the behavioural characterization of WM function in ASD suggests both intact [33][34][35][36][37] WM performance on simple memory tasks, and impaired [33][34][35][36]38] performance on more complex tasks, relative to TD individuals. Discrepancies in the literature may be due partly to methodological inconsistencies regarding task choice, comparison groups or population age. Basic WM abilities are intact in high-functioning ASD (see [39] for review). The majority of studies that found significant group differences assessed WM using tasks with increased complexity and/or cognitive load and thereby imposed heavier demands on WM and executive functions [9,34,38]. Further, there is evidence of primarily visuo-spatial WM impairment, whereas verbal WM appears relatively intact in individuals with ASD [35][36][37].
Only a few studies have examined the neural correlates of WM function in ASD, and no neuroimaging studies of pre-adolescent children exist. Overall, WM processes are largely subserved by the prefrontal and parietal cortices [4,[40][41][42][43][44]. Converging literature identifies a broad system of prefrontal, premotor, dorsal cingulate and posterior parietal activation in visual WM tasks see [44], and neuroimaging studies in ASD provide evidence for atypical activity in these regions [10,12,14]. Using an oculomotor visuo-spatial WM task, Luna and colleagues [12] found behavioural impairments in WM, as well as reduced activation in the dorsolateral prefrontal cortex (DLPFC) and posterior cingulate cortex in adults with ASD compared to controls. Interestingly, there was no evidence of impaired activation in other areas known to support WM neural circuitry, such as the anterior cingulate cortex. A study by Koshino et al. [10] demonstrated bilateral activation in the DLPFC in adult controls, whereas adults with ASD showed limited activation in the left, and more right hemisphere recruitment of prefrontal regions, despite the absence of differences in performance accuracy on an n-back letter task. The authors suggested that while typical adults processed letter stimuli using verbal codes, those with ASD employed visual strategies [10], supporting the idea that WM deficits in ASD may be attributed to less efficient processing strategies. In the only study of adolescents, Silk et al. [14] used a visuo-spatial WM mental rotation task, where behavioural performance was similar for ASD and control groups, and observed impaired cortical activation in the frontal lobes in the ASD group, including the anterior cingulate, DLPFC and caudate nucleus, but normal activation in the parietal cortices relative to controls. These findings suggest dysfunctional frontostriatal networks in ASD. More work is needed to understand the neural correlates and developmental trajectory of WM in childhood in ASD, as limited research is available reporting on pre-adolescent children.
One of the most common experimental paradigms used to manipulate cognitive load in studying WM is the 'n-back' protocol [10,41,[43][44][45][46][47]. In a typical n-back task, participants view a series of stimuli and indicate whether the currently presented stimulus matches one presented 'n' (e.g., 0, 1, 2 or 3) trials prior. As difficulty level increases, the number of interfering stimuli between the target and relevant stimulus increases, requiring the utilization of different mental strategies at each level (e.g. 0-back, recognition; 1-back, maintenance; 2-back, maintenance and monitoring). This manipulation consequently increases both memory load and executive function demand (i.e. strategy needed to solve the task) in a non-linear fashion from one level to the next, making function-specific changes difficult to quantify and link with brain areas. Therefore, we used a 1-back colour matching task (CMT) [48] which systematically manipulated memory load while keeping executive function constant across all difficulty levels, allowing a direct investigation of the influence of cognitive load on WM. Specifically, executive schemes (i.e. procedural strategies for solving the task) are constant across levels of items in CMT; what varies with each level is the number of relevant items (colours) to be remembered. Difficulty level was parametrically graded based on behavioural agedependent growth patterns observed using this task in previous work [48]. These observations point to a linear pattern of WM development, which may also be evident in the neural processes across increasing load. Our task is novel, and it captures the brain regions associated with this linear pattern of function across cognitive load in preadolescent children. In a previous study from our group, typical adults showed positive linear relations between cortical activity and CMT task difficulty level in areas involved in WM function [49]. Negative linear relations were found in areas typically associated with the default mode network (DMN). The DMN, which has been found in a wide range of neuroimaging studies, is a network of brain regions, including the medial prefrontal cortex, posterior cingulate and inferior parietal lobules, characterized by decreased activation during goal-oriented or attentiondemanding tasks [50]. Neurodevelopmental disorders, including ASD, have been associated with abnormal function [51,52] and structure [53] of the DMN that may interfere with cognitive function.
The present study identified and compared the neural activity underlying visuo-spatial WM capacity in children with and without ASD using fMRI, with manipulation of cognitive load. Given the increasing difficulty levels of the current task, we hypothesized that children with ASD would show poorer performance than matched TD controls on behavioural measures of WM. Further, frontal and posterior parietal cortical areas related to visuo-spatial WM capacity would be under-recruited in children with ASD relative to TD controls, and this difference would increase with load. We expected activity to be linearly modulated (positively for WM areas, negatively in default-mode areas) by task difficulty in all children (comparable to previous work with adults), but these trends would be less marked in children with ASD, particularly at higher cognitive loads.
With the protracted maturation of the frontal lobes, and their susceptibility to developmental anomalies, understanding development in these regions is crucial, particularly in populations with frontal lobe abnormalities, as seen in ASD. As previous research has focused on studying WM in adults and adolescents with ASD, there is a critical need to study WM in pre-adolescents with ASD, who exhibit both cognitive and neurological differences. Given the links between social function, school success and executive function ability, investigating the neural bases of WM deficits in children with ASD will contribute to our knowledge of the underlying causes of ASD-related behaviour. Previous research highlights the importance of investigating cognitive impairments in ASD that may arise from brain abnormalities and drive behavioural symptoms [54]. Further, exploring WM correlates will supplement behavioural phenotypes of ASD. This work will help identify the nature of atypical brain development, with future expectations of establishing age-appropriate interventions that can effectively target WM function and, in turn, other symptoms of ASD.
Exclusion criteria for all participants were the presence of any current significant axis I psychiatric comorbidities [2], neurological disorders, medical illnesses, prematurity, uncorrected vision, colour blindness, IQ < 70 as well as standard MRI contraindicators (e.g. ferromagnetic implants). A history of developmental delay, learning disability and attention deficit hyperactivity disorder (ADHD) was used to exclude TD children only; however, these factors were also not current primary diagnoses in the ASD group. Six children with ASD were each on one psychotropic medication (Strattera, Biphentin, Fluoxetine, Concerta, Abilify and Atomoxetine). Their fMRI data were examined in comparison to children with ASD who were not taking medication, and the data did not differ between these subgroups (see Additional file 1). Children were recruited through community support centres, parent support groups, email lists, hospital ads and private schools. Informed consent, clinical and cognitive testing and MRI scanning were performed at the Hospital for Sick Children in Toronto. Experimental procedures were approved by the Research Ethics Board at the Hospital for Sick Children. All children gave informed assent, and the parents provided informed written consent.
Clinical diagnoses of ASD were confirmed in all cases with a combination of expert clinical judgement, clinical records and the Autism Diagnostic Observation Schedule (ADOS) [56], which was administered by a trained individual who maintains inter-rater research reliability. All children completed the backwards digit recall, listening recall, digit recall, mazes memory and block recall subscales of the Working Memory Test Battery for Children (WMTB-C) [57] to supplement behavioural data collected during fMRI tasks. See Table 1 for demographic, neuropsychological and clinical test characteristics.

The colour matching task
Children were required to attend to coloured figures of a clown presented one at a time in sequence. Children were taught to ignore the clown's face and irrelevant colours (blue and green) and focus on other, relevant colours (yellow, purple, pink, red, orange, brown and grey). Tasks that contain misleading or irrelevant factors that evoke interference have been shown to be a more suitable measure of WM capacity, yielding reliable estimates of developmental growth in adults [48] and children (Powell, Arsalidou, Vogan, and Taylor: Controlled inference and assessments of developmental working memory capacity: Evidence from letter and colour matching tasks, submitted). The number of 'n' relevant colours (capacity) in the figures was increased by one for each increase in difficulty level. CMT has two integral functions that require mental attention, in which participants must first actively extract the relevant colours embedded in the clown figure and second check for a possible match with colours of the criterion set. As such, items with n (e.g. 1) relevant colours will have difficulty level of n + 2 (e.g. 3; Figure 1A). Children indicated after each item whether the relevant colours of the current figure matched those from the immediately preceding figure (i.e. 1-back), disregarding colour location and repetition. Using a keypad with the right hand, children responded in the scanner by pushing a button for 'same' when the clown figure was wearing the same colours as the previous clown and 'different' when the clown figure was wearing different colours. All children successfully completed practice trials on a computer outside the scanner with accuracy of 80% or greater. A total of 24 task blocks (168 task trials) and 24 baseline blocks were presented across four runs. Each run consisted of six 32-s blocks, one for each difficulty level (six levels in total); within each block, there was a constant difficulty level, and all difficulty levels were randomized within each run ( Figure 1B). The same four runs were presented to all children in the same order. Task blocks consisted of eight stimuli of the same difficulty and alternated with 20-s baseline blocks ( Figure 1C), where clowns were coloured only in blue and green, and children were instructed to look at the figures but not respond. Participants had 3 s to view a stimulus and respond, followed by a 1-s inter-stimulus interval where a fixation cross was presented ( Figure 1D). The fMRI task took approximately 22 min of scan time.
Performance data were recorded for both accuracy (proportion correct) and reaction time; items were correct if responded correctly within 3 s of stimulus presentation. To ensure adequate task completion, children were excluded if they did not achieve at least 60% accuracy (averaged across four runs) on the easiest two difficulty levels,  and also excluded if they did not have at least two out of four runs where 50% or more of the blocks were acceptable in terms of performance (60% accuracy) and motion. Motion was deemed acceptable if children moved less than 1.5 mm from their median head position in at least 60% of the volumes within a task block. See the fMRI preprocessing section below for a description of displacement calculations. A 60% accuracy criterion was chosen because while we could be sure that children were performing better than chance (50%), it was not too stringent. Motion parameters were also entered into the fMRI preprocessing pipeline.

Image acquisition
All imaging data were acquired using a 3 T Siemens Trio MRI scanner with a 12-channel head coil. Head stabilization and motion restriction were achieved with foam padding. The structural scan was a high-resolution

WMTB-C
Standardized scores on the subscales were compared across group using t tests to determine if there were differences between ASD and TD children on these neuropsychological measures of WM.

fMRI data analyses
Image preprocessing of functional data was performed using a combination of standard AFNI [58] and FMRIB's Software Library (FSL) [59] tools. The first three volumes of each run were discarded for scanner stabilization. After slice timing and motion correction, data were smoothed using a 6-mm FWHM Gaussian kernel, temporally filtered (lower and upper cutoff frequencies of 0.01 and 0.2 Hz, respectively) and converted to percent signal change from the baseline volumes. Before group-level analyses, images were registered to the Montreal Neurological Institute (MNI) 152 template. The maximum Euclidean displacement (MD) travelled by any voxel within the brain was calculated from the six rigid body transformation parameters for each volume. This MD metric was used to flag volumes with unacceptable motion, as described above. The average MD for each subject was used to explore group differences in head motion. Although more motion was found in children with ASD (M = 0.58 mm, SD = 0.54 mm) than TD children (M = 0.29 mm, SD = 0.22 mm), t (34) = 2.07, p = 0.05, both groups had minimal average motion of under 0.60 mm. To control for motion, MD was also included as a covariate of no interest in the GLM. Data were analyzed with the FSL fMRI Expert Analysis Tool (FEAT) [60]. Data were fit first to a block-design general linear model convolved with a gamma function to model haemodynamic response, using the task parameters (D3 to D6). To examine areas that linearly modulated as a function of difficulty, linear trend analyses were conducted from D3 to D6 using fixed-effects higher level modelling. Individual results were then averaged across runs for each subject in a second-level analysis. Betweengroup comparisons were carried out using FMRIB's Local Analysis of Mixed Effects-1 (FLAME 1) [59] to obtain an (See figure on previous page.) Figure 1 Protocol description of the colour matching task (CMT). (A) There were six levels of difficulty where the number of relevant colours (yellow, purple, pink, red, orange, brown and grey) increased by one to increase the difficulty level. Difficulty = (# of colours) + 2. Children were taught to ignore the clown's face, colour location, colour repetition and irrelevant colours (blue and green). (B) It was a block design task, where each run consisted of six 32-s task blocks (for each difficulty) followed by 20-s baseline blocks where clowns are presented in only blue and green (ignore). Task blocks were presented pseudo-randomly within each run. (C) Example of part of a sequence in a baseline block. Stimuli were presented for 3 s followed by 1-s inter-stimulus cross. Children were instructed to not respond. (D) Example of part of a sequence in a task block; children indicated if the current clown was wearing the same or different colours as the previous clown. Stimuli were presented for 3 s followed by a 1-s inter-stimulus fixation cross.
accurate between-subject variance estimation, which increased our ability to detect real activation [60]. Significant activations were reported using cluster-based thresholding determined by Z > |2.3| and a corrected cluster significance threshold of p < 0.05. Regions of interest (ROIs) were selected from the local maxima of areas showing significant group differences between TD and ASD groups in the linear trend analyses. Average percent signal change and standard error scores were extracted from spherical ROIs (6 mm radius) centred about the local maxima of group difference maps, and the mean peak cluster signal change for each group was plotted as a function of difficulty to further examine activation patterns across cognitive load.

Behavioural data
There was a significant effect of group on accuracy (Figure 2A), F (1, 34) = 5.15, p = 0.03, η 2 = 0.13, which was driven by TD children (M = 0.80, SD = 0.12) performing more accurately than ASD children (M = .71, SD = 0.12) only at D5 (t (34) = 2.30, p = 0.03). As such, comparisons of brain activity between control and ASD groups were made under comparable accuracy scores across most levels. There was a significant main effect of difficulty on accuracy, F (3, 34) = 53.87, p < 0.001, η 2 = 0.61, with performance accuracy decreasing as a function of difficulty level in both groups. Post hoc pairwise comparisons, adjusted for multiple comparisons using Bonferroni, revealed that accuracy on the most difficult level was significantly different from each other at p < 0.05, the exceptions being between D5 and D6 in the ASD group (Table 2), and between D3 and D4, and between D4 and D5 in TD children (Table 2). We also ran a supplementary analysis, inputting performance as a covariate. There were no areas of activation that correlated with performance in either group, suggesting that findings were not confounded by performance per se. Also, there were no significant differences in performance across runs in the TD (F (3,48) =1.21, p = 0.32) or ASD (F (3,45) = 0.50, p = 0.69) group, and therefore patterns of performance did not reflect fatigue across time.
Overall, response times increased with increasing difficulty in TD children but only increased up until D5 in children with ASD ( Figure 2B). There was no main effect of group on response times (F (1, 34) = 2.33, p = 0.14), but there was a main effect of difficulty level, F (3, 34) = 54.60, p < 0.001, η 2 = 0.62. Post hoc comparisons in the ASD group showed that response times differed between difficulty levels, except between D4 and D6, and between D5 and D6 (Table 3). In TD children, response times significantly differed between difficulty levels, except between D4 and D6 and D5 and D6 (Table 3).
For the WMTB-C, TD children had significantly higher scores on the listening recall (t (34) = 2.48, p < 0.05) and backward digit recall (t (34) = 2.28, p < 0.05) subtests than children with ASD. The groups did not differ on any other subtests (see Table 1 for scores on the WMTB-C).

fMRI Data Task-related activation within groups
The primary objectives of the fMRI analyses were to investigate the pattern of brain activity exhibited as a function of cognitive load (i.e. difficulty level) and determine if this pattern differed in children with and without ASD. Linear trend analyses (D3 to D6) showed that while some brain areas increased in activity as a function of difficulty level, others decreased. 'Increasing activation' refers to an increase in BOLD signal with increasing load (i.e. positive linear relations between cortical activity and task difficulty level) and 'decreasing activation' refers to a decrease in BOLD signal with increasing load (i.e. negative linear relations between cortical activity and task difficulty level). As shown in Figure 3, magnitude of the signal change increased with difficulty in TD children, particularly between D5 and D6, in bilateral fusiform (BA37), precuneus (BA7), inferior frontal gyri (BA45), right DLPFC (BA9) and bilateral dorsal cingulate/dorsal medial prefrontal cortex (B32/8) extending to the anterior cingulate (B33/24). In the bilateral middle occipital gyrus (BA19), a positive linear change in activation with increasing task difficulty was seen up to D6 (Table 4). Activation found in the bilateral posterior cingulate (BA 23/31) and anterior medial prefrontal gyrus (BA10) decreased (i.e. showed a negative linear relation) as a function of difficulty level (Table 4).
Children with ASD did not show the same trends in activation, particularly in the frontal brain regions. A positive linear relation between cortical activity and task difficulty was seen only in the posterior brain regions, including the bilateral middle occipital (BA19) and fusiform gyri (BA37, Table 5). Similar to TD children, activation in the posterior cingulate (BA23/31) and medial prefrontal gyrus (BA10/32) decreased (i.e., showed a negative linear relation) as a function of difficulty level in children with ASD (Table 5).

Between-group differences in task-related activation
Between-group analyses of each difficulty level separately showed no significant differences (see Additional file 2). However, there were significant group effects in linear activation trends as a function of difficulty level. As shown in Figure 4, three regions demonstrated group differences in activation. These areas included the bilateral precuneus (BA7), right DLPFC (BA9) and left dorsal medial premotor cortex (BA8) ( Table 6). In these regions, TD children showed significant positive linear relations between cortical activity and difficulty level, particularly between D5 and D6, whereas children with ASD did not show these trends between D5 and D6. The magnitude of the signal change in frontal and parietal areas increased largely with the most difficult level when TD children tried to meet WM demand, but children with ASD failed to show this increase. See Figure 5 for graphs of percent signal change with standard error for ROIs of brain regions that showed significantly different linear patterns between TD and ASD groups.

Discussion
This is the first study to examine the neural correlates of visuo-spatial WM in pre-adolescent children with ASD relative to TD children. Whereas children with and without ASD did not differ in brain activation during WM function irrespective of load, between-group differences were observed in the linear activation trends across difficulty level. In other words, children with and without ASD differ in how they modulate WM processes during tasks that increase in difficulty level. Using a task that isolated cognitive load on WM, we observed positive linear relations between cortical activity and task difficulty level in prefrontal and parietal regions, particularly between the highest difficulty levels, in TD children that were significantly different from activation trends seen in children with ASD. These areas included bilateral precuneus, right DLPFC and left medial premotor cortices. Other areas known to underlie WM function, including the anterior cingulate and inferior frontal gyri, followed this same pattern but the linear activation trends did not differ significantly between TD and ASD subjects, likely due to subject variability. Overall, TD children demonstrated an opposing system of cognitive processes where areas related to task difficulty (frontal regions) increased in activity and areas associated with the brain's DMN (posterior cingulate and anterior medial frontal gyrus) decreased in activity   with increasing cognitive demand; this pattern of activation was absent in children with ASD. There were significant differences in performance on only two verbal WM subtests of the WMTB-C. Further, TD and ASD groups performed similarly on the CMT (except for a small difference at D5), suggesting that the findings are not confounded by behavioural differences. TD children showed increased recruitment of both frontal and posterior parietal regions as a function of cognitive load, whereas the ASD group showed increased recruitment of only the posterior visual regions, including the left and right fusiform and middle occipital gyri. Our observed pattern of findings in children with ASD are consistent with findings in adults with ASD who also relied on posterior regions related to low-level cognition during WM tasks, rather than areas of high level-cognition, such as in the prefrontal cortices [10]. More activation in occipital-temporal areas in ASD has been proposed to reflect visually based processing styles and a tendency to rely primarily on visual features and details of objects, rather than on WM [61]. These findings are in line with the enhanced perceptual functioning model [62] that suggests individuals with ASD display superior activation of visuo-perceptual regions in association with a reduced activation in higher order frontal areas. The present study extends this literature in adults by highlighting that this processing style appears in young children with ASD as well. The same pattern of activation has also been found in socially relevant tasks [63,64]. Although greater posterior activation may be adequate for visuo-spatial processing, it is likely inefficient for more complex cognitive operations, such as language comprehension [19] and WM [10], which would have implications for the social and executive deficits typically observed in ASD. Less extensive use of prefrontal areas may be the result of early abnormal prefrontal development [15,16] and connectivity [10,19] reported in individuals with ASD.
The finding of stable activation with increasing cognitive load within the prefrontal cortex may be most significant for ASD due to the crucial role of frontal regions in WM and executive function [10,43,44,[65][66][67]. Similar to adults [10,12] and adolescents [14] with ASD, we found that pre-adolescent children with ASD demonstrated different activation patterns in the DLPFC (BA 9) during WM function. Specifically, in the DLPFC, TD children showed a graded increase in activity from D3 to D6, whereas activity within this region did not increase at the same magnitude in children with ASD. The DLPFC is believed to play a critical role in holding information 'online' [26] and mediating strategic organization and data compression processes [68,69], hence its sensitivity to increasing cognitive demand in ours and other studies [68][69][70][71]. Stable activation across difficulty level in this area suggests that young children with ASD may fail to use appropriate organizational strategies, such as Results from linear trend analyses from D3 to D6 for TD children. Areas that increased as a function of difficulty level (A) are associated with WM and visuo-spatial processing, whereas areas that decreased as a function of difficulty level (B) are associated with the default mode network. MNI coordinates represent the peak Z value of the cluster; X peak local maximas within cluster.
'chunking' methods, which facilitate WM by simplifying cognitive load/demand. Furthermore, normative studies have demonstrated that with increasing age, individuals rely more on DLPFC in WM processes [43,66,72], suggesting that specialization of this region for WM coincides with structural maturation across development. It is possible that atypical abnormal growth patterns in ASD in the frontal cortex, the DLPFC in particular [15], adversely impacts its functional integrity; these speculations require further research.
Failure to increase recruitment of the precuneus across difficulty level in ASD subjects is also of interest, given evidence supporting parietal involvement in addition to the prefrontal cortex during spatial WM processing [44,66,70,73]. The precuneus forms part of the proposed occipito-parietal network, or alternatively the visual 'dorsal pathway' , responsible for spatial visual processing (e.g. object location) [74]. Conversely, the occipito-temporal, or visual 'ventral pathway' , includes the fusiform gyri and is critical to object identification (e.g. colour and shape) Results from linear trend analyses from D3 to D6 for children with ASD. Areas that increased as a function of difficulty level (A) are associated with WM and visuo-spatial processing, whereas areas that decreased as a function of difficulty level (B) are associated with the default mode network. MNI coordinates represent the peak Z value of the cluster; X peak local maximas within cluster. Figure 4 Results from between-group comparisons. Significant activations using cluster-based thresholding determined by Z > 2.3 and a corrected cluster significance threshold of p = 0.05. Areas in red depict regions where the control children showed greater positive linear relations between cortical activity and task difficulty level than the ASD group. There were no areas where the ASD group showed greater linear activation trends across difficulty level in the negative or positive direction than controls. Prec precuneus, DLPFC dorsolateral prefrontal cortex, dmPMC dorsal medial premotor cortex. [74]. Given that CMT was a visuo-spatial task with gradual increases in WM that needed to be processed, areas responsive to both spatial search (occipito-parietal regions) and categorization (occipito-temporal regions) should show increased recruitment across cognitive load. Although group differences in fusiform activity were absent, we observed a significant difference in the precuneus; TD children showed positive linear relations between cortical activity and difficulty level, particularly between D5 and D6, whereas the ASD group did not show this linear change. Impaired parietal activity in pre-adolescent children with ASD lends support to the proposed dorsal stream deficits within this population as well as in other developmental disorders [75][76][77][78], while ventral processing is relatively intact. Furthermore, normative developmental studies suggest that the dorsal pathway has a more protracted maturational trajectory than the ventral stream [66], and WM fMRI studies demonstrate greater parietal [4,73] but less fusiform recruitment [66] with age. Thus, spared fusiform/ occipital activity and abnormal parietal functioning in our ASD group may reflect immature WM processing typically seen in very young TD children. However, it is not clear whether the recruitment of mature neural substrates underlying WM processes is delayed, persistently weak or arrested. Future longitudinal work will help clarify the developmental path of WM neural circuitry in ASD. Results from between group comparisons of the linear trend analyses from D3 to D6. All regions reported are areas in which TD children showed greater positive linear relations between cortical activity and difficulty level (i.e. increasing BOLD signal with increasing difficulty) than children with ASD. There were no areas where children with ASD showed greater linear relations between cortical activity and difficulty level than TD children. MNI coordinates represent the peak Z value of the cluster; X peak local maximas within cluster. Figure 5 Mean peak cluster percent signal changes and standard error. As a function of difficulty between task difficulty and baseline conditions in areas where children with ASD significantly differed from TD children in the linear trend analyses. DLPFC dorsolateral prefrontal cortex.

Conclusions
Overall, the current study fills a significant gap in our knowledge of neural substrates implicated in visuo-spatial WM functioning among pre-adolescent children and how they differ in children with ASD. Similar to adults, we found that children with ASD show prefrontal and parietal system abnormalities and tend to rely on posterior brain regions associated with lower level cognitive processing. The present study also converges on the growing body of literature proposing dorsal stream deficits in visuo-spatial processing in this population. This study contributes new information on WM differences in school-aged children with and without ASD, demonstrating that even in preadolescents, there are significant differences in brain activation patterns during WM processing with increasing cognitive load that differentiate the groups. In light of similar task performance between groups, it is also important to entertain the idea that higher order processing may not be mandatory in ASD with tasks that can be processed efficiently using a perceptual processing approach because of their proposed enhanced perceptual functioning [62]. Conversely, higher order control may be necessary in TD children who do not have such perceptual expertise. Although data may support this argument, ASD is generally considered a network/brain disorder, and what is interesting is how networks in children with ASD differ from TD children with complex tasks. Group effects were driven by the most difficult level analyzed, which may be due to strategies used by TD children necessary to encode the increased number of colours in the highest task load (i.e. meet high WM demand) that children with ASD are unable to employ. Similarly, previous research reports behavioural differences with tasks that have more complex WM demands e.g. [9,34,39]. It is important to consider the limitations of the current study when interpreting results. With comparable CMT behavioural performance between ASD and TD groups, we eliminated performance as a confounding factor. Consequently, our sample was less representative of lowfunctioning individuals with ASD, and thus, results are generalizable to higher functioning individuals only. Future fMRI research is required to understand WM function across various levels of functioning and a range of symptoms. In addition, given our choice of control subjects (TD children), findings can only provide information about differences from the norm. Comparisons to other atypical populations who share similar cognitive but different clinical profiles (e.g. ADHD) [5] will further our understanding about the neural patterns that are unique to ASD, potentially explaining characteristic behaviour in this complex group. Lastly, results are also limited by the relatively small sample size, due to assessing a complex cognitive ability while scanning, which greatly increased the amount of children's data that was not usable.
The findings from the current study have a number of significant implications despite these limitations. Several researchers have stressed the importance of executive cognitive skills for social function [5,79], a core deficit of ASD. However, there remains a gap in knowledge regarding the link between neuropsychopathology and clinical symptoms of autism. With respect to our findings of atypical neural activity underlying WM processes in ASD, this may impair the ability to hold information 'online' that may affect one's ability to evaluate and select appropriate responses during peer interactions, translating into socially inappropriate behaviour. Social behaviour requires complex cognitive processing, further highlighting the importance of our results of increasing group differences in brain activity with increases in complexity. Future neuroimaging studies could work towards understanding the relation between ASD symptomology and neural activation patterns associated with WM processing. Overall, our findings will help guide future longitudinal work by localizing areas of vulnerability to developmental disturbances and allow health care providers to carefully monitor their development. Developmental information will allow us to identify the nature and timing of atypical development, which is critical in establishing age-appropriate cognitive or pharmacological remediation for WM function and behavioural deficits in ASD.

Additional files
Additional file 1: fMRI data for medicated and non-medicated children with ASD. Percent signal change as a function of difficulty between task difficulty and baseline conditions in children with ASD who were on medication versus those who were not. Areas of the brain shown are from regions where children with ASD significantly differed from TD children in the linear trend analyses. Potential differences between medicated and non-medicated children with ASD were also examined statistically using the FSL FEAT, and no significant differences were found between children with ASD who were and were not on medication. However, due to low N (only six subjects on medication), this statistical test may not be reliable. Therefore scatter plots were created to visually examine the data for significant group differences, and this reaffirmed that medication does not appear to affect the findings.
Additional file 2: Individual group activation maps for children with ASD and typically developing children at all levels (D3 to D6) of CMT. Significant activations using cluster-based thresholding determined by Z > |2.3| and a corrected cluster significance threshold of p = 0.05. Areas in red and blue depict regions with significantly higher and lower BOLD signal than baseline, respectively. Between-group comparisons showed no areas of significant difference between children with and without ASD at any single difficulty level. participated in design, analyses and revising manuscript. All authors read and approved the final manuscript.