- Open Access
Gray matter covariations and core symptoms of autism: the EU-AIMS Longitudinal European Autism Project
Molecular Autism volume 11, Article number: 86 (2020)
Voxel-based morphometry (VBM) studies in autism spectrum disorder (autism) have yielded diverging results. This might partly be attributed to structural alterations being associating with the combined influence of several regions rather than with a single region. Further, these structural covariation differences may relate to continuous measures of autism rather than with categorical case–control contrasts. The current study aimed to identify structural covariation alterations in autism, and assessed canonical correlations between brain covariation patterns and core autism symptoms.
We studied 347 individuals with autism and 252 typically developing individuals, aged between 6 and 30 years, who have been deeply phenotyped in the Longitudinal European Autism Project. All participants’ VBM maps were decomposed into spatially independent components using independent component analysis. A generalized linear model (GLM) was used to examine case–control differences. Next, canonical correlation analysis (CCA) was performed to separately explore the integrated effects between all the brain sources of gray matter variation and two sets of core autism symptoms.
GLM analyses showed significant case–control differences for two independent components. The first component was primarily associated with decreased density of bilateral insula, inferior frontal gyrus, orbitofrontal cortex, and increased density of caudate nucleus in the autism group relative to typically developing individuals. The second component was related to decreased densities of the bilateral amygdala, hippocampus, and parahippocampal gyrus in the autism group relative to typically developing individuals. The CCA results showed significant correlations between components that involved variation of thalamus, putamen, precentral gyrus, frontal, parietal, and occipital lobes, and the cerebellum, and repetitive, rigid and stereotyped behaviors and abnormal sensory behaviors in autism individuals.
Only 55.9% of the participants with autism had complete questionnaire data on continuous parent-reported symptom measures.
Covaried areas associated with autism diagnosis and/or symptoms are scattered across the whole brain and include the limbic system, basal ganglia, thalamus, cerebellum, precentral gyrus, and parts of the frontal, parietal, and occipital lobes. Some of these areas potentially subserve social-communicative behavior, whereas others may underpin sensory processing and integration, and motor behavior.
Autism spectrum disorder (henceforth autism) is an early onset neurodevelopmental condition characterized by core deficits in social interaction and communication, along with restrictive interests and behavior, and sensory abnormalities . Magnetic resonance imaging (MRI) studies have increased our understanding of the neuroanatomical underpinnings of autism and show that autism is associated, at the group level, with brain structural changes . However, many results are not robust across different studies. For example, two studies using the same large-scale open access Autism Brain Imaging Data Exchange (ABIDE) dataset  came to different conclusions with regard to the volume of the pallidum [4, 5]. Also, across whole brain approaches investigating cortical (i.e., cortical thickness and surface area) and subcortical (i.e., volume) features have been inconsistent; two large-scale pooled estimate analytical studies observed diverging results of cortical changes in autism [6, 7]. Similarly, autism studies quantifying voxel-wise gray matter (GM) density also found divergent results of GM differences between autism diagnosed and control individuals; for instance, meta-analyses reported diverse changes of GM morphometry in autistic individuals on average, reporting either increased or decreased density of right inferior temporal gyrus in autism [8, 9]. Even when taking age into account, studies still observed different structural brain alterations in children and adolescents with autism [10, 11].
A commonality to all these studies is their reliance on mass-univariate statistics. This approach identifies alterations in isolated regions or voxels but ignores possible relationships between them. The brain is a complex system of interconnected networks, and research into the neural basis of autism has moved away from focusing on local abnormalities into conceptualizing autism as a disorder of alterations in structural and functional brain connectivity . This implies that structural brain alterations in autism likely reflect the combined influence of several regions and are not confined to one specific region [13, 14]. The present paper aims to advance prior work on brain structural neural correlates of autism in two ways. First, we aim to move away from the standard univariate approach and incorporate an alternative that adheres more closely to the hypothesis of autism as a disconnection syndrome , thus providing greater sensitivity for between-group effects. For this purpose, we identify inter-regional sources of structural covariation using independent component analysis (ICA) , a data-driven unsupervised approach that allows the identification of interconnected brain regions across the whole brain. It has previously been applied successfully to identify covariance of brain morphometry in patients with psychiatric disorders [16,17,18]. Second, we move beyond the categorical autism case–control comparison towards exploring associations between brain structure and symptom dimensions or profiles of autism. Although former studies have used univariate approaches to explore the relationship between brain substrates and clinical phenotypes [6, 19], such associations are potentially the consequence of integrated effects across multiple symptoms dimensions and brain regions, rather than simple associations between a specific brain region and a specific symptom dimension. To study such multidimensional associations multivariate methods are effective [20, 21] and here we achieve this integration using Canonical Correlation Analysis (CCA) .
In summary, we investigate alterations in GM morphometric covariations in a deeply phenotyped large European autism case–control sample [23, 24] that allows us to improve our understanding of correlated structural brain alterations in autism. Subsequently, we focus on the covariation between the identified structural features and symptom behavior profiles among individuals with autism.
The participants were selected from the first wave of the European Autism Interventions—A Multicentre Study for Developing New Medications (EU-AIMS) Longitudinal European Autism Project (LEAP) dataset, which is a large multicenter study that aims to identify and validate biomarkers for autism . In total, six centers are involved: Institute of Psychiatry, Psychology and Neuroscience, King’s College London, United Kingdom; Autism Research Centre, University of Cambridge, United Kingdom; Radboud University Medical Centre, Nijmegen, the Netherlands; University Medical Centre Utrecht, the Netherlands; Central Institute of Mental Health, Mannheim, Germany; and University Campus Bio-Medico, Rome, Italy. Each participant underwent clinical, cognitive, and MRI assessment. Autism diagnoses were confirmed by clinicians according to the Diagnostic and Statistical Manual-IV (DSM-IV), International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10), or DSM-5. The study was approved by local ethical committees in each participating center, and written informed consent was provided by all participants and/or their legal guardians (for those < 18 years old). For further details on experimental design and clinical assessments, see [23, 24].
In the present study, we selected participants with available structural MRI data. All images were inspected visually and participants were excluded in cases of brain injury or structural abnormalities (e.g., enlarged ventricles or cysts), excessive head motion, or preprocessing failure (n = 29). We excluded the participants from the Rome site due to low sample size (n = 1). We also excluded participants without full-scale intelligence quotient (FSIQ) data in the further statistical analyses (n = 5). This resulted in a sample of 599 participants from 5 sites, including 347 individuals with autism and 252 typically developing (TD) controls. Demographic and clinical information is shown in Table 1.
We used the Autism Diagnostic Interview-Revised (ADI)  and the Autism Diagnostic Observational Schedule 2 (ADOS)  to quantify past (ever and previous 4-to-5 years) and current autism symptoms of the social interaction, communication, and restricted repetitive behaviors (RRB) domains. We used T-scores (age- and sex-adjusted) of the Social Responsiveness Scale 2nd Edition (SRS)  in the autism group to assess severity of autistic traits/symptoms and the Repetitive Behavior Scale-Revised (RBS)  to measure repetitive and rigid behaviors associated with autism. Moreover, sensory processing abnormalities of autism were assessed with the Short Sensory Profile (SSP) . To examine associations between clinical features in autism and brain measures, we created two sets of clinical measures: (1) the subscale scores of ADI-R and ADOS, both instruments were rated by qualified examiners, and (2) the total scores of SRS, RBS, and SSP, we included parent-rated reports only and limited the analyses to within the autism group. Further, concerning the potential effect of comorbidity with Attention Deficit Hyperactivity Disorder (ADHD), we included comorbidity with ADHD as an additional covariate in analyses. ADHD symptoms were assessed with the ADHD DSM-5 rating scale that includes symptom scales of inattention and hyperactivity/impulsivity symptoms. The ADHD DSM-5 rating scale was based on either parent-report or self-report scores; self-report scores were only used when parent-reports were unavailable. The categorical output of the ADHD rating scale was used in this study. The summary for each of these clinical measures can be found in Table 1.
MRI data acquisition
All participants were scanned on 3 T MRI scanners (University of Cambridge: Siemens Verio; King’s College London: GE Medical Systems Discovery MR 750; Mannheim University: Siemens TimTrio; Radboud University: Siemens Skyra; Rome University: GE Medical Systems Sigma HDxTt; Utrecht University: Philips Medical Systems Achieva/Ingenia CX). High-resolution structural T1-weighted images were acquired with full head coverage, at 1.2 mm thickness with 1.1 × 1.1 mm in-plane resolution. For all other scanning parameters, please see Additional file 1: Subsection 1.
GM density estimation
Voxel-based morphometry (VBM) is a spatially unbiased whole-brain approach that extracts voxel-wise GM density (the amount of GM at a voxel) estimations. We performed VBM analyses using the CAT12 toolbox  in SPM12 (Wellcome Department of Imaging Neuroscience, London, UK). T1-weighted images were automatically segmented into GM, white matter, and cerebrospinal fluid and affine registered to the MNI template to improve segmentation. All resulting segmented GM maps were then used to generate a study-specific template and registered to MNI space via a high-dimensional, nonlinear diffeomorphic registration algorithm (DARTEL) . A Jacobian modulation step was included using the flow fields to preserve voxel-wise information on local tissue volume. Images were smoothed with a 4 mm full-width half-max (FWHM) isotropic Gaussian kernel.
A full quality control report was generated by the CAT SPM pipeline for each participant that included visualizations of the segmentation, which were evaluated visually, and quantitative quality measures including mean correlation from sample homogeneity module and weighted overall image quality ratings that were additionally used to detect and exclude images of insufficient quality for inclusion in analysis. We visually checked the images with the mean correlation smaller than three standard deviations from the sample mean. Accordingly, 5 participants required visual inspection after preprocessing, whereas none of them were observed obvious artifacts, and they were consequently included in the final analyses.
Structural ICA decomposition
All participants’ VBM data were simultaneously decomposed into 100 spatially independent sources of spatial variation using MELODIC-ICA . Such ICA decomposition provides, at each independent component, a brain map reflecting a pattern of GM density covariation across participants, and a participant’s loading vector reflecting the contribution of each participant to each component. Higher model orders (i.e., number of components) tend to extract spatially sparser maps, and hence participant loading vectors variance would be explained by smaller subsets of participants. Therefore, 100 components were chosen to capture as much variation as possible while remaining statistically powered. However, the choice of model order is heuristic, we thereby examined the dependence on model order using different dimensionalities; more precisely, in addition to the 100 dimensional factorization, we also considered an automatic dimension estimation approach as implemented in MELODIC-ICA and a 50 dimensional independent component factorization.
For completeness we first performed a standard mass-univariate statistical analyses directly on the GM densities. To that end we used a generalized linear model (GLM) to detect group differences (autism vs. TD) on GM densities using the FMRIB Software Library v6.0 (FSL) . Participants’ VBM maps were considered as the dependent variable and diagnostic group as the independent factor, with age, sex, FSIQ, and scan site as covariates. Significance was assessed using permutation testing (5000 permutations) and correction for multiple comparisons was achieved using Threshold-Free Cluster Enhancement (TFCE, two-tailed, threshold at p < 0.05) [34, 35].
Next, we considered the results of the ICA factorization of the VBM data. A GLM was used to examine differences between autistic and TD individuals, by using each participant's loading to each component as a dependent variable, diagnostic group as independent variable, and age, sex, FSIQ, and scan site as regressors. To avoid the results of case–control differences being biased by IQ, we repeated the same procedure by excluding participants in Schedule D (FSIQ < 75, more details see Additional file 1: Subsection 2). Considering that autism is highly comorbid with ADHD , we also controlled for comorbidity with ADHD by adding a dummy-coded variable (with/without ADHD) to the original GLM analyses (autism vs. TD) in the case–control ICA analysis. Additionally, we also conducted analyses that included a dimensional score (instead of categorical code) of the two ADHD subscales as an additional covariate in the group effect analyses (see Additional file 1: Subsection 3). Furthermore, we separately checked the effects of age squared, age-by-group, age squared-by-group, sex-by-group interactions, anxiety, depression comorbidity, medication use and image quality on case–control differences of structural covariance. The results can be found in Additional file 1: Subsection 3. Multiple comparison correction was implemented using false discovery rate (FDR) (p < 0.05) .
We also explored independently the relationships between each estimated brain component and subscales of ADI and ADOS, SRS, RBS, and SSP in the autism group using GLM analyses and again correction for multiple comparisons was implemented with the FDR method (p < 0.05). Then, to simultaneously explore the relationship between all the brain structural phenotypes estimated through ICA and all symptom phenotypes in the autism group, we used CCA. In the considered scenario, CCA is able simultaneously to learn linear projections (via canonical coefficients) of the brain structural sources and the behavioral measures that maximize the correlation between them at the participant level. Canonical variates are generated independently for the brain and the behavioral datasets according to the product of the canonical coefficients (learned through CCA) and the original datasets. We referred to each pair of canonical variates as CCA mode . The canonical coefficients separately represent weights for each brain variable and each symptom measure. For interpretation of the contribution of each independent source and each clinical measure to the CCA mode, their canonical coefficients were corrected as noted in . In this work we focused on the first canonical mode since it provides the strongest associations and we denoted it as the main CCA mode. Furthermore, for each CCA analysis the statistical significance of the main CCA mode was determined by permutation testing (10,000 permutations, Bonferroni corrected p < 0.05/number of CCA modes). Here, we performed two separate CCA analyses to link the independent components participants’ contributions to subsets of behavioral measures; in the first CCA analyses (CCA1) we included the subscales of ADI and ADOS as clinical measures and in the second (CCA2) we used total scores of SRS, RBS, and SSP. The reliability of the CCA results presented as well as its dependence on the number of participants were tested using a leave-one-out cross-validation approach (Additional file 1: Subsection 4).
The standard mass-univariate GLM analysis of the VBM data comparing cases and controls did not show significant group differences for voxel-wise GM volumes after multiple comparison correction. We provide the uncorrected statistical results (p < 0.05) in Additional file 1: Subsection 5.
Group effect on ICA decomposition
The structural data ICA decomposition provided a set of 100 independent spatial sources, each of which is connected to a vector that depicts the degree of each participant’s contribution to the corresponding components. For clarity, we further refer to these vectors as the participant loadings. Post-hoc GLM analyses of these participant loadings showed case–control differences at nine independent components (ICs) (p < 0.05, i.e., IC10, IC13, IC14, IC15, IC23, IC28, IC31, IC48, and IC 99, see Additional file 1: Subsection 6). Of these, two components, IC10 (\(\beta = - 0.147\), p = 8.850 × 10–5, effect size [Cohen’s d] d = − 0.358) and IC14 (\(\beta = - 0.132\), p = 5.450 × 10–4, d = − 0.321), survived multiple comparison correction (FDR corrected, p < 8.072 × 10–4). These results were not driven by age, sex, or scan site.
In Fig. 1, we present summary images reflecting the brain areas involved in the structural variances occurring at these two components. The top row of Fig. 1 shows that IC10 primarily relates to structural variation in the bilateral insula, inferior frontal gyrus (IFG), orbitofrontal cortex (OFC), and caudate nuclei. Among these brain regions, the bilateral caudate exhibits alterations in the opposite direction to the others. Given the negative beta coefficient obtained from the GLM analysis between participant loadings at IC10 and the diagnosis group labels, individuals with autism demonstrate increased GM densities in the bilateral caudate and decreased densities in the bilateral insula, IFG, and OFC. The bottom row of Fig. 1 shows that IC14 mainly involves variations in the bilateral amygdala, hippocampus, and parahippocampal gyrus (PHG). Similarly, according to the sign of the beta values obtained through the GLM, the autism group shows decreased densities in the areas involved in IC14.
The robustness of the ICA results to the model order choice was evaluated by considering, in addition to the original 100-dimensional factorization, an automatic dimensionality estimation procedure resulting in a 91-dimensional factorization, and a 50-dimensional factorization. We observed that the main components reported (IC10 and IC14) are highly reproducible independent of the model order choice. For details, see Additional file 1: Subsection 7.
To validate the ICA results not being biased by low IQ participants, an additional validation was performed by taking FSIQ into account to exclude the participants in Schedule D from the ICA factorization. This showed that a unique IC, corresponding to the original IC10, survived FDR correction (see Additional file 1: Subsection 3). Further, the effect of comorbidity with ADHD on brain structural variations was determined using data from 500 participants (for detailed demographic information, see Additional file 1: Subsection 3). This analysis showed that IC14 remained significantly associated with the autism group (FDR corrected, p = 9.669 × 10–4). However, IC10 was no longer associated with autism (p = 0.004).
Further, post-hoc GLM analyses of the relationships between brain ICs and symptom ratings did not provide any significant associations (Additional file 1: Subsection 8).
Relating gray matter spatial variation patterns to symptoms profiles
As a final step, we applied CCA to examine the associations between the 100 components and the two sets of clinical measures among the autism cases only. The CCA1 (linking ADI and ADOS subscale scores to brain measures), involved 325 autism participants and showed a Bonferroni corrected (p = 0.05/5 = 0.010) significant relationship (Fig. 2a,c, r = 0.701, permutation p = 0.008). In this main CCA mode, IC16, IC61, IC89, and IC14 were the highest contributors to the correlation with autism symptoms, and the ADOS RRB subscale loaded most on the association with the brain measures (Fig. 3a, b). Among the four components, IC16 mainly involved density variations in bilateral thalamus and putamen (canonical weight: 0.447), IC61 in right lateral occipital and left superior parietal lobe (canonical weight: − 0.366), and IC89 in the right precentral gyrus (canonical weight: − 0.333). For details, see Additional file 1: Subsection 9. Note that IC14 is among the components previously reported showing linear significant case–control group effects. The regions involved in IC14 were mentioned above (canonical weight: − 0.312). Since higher scores of the ADI and ADOS reflect more severe autism symptoms, positive values of IC16 suggest that higher loading on this component is related to more severe symptoms in autism, and negative values of IC61, IC89, and IC14 meant that lower loadings on these three ICs are associated with more severe symptoms. In Fig. 2a, participants were color coded according to their ADOS-RRB scores to illustrate how the ADOS-RRB score drives the canonical correlation.
In CCA2 we linked SRS, RBS, and SSP scores to the brain measures of 194 individuals with autism, which is 55.9% of all participants with autism (lower number due to missing questionnaire data). We found a Bonferroni corrected (p = 0.05/3 = 0.017) significant relationship (Fig. 2b, r = 0.840, permutation p = 0.002, Fig. 2d). In this main CCA mode, IC82, IC99, and IC100 were the highest contributors to the correlation with behavior profiles, and SSP score loaded most on the association with the brain measures in the autism group (Fig. 3c, d). IC82 mainly involved variations in the bilateral cerebellum (canonical weight: 0.414), IC99 in the left lateral occipital and parietal lobe, and bilateral precentral gyrus (canonical weight: 0.277), and IC100 in the left inferior frontal gyrus and right middle frontal lobe (canonical weight: 0.262). For details, see Additional file 1: Subsection 9. Similarly, lower loadings on these three ICs were related to more severe symptoms. In Fig. 2b, each participant was color coded according to their SSP score, and it shows how SSP score drives the correlation. In this case, both IC10 and IC14 ranked outside the top 20 of the 100 components, suggesting that these two components with significant case–control difference have no strong contribution to the CCA2 correlation. However, for completeness, direct interpretation (referring to uncorrected coefficients) of the CCA2 weights ranks IC14 as the third strongest contributor to this canonical correlation (Additional file 1: Subsection 10).
The CCA robustness analyses indicated that the main CCA modes of both CCA analyses were reliably estimated in a leave-one-subject out setting (Additional file 1: Subsection 4). In CCA1, the weights of the main CCA mode of each leave-one-out analysis correlated on average above 0.94 with the weights of original main CCA mode in brain loadings and above 0.95 in behavior phenotypes when the sample was bigger than 122 participants. In CCA2, the weights of the main CCA mode related on average above 0.92 in brain loadings and above 0.96 in behavior profiles when the sample was bigger than 111 participants. Both CCA analyses are no reproducible for sample sizes smaller than (approximately) 100 participants.
The present study used a reliable approach to quantify inter-individual differences in GM morphometry covariations in a deeply phenotyped large sample of individuals with and without autism. The standard, univariate VBM analysis did not show significant case–control differences. We then utilized an ICA decomposition of all participants GM density images, and similarly performed a case–control post-hoc statistical analyses. This analysis showed that autism was significantly associated with alterations in two independent sources of GM density covariations. These findings corroborated our hypothesis that alterations in brain morphometry in autism are associated with the combined influence of several regions rather than with a single region. In a following step, we applied CCA to explore multivariate associations between sets of continuous measures of core symptoms and sets of ICA-derived morphometry measures within the autism group, and were able to identify significant relationships between brain components and symptom profiles. Notably, one of the components which showed significant case–control differences was also among the highest loading components in the CCA.
We found nonsignificant case–control difference in univariate analysis in current study, which is inconsistent with previous studies (e.g., [8, 9]). In light of the substantial biological heterogeneity among autistic individuals, it is expected that findings differ across datasets, especially those of smaller size. Specifically, the current study included the individuals with a relative wide age and IQ range, and did not exclude co-occurring psychiatric symptoms. Additionally, integrating the findings from multivariate analyses in our study, the absence of diagnostic differences at single region/voxel underscores the detection sensitivity of group effect of a multivariate approach which evidently verified our hypothesis.
Our findings showed two covarying sets of brain areas that structurally differed between cases and controls. While one source of GM density covariation, IC10, mainly related to the bilateral insula, IFG, OFG, and caudate, another source, IC14, primarily involved the bilateral amygdala, hippocampus, and PHG. The brain regions within each component are anatomically clustered and symmetrical, which indicates that the independent structural covariation alteration in the GM of individuals with autism is concentrated in nearby brain areas. This is in line with a previous study that used a similar approach . It is further in line with organizing principles of the brain that regions tend to be more interconnected when they are located close to each other [41, 42]. However, when we compared the regions loading on the two components, the covarying regions of each component distribute in different brain locations. This suggests that neuroanatomic alterations underlying autism are more widely distributed at the whole brain level. It is of note that, when accounting for ADHD comorbidity, IC14 remained significant but IC10 did not. This suggests that IC14 is more specifically related to autism associated structural variations, even after linearly accounting for ADHD effects, while IC10 might reflect variations associated with both autism and ADHD.
The brain regions with high loadings on either of these two components, i.e., insula, amygdala, hippocampus and PHG have lower densities in autism and have earlier been associated with autism [9, 43]. The opposite direction of the alteration of the caudate nucleus in autism has also previously been found . This is not the case for the IFG and OFG, which showed lower densities in autism in our study, where prior studies found mixed results [45, 46]. Importantly, the brain regions identified by our analyses have earlier been implicated in the neurobiology and/or neurocognition of autism. In IC10, structural and/or functional alterations of the insula, IFG, and OFC have been associated with social and non-social cognitive impairments in autism [46,47,48,49]. A meta-analysis reported abnormal functional activations of the insula, IFG, and OFG during social cognition tasks in autism . Additionally, variance of the caudate nucleus volume was found to correlate with the severity of RRB symptoms in autism . Together with deviant structural and functional connectivity between frontal cortical areas and striatum in autism [47, 51, 52], structural covariation in striatum and frontal areas may underlie atypical functional fronto-striatal connectivity, and this has been associated with repetitive behavior and executive functioning impairments in autism [3, 45]. In the present study, the density of caudate nuclei increase simultaneously with densities decreasing on other areas in autism, which fits with the results of a few functional studies that indicate inverse functional changes of these areas . Particularly, the special pattern of GM densities changes in frontal and striatal areas might serve an important role in autism-related symptoms.
In IC14, we found decreased densities of amygdala and hippocampus, where the structural alterations have previously been related to social deficits in autism [6, 54]. The amygdala, hippocampus, and PHG subserve cognitive and emotional functions that were found abnormal in individuals with autism [50, 55, 56]. In addition to being involved in emotion and face processing, the three areas have been proposed as structures critical for working memory in autism . Furthermore, these cognitive domains exert bidirectional effects on each other, with atypical social-emotional processing influencing memory performance in individuals with autism, and memory being involved in complex information processing and executive functioning, which in turn affects social cognition [57, 58]. In sum, given the potential functional interactions between these three brain areas, the structural covariance alterations of the amygdala, hippocampus, and PHG found in our study, may underlie or contribute to abnormal functional connections of these areas, and thus underlie poor performance on social cognition and memory tasks in individuals with autism.
Our multivariate correlation analyses moved from the case–control comparison to the use of continuous symptoms among individuals with autism and identified two prominent relationships between all structural brain covariances and symptoms in autism. Three of the four brain components that ranked top in this analysis did not show case–control differences, while there was one component (IC14) that differed between cases and controls and also significantly correlated with the severity of autism symptoms assessed by ADI and ADOS. The brain areas loading high on the brain components identified in the CCA are somewhat different from those implicated in the case–control analyses. These former brain areas are the thalamus, putamen, precentral gyrus, and parietal and occipital lobes in CCA1, and the cerebellum, frontal lobe, and again precentral gyrus, and parietal and occipital lobes in CCA2. These are foremost areas of the brain implicated in the processing and higher order integration of sensory information and motor functions. This makes sense since repetitive, rigid and stereotyped behaviors and abnormal sensory behaviors seem to drive the brain-behavior associations much more than the measures on social-communication symptoms. Note that variance within the different autism symptom domains (social-communication, repetitive behaviors and sensory abnormalities) was similar and cannot explain the differential symptom-brain associations.
Overall, the results of our multivariate analyses on case–control differences and on continuous measures of symptom severity among those with autism demonstrate the complexity of the brain morphometry correlates of autism. Brain areas involved are scattered across the whole brain and include the limbic system, basal ganglia, thalamus, cerebellum, precentral (motor) gyrus, and parts of the frontal, parietal, and occipital lobes.
Strengths and limitations
The strengths of our study are the use of a large deeply phenotyped sample, bottom-up data-driven analyses, a multivariate approach for examining brain-symptom associations, and a large set of continuous measures of core autism symptoms. There are several limitations in our study. First, only 55.9% of the autism group had complete questionnaire data on continuous parent-reported symptom measures, which may have lowered statistical power for this analysis. Second, due to the clinical characteristic of autism, the participants showed significant differences in the proportion of sex distribution. Despite this we regressed out the sex effect in the group difference analyses, future studies may consider balancing the sex distributions in case and control groups; however this may be difficult to achieve. Third, albeit we ran a series of sensitivity and control analyses, the effect of other potential sources of variance, such as the complex effects of medication use that may influence our findings.
We demonstrate brain morphometry differences between individuals with autism and typical controls in the inter-regional covariation of the insula, frontal area, caudate, amygdala, hippocampus, and PHG. Further, we highlight associations between covariation in density of the thalamus, putamen, precentral gyrus, frontal, parietal, and occipital lobes, and the cerebellum, and core autism symptoms, in particular repetitive behaviors and abnormal sensory behavior. Future studies may link our morphometry findings with data on brain function obtained from cognitive tests and/or functional and resting-state MRI, and with genomics data.
Availability of data and materials
Data collected in EU-AIMS LEAP are stored and curated at the central EU-AIMS database at the Pasteur Institute in Paris. The database is open to members of the wider scientific community upon request and submission of a paper and data analytic proposal.
Attention Deficit Hyperactivity Disorder
Autism Diagnostic Interview-Revised
Autism Diagnostic Observational Schedule 2
Canonical correlation analysis
Diagnostic and Statistical Manual-IV
European Autism Interventions—A Multicentre Study for Developing New Medications
False discovery rate
Full-scale intelligence quotient
Generalized linear model
Independent component analysis
International Statistical Classification of Diseases and Related Health Problems 10th Revision
Inferior frontal gyrus
Magnetic resonance imaging
Repetitive Behavior Scale-Revised
Restricted and repetitive behaviors
Social Responsiveness Scale 2nd Edition
Short Sensory Profile
Threshold-Free Cluster Enhancement
American Psychiatric Association. Diagnostic and statistical manual of mental disorders. 5th ed. Arlington: Author; 2013.
Ecker C, Bookheimer SY, Murphy DG. Neuroimaging in autism spectrum disorder: brain structure and function across the lifespan. Lancet Neurol. 2015;14(11):1121–34.
Di Martino A, Yan CG, Li Q, Denio E, Castellanos FX, Alaerts K, et al. The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. Mol Psychiatry. 2014;19(6):659–67.
Turner AH, Greenspan KS, van Erp TGM. Pallidum and lateral ventricle volume enlargement in autism spectrum disorder. Psychiatry Res Neuroimaging. 2016;252:40–5.
Zhang W, Groen W, Mennes M, Greven C, Buitelaar J, Rommelse N. Revisiting subcortical brain volume correlates of autism in the ABIDE dataset: effects of age and sex. Psychol Med. 2018;48(4):654–68.
van Rooij D, Anagnostou E, Arango C, Auzias G, Behrmann M, Busatto GF, et al. Cortical and subcortical brain morphometry differences between patients with autism spectrum disorder and healthy individuals across the lifespan: results from the ENIGMA ASD working group. Am J Psychiatry. 2018;175(4):359–69.
Bedford SA, Park MTM, Devenyi GA, Tullo S, Germann J, Patel R, et al. Large-scale analyses of the relationship between sex, age and intelligence quotient heterogeneity and cortical morphometry in autism spectrum disorder. Mol Psychiatry. 2019;25(3):614–28.
DeRamus TP, Kana RK. Anatomical likelihood estimation meta-analysis of grey and white matter anomalies in autism spectrum disorders. Neuroimage Clin. 2015;7:525–36.
Carlisi CO, Norman LJ, Lukito SS, Radua J, Mataix-Cols D, Rubia K. Comparative multimodal meta-analysis of structural and functional brain abnormalities in autism spectrum disorder and obsessive-compulsive disorder. Biol Psychiatry. 2017;82(2):83–102.
Liu J, Yao L, Zhang W, Xiao Y, Liu L, Gao X, et al. Gray matter abnormalities in pediatric autism spectrum disorder: a meta-analysis with signed differential mapping. Eur Child Adolesc Psychiatry. 2017;26(8):933–45.
Cai J, Hu X, Guo K, Yang P, Situ M, Huang Y. Increased left inferior temporal gyrus was found in both low function autism and high function autism. Front Psychiatry. 2018;9:542.
Vasa RA, Mostofsky SH, Ewen JB. The disrupted connectivity hypothesis of autism spectrum disorders: time for the next phase in research. Biol Psychiatry Cogn Neurosci Neuroimaging. 2016;1(3):245–52.
Geschwind DH, Levitt P. Autism spectrum disorders: developmental disconnection syndromes. Curr Opin Neurobiol. 2007;17(1):103–11.
Picci G, Gotts SJ, Scherf KS. A theoretical rut: revisiting and critically evaluating the generalized under/over-connectivity hypothesis of autism. Dev Sci. 2016;19(4):524–49.
Beckmann CF, Smith SM. Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging. 2004;23(2):137–52.
Xu L, Groth KM, Pearlson G, Schretlen DJ, Calhoun VD. Source-based morphometry: the use of independent component analysis to identify gray matter differences with application to schizophrenia. Hum Brain Mapp. 2009;30(3):711–24.
Hafkemeijer A, Moller C, Dopper EG, Jiskoot LC, van den Berg-Huysmans AA, van Swieten JC, et al. Differences in structural covariance brain networks between behavioral variant frontotemporal dementia and Alzheimer’s disease. Hum Brain Mapp. 2016;37(3):978–88.
Vanasse TJ, Fox PM, Barron DS, Robertson M, Eickhoff SB, Lancaster JL, et al. BrainMap VBM: An environment for structural meta-analysis. Hum Brain Mapp. 2018;39(8):3308–25.
Ecker C, Suckling J, Deoni SC, Lombardo MV, Bullmore ET, Baron-Cohen S, et al. Brain anatomy and its relationship to behavior in adults with autism spectrum disorder: a multicenter magnetic resonance imaging study. Arch Gen Psychiatry. 2012;69(2):195–209.
Smith SM, Nichols TE. Statistical challenges in “big data” human neuroimaging. Neuron. 2018;97(2):263–8.
Mulders P, Llera A, Tendolkar I, van Eijndhoven P, Beckmann C. Personality profiles are associated with functional brain networks related to cognition and emotion. Sci Rep. 2018;8(1):13874.
Hotelling H. Relations between two sets of variates. Biometrika. 1936;28:312–77.
Loth E, Charman T, Mason L, Tillmann J, Jones EJH, Wooldridge C, et al. The EU-AIMS Longitudinal European Autism Project (LEAP): design and methodologies to identify and validate stratification biomarkers for autism spectrum disorders. Mol Autism. 2017;8:24.
Charman T, Loth E, Tillmann J, Crawley D, Wooldridge C, Goyard D, et al. The EU-AIMS Longitudinal European Autism Project (LEAP): clinical characterisation. Mol Autism. 2017;8:27.
Rutter M, Le Couteur A, Lord C. Autism diagnostic interview-revised. Los Angeles: Western Psychological Services; 2003.
Lord C, Rutter M, DiLavore PC, Risi S, Gotham K, Bishop S. Autism diagnosticobservation schedule, Second Edition (ADOS-2) manual (part I): modules 1–4. Torrance: Western Psychological Services; 2012.
Constantino JN, Gruber CP. Social responsiveness scale. 2nd ed. Los Angeles: Western Psychological Services; 2012.
Bodfish JW, Symons FJ, Parker DE, Lewis MH. Varieties of repetitive behavior in autism: comparisons to mental retardation. J Autism Dev Disord. 2000;30(3):237–43.
Tomchek SD, Dunn W. Sensory processing in children with and without autism: a comparative study using the short sensory profile. Am J Occup Ther. 2007;61(2):190–200.
Gaser C, Dahnke R. https://dbm.neuro.uni-jena.de/cat/. Accessed 2 Sept 2018.
Ashburner J. A fast diffeomorphic image registration algorithm. Neuroimage. 2007;38(1):95–113.
Groves AR, Smith SM, Fjell AM, Tamnes CK, Walhovd KB, Douaud G, et al. Benefits of multi-modal fusion analysis on a large-scale dataset: life-span patterns of inter-subject variability in cortical morphometry and white matter microstructure. Neuroimage. 2012;63(1):365–80.
Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. Fsl Neuroimage. 2012;62(2):782–90.
Eklund A, Nichols TE, Knutsson H. Cluster failure: why fMRI inferences for spatial extent have inflated false-positive rates. Proc Natl Acad Sci USA. 2016;113(28):7900–5.
Winkler AM, Ridgway GR, Douaud G, Nichols TE, Smith SM. Faster permutation inference in brain imaging. Neuroimage. 2016;141:502–16.
Simonoff E, Pickles A, Charman T, Chandler S, Loucas T, Baird G. Psychiatric disorders in children with autism spectrum disorders: prevalence, comorbidity, and associated factors in a population-derived sample. J Am Acad Child Adolesc Psychiatry. 2008;47(8):921–9.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodolo). 1995;57(1):289–300.
Smith SM, Nichols TE, Vidaurre D, Winkler AM, Behrens TE, Glasser MF, et al. A positive-negative mode of population covariation links brain connectivity, demographics and behavior. Nat Neurosci. 2015;18(11):1565–7.
Haufe S, Meinecke F, Gorgen K, Dahne S, Haynes JD, Blankertz B, et al. On the interpretation of weight vectors of linear models in multivariate neuroimaging. Neuroimage. 2014;87:96–110.
Di X, Biswal BB. Similarly expanded bilateral temporal lobe volumes in female and male children with autism spectrum disorder. Biol Psychiatry Cogn Neurosci Neuroimaging. 2016;1(2):178–85.
Beul SF, Grant S, Hilgetag CC. A predictive model of the cat cortical connectome based on cytoarchitecture and distance. Brain Struct Funct. 2015;220(6):3167–84.
Oligschlager S, Xu T, Baczkowski BM, Falkiewicz M, Falchier A, Linn G, et al. Gradients of connectivity distance in the cerebral cortex of the macaque monkey. Brain Struct Funct. 2019;224(2):925–35.
Nickl-Jockschat T, Habel U, Michel TM, Manning J, Laird AR, Fox PT, et al. Brain structure anomalies in autism spectrum disorder–a meta-analysis of VBM studies using anatomic likelihood estimation. Hum Brain Mapp. 2012;33(6):1470–89.
Langen M, Bos D, Noordermeer SD, Nederveen H, van Engeland H, Durston S. Changes in the development of striatum are involved in repetitive behavior in autism. Biol Psychiatry. 2014;76(5):405–11.
Foster NE, Doyle-Thomas KA, Tryfon A, Ouimet T, Anagnostou E, Evans AC, et al. Structural gray matter differences during childhood development in autism spectrum disorder: a multimetric approach. Pediatr Neurol. 2015;53(4):350–9.
Sato W, Kochiyama T, Uono S, Yoshimura S, Kubota Y, Sawada R, et al. Reduced gray matter volume in the social brain network in adults with autism spectrum disorder. Front Hum Neurosci. 2017;11:395.
Kana RK, Keller TA, Minshew NJ, Just MA. Inhibitory control in high-functioning autism: decreased activation and underconnectivity in inhibition networks. Biol Psychiatry. 2007;62(3):198–206.
Wadsworth HM, Maximo JO, Donnelly RJ, Kana RK. Action simulation and mirroring in children with autism spectrum disorders. Behav Brain Res. 2018;341:1–8.
Nomi JS, Molnar-Szakacs I, Uddin LQ. Insular function in autism: Update and future directions in neuroimaging and interventions. Prog Neuropsychopharmacol Biol Psychiatry. 2019;89:412–26.
Patriquin MA, DeRamus T, Libero LE, Laird A, Kana RK. Neuroanatomical and neurofunctional markers of social cognition in autism spectrum disorder. Hum Brain Mapp. 2016;37(11):3957–78.
Langen M, Leemans A, Johnston P, Ecker C, Daly E, Murphy CM, et al. Fronto-striatal circuitry and inhibitory control in autism: findings from diffusion tensor imaging tractography. Cortex. 2012;48(2):183–93.
Di Martino A, Kelly C, Grzadzinski R, Zuo XN, Mennes M, Mairena MA, et al. Aberrant striatal functional connectivity in children with autism. Biol Psychiatry. 2011;69(9):847–56.
Clements CC, Zoltowski AR, Yankowitz LD, Yerys BE, Schultz RT, Herrington JD. Evaluation of the social motivation hypothesis of autism: a systematic review and meta-analysis. JAMA Psychiatry. 2018;75(8):797–808.
Baribeau DA, Dupuis A, Paton TA, Hammill C, Scherer SW, Schachar RJ, et al. Structural neuroimaging correlates of social deficits are similar in autism spectrum disorder and attention-deficit/hyperactivity disorder: analysis from the POND Network. Transl Psychiatry. 2019;9(1):72.
Aoki Y, Cortese S, Tansella M. Neural bases of atypical emotional face processing in autism: A meta-analysis of fMRI studies. World J Biol Psychiatry. 2015;16(5):291–300.
Klapwijk ET, Aghajani M, Colins OF, Marijnissen GM, Popma A, van Lang ND, et al. Different brain responses during empathy in autism spectrum disorders versus conduct disorder and callous-unemotional traits. J Child Psychol Psychiatry. 2016;57(6):737–47.
Boucher J, Mayes A, Bigham S. Memory in autistic spectrum disorder. Psychol Bull. 2012;138(3):458–96.
Barendse EM, Hendriks MP, Jansen JF, Backes WH, Hofman PA, Thoonen G, et al. Working memory deficits in high-functioning adolescents with autism spectrum disorders: neuropsychological and neuroimaging correlates. J Neurodev Disord. 2013;5(1):14.
We thank all participants and their families for participating in this study. We gratefully acknowledge the contributions of all members of the EU-AIMS LEAP group: Jumana Ahmad, Sara Ambrosino, Bonnie Auyeung, Tobias Banaschewski, Simon Baron-Cohen, Sarah Baumeister, Christian F. Beckmann, Sven Bölte, Thomas Bourgeron, Carsten Bours, Michael Brammer, Daniel Brandeis, Claudia Brogna, Yvette de Bruijn, Jan K. Buitelaar, Bhismadev Chakrabarti, Tony Charman, Ineke Cornelissen, Daisy Crawley, Flavio Dell’Acqua, Guillaume Dumas, Sarah Durston, Christine Ecker, Jessica Faulkner, Vincent Frouin, Pilar Garcés, David Goyard, Lindsay Ham, Hannah Hayward, Joerg Hipp, Rosemary Holt, Mark H. Johnson, Emily J.H. Jones, Prantik Kundu, Meng-Chuan Lai, Xavier Liogier D’ardhuy, Michael V. Lombardo, Eva Loth, David J. Lythgoe, René Mandl, Andre Marquand, Luke Mason, Maarten Mennes, Andreas Meyer-Lindenberg, Carolin Moessnang, Nico Mueller, Declan G.M. Murphy, Bethany Oakley, Laurence O’Dwyer, Marianne Oldehinkel, Bob Oranje, Gahan Pandina, Antonio M. Persico, Annika Rausch, Barbara Ruggeri, Amber Ruigrok, Jessica Sabet, Roberto Sacco, Antonia San José Cáceres, Emily Simonoff, Will Spooren, Julian Tillmann, Roberto Toro, Heike Tost, Jack Waldman, Steve C.R. Williams, Caroline Wooldridge, and Marcel P. Zwiers.
This work is primarily supported by the EU-AIMS consortium (European Autism Interventions), which receives support from Innovative Medicines Initiative Joint Undertaking Grant No.115300, the resources of which are composed of financial contributions from the European Union’s Seventh Framework Programme (Grant No. FP7/2007-2013), from the European Federation of Pharmaceutical Industries and Associations companies’ in-kind contributions; and by the AIMS-2-TRIALS consortium (Autism Innovative Medicine Studies-2-Trials), which has received funding from the Innovative Medicines Initiative 2 Joint Undertaking under Grant Agreement No. 777394, and this Joint Undertaking receives support from the European Union's Horizon 2020 research and innovation programme and EFPIA and AUTISM SPEAKS, Autistica, SFARI. TM is supported by a China Scholarship Council Grant (No 201806010408). This work has been further supported by the European Union Seventh Framework Programme Grant Nos. 602805 (AGGRESSOTYPE) (to JKB), 603016 (MATRICS) (to JKB), and 278948 (TACTICS) (to JKB); European Community’s Horizon 2020 Programme (H2020/2014–2020) Grant Nos. 643051 (MiND) (to JKB), 642996 (BRAINVIEW) (to JKB) and 847818 (CANDY) (to JKB and CFB); the Netherlands Organization for Scientific Research VICI Grant No. 2020/TTW/00836465 (to CFB); Wellcome Trust Collaborative Award Grant No. 215573/Z/19/Z (to CFB); the Autism Research Trust (to SBC).
Ethics approval and consent to participate
Ethical approval for this study was obtained through ethics committees at each site.
Consent for publication
Consent for publication was obtained from all participants prior to the study.
JKB has been a consultant to, advisory board member of, and a speaker for Janssen Cilag BV, Eli Lilly, Shire, Lundbeck, Roche, and Servier. He is not an employee of any of these companies, and not a stock shareholder of any of these companies. He has no other financial or material support, including expert testimony, patents or royalties. CFB is director and shareholder in SBGNeuro Ltd. TB served in an advisory or consultancy role for Lundbeck, Medice, Neurim Pharmaceuticals, Oberberg GmbH, Shire, and Infectopharm. He received conference support or speaker’s fee by Lilly, Medice, and Shire. He received royalties from Hogrefe, Kohlhammer, CIP Medien, and Oxford University Press. TC has received consultancy from Roche and received book royalties from Guildford Press and Sage. DGM has been a consultant to, and advisory board member, for Roche and Servier. He is not an employee of any of these companies, and not a stock shareholder of any of these companies. The present work is unrelated to the above grants and relationships. The other authors report no biomedical financial interests or potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Subsection 1. Acquisition parameters. Subsection 2. Demographic information of each schedule. Subsection 3. Sensitivity analyses. Subsection 4. Leave-one-out (LOO) validation of CCA results. Subsection 5. Group differences at voxel-wise gray matter volumes. Subsection 6. Independent components with case-control differences. Subsection 7. Robustness assessment of ICA model orders. Subsection 8. GLM results of the association between brain components and symptom profiles. Subsection 9. Components with highest loadings in CCA. Subsection 10. Uncorrected main CCA mode loadings of each component.
About this article
Cite this article
Mei, T., Llera, A., Floris, D.L. et al. Gray matter covariations and core symptoms of autism: the EU-AIMS Longitudinal European Autism Project. Molecular Autism 11, 86 (2020). https://doi.org/10.1186/s13229-020-00389-4
- Magnetic resonance imaging
- Voxel-based morphometry
- Independent component analysis
- Canonical correlation analysis