Figure 1. Random-field theory-based cluster-corrected (P < .05) maps for cortical thickness (A) and cortical volume (B). Relative deficits in adults with autism spectrum disorder compared with controls are displayed in red/yellow, while excesses are displayed in blue/cyan. There were no clusters of significant differences for surface area.
Figure 2. Spatially distributed differences in cortical thickness (A), surface area (B), and cortical volume (C) between individuals with autism spectrum disorder and controls (P < .05, uncorrected). Relative deficits in adults with autism spectrum disorder compared with controls are displayed in red/yellow, while excesses are displayed in blue/cyan. Surfaces are presented in lateral, medial, and frontal views for the left and right pial (outer) surface.
Figure 3. Significant correlations between cortical thickness and surface area in controls (P < .01, uncorrected). White lines indicate regions of between-group differences in surface area or cortical thickness (P < .05, uncorrected). Positive correlations are displayed in red/yellow, while negative correlations are displayed in blue/cyan.
Figure 4. Overlap between patterns of differences in cortical thickness and surface area (A), with letters indicating the regions with uncorrelated measures of cortical thickness and surface area as shown in the scatterplots (B-E). ASD indicates autism spectrum disorder.
Ecker C, Ginestet C, Feng Y, Johnston P, Lombardo MV, Lai MC, Suckling J, Palaniyappan L, Daly E, Murphy CM, Williams SC, Bullmore ET, Baron-Cohen S, Brammer M, Murphy DGM; MRC AIMS Consortium. Brain surface anatomy in adults with autism: the relationship between surface area, cortical thickness, and autistic symptoms. JAMA Psychiatry. Published online November 26, 2012. doi:10.1001/jamapsychiatry.2013.265.
eAppendix. MRI data acquisition and measures of areal expansion/compression.
eFigure. Brain map (voxel-based morphometry) of significant gray matter differences in ASD (N = 178; P = .0041) rendered to the cortical surface.
Ecker C, Ginestet C, Feng Y, Johnston P, Lombardo MV, Lai M, Suckling J, Palaniyappan L, Daly E, Murphy CM, Williams SC, Bullmore ET, Baron-Cohen S, Brammer M, Murphy DGM, MRC AIMS Consortium FT. Brain Surface Anatomy in Adults With AutismThe Relationship Between Surface Area, Cortical Thickness, and Autistic Symptoms. JAMA Psychiatry. 2013;70(1):59-70. doi:10.1001/jamapsychiatry.2013.265
Author Affiliations: Department of Forensic and Neurodevelopmental Sciences, Institute of Psychiatry, King's College London, London (Drs Ecker, Ginestet, Feng, Johnston, C. M. Murphy, Williams, Brammer, and D. G. M. Murphy and Ms Daly), Autism Research Centre (Drs Lombardo, Lai, and Baron-Cohen) and Brain Mapping Unit (Drs Suckling and Bullmore), Department of Psychiatry, University of Cambridge, Cambridge, and Division of Psychiatry, University of Nottingham, Nottingham (Dr Palaniyappan), England.
Context Neuroimaging studies of brain anatomy in autism spectrum disorder (ASD) have mostly been based on measures of cortical volume (CV). However, CV is a product of 2 distinct parameters, cortical thickness (CT) and surface area (SA), that in turn have distinct genetic and developmental origins.
Objective To investigate regional differences in CV, SA, and CT as well as their relationship in a large and well-characterized sample of men with ASD and matched controls.
Design Multicenter case-control design using quantitative magnetic resonance imaging.
Setting Medical Research Council UK Autism Imaging Multicentre Study.
Participants A total of 168 men, 84 diagnosed as having ASD and 84 controls who did not differ significantly in mean (SD) age (26  years vs 28  years, respectively) or full-scale IQ (110  vs 114 , respectively).
Main Outcome Measures Between-group differences in CV, SA, and CT investigated using a spatially unbiased vertex-based approach; the degree of spatial overlap between the differences in CT and SA; and their relative contribution to differences in regional CV.
Results Individuals with ASD differed from controls in all 3 parameters. These mainly consisted of significantly increased CT within frontal lobe regions and reduced SA in the orbitofrontal cortex and posterior cingulum. These differences in CT and SA were paralleled by commensurate differences in CV. The spatially distributed patterns for CT and SA were largely nonoverlapping and shared only about 3% of all significantly different locations on the cerebral surface.
Conclusions Individuals with ASD have significant differences in CV, but these may be underpinned by (separable) variations in its 2 components, CT and SA. This is of importance because both measures result from distinct developmental pathways that are likely modulated by different neurobiological mechanisms. This finding may provide novel targets for future studies into the etiology of the condition and a new way to fractionate the disorder.
Autism spectrum disorder (ASD) is characterized by a triad of symptoms, namely impaired social communication, deficits in social reciprocity, and repetitive and stereotyped behavior.1,2 There is consensus that ASD is a highly genetic neurodevelopmental condition that is accompanied by differences in brain anatomy.
Evidence for neuroanatomical differences in ASD comes from a variety of postmortem and structural neuroimaging studies (reviewed by Amaral et al3 and Toal et al4). For example, differences have been described in the cerebellum,5 amygdala-hippocampal complex,6- 9 frontotemporal regions,6,7,10 and caudate nucleus.11,12 Most of these prior anatomical studies were based on volumetric analyses. However, in particular for the cortical volume (CV)—owing to its specific anatomical characteristics—it would be informative for the understanding of the biological basis of ASD to determine whether differences in CV are driven by differences in cortical thickness (CT) or surface area (SA), or a combination of both.
The existence of such independent variations in CT and SA is potentially of great importance as there is evidence to suggest that they have distinct genetic determinants,13 contrasting phylogeny,14 and differing developmental trajectories.15 A recent twin study suggests that although both CT and SA are highly heritable, they are unrelated genetically and follow discrete genetic mechanisms.13 Also, it has been proposed that measures of CT and SA reflect different aspects of the underlying neural architecture. For instance, it is now widely recognized that the cerebral cortex is organized into ontogenetic columns16 and that cells within a column share a common origin before migrating to their location within the cortex during development (radial unit hypothesis).14,17,18 In this model, the number of cells within cortical columns mediates CT, whereas SA primarily reflects the number of columns within a cortical region.14 Measurements of CV may therefore reflect structural properties that are unique to cortical SA or unique to CT. Thus, it is necessary to explore CT and SA separately to better understand the neurobiological mechanisms associated with brain abnormalities in ASD and to refine the autism endophenotype for future etiological studies.
To date, relatively few studies have investigated region-specific differences in CT in individuals with ASD, and those that are available reported inconsistent results. In children, some studies report increases in CT of frontotemporal regions,19,20 while others note the opposite trend in similar regions.21 In adults, most studies report cortical thickening of the frontal cortex,22,23 while the temporal lobe may be thicker22 or thinner24 in ASD. To our knowledge, only 2 prior studies on ASD examined CV, SA, and CT in the same sample of individuals. First, a recent study by Hazlett et al25 suggests that increased CV may be associated with increased cortical SA rather than CT. Second, Raznahan et al26 reported altered neurodevelopmental trajectories for CV and CT (also see the articles by Mak-Fan et al19 and Scheel et al22), but not SA, in a cross-sectional study of children and adults with ASD. Notably, a differential growth trajectory for CT and SA has also been reported in neurotypical subjects, in whom frontal SA and CV decrease with increasing age (ages 12-16 years), while frontal CT increases during maturation.27These prior studies were important first steps toward establishing the specific determinants of volumetric differences in ASD, and they further highlight the need for investigating CT and SA in isolation.
Our aim was to investigate regional differences in CV, SA, and CT as well as their relationship in a large and well-characterized sample of men with ASD and matched controls. We used a spatially unbiased vertex-based approach that provides measures of CV, SA, and CT at several thousand points across the cortical sheet. This allowed us to investigate the following: (1) the spatially distributed networks of differences in CT and SA; (2) the degree of spatial overlap between them; and (3) their relative contribution to observed differences in regional CV. It was hypothesized that individuals with ASD show neuroanatomical differences in predominantly frontotemporal regions and that these differences are associated with the severity of autistic symptoms.
Eighty-four right-handed men with ASD and 84 matched controls aged 18 to 42 years were recruited by advertisement and subsequently assessed at 1 of 3 centers: the Institute of Psychiatry, King's College London; the Autism Research Centre, University of Cambridge; and the Autism Research Group, University of Oxford. Approximately equal ratios of cases to controls were recruited at each site: London, 38:38; Cambridge, 31:29; and Oxford, 15:17.
Exclusion criteria for all participants included a history of major psychiatric disorder, head injury, genetic disorder associated with autism, or any other medical condition affecting brain function. We excluded participants with substance abuse and participants on antipsychotic medications, mood stabilizers, or benzodiazepines. All participants with ASD were diagnosed according to International Statistical Classification of Diseases, 10th Revision research criteria and confirmed using the Autism Diagnostic Interview–Revised28 (ADI-R) to ensure that all participants with ASD met the criteria for childhood autism. All cases of ASD reached ADI-R algorithm cutoffs in the 3 domains of impaired reciprocal social interaction, communication, and repetitive behaviors and stereotyped patterns, although failure to reach cutoff in 1 of the domains by 1 point was permitted (Table 1).
Current symptoms were assessed using the Autism Diagnostic Observation Schedule29and were not used as an inclusion criterion. Overall intellectual ability was assessed using the Wechsler Abbreviated Scale of Intelligence.30 All participants fell within the high-functioning range on the spectrum defined by a full-scale IQ greater than 70. All participants gave informed written consent in accordance with ethics approval by the National Research Ethics Committee, Suffolk, England.
All participants were scanned with contemporary magnetic resonance imaging (MRI) scanners operating at 3T and fitted with an 8-channel receive-only head coil (GE Medical Systems HDx at the Autism Research Centre, University of Cambridge and the Institute of Psychiatry, King's College London; Siemens Medical Systems Trim Trio at the Autism Research Group, University of Oxford). A specialized acquisition protocol with quantitative imaging (driven-equilibrium single-pulse estimation of T1) was used to ensure standardization of structural MRI scans across the 3 scanner platforms. This protocol has been validated and is described elsewhere31 (eAppendix, eReferences, and eFigure).
Scans were initially screened by a radiologist to exclude clinically significant abnormalities and to assess the existence of movement. Scans of insufficient quality were excluded from the analysis. Brain anatomy of the data has already been investigated with another approach in a previous study.32 As the present analysis requires an advanced quality standard, 9% of the recently used MRI scans had to be rejected as being not sufficiently applicable for the procedure.
The FreeSurfer analysis suite (http://surfer.nmr.mgh.harvard.edu/) was used to derive models of the cortical surface in each T1-weighted image. These well-validated and fully automated procedures have been extensively described elsewhere.33- 39 In brief, a single filled white matter volume was generated for each hemisphere after intensity normalization, skull stripping, and image segmentation using a connected components algorithm.33 Then, a surface tessellation was generated for each white matter volume by fitting a deformable template. This resulted in a triangular cortical mesh for gray and white matter surfaces consisting of approximately 150 000 vertices (ie, points of triangles) per hemisphere.
Measures of CT are the closest distance from the gray and white matter boundary to the gray matter and cerebrospinal fluid boundary at each vertex on the tessellated surface.35 Vertex-based estimates of SA were obtained by computing the average of the area of the triangles incident to that vertex (ie, sharing that vertex) in a standardized, spherical atlas-space surface tessellation when mapped into the individual subject space (eAppendix, eReferences, and eFigure). This provides point-by-point estimates of the relative areal expansion or compression of each location in atlas space.40,41 Here, the local SA was used interchangeably with areal expansion or compression. Estimates of regional CV were derived by multiplying CT measures by their areal expansion or compression at each vertex. These measures are thus different from conventional measures of brain volume resulting from the standard FreeSurfer pipeline, and they indicate the degree of volumetric expansion or compression at each vertex. To improve the ability to detect population changes, each parameter was smoothed using a 15-mm surface-based smoothing kernel.
Group differences in total brain volume, gray matter volume, mean CT, and SA as estimated by FreeSurfer42 were assessed using t tests for independent samples.
Statistical analysis was conducted using the SurfStat toolbox (http://www.math.mcgill.ca/keith/surfstat/) for Matlab (R2010b; MathWorks). Parameter estimates for each measure (CV, CT, and SA) and the main effect of group (Gi) were estimated by regression of a general linear model at each vertex i and subject j, with center (Ci) as a categorical fixed-effects factor and age, IQ, and a total brain measure (indicated by Bi: total brain volume for CV, mean CT for CT, and total SA for SA) as continuous covariates: yij = β0 + β1Gj + β2Cj + β3Age j + β4IQ j + β5 Bj + ϵ ij, where ϵ is the residual error. Between-group differences were estimated from the fixed-effect coefficient β1 normalized by the corresponding standard error. Corrections for multiple comparisons across the whole brain were performed using random-field theory-based cluster-corrected analysis43 for nonisotropic images using a P < .05 (2-tailed) cluster significance threshold.
To explore the wider neural systems underlying ASD and to investigate the degree of spatial overlap for the 3 different measures, we reexamined between-group differences heuristically using an uncorrected threshold of P < .05 with a cluster threshold of 50 vertices using the same model(s) as described earlier. This resulted in spatially distributed binary patterns of differences unique to CT and/or SA, as well as their overlap, regardless of the sign (ie, based on their statistical threshold). First, a χ2 test was used to compare frequencies of unique or overlapping differences in each morphometric parameter (ie, contingency tables), testing the null hypothesis that differences in CT and SA are equally distributed. Second, a simulation strategy was used to assess whether the observed degree of overlap between differences in CT and SA is consistent with the idea of 2 spatially independent patterns. This hypothesis was tested on the basis of N = 10 000 randomly generated difference maps (ie, maps containing random t values, thresholded at P < .05) for CT and SA. The extent of overlap (ie, number of vertices with differences in CT and SA) was then assessed in each of the 10 000 overlapping patterns to derive a probability value of obtaining a given percentage of overlap on the basis of randomly varying patterns of differences.
The relationships between regional anatomical abnormalities and domains of symptom severity were explored using Pearson correlation coefficients. Within the ASD group, we examined correlations between volumetric features (ie, mean CT and CV) within and of clusters showing a significant between-group difference and the 3 domains of the ADI-R measuring symptoms at ages 4 to 5 years as well as the total Autism Diagnostic Observation Schedule scores of current symptoms.
Participant demographic characteristics and global brain measures are shown in Table 1. There were no differences (2-tailed) between people with ASD and controls in age (t166 = −1.82; P = .07), full-scale IQ (t166 = 1.78; P = .07), or verbal IQ (t166 = −0.22; P = .82). However, individuals with ASD had significantly lower performance IQs (t166 = −3.54; P = .001). There were no significant group differences in total brain volume (t166 = −0.19; P = .84), total SA (t166 = −1.20; P = .23), or mean CT (t166 = −0.51; P = .61). Using Levene's test for equality of variances,44 we found no evidence of heterogeneity of variance in total brain measures.
Following correlation for multiple comparisons across the brain, individuals with ASD had significantly increased CT in a large cluster located in the left lateral prefrontal cortex (P < .003). This cluster had 2 peaks in the ventrolateral and rostrolateral prefrontal cortex (peak Talairach coordinates: −36, 42, 4 and −36, 20, 22, respectively). In addition, individuals with ASD had significantly decreased CT of the right anterior temporal lobe (P < .04) including the superior, middle, and inferior temporal gyrus (Figure 1A).
Individuals with ASD had significantly decreased CV in a large cluster in the bilateral orbitofrontal cortex (peak Talairach coordinates, left: −19, 44, −13; right: 8, 36, −19; Figure 1B) relative to controls (P < .005). We did not observe any significant cluster between-group differences in SA.
To examine the underlying contribution of variations in CT and SA to the observed differences in CV and to investigate the spatial overlap between CT and SA, we reexamined the difference maps heuristically using a more lenient statistical threshold of P < .05 (uncorrected). This resulted in spatially distributed patterns of differences indicating the wider neural systems implemented in ASD.
Overall, individuals with ASD had increased CT across the frontal lobe including the medial prefrontal cortex, the superior frontal gyrus, the dorsal and rostral lateral prefrontal cortex, and the left pars opercularis and pars triangularis area of the inferior frontal gyrus—overlapping with the Broca area (Figure 2A). Individuals with ASD also had thicker temporal cortices of the right posterior temporal lobe (eg, middle and superior temporal gyrus) and left angular gyrus. In contrast, individuals with ASD had lower CT bilaterally in the anterior part of the superior temporal lobe and in the parahippocampal cortex, in the right anterior fusiform gyrus and occipital cortex, and in the left subgenual anterior cingulate cortex (Table 2).
Individuals with ASD showed decreased SA (Figure 2B) of the following: (1) right dorsolateral prefrontal cortex and bilateral orbitofrontal cortex; (2) left inferior temporal gyrus and anterior cingulate; and (3) supplementary and presupplementary motor cortices. In contrast, the ASD group had increased SA of the temporoparietal junction including the supramarginal gyrus, the parahippocampal gyrus, and the medial superior parietal cortex (Table 3).
Across both hemispheres, we found that the number of vertices with a significant difference in CT only was approximately equal to those with a significant difference in SA only (51% vs 45%, respectively; χ22 = 0.37; P = .83) (Table 4).
The patterns of significant differences in CT and SA were largely nonoverlapping and shared only about 3% of all significantly different spatial locations on the cerebral surface, ie, the probability that any 1 vertex has an overlapping difference in both CT and SA is very low. Simulations revealed that the probability of obtaining the same degree of overlap (ie, 3%) between 2 randomly generated difference maps is as high as 85%. The observed percentage of overlap is hence consistent with the idea that differences in SA and CT are spatially independent and that SA and CT contribute in a unique way to the group differences in CV.
Regions with a significant group difference in SA or CT also did not overlap with regions displaying a significant association between SA and CT in general (Figure 3). The observed difference maps therefore coincided with regions displaying generally low coupling between SA and CT.
Nevertheless, some regions showed significant abnormalities in both CT and SA (Figure 4A). These included the left hemisphere superior frontal gyrus, right medial orbitofrontal cortex, left medial frontal cortex, and right hemisphere ventrolateral prefrontal cortex. Within these regions, measures of CT and SA were also not correlated (all P > .20) and hence uniquely contributed to measures of CV (see scatterplots in Figure 4B-E).
Of all the underlying differences in CV (P < .05, uncorrected) (Figure 2C), a total of 67% were modulated by (ie, overlapped with) significant differences in CT only (8%), SA only (56%), or both CT and SA (5%) (Table 5). Thus, differences in SA explained a significantly larger proportion of differences in CV than did differences in CT (χ22 = 46; P < .001). The remaining 31% of all significant differences in CV could not be explained by significant differences in either SA or CT and must therefore be due to a combination of subthreshold differences in CT or SA (ie, P > .05, uncorrected).
Within the ASD group, there were significant positive correlations between the left dorsolateral frontal cluster, where individuals with ASD displayed a significant increase in CT, and higher scores on the ADI-R communication (r = 0.23; P = .02) and repetitive (r = 0.26; P = .009) domains as well as the total ADI-R scores (r = 0.25; P = .01). We also observed a marginally significant correlation between differences in CT of the left temporal lobe and ADI-R symptoms in the repetitive domain (r = 0.18; P = .47) (Table 6).
In the orbitofrontal cortex, where individuals with ASD showed significant reductions in CV, we found a significant negative correlation between symptom severity on the ADI-R social domain in the left (r = −0.23; P = .03) and right (r = −0.22; P = .04) hemispheres. Thus, individuals with more severe social autistic symptoms at ages 4 to 5 years displayed significantly smaller CV of the orbitofrontal lobes.
To our knowledge, this is the first study to examine regional differences in CV on the basis of its 2 components, CT and SA, in a large and well-characterized sample of men with ASD. We found that individuals with ASD had significant differences in CT and CV, which primarily affected the frontal and temporal lobes. These differences were caused by subtle (ie, subthreshold) variations in CT and SA that were largely nonoverlapping. Furthermore, these spatially independent patterns of differences did not coincide with regions displaying significant associations between SA and CT, but instead overlapped with areas of generally low coupling between them. Differences in CV observed in adults with ASD may therefore be underpinned by separable variations in its 2 components, CT and SA, which reflect variations in 2 independent neuroanatomical measures. This is of importance because CT and SA result from distinct developmental pathways that are likely modulated by different neurobiological mechanisms. This finding may provide novel targets for future studies into the etiology of the condition and a new way to fractionate the disorder.
Research agrees that ASD is a highly heterogeneous condition with multiple causes (ie, complex genetics45) and various phenotypes, which makes the neuroanatomy of ASD inherently difficult to describe. Although several autistic core structures have been highlighted by previous studies, growing evidence suggests that ASD is a neural systems condition characterized by subtle differences in large-scale neural systems rather than isolated regions with large effects.32 Herein, we first demonstrated that individuals with ASD have significant differences in CT and CV—mainly located in frontal and temporal regions—that significantly correlated with measures of symptom severity. However, these clusters of significant between-group differences were caused by subtle and spatially distributed variations in CT and SA. Our results therefore agree with the notion that the neuroanatomy of ASD is not confined to individual brain regions but rather affects wider neural systems, which may be difficult to detect using mass-univariate techniques (also see the article by Ecker et al46).
Second, we found that the spatially distributed patterns of differences in CT and SA displayed statistically independent sources of variability (ie, were largely nonoverlapping) and thus potentially reflect different neuropathological processes. In terms of phylogeny, it is now widely believed that CT and SA originate from different types of progenitor cells, which divide in the ventricular zone to produce glial cells and neurons. Cortical thickness has been related primarily to intermediate progenitor cells (neurogenic transient amplifying cells in the developing cerebral cortex),47 which divide symmetrically at basal (nonsurface) positions of the ventricular surface and only produce neurons.48,49 These neurons then migrate along radial glial fibers to form ontogenetic columns (ie, radial units). According to the radial unit hypothesis,14 CT depends on the neuronal output from each radial unit—amplified by intermediate progenitor cells—and therefore reflects the number of neurons produced in each unit. On the other hand, SA has mainly been related to radial unit progenitor cells, which divide at the apical (ventricular) surface. The early proliferation of radial unit progenitor cells leads to an increase in the number of proliferation units, which in turn results in an increase in SA.47 In other words, SA is related to the number of ontogenetic columns. Therefore, the spatially distributed nonoverlapping patterns of differences in CT and SA in adults with ASD most likely reflect the end result of different phylogenetic processes and particularly affect frontal and temporal regions.
The existence of distinct phylogenetic processes for CT and SA also implies that at least 2 different genetic (or other) mechanisms may be involved in their etiology and regulation. The genetic dissociation of CT and SA has previously been explored in both animal and human studies. For instance, in the mouse, mutation of the genes PAX6, LRP6, and NGN1/2 modifies the abundance of intermediate progenitor cells and results in parallel increases in CT but not SA.50 In humans, mutations in PAX651 or TBR252 are associated with a reduction in CT relative to SA. In contrast, SA but not CT is modulated by variations in MECP2 within specific cortical regions (eg, cuneus, fusiform gyrus).50 A direct link between these genes and the neuroanatomical phenotype of ASD remains to be established; however, there is preliminary evidence that some of these genes are associated with ASD.53- 55 The spatially distributed pattern of differences for CT and SA may therefore reflect distinct genetic etiologies—as ASD is a neurodevelopmental disorder—and offer new targets for further exploration of the genetic mechanisms leading to ASD. Future studies of cortical anatomy in ASD will therefore need to measure more than CV alone to identify the potentially different etiological factors leading to ASD. In addition, this study highlights the need for future histological studies to identify the specific cytoarchitectonic correlates of ASD as these are not directly accessible using MRI.
We further observed that most differences in CV (67%) were driven by differences in SA rather than CT. This finding agrees with previous studies demonstrating that interindividual variation in brain volume in healthy adults is driven predominantly by differences in SA.56- 58 Further, Hazlett et al25 recently reported that brain enlargement observed in toddlers with ASD may be associated with increased SA but not CT. However, as SA is a 2-dimensional measure (ie, millimeters squared) and CT is a 1-dimensional measure (ie, millimeters), differences in SA will naturally have a stronger impact on measures of CV. This implies that measures of brain volume are generally more similar to measures of SA and that CT should be explored in isolation. Thus, while our study suggests that differences in brain volume in ASD may primarily be caused by differences in SA, individuals with ASD also have significant differences in CT. These may be more difficult to measure in individuals with ASD as the delineation between the gray and white matter boundary is more difficult owing to the presence of supernumerary neurons beneath the cortical plate.59 Another factor that may influence differences in CT and SA is gyrification (ie, cortical folding). In heavily folded brain areas, vertex-based measures of cortical expansion may be decreased as more vertices are needed to describe the more complex regional geometry. In addition, increased gyrification has been associated with increased CT.60,61 Thus, future studies are needed to establish the relationship between volumetric measures and gyrification in ASD.
Regional differences in CT and CV were predominantly observed in frontotemporal regions. This is consistent with prior reports suggesting that people with ASD have differences in frontal lobe neuronal integrity,62 function,63- 65 anatomy,66- 69 and connectivity.7 Furthermore, it has been suggested by some that individuals with ASD have a delay in frontal lobe maturation70 and that abnormalities in frontal lobe development71 may underpin some of the social impairments reported in people with ASD (eg, deficits in social cognition72,73). Our findings and the work of others thus suggest a high degree of neuroanatomical vulnerability in frontal regions that may be modulated by (mainly) nonoverlapping differences in both SA and CT. If we accept that CT and SA reflect different aspects of the neural architecture, our finding suggests that abnormalities in surface anatomy of specific frontal regions in ASD may be caused by different neuropathological processes.
The same applies for the temporal lobe, where individuals with ASD had significantly thinner cortices than controls, mainly in the middle and superior regions of the anterior temporal lobe. Gray matter differences in temporal regions have previously been reported in ASD during childhood,74,75 adolescence,68 and young adulthood6,69 and have functionally been linked to processing of biological movement (ie, movement characteristic of living organisms), social perceptions, and theory of mind tasks,76- 79 which have all been reported as affected in ASD.80 Moreover, based on the radial unit hypothesis, an increase in frontal CT in ASD would imply supernumerary neurons within minicolumns, whereas a regional decrease of SA (eg, in the orbitofrontal cortex) may be related to fewer or narrower units. However, these conclusions remain to be validated by future histological studies as existing postmortem studies in ASD are rare and available samples are small.81,82 Furthermore, most postmortem studies are restricted to few and small regions of the cortex, so it is difficult to generalize findings across the brain.
Our study raises a number of methodological issues. First, we investigated neuroanatomy in a sample of high-functioning men diagnosed with the ADI-R as having ASD. As ASD is a spectrum condition, our sample therefore represents a specific subpopulation of the autistic phenotype and results might not generalize to other groups on the autism spectrum (eg, individuals with intellectual disability) or to females with the condition. Second, a multicenter design was used for MRI data acquisition to overcome single-site recruitment limitations. However, we used a recently developed acquisition protocol that standardizes structural MRI data across multiple platforms31 and also accounted for intersite effects in the statistical model.83,84 Therefore, the detected between-group differences cannot be fully explained by these limitations. Third, the data build on a previous study published recently by our group using voxel-based morphometry.32 Although both techniques lead to similar results in frontal regions, there is little overlap between results in the temporal and occipital lobes. Such discrepancies have been noted previously85,86 and most likely reflect different approaches to image normalization and registration, which may lead to differential spatial patterns of variance and thus to differential spatial sensitivity to a fixed effect size. Between-group tests on the basis of different volume metrics are hence conditioned by the variance and consequently yield to different patterns of results.83 Also, the volumetric measures investigated in the present study were based on measures of areal expansion or compression (relative to a template) and hence do not provide absolute measures of regional brain volume. Last, to our knowledge, this is the first study investigating regional measures of cortical SA in individuals with ASD using a spatially unbiased vertex-based approach and their relationship to differences in CT. Hence, the results should be considered exploratory (ie, proof of concept), requiring future replication.
To conclude, differences in CV observed in individuals with ASD may be underpinned by CT and SA, which exhibited statistically independent sources of variability. These 2 factors are likely to be the result of distinct developmental pathways that are modulated by different neurobiological mechanisms. Both CT and SA would thus benefit from being explored in isolation to elucidate the etiology and neurobiology of ASD.
Correspondence: Christine Ecker, MSc, PhD, Department of Forensic and Neurodevelopmental Sciences, PO50, Institute of Psychiatry, King's College London, De Crespigny Park, London SE5 8AF, England (email@example.com).
Submitted for Publication: October 11, 2011; final revision received February 14, 2012; accepted March 30, 2012.
Published Online: November 26, 2012. doi:10.1001/jamapsychiatry.2013.265
MRC AIMS Consortium Members: Anthony J. Bailey, MD, PhD, Simon Baron-Cohen, BA, PhD, FBPsC, MPhil, Patrick F. Bolton, PhD, FRCPsych, Edward T. Bullmore, MD, PhD, FRCPsych, FMedSci, Sarah Carrington, PhD, Bishmadev Chakrabarti, PhD, Eileen Daly, BA, Sean C. Deoni, PhD, Christine Ecker, BSc, MSc, PhD, Francesca Happe, PhD, Julian Henty, PhD, Peter Jezzard, PhD, Patrick Johnston, BSc, PhD, Derek K. Jones, PhD, Meng-Chuan Lai, MD, PhD, Michael V. Lombardo, BA, PhD, Anya Madden, BSc, Diane Mullins, MD, Clodagh M. Murphy, FRCPsych, Declan G. M. Murphy, MBBS, FRCPsych, MD, Greg Pasco, PhD, Susan Sadek, MD, Debbie Spain, BSc, Rose Steward, BSc, John Suckling, PhD, Sally Wheelwright, BSc, and Steven C. Williams, BSc, PhD.
Conflict of Interest Disclosures: Dr Bullmore is employed half time by GlaxoSmithKline and holds shares in GlaxoSmithKline.
Funding/Support: The MRC AIMS Consortium is a collaboration of autism research centers in the United Kingdom including the Institute of Psychiatry, King's College London; the Autism Research Centre, University of Cambridge; and the Autism Research Group, University of Oxford. It is funded by the Medical Research Council and headed by the Department of Forensic and Neurodevelopmental Sciences, Institute of Psychiatry. This work was supported by the AIMS Consortium through grant G0400061 from the Medical Research Council and by the NIHR Biomedical Research Centre for Mental Health at King's College London, Institute of Psychiatry, and South London & Maudsley NHS Foundation Trust. Scanning took place in 3 centers (London, Oxford, and Cambridge; local principal investigators were Drs Bailey, Baron-Cohen, Bullmore, and D. G. M. Murphy).
Additional Contributions: We are grateful to those who agreed to be scanned and who gave their time so generously to this study.