Skip to main content

The neuroanatomical substrates of autism and ADHD and their link to putative genomic underpinnings

Abstract

Background

Autism spectrum disorders (ASD) are neurodevelopmental conditions accompanied by differences in brain development. Neuroanatomical differences in autism are variable across individuals and likely underpin distinct clinical phenotypes. To parse heterogeneity, it is essential to establish how the neurobiology of ASD is modulated by differences associated with co-occurring conditions, such as attention-deficit/hyperactivity disorder (ADHD). This study aimed to (1) investigate between-group differences in autistic individuals with and without co-occurring ADHD, and to (2) link these variances to putative genomic underpinnings.

Methods

We examined differences in cortical thickness (CT) and surface area (SA) and their genomic associations in a sample of 533 individuals from the Longitudinal European Autism Project. Using a general linear model including main effects of autism and ADHD, and an ASD-by-ADHD interaction, we examined to which degree ADHD modulates the autism-related neuroanatomy. Further, leveraging the spatial gene expression data of the Allen Human Brain Atlas, we identified genes whose spatial expression patterns resemble our neuroimaging findings.

Results

In addition to significant main effects for ASD and ADHD in fronto-temporal, limbic, and occipital regions, we observed a significant ASD-by-ADHD interaction in the left precentral gyrus and the right frontal gyrus for measures of CT and SA, respectively. Moreover, individuals with ASD + ADHD differed in CT to those without. Both main effects and the interaction were enriched for ASD—but not for ADHD-related genes.

Limitations

Although we employed a multicenter design to overcome single-site recruitment limitations, our sample size of N = 25 individuals in the ADHD only group is relatively small compared to the other subgroups, which limits the generalizability of the results. Also, we assigned subjects into ADHD positive groupings according to the DSM-5 rating scale. While this is sufficient for obtaining a research diagnosis of ADHD, our approach did not take into account for how long the symptoms have been present, which is typically considered when assessing ADHD in the clinical setting.

Conclusion

Thus, our findings suggest that the neuroanatomy of ASD is significantly modulated by ADHD, and that autistic individuals with co-occurring ADHD may have specific neuroanatomical underpinnings potentially mediated by atypical gene expression.

Introduction

Autism spectrum disorders (ASD) are highly heterogeneous neurodevelopmental conditions characterized by differences in social interaction and communication, alongside restricted, repetitive behaviors and interests [1]. These features are associated with a different development of the brain [2] and regional differences in neuroanatomy [3, 4]. Yet, neuroanatomical differences in autism are highly variable across individuals and likely underpin different clinical phenotypes.

Autism is a condition with a high prevalence of co-occurring conditions. More specifically, 72% autistic adolescents have at least one co-occurring psychiatric disorder, and often meet diagnostic criteria for two or more co-occurring conditions [5]. Among these, attention-deficit/hyperactivity disorder (ADHD) is the most prevalent, with estimates ranging between 16 and 80% of autistic individuals also having a clinical diagnosis of ADHD [5,6,7,8]. The large range in prevalence estimates has mainly been attributed to variability in diagnostic methods employed [5, 6, 8]. Given the large degree of co-occurrence between autism and ADHD, it is likely that inter-individual differences in the neuroanatomy of autism are confounded by the severity and nature of ADHD features and so contribute to the highly complex clinical and neurobiological phenotype that is characteristic for autism. Overall, autism is characterized by neuroanatomical differences in several large-scale neural systems that include fronto-temporal and fronto-parietal regions, the amygdala-hippocampal complex, the cerebellum, basal ganglia, and anterior and posterior cingulate regions [2,3,4]. Many of these brain regions overlap with core components of the so-called social brain network, which comprises a set of brain regions that mediate functions related to social cognition and/or emotional processing [9]. Similarly, ADHD has been linked to neuroanatomical differences in frontal regions (e.g., dorsolateral and ventromedial prefrontal cortex), the parietal cortex, cingulate cortices, basal ganglia, the limbic lobe (e.g., amygdala and hippocampus) and the cerebellum [10,11,12]. There is thus a large degree of overlap between the neural networks underpinning autism and ADHD. Even though the neural systems subserving autism and ADHD seem to overlap to a large degree, neuroanatomical differences between both conditions have also been reported. For example, the neuroanatomical differences in the temporal lobe are less commonly reported in ADHD than in autism, while differences in cingulate brain regions seem to be more common in ADHD than autism. Moreover, it has been reported that individuals with ADHD (without co-occurring autism) have significantly reduced gray matter (GM) volume in fronto-temporal areas relative to autistic individuals, as well as significantly increased volume in parietal areas [13]. Another study reported that, compared controls, autistic individuals show increased cortical thickness (CT), while individuals with ADHD show for example decreased CT across the brain [14]. Taken together, these studies suggest that ASD and ADHD may have separable neuroanatomical underpinnings.

So far, only few studies have explored the neuroanatomy of autism with co-occurring ADHD. A study by Mizuno et al. [15] reported that autistic individuals with co-occurring ADHD have significantly decreased volume in the left postcentral gyrus compared to typically developing controls [15]. However, as this study only examined individuals with a diagnosis of both autism and ADHD relative to neurotypical controls, it was not possible to disentangle the relative impact of autism and ADHD on the level of brain anatomy. Another study observed significant differences in gray matter volume and surface area (SA) in the post- and precentral gyrus when comparing individuals with autism to individuals with autism and ADHD [16]. Further, there are reports of significant ASD-by-ADHD interactions in the parietal, temporal and limbic lobe for measures of CT, and in the right cingulate cortex for measures of cortical volume [17]. However, in this study, the interaction terms did not survive multiple comparison correction [17]. Neuroanatomical variability may therefore exist not only between neurodevelopmental conditions, but also between distinct measures of neuroanatomy.

Traditional investigations of brain structure have mainly focused on cortical volume, measured on the whole-brain, regional, and/or vertex-level. However, cortical volume is, by definition, the product of two different cortical features, namely CT and SA. Previous research by our group has shown that measures of brain volume, therefore, do not characterize a specific (i.e., unique) aspect of the neural architecture, but can be explained by separable variations in CT and SA (e.g., [18]). These two morphometric features also have distinct genetic determinants and developmental trajectories [19,20,21,22]. Thus, to characterize neuroanatomical variability associated with ASD and/or ADHD, it is important to investigate CT and SA separately. Using data provided by the EU-AIMS Longitudinal European Autism Project (LEAP; [23]), a large-scale European research collaboration on autism, we therefore examined the neuroanatomical underpinnings of autism with and without co-occurring ADHD symptomatology based on measures of CT and SA. Given the conceptual complexity of both ASD and ADHD [24,25,26], we utilized a categorical rather than a dimensional design that allowed us to establish the extent to which ADHD symptoms modulate the neuroanatomy of autism.

Moreover, leveraging the spatial gene expression data provided by the Allen Human Brain Atlas (AHBA; [27]) we aimed at linking these neuroanatomical differences to putative molecular underpinnings. Overall, neuropsychiatric conditions are highly heritable with an estimated heritability of 83% for autism [28], and 75% for ADHD [29]. However, the genetic architecture of autism is complex, involving hundreds (or more) genetic variants that mediate wider autism traits [30]. Many of these genes map onto biological pathways underpinning neural development [30,31,32]. The genetic architecture of ADHD is equally complex, and implicates mechanisms underpinning early embryonic development and cognitive abilities [33]. Previous studies also suggest that there is a significant genetic overlap between autism and ADHD. Overall, the genetic correlation between both conditions is 0.37 [34], which is higher than the genetic correlation between autism and other psychiatric disorders (e.g., r = 0.22 for autism and schizophrenia (SCZ), r = 0.14 for autism and bipolar disorder; [34]). Genome-wide association studies (GWAS) report several single nucleotide polymorphisms (SNPs) that are significantly associated with autism or ADHD, and a significant overlap between the SNPs associated with autism and those associated with ADHD has been reported [35]. These SNPs map onto genes that are predominantly expressed in the brain, and play a crucial role in neuronal migration, neuronal development, neuromodulation of neurotransmission, and in general brain development [35]. In the present study, we therefore also examined whether brain regions where the neuroanatomy of autism is significantly modulated by ADHD are enriched for genes that have previously been linked to autism and/or ADHD on the genetic and transcriptomic level.

Methods and materials

Participants

This study utilized data provided by the EU-AIMS LEAP project [23], which is a European multicenter study on stratification biomarkers for autism (www.aims-2-trials.eu). A comprehensive description of the sample has been published elsewhere [23]. In brief, the total sample for which usable structural MRI data were available included N = 638 individuals (N = 359 with ASD and N = 279 non-autistic controls). For the present study, a subset of N = 533 individuals aged between 6 and 30 years was selected that included all individuals with available data from the ADHD DSM-5 rating scale [1] (Table 1). Based on autism diagnostic criteria according to DSM-IV [36], DSM-IV-TR [37], DSM-5 [1], or ICD-10 [38] and whether participants met criteria for ADHD according to the ADHD DSM-5 rating scale [1], we subdivided the autistic individuals into two subgroups: ASD only (N = 170; 112 males, 58 females) and ASD with co-occurring ADHD [further referred to as ASD + ADHD group (N = 142; 107 males, 35 females)]. From the neurodiverse control group without an ASD diagnosis, we collected those individuals who met criteria for ADHD as a third subgroup [further referred to as ADHD only group (N = 25; 14 males, 11 females)] and those who did not meet ADHD criteria as our fourth subgroup [further referred to as typically developing (TD) controls (N = 196; 124 males, 72 females)]. The ADHD DSM-5 rating scale was based on either parent- or self-report scores depending on participants’ age. For this study we used the categorical output of the ADHD rating scale, which measures the presence of symptoms, evaluated on a 0–3 scale (0 = not at all to 3 = very often). To be evaluated with a clinical concern, a participant must at least score six responses with a 2 or 3 (“Often”, “Very Often”) (for more details see [23]). Therefore, the rating scale output is not a clinical diagnosis per se but the symptom count is a proxy for one. Notably, subsets were not matched for age, IQ, and sex to maximize the available sample size. We therefore controlled for these measures in all subsequent analyses (see below). Given the wide range in full scale IQ (FSIQ) across individuals, we also repeated our analyses in a smaller subset of individuals, which excluded participants with a mild intellectual disability (ID, i.e., FSIQ < 70). The results of these analyses are presented in the Additional file 1: Effects of Intellectual Disability. For further details on demographics and exclusion criteria see Additional file 1: Sample Description. An independent ethics committee at each site approved the study, and written informed consent was obtained for all participants.

Table 1 Characteristics of participants with ASD, ADHD, ASD + ADHD and control subjects

MRI data acquisition

All participants underwent MR imaging in 3-T scanners, at six different sites (University of Cambridge and King’s College London, U.K.; Central Institute for Mental Health, Mannheim, Germany; Radboud University Medical Centre and University Medical Centre Utrecht, the Netherlands; Rome University, Italy). High-resolution structural T1-weighted volumetric images were acquired with full head coverage, at 1.2 mm thickness with 1.2 × 1.2-mm in-plane resolution (see Additional file 1: Table S1 for details).

Cortical surface reconstructions using FreeSurfer

Usable structural MRI data were initially available for 708 individuals in the LEAP sample. FreeSurfer v6.0.0 software (http://surfer.nmr.mgh.harvard.edu/) was used to derive models of the cortical surface for each T1-weighted image. These well-validated and fully automated procedures have been previously described elsewhere [39,40,41,42,43]. Each reconstructed surface underwent strict quality assessments (see Additional file 1: MRI Data Quality Assessments), resulting in a final sample of N = 638 individuals [44]. In brief, three independent raters judged the quality of each scan with three possible decisions: include, exclude, edit [44]. In our study, we examined measures of CT and SA. Measures of CT were computed as the closest distance from the outer (i.e., pial) to the inner (i.e., white matter) boundary at each vertex on the tessellated surface [43]. Measures of surface area were quantified as the area of the cortex at a given point on the cortical surface (i.e., the sum of faces in the polygon mesh representation of the cortex at a particular vertex) as previously described by Winkler et al. [45]. For each participant we also computed mean CT across the entire brain (CT0), as well as total SA (SAtotal). To improve detection of population changes, each parameter was smoothed using a 15-mm surface based smoothing kernel.

Surface-based statistical analyses of cortical thickness and surface area

Statistical analyses were performed using the SurfStat toolbox (https://www.math.mcgill.ca/keith/surfstat/) for MATLAB version R2021a (https://www.mathworks.com), and R for Statistical Computing (www.r-project.org). Vertex-wise differences in neuroanatomy (Y) were quantified using a general linear model (GLM) with (1) ASD, ADHD, sex, and site (see Additional file 1: Site effects for further details, also using ComBat batch correction) as fixed-effect factors, (2) an ASD-by-ADHD interaction term, and (3) linear and quadratic age, FSIQ, and a total brain measure (i.e., CT0 or SAtotal, respectively) as continuous covariates:

$$\begin{aligned} Y_{i} = & \beta_{0} + \beta_{1} {\text{ASD}} + \beta_{2} {\text{ADHD}} + \beta_{3} {\text{ASD}}*{\text{ADHD}} + \beta_{4} {\text{Sex}} + \beta_{5} {\text{Age}} + \beta_{6} {\text{Age}}^{2} \\ & + \beta_{7} {\text{FSIQ}} + \beta_{8} {\text{Site}} + \beta_{9} {\text{total}}\;{\text{Brain}} + \varepsilon_{i} \\ \end{aligned}$$

where εi is the residual error at vertex i. Continuous covariates were mean centered across groups to improve interpretability of the coefficients. Corrections for multiple comparisons across the whole brain were performed using random field theory (RFT) based cluster analysis with a cluster-forming and cluster-based threshold of p < 0.05 (two-tailed). Effect sizes associated with each model term were assessed using Cohen’s f, where values of 0.25, 0.5, and 0.75 indicate small, medium, and large effects, respectively.

To specifically identify neuroanatomical differences between ASD individuals with and without a co-occurring diagnosis of ADHD, we fitted an additional GLM on a subset of the total sample that only included participants with ASD and examined the main effect of ADHD (N = 312), i.e.,

$$Y_{i} = \beta_{0} + \beta_{1} {\text{Group}} + \beta_{2} {\text{Sex}} + \beta_{3} {\text{Age}} + \beta_{4} {\text{Age}}^{2} + \beta_{5} {\text{FSIQ}} + \beta_{6} {\text{Site}} + \beta_{7} {\text{total}}\;{\text{Brain}} + \varepsilon_{i }$$

Gene expression decoding analysis

To link our neuroanatomical findings to putative genomic underpinnings, we leveraged the spatial gene expression data from the Allen Human Brain Atlas (AHBA; [27]) to identify a list of genes with a spatial pattern of expression that resembles the neuroanatomical patterns highlighted by our statistical neuroimaging analyses. To this aim, we initially uploaded the statistical t-maps associated with the main effect of ASD, the main effect of ADHD, as well as the ASD-by-ADHD interaction term for CT and SA (Fig. 4A, C) to the Neurovault server (https://neurovault.org). Next, using python code embedded within Neurovault and Neurosynth (https://neurosynth.org), we performed a gene expression decoding analysis that statistically assesses the spatial correlation between our statistical maps and the pattern of gene expression for each of a total of 20,787 protein coding genes [46]. To do so, the six AHBA donor brains are initially co-registered with the MNI atlas (also used by FreeSurfer) using nonlinear registration (transcriptomic alignment). At each sampling site (i.e., probe), a spherical region-of-interest (ROI) is drawn (default radius r = 4 mm), and the statistical test parameter in each FreeSurfer overlay was averaged within each ROI. This resulted in a spatial vector of values for each donor, which was subsequently correlated with the normalized gene expression data. Here, the analysis constructs a linear model for each donor brain, where the slopes encode the spatial correlation between each gene’s expression pattern to the statistical neuroimaging map (random effects model). In line with the input maps, these analyses were restricted to cortical tissue. The slopes are then subjected to a one-sample t test to identify genes whose spatial expression patterns are consistently highly similar to the imaging maps (i.e., across donor brains). The derived list of genes was thresholded at p < 0.01. We chose this ‘liberal’ threshold as this analysis did not constitute a hypothesis test per se, but rather a selection step aimed at yielding an initial list of genes for the subsequent enrichment analyses. Given that both sides of our imaging contrasts were of equal relevance, we considered both positive and negative t-statistic values.

In addition to the Neurosynth decoding approach, we performed gene expression decoding using a General Least Squares (GLS) approach that also accounted for spatial autocorrelations in embedded transcriptomic maps, to assess the robustness of our findings. The GLS approach is described in detail in Additional file 1: General Least Square (GLS)-decoding.

Gene enrichment analyses

Next, we performed several gene enrichment analyses to establish the biological relevance and functional role of decoded genes. All enrichment testing was performed using the GeneOverlap package in R (https://doi.org/10.18129/B9.bioc.GeneOverlap). Specifically, we tested the decoded gene lists for an enrichment with different gene sets known to be associated with ASD at the genetic and transcriptomic level. At the genetic level, this included the 102 rare and de novo protein truncating variants identified in the largest exome sequencing study of autism worldwide [47]. We also included an ASD-related gene list compiled by SFARI (categories S, 1, 2, and 3 downloaded in November 2020 from https://gene.sfari.org/). At the transcriptomic level, we included a list of differentially expressed genes (DEGs) (upregulated/downregulated) in post-mortem cortex tissue in ASD [48], and genes that are differentially expressed in specific cell types in ASD [49]. Moreover, we included genes from differentially expressed co-expression modules in ASD that map onto specific biological processes [50]. ADHD-related genes were derived based on a GWAS study published by Demontis and colleagues [33]. Here, we used the MAGMA plug-in on the FUMA GWAS annotation platform (https://fuma.ctglab.nl) to perform a GWAS analysis [51]. Here, variants were mapped onto genes based on their exact position, and aggregated association p values were calculated for each gene. Taking into account the sample composition, a European ancestry reference from 1000 Genomes phase 3 was used as reference panel. Bonferroni correction was used to set the significance threshold (correcting for all N = 10,894 gene sets tested from MsigdB v5.2; [52]) The resulting gene list consisted of N = 22 ADHD-related genes, which was used for further gene enrichment analyses. Our enrichment tests generated enrichment odds ratios, hypergeometric p values, and FDR-corrected p values using a background total of the 20,787 Neurosynth genes. Last, we examined the percentage overlap between genes associated with both main effects and the interaction term relative to the total number of genes significantly associated with the imaging phenotypes (see Additional file 1: Figure S7 for details).

Results

Subject demographics

There was no significant difference in the ratio of males to females between any of the four subgroups (p = 0.06642) (Table 1). There was, however, a significant difference in age between the ASD group and the co-occurring group (ASD only: 18.45 ± 5.6; ASD + ADHD: 16.38 ± 5.30; p = 0.0061) (Table 1). Subgroups also differed significantly in FSIQ (p < 2e−16) [with TD having the highest IQ (108.44 ± 14.1) and the ADHD only group the lowest (79.8 ± 21.1)] (Table 1). We therefore covaried for these potential confounds in all subsequent analyses.

Main effect of ASD and ADHD on CT and SA

Following RFT-based cluster-correction (p < 0.05, two-tailed), the main effect of autism was associated with increased CT in the anterior-cingulate cortex (ACC; approximate Brodmann areas (BA) 24/33), in the left superior and middle temporal gyrus (BA 21/22), and the right precuneus cortex (BA 31). By contrast, autism was associated with decreased CT in the parietal cortex (BA 7) and the middle frontal gyrus (BA 6/8/9) (Fig. 1A). We also found that the main effect of autism was associated with decreased SA in the right anterior-cingulate cortex (BA 24/33) relative to non-autistic individuals (Fig. 2A; Additional file 1: Table S2). For the main effect of ADHD, we observed decreases in CT in right superior frontal gyrus (BA 6/8/9), right cingulate cortex (BA 23/24/33), and the right precuneus cortex (BA 31) (Fig. 1B). ADHD was also associated with increased SA in left parahippocampal gyrus (BA 27/28) only (Fig. 2B; Additional file 1: Table S3; see Additional file 1: Figure S3 for unthresholded t-maps). Vertex-level effect sizes (Cohen’s f) for the main effects of autism and ADHD, and all other model terms are displayed in the Additional file 1: Figure S4.

Fig. 1
figure 1

Differences of cortical thickness for main effects of ASD, ADHD, and ASD-by-ADHD interaction. AC Random field theory (RFT)-based cluster corrected t-maps (p < 0.05, 2-tailed) for CT. A Significant decreases in CT in ASD compared to non-ASD are displayed in blue and significant increases are displayed in orange. B Significant decrease in CT in ADHD compared to non-ADHD is displayed in blue and significant increase is displayed in orange, C the ASD × ADHD interaction effect for CT, D mean CT at right cingulate cortex for the interaction effect, E mean CT at right temporal lobe for the interaction effect and F mean CT at left parietal cortex for the interaction effect. ASD autism spectrum disorder, ADHD attention-deficit/hyperactivity disorder, TD typically developing, t t-statistic, CT cortical thickness, L left, R right

Fig. 2
figure 2

Differences of surface area for main effects of ASD, ADHD, and ASD-by-ADHD interaction. AC Random field theory (RFT)-based cluster corrected t-maps (p < 0.05, 2-tailed) for SA. A Significant decrease in SA in ASD compared to non-ASD is displayed in blue and significant increase is displayed in orange. B Significant decrease in SA in ADHD compared to non-ADHD is displayed in blue and significant increase is displayed in orange, C the ASD × ADHD interaction effect for SA. D Mean SA at right frontal gyrus for the interaction effect. ASD autism spectrum disorder, ADHD attention-deficit/hyperactivity disorder, TD typically developing, t t-statistic, SA surface area, L left, R right

Significant interactions between ASD and ADHD

In addition to the main effects of group, we examined an ASD-by-ADHD interaction term. We identified several cortical regions, where the neuroanatomy of autism (or ADHD) was significantly modulated by co-occurring ADHD (or ASD) including the left parietal cortex (BA 7), precentral and superior frontal gyrus (BA 6/8), the right temporal gyrus (BA 20/21/22), cingulate cortex (BA 23/24/33), and precuneus cortex (BA 31) for measures of CT (Fig. 1C). In these regions, we observed different types of interactions. For example, in the cingulate cortex CT in the ASD + ADHD group was significantly increased compared to the other groups (Fig. 1D), whereas in the precentral gyrus CT in the ASD + ADHD group was significantly decreased compared to all other groups (Fig. 1F). We also found a significant cluster in the right temporal gyrus with no significant difference between both ASD groups. Here, the interaction was mainly driven by the difference between individuals with ADHD and TD controls, which was significantly increased for measures of CT (Fig. 1E; see Additional file 1: Figure S5A for interaction plots of all significant clusters). For measures of SA, significant differences in the right dorsolateral prefrontal gyrus (BA 46) (Fig. 2C) were found. Here, the interaction was mainly driven by the difference between individuals with ADHD and TD controls, which was significantly decreased, while small or no differences were observed between controls and both autism subgroups (Fig. 2D; see Additional file 1: Figure S5B for interaction plots of all significant clusters; Additional file 1: Table S4).

Significant differences between ASD + ADHD and ASD only

To specifically identify differences between the ASD + ADHD and the ASD only group, we examined a subsample of individuals with ASD. We observed that individuals with ASD + ADHD had significantly decreased CT relative to individuals with ASD only in the left precentral, caudal middle frontal, and postcentral gyrus (BA 6/8/9) (Fig. 3; Additional file 1: Table S4). There were no significant differences for measures of SA.

Fig. 3
figure 3

Differences of cortical thickness for ASD + ADHD versus ASD only. Random-field-theory-based cluster-corrected t-maps (p < 0.05, 2-tailed). A Significant decrease in CT in ASD + ADHD compared to ASD only is displayed in blue and significant increase is displayed in orange, B mean CT at left parietal cortex for the ASD + ADHD subgroup (blue) and the ASD only subgroup (yellow). ASD autism spectrum disorder, ADHD attention-deficit/hyperactivity disorder, t t-statistic, CT cortical thickness, L left, R right. **p < 0.01

Gene set enrichment analyses

Here, we decoded the t-maps for the main effects of autism and ADHD, as well as the CT and SA maps for the ASD-by-ADHD interaction terms (i.e., six maps in total). This resulted in sets of (1) NCT = 598 and NSA = 259 significant genes associated with the main effect of autism for CT and SA, respectively, (2) NCT = 272 and NSA = 1005 significant genes associated with the main effect of ADHD, and (3) NCT = 1258 and NSA = 233 significant genes associated with the ASD-by-ADHD interaction term. Within these gene sets, we observed a significant enrichment of genes known to be associated with autism, and particularly for genes that are known to be differentially expressed in autism during childhood and adolescence (see Fig. 4 for details). More specifically, in all decoded CT maps (i.e., for main effect of autism, ADHD, and the interaction term), we observed a significant enrichment of gene co-expression module CTX.M20, which is known to be upregulated in the autism cortex and has been linked to Gene Ontology terms representing developmental processes and the regulation of cell differentiation [50]. Moreover, we observed a significant enrichment of autism genes listed in the SFARI database. Notably, there was no significant enrichment of autism- or ADHD-related risk genes resulting from GWAS studies. However, these were also among the smallest gene sets tested with only N = 22 genes for ADHD (Fig. 4A, B).

Fig. 4
figure 4

Genomic underpinnings of neurodevelopmental deviations in cortical thickness and surface area in ASD, ADHD, and ASD + ADHD. Panel A and C show the t-maps of the main effect of ASD, the main effect of ADHD and the ASD-by-ADHD interaction term for measures of CT (A) and SA (C). Panel B (CT) and D (SA) show significant odds ratios at a false discovery rate (FDR) corrected p threshold of 0.05 resulting from the gene set enrichment analyses for genes expressed in the different output maps. Gene sets were subdivided into sets with differential gene expression in ASD, sets representing ASD risk genes, and a set representing ADHD risk genes. Gene sets are annotated and labeled based on their original publication. ASD autism spectrum disorder, ADHD attention-deficit/hyperactivity-disorder, CTX cortex, DEG differentially expressed gene, down downregulated expression in ASD, up upregulated expression in ASD, NGenes number of genes in each gene set, CT cortical thickness, SA surface area. *p < 0.05 (FDR-corrected), **p < 0.01 (FDR-corrected); red squares indicate consistently significant enrichments across both approaches (see Additional file11: Figure S6)

Different gene sets were found to be enriched in the SA maps (Fig. 4C, D). Here, in all decoded maps, we observed an enrichment of gene co-expression module CTX.M16, which is associated with neuronal markers and synaptic genes [53], and is known to be downregulated in the autism cortex. In contrast to measures of CT, however, where most gene sets were enriched in both main effects for autism and ADHD, gene set enrichments were more prominent in the main effect of ADHD for measures of SA. Autism-related genes hence seem to have a pattern of expression that more closely resembles the differences in CT in autism rather than differences in SA. Moreover, autism-related genes are equally important to the neuroanatomy underpinning ADHD symptomatology, particularly when examining measures of SA. Gene sets resulting from the GLS-decoding approach are shown in the Additional file 1: Figure S6.

Overall (i.e., across gene sets), there was little overlap between genes significantly associated with the main effect of ASD, ADHD, and/or the ASD-by-ADHD interaction terms on the transcriptomic level (see Additional file 1: Figure S7). More specifically, a total of 1,3% of all genes identified as being significantly enriched in any of the imaging phenotypes was shared between the main effect of ASD and ADHD for measures of CT, and a total of 8.1% of genes associated with these main effects for measures of SA.

Discussion

This study examined neuroanatomical differences between autistic individuals with and without co-occurring ADHD relative to individuals with ADHD only and non-autistic controls. Moreover, to bridge the gap between macroscopic and microscopic differences, we examined whether the spatial patterns of neuroanatomical differences in CT and SA are enriched for genes implemented in the etiology of autism and ADHD, and genes known to be differentially expressed in autism. We established that it is possible to separate the effect of autism from the effect of ADHD on the level of neuroanatomy based on measures of CT and SA. Notably, we also observed significant ASD-by-ADHD interactions, and significant differences in CT between ASD individuals with ADHD and those without. This suggests that the neuroanatomy of ASD is significantly modulated by co-occurring ADHD. Further we showed a significant association of autism and ADHD-related patterns of neuroanatomical variability with autism—but not ADHD-related genes. Our study thus provides important novel insights into the neurobiological and putative genomics underpinning the complex clinical phenotype of autism.

Main effect of ASD and ADHD on CT and SA

Neuroanatomical differences in autism and ADHD are well documented and primarily affect fronto-temporal and fronto-parietal regions for autism [2], as well as cingulate cortices for ADHD [11]. In the present study, we examined the neuroanatomical underpinnings of autism and ADHD within a 2 × 2 factorial design that included (1) ASD individuals with co-occurring ADHD, (2) ASD only individuals, (3) ADHD only individuals, and (4) non-autistic controls. While it remains a topic of debate whether autism should be encoded as a categorical variable [54], this design allowed us to identify a set of brain regions where neuroanatomical variability in CT and SA are uniquely attributable to either autism or ADHD. By examining both main effects, we were able to separate autism from ADHD based on patterns of neuroanatomical variability in fronto-temporal regions (e.g., left superior and middle temporal gyrus), and in limbic and prefrontal regions (i.e., cingulate cortex, parahippocampal gyrus, and the superior frontal gyrus). Many of the brain regions associated with the main effect of ASD have previously been linked to the wider neural systems that mediate functions related to social cognition and/or emotional processing, i.e., core autism traits [9]. Prefrontal areas such as the superior frontal gyrus, which showed significant differences in CT for the main effect of ADHD, are important for executive functioning, attention and motor planning, which are all functions known to be impaired in ADHD [55, 56]. Our findings of significant neuroanatomical differences associated with ASD and ADHD are thus in line with previous reports using a categorical approach; however, it does not rule out the necessity of future studies using dimensional approaches. While a categorical approach is particularly well suited to examine neuroanatomical differences associated with both ASD and ADHD (as well as their interaction), future dimensional research would complement our study by accounting for the inter-individual differences in the broader phenotypes of ASD and ADHD. To date, autism diagnostic is still based on a categorical cutoff that remains largely unchanged across the human lifespan.

In addition, as studies examining the neuroanatomical underpinnings of autism and co-morbid ADHD are still rare, and there is currently little knowledge on basic case–control differences, elicited in adequately powered samples. In terms of ASD, several lines of ASD research also converge in suggesting that the condition might be most effectively understood as a categorical fixed effect rather than a dimensional construct (for review see [57]). For example, Frazier et al. [58] examined several indicators of ASD (e.g., eye gaze metrics from social stimulus paradigms), establishing a categorical data structure that closely corresponded to a diagnosis of ASD across measures [58]. ADHD, on the other hand, is widely considered a dimensional construct with fluctuations in symptom severity and profiles across the human lifespan [11, 59, 60]. For example, only 15% of children with ADHD still show criteria for a diagnosis during adulthood [60], and symptoms often also change from a more externalizing nature during childhood to a more internalizing nature during adulthood [11, 59]. Future studies might therefore benefit from employing both a categorical and a dimensional approach to examine the neuroanatomy associated with ASD and ADHD.

The effect sizes associated with the main effect of ASD and ADHD were relatively small, overall, and did not exceed a value of 0.2 on the vertex level. Small-to-medium effects associated with between-group differences in neuroanatomy have previously been observed even in large-scale neuroimaging studies with N > 500 that compare ASD or ADHD individuals with neurotypical controls [12, 14]. The low effects reported by neuroimaging studies seem to be driven by small differences in mean, and larger phenotypic variability in neurodiverse study population [44], thus highlighting the significant inter-individual heterogeneity that is characteristic for neurodevelopmental conditions on the level of etiology, neurobiology, and symptomatology. Taken together, our findings indicate that autism and ADHD may have separable neuroanatomical underpinnings that may mediate differences in clinical phenotypes. The main effects of autism and ADHD are, however, only interpretable in the absence of a significant ASD-by-ADHD interaction.

Significant interactions between ASD and ADHD

While previous studies have focused on the effects of ASD and ADHD in separate case–control designs [14, 61], we examined to what extent, and where in the brain, the neuroanatomy of ASD is modulated by ADHD symptomatology (i.e., the ASD-by-ADHD interaction). We observed a significant ASD-by-ADHD interaction term in the right temporal gyrus, the right cingulate gyrus and the left precentral gyrus for measures of CT, and in the right dorsolateral prefrontal gyrus for measures of SA. Some of these clusters (right temporal gyrus and right dorsolateral prefrontal gyrus) were mainly driven by neuroanatomical differences between individuals with ADHD only relative to controls, while small or no differences were observed between non-autistic controls and autistic individuals (including those with co-occurring ADHD). ADHD-related neuroanatomical variability in these brain regions thus seems to be modulated, or even masked, by also having a diagnosis of autism. However, other clusters associated with the ASD-by-ADHD interaction term, e.g., in the right cingulate gyrus and left precentral gyrus, also showed significant differences between the ASD participants with co-occurring ADHD and those with ASD only. A significant interaction between autism and ADHD indicates that the neuroanatomical underpinnings of ADHD in autistic individuals cannot be explained by either a diagnosis of autism or a diagnosis of ADHD alone.

Significant differences between ASD + ADHD and ASD only

Our additional analysis within the ASD ± ADHD subsample showed significant differences between both groups in the left precentral gyrus for measures of CT. Here, individuals with co-occurring autism and ADHD had a significantly thinner cortex compared to those individuals with autism alone. For measures of SA, no significant clusters for this contrast were observed. Taken together, our findings suggest that (1) the neuroanatomy of autism is significantly modulated by co-occurring ADHD symptomatology, and (2) although ADHD is highly co-occurring in autism, both conditions seem to have separable neuroanatomical underpinnings.

Gene set enrichment analyses

In a next analysis step, we examined whether the patterns of CT and SA attributed to the main effects and the interaction term are linked to the genetic etiology of autism or ADHD. Similar to Romero-Garcia et al. [62], we found that all decoded CT maps were enriched for genes and co-expression modules that have previously been implicated in the etiology of autism, and particularly for genes associated with brain development and the regulation of cell differentiation [50]. For measures of SA, an enrichment of genes and co-expression modules that are associated with neuronal markers and synaptic genes was found [53]. Notably, the pattern of neuroanatomical differences in SA associated with the main effect of ADHD also showed a robust enrichment of autism-related genes, implying that these genes may not be specific to autism, but may also affect the neuroanatomy of other neurodevelopmental conditions such as ADHD.

This finding agrees with previous genetic studies demonstrating that ASD and ADHD are genetically correlated, with a SNP-based genetic correlation of 0.37 [34]. Many of the genetic loci implicated in ASD and ADHD are therefore ‘pleiotropic’, i.e., can influence two or more clinical phenotypes. Moreover, the effects of ASD susceptibility genes on the brain are known to be pleiotropic and may exert their influence via gene regulatory mechanisms during childhood and adolescence [30, 63, 64]. This could also account for highly individualized patterns of neuroanatomical variations or ‘fingerprints’ that are commonly observed in ASD [44]. It is therefore possible that similar genotypes underpin distinct phenotypes, which could explain why we observed an enrichment of ASD-related genes both in the ASD and ADHD imaging phenotype, even though the neuroanatomical patterns characteristic for each phenotype were different.

Moreover, while ASD and ADHD are known to be genetically related, little is currently known about the functional involvement of genes that are specific to ASD and/or ADHD. Moreover, distinct imaging phenotypes could be caused by more complex genetic interactions between ASD- and ADHD-related genes, which cannot be identified based on the gene set enrichment analyses performed in the present study. This highlights the need for conducting future large-scale genome wide association studies (GWAS) to link genetic variation associated with neurodevelopmental conditions to differences in neuroanatomical phenotypes.

Also, while we observed a significant enrichment of autism-related genes, there was no significant enrichment of ADHD risk genes in the main effect of ASD or ADHD. However, the number of recognized ADHD-susceptibility genes remains as yet small, which limits the statistical power of enrichment tests. Moreover, little is currently known about the functional role of ADHD genes and their impact on cortical gene expression. Future genetic and/or transcriptomic studies are therefore needed to provide further insights into the genomic underpinnings of ADHD, and to expand the set of genes that might be utilized to link ADHD-related differences in brain anatomy to putative underlying mechanisms.

Limitations

The present study needs to be interpreted in light of several limitations. Although we employed a multicenter design to overcome single-site recruitment limitations, our sample size of N = 25 individuals in the ADHD only group is relatively small compared to the other subgroups, which limits the generalizability of the results. Also, we assigned subjects into ADHD positive groupings according to the DSM-5 rating scale. While this is sufficient for obtaining a research diagnosis of ADHD, our approach did not take into account for how long the symptoms have been present, which is typically considered when assessing ADHD in the clinical setting. It will therefore be important to repeat our analyses in a bigger sample of pure ADHD cohorts that were diagnosed according to clinical gold standards to make our study comparable to previous findings in these cohorts. Furthermore, future studies employing a dimensional approach are also required to provide a more nuanced picture of the heterogeneity associated with autism. Last, our gene expression decoding analysis was based on the Allen Human Brain Atlas [27], which is the most comprehensive gene expression atlas to date. However, the Allen atlas is based on adult donors exclusively, and provides a coverage that is significantly lower than the spatial resolution of our neuroimaging data. We therefore acknowledge the importance of repeating the analyses in high-resolution age-specific gene expression data sets once these become available, to corroborate the important link between molecular and macroscopic pathology in autism.

Conclusion

Our findings indicate that the neuroanatomy of ASD is significantly modulated by ADHD, and that autistic individuals with co-occurring ADHD may have specific neuroanatomical underpinnings potentially mediated by atypical gene expression.

Availability of data and materials

Due to ethical restrictions, we are currently not able to make the data publicly available. However, it is planned to release the data for perusal of the wider scientific community within the next couple of years once legal requirements are in place.

Abbreviations

ASD:

Autism spectrum disorders

ADHD:

Attention-deficit/hyperactivity disorder

GM:

Gray matter

CT:

Cortical thickness

LEAP:

Longitudinal European Project

SA:

Surface area

AHBA:

Allen Human Brain Atlas

SCZ:

Schizophrenia

GWAS:

Genome-wide association studies

SNP:

Single nucleotide polymorphisms

TD:

Typically developing

GLM:

General linear model

FSIQ:

Full-scale IQ

ID:

Intellectual disability

RFT:

Random field theory

BA:

Brodmann areas

References

  1. American Psychiatric Association (APA). Diagnostic and statistical manual of mental disorders (DSM-5®). 2013.

  2. Ecker C, Bookheimer SY, Murphy DGM. Neuroimaging in autism spectrum disorder: brain structure and function across the lifespan. Lancet Neurol. 2015;14(11):1121–34.

    PubMed  Google Scholar 

  3. 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.

    PubMed  Google Scholar 

  4. Ecker C. The neuroanatomy of autism spectrum disorder: an overview of structural neuroimaging findings and their translatability to the clinical setting. Autism. 2017;21:18–28.

    PubMed  Google Scholar 

  5. Leyfer OT, Folstein SE, Bacalman S, Davis NO, Dinh E, Morgan J, et al. Comorbid psychiatric disorders in children with autism: interview development and rates of disorders. J Autism Dev Disord. 2006;36(7):849–61.

    PubMed  Google Scholar 

  6. Hanson E, Cerban BM, Slater CM, Caccamo LM, Bacic J, Chan E. Brief report: prevalence of attention deficit/hyperactivity disorder among individuals with an autism spectrum disorder. J Autism Dev Disord. 2013;43(6):1459–64.

    PubMed  Google Scholar 

  7. Baker BL, Blacher J. Disruptive behavior disorders in adolescents with ASD: comparisons to youth with intellectual disability or typical cognitive development. J Ment Health Res Intellect Disabil. 2015;8(2):98–116.

    Google Scholar 

  8. de Bruin EI, Ferdinand RF, Meester S, de Nijs PFA, Verheij F. High rates of psychiatric co-morbidity in PDD-NOS. J Autism Dev Disord. 2007;37(5):877–86.

    PubMed  Google Scholar 

  9. Blakemore SJ. The social brain in adolescence. Nat Rev Neurosci. 2008;9(4):267–77.

    CAS  PubMed  Google Scholar 

  10. Hoogman M, Bralten J, Hibar DP, Mennes M, Zwiers MP, Schweren LSJ, et al. Subcortical brain volume differences in participants with attention deficit hyperactivity disorder in children and adults: a cross-sectional mega-analysis. Lancet Psychiatry. 2017;4(4):310–9.

    PubMed  PubMed Central  Google Scholar 

  11. Faraone SV, Asherson P, Banaschewski T, Biederman J, Buitelaar JK, Ramos-Quiroga JA, et al. Attention-deficit/hyperactivity disorder. Nat Rev Dis Primers. 2015;1.

  12. Hoogman M, Muetzel R, Guimaraes JP, Shumskaya E, Mennes M, Zwiers MP, et al. Brain imaging of the cortex in ADHD: a coordinated analysis of large-scale clinical and population-based samples. Am J Psychiatry. 2019;176(7):531–42.

    PubMed  PubMed Central  Google Scholar 

  13. Brieber S, Neufang S, Bruning N, Kamp-Becker I, Remschmidt H, Herpertz-Dahlmann B, et al. Structural brain abnormalities in adolescents with autism spectrum disorder and patients with attention deficit/hyperactivity disorder. J Child Psychol Psychiatry. 2007;48(12):1251–8.

    PubMed  Google Scholar 

  14. Park MTM, Raznahan A, Shaw P, Gogtay N, Lerch JP, Chakravarty MM. Neuroanatomical phenotypes in mental illness: identifying convergent and divergent cortical phenotypes across autism, ADHD and schizophrenia. J Psychiatry Neurosci. 2018;43(3):201–12.

    PubMed  PubMed Central  Google Scholar 

  15. Mizuno Y, Kagitani-Shimono K, Jung M, Makita K, Takiguchi S, Fujisawa TX, et al. Structural brain abnormalities in children and adolescents with comorbid autism spectrum disorder and attention-deficit/hyperactivity disorder. Transl Psychiatry. 2019;9(1):332.

    PubMed  PubMed Central  Google Scholar 

  16. Mahajan R, Dirlikov B, Crocetti D, Mostofsky SH. Motor circuit anatomy in children with autism spectrum disorder with or without attention deficit hyperactivity disorder. Autism Res. 2016;9(1):67–81.

    PubMed  Google Scholar 

  17. Nickel K, Tebartz van Elst L, Manko J, Unterrainer J, Rauh R, Klein C, et al. Inferior frontal gyrus volume loss distinguishes between autism and (comorbid) attention-deficit/hyperactivity disorder—a FreeSurfer analysis in children. Front Psychiatry. 2018;9(October):1–11.

    Google Scholar 

  18. Ecker C, Murphy DGM. Neuroimaging in autism—from basic science to translational research. Nat Rev Neurol. 2014;10(2):82–91.

    PubMed  Google Scholar 

  19. Huttenlocher PR. Synaptic density in human frontal cortex—developmental changes and effects of aging. Brain Res. 1979;163(2):195–205.

    CAS  PubMed  Google Scholar 

  20. Nemani JK, Kostovi I, Nemani O. Prenatal and perinatal development of radial cell columns in the human auditory cortex. Acta Otolaryngol. 1984;97(5–6):489–95. https://doi.org/10.3109/00016488409132926.

    Article  Google Scholar 

  21. Rakic P. Neuronal migration and contact guidance in the primate telencephalon. Postgrad Med J. 1978;54(Suppl 1):25–40.

    PubMed  Google Scholar 

  22. Pretzsch CM, Floris DL, Schäfer T, Bletsch A, Gurr C, Lombardo MV, et al. Cross-sectional and longitudinal neuroanatomical profiles of distinct clinical (adaptive) outcomes in autism. Mol Psychiatry. 2023. https://www.nature.com/articles/s41380-023-02016-z.

  23. 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(1):1–21.

    Google Scholar 

  24. Wozniak RH, Leezenbaum NB, Northrup JB, West KL, Iverson JM. The development of autism spectrum disorders: variability and causal complexity. Wiley Interdiscip Rev Cogn Sci. 2017;8(1–2): e1426.

    Google Scholar 

  25. Joon P, Kumar A, Parle M. What is autism? Pharmacol Rep. 2021;73(5):1255–64.

    PubMed  Google Scholar 

  26. Luo Y, Weibman D, Halperin JM, Li X. A review of heterogeneity in attention deficit/hyperactivity disorder (ADHD). Front Hum Neurosci. 2019;11:13.

    Google Scholar 

  27. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489(7416):391–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Sandin S, Lichtenstein P, Kuja-Halkola R, Hultman CM, Larsson H, Reichenberg A. The heritability of autism spectrum disorder. JAMA. 2017;318(12):1182.

    PubMed  PubMed Central  Google Scholar 

  29. Cross-Disorder Group of the Psychiatric Genomics Consortium. Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nat Genet. 2013;45(9):984–94.

    PubMed Central  Google Scholar 

  30. Geschwind DH. Genetics of autism spectrum disorders. Trends Cogn Sci. 2011;15(9):409–16.

    PubMed  PubMed Central  Google Scholar 

  31. Gilman SR, Iossifov I, Levy D, Ronemus M, Wigler M, Vitkup D. Rare de novo variants associated with autism implicate a large functional network of genes involved in formation and function of synapses. Neuron. 2011;70(5):898–907.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Huguet G, Benabou M, Bourgeron T. The genetics of autism spectrum disorders. In: Sassone-Corsi P, Christen Y, editors. A time for metabolism and hormones. Springer; 2016. p. 101–29.

    Google Scholar 

  33. Demontis D, Walters GB, Athanasiadis G, Walters R, Therrien K, Farajzadeh L, et al. Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains. medRxiv. 2022;15:2022.02.14.22270780.

  34. Cross-Disorder Group of the Psychiatric Genomics Consortium. Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell. 2019;179(7):1469-1482.e11.

    PubMed Central  Google Scholar 

  35. Rommelse NNJ, Franke B, Geurts HM, Hartman CA, Buitelaar JK. Shared heritability of attention-deficit/hyperactivity disorder and autism spectrum disorder. Eur Child Adolesc Psychiatry. 2010;19(3):281–95.

    PubMed  PubMed Central  Google Scholar 

  36. American Psychiatric Association (APA). Diagnostic and statistical manual of mental disorders. 1994.

  37. Association AP. Diagnostic and statistical manual of mental disorders. 2000.

  38. WH. O. The ICD-10 Classification of mental and behavioural disorders: diagnostic criteria for research. 1993;2013.

  39. Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis. Neuroimage. 1999;9(2):179–94.

    CAS  PubMed  Google Scholar 

  40. Fischl B, Sereno MI, Dale AM. Cortical surface-based analysis. Neuroimage. 1999;9(2):195–207.

    CAS  PubMed  Google Scholar 

  41. Ségonne F, Dale AM, Busa E, Glessner M, Salat D, Hahn HK, et al. A hybrid approach to the skull stripping problem in MRI. Neuroimage. 2004;22(3):1060–75.

    PubMed  Google Scholar 

  42. Jovicich J, Czanner S, Greve D, Haley E, van der Kouwe A, Gollub R, et al. Reliability in multi-site structural MRI studies: effects of gradient non-linearity correction on phantom and human data. Neuroimage. 2006;30(2):436–43.

    PubMed  Google Scholar 

  43. Fischl B, Dale AM. Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc Natl Acad Sci. 2000;97(20):11050–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Ecker C, Pretzsch CM, Bletsch A, Mann C, Schaefer T, Ambrosino S, et al. Interindividual differences in cortical thickness and their genomic underpinnings in autism spectrum disorder. Am J Psychiatry. 2021;(7):appi.ajp.2021.2.

  45. Winkler AM, Sabuncu MR, Yeo BTT, Fischl B, Greve DN, Kochunov P, et al. Measuring and comparing brain cortical surface area and other areal quantities. Neuroimage. 2012;61(4):1428–43.

    PubMed  Google Scholar 

  46. Gorgolewski1 KJ, Fox AS, Chang L et al. Tight fitting genes: finding relations between statistical maps and gene expression patterns [Internet]Hamburg. In: 2014. https://f1000research.com/posters/1097120.

  47. Satterstrom FK, Kosmicki JA, Wang J, Breen MS, De Rubeis S, An JY, et al. Large-scale exome sequencing study implicates both developmental and functional changes in the neurobiology of autism. Cell. 2020;180(3):568-584.e23.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Gandal MJ, Zhang P, Hadjimichael E, Walker RL, Chen C, Liu S, et al. Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science (1979). 2018. https://doi.org/10.1126/science.aat8127.

    Article  Google Scholar 

  49. Velmeshev D, Schirmer L, Jung D, Haeussler M, Perez Y, Mayer S, et al. Single-cell genomics identifies cell type–specific molecular changes in autism. Science (1979). 2019;364(6441):685–9. https://doi.org/10.1126/science.aav8130.

    Article  CAS  Google Scholar 

  50. Parikshak NN, Swarup V, Belgard TG, Irimia M, Ramaswami G, Gandal MJ, et al. Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism. Nature. 2016;540(7633):423–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1826.

    PubMed  PubMed Central  Google Scholar 

  52. Watanabe K. FUMA SNP2GENE Tutorial [Internet]. https://fuma.ctglab.nl/tutorial#magma.

  53. Kissel LT, Werling DM. Neural transcriptomic analysis of sex differences in autism spectrum disorder: current insights and future directions. Biol Psychiatry. 2022;91(1):53–60.

    CAS  PubMed  Google Scholar 

  54. Gudbrandsen M, Bletsch A, Mann C, Daly E, Murphy CM, Stoencheva V, et al. Neuroanatomical underpinnings of autism symptomatology in carriers and non-carriers of the 22q11.2 microdeletion. Mol Autism. 2020;11(1):1–15.

    Google Scholar 

  55. Shaw P, Eckstrand K, Sharp W, Blumenthal J, Lerch JP, Greenstein D, et al. Attention-deficit/hyperactivity disorder is characterized by a delay in cortical maturation. Proc Natl Acad Sci. 2007;104(49):19649–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Castellanos FX, Tannock R. Neuroscience of attention-deficit/hyperactivity disorder: the search for endophenotypes. Nat Rev Neurosci. 2002;3(8):617–28.

    CAS  PubMed  Google Scholar 

  57. Coghill D, Sonuga-Barke EJS. Annual research review: categories versus dimensions in the classification and conceptualisation of child and adolescent mental disorders—implications of recent empirical study. J Child Psychol Psychiatry. 2012;53(5):469–89.

    PubMed  Google Scholar 

  58. Frazier TW, Chetcuti L, Al‐Shaban FA, Haslam N, Ghazal I, Klingemier EW, et al. Categorical versus dimensional structure of autism spectrum disorder: a multi‐method investigation. JCPP Adv. 2023.

  59. Matte B, Rohde LA, Grevet EH. ADHD in adults: a concept in evolution. ADHD Atten Deficit Hyperact Disord. 2012;4(2):53–62.

    Google Scholar 

  60. Faraone SV, Biederman J, Mick E. The age-dependent decline of attention deficit hyperactivity disorder: a meta-analysis of follow-up studies. Psychol Med. 2006;36(2):159–65.

    PubMed  Google Scholar 

  61. Dougherty CC, Evans DW, Myers SM, Moore GJ, Michael AM. A comparison of structural brain imaging findings in autism spectrum disorder and attention-deficit hyperactivity disorder. Neuropsychol Rev. 2016;26(1):25–43.

    PubMed  Google Scholar 

  62. Romero-Garcia R, Warrier V, Bullmore ET, Baron-Cohen S, Bethlehem RAI. Synaptic and transcriptionally downregulated genes are associated with cortical thickness differences in autism. Mol Psychiatry. 2019;24(7):1053–64.

    CAS  PubMed  Google Scholar 

  63. Courchesne E, Pramparo T, Gazestani VH, Lombardo MV, Pierce K, Lewis NE. The ASD living biology: from cell proliferation to clinical phenotype. Mol Psychiatry. 2019;24(1):88–107.

    PubMed  Google Scholar 

  64. Parikshak NN, Luo R, Zhang A, Won H, Lowe JK, Chandran V, et al. Integrative functional genomic analyses implicate specific molecular pathways and circuits in autism. Cell. 2013;155(5):1008–21.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge the contributions 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 Bast, Declan G.M. Murphy, Bethany Oakley, Laurence O’Dwyer, Marianne Oldehinkel, Bob Oranje, Gahan Pandina, Antonio M. Persico, 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. The primary contact for the EU-AIMS LEAP Group is Declan G. Murphy (Email: pa-dmurphy@kcl.ac.uk).

Funding

Open Access funding enabled and organized by Projekt DEAL. This project has received funding from the Innovative Medicines Initiative 2 Joint Undertaking (IMI-JU2) under grant agreement No 115300 for the project EU-AUIMS and No 777394 for the project AIMS-2- TRIALS. This Joint Undertaking receives support from the European Union's Horizon 2020 research and innovation programme and EFPIA and AUTISM SPEAKS, Autistica, SFARI. This work has received funding from Horizon Europe (Grant Agreement No. 101057385) and from UK Research and Innovation (UKRI) under the UK government’s Horizon Europe funding guarantee (Grant No. 10039383) (R2D2-MH). CE gratefully acknowledges support from the German Research Foundation (DFG) under the Heisenberg Programme (EC480/1-1 and EC480/2-1). DGM also acknowledges support from the NIHR Maudsley Biomedical Research Centre. DLF is supported by funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement No 101025785.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

The “European Autism Interventions—A Multicenter Study for Developing New Medications (EU-AIMS) Longitudinal European Autism Project (LEAP)” group collaboratively collected the data. LMB conceptualized the study and conceived the analysis plan with support from TS, CG, and CE. LMB conducted the data analysis and wrote the first draft of manuscript. TS and CE provided the software used for the analyses. JL, HS, AB, CMP, BO, EL, DLF, JKB, CFB, ToB, TC, EJHJ, JT, CHC, TB, DGM, and all other authors provided feedback on the manuscript versions and agreed to the final draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lisa M. Berg.

Ethics declarations

Ethics approval and consent to participate

An independent ethics committee at each site approved the study, and written informed consent was obtained for all participants.

Site

Ethics committee

ID/reference no

KCL

London-Central and Queen Square Health Research Authority Research Ethics Committee

13/LO/1156

UCAM

RUNMC

Radboud universitair medisch centrum Instituut Waarborging Kwaliteit en Veiligheid Commissie Mensgebonden Onderzoek Regio Arnhem-Nijmegen (Radboud University Medical Centre Institute Ensuring Quality and Safety Committee on Research Involving Human Subjects Arnhem-Nijmegen)

2013/455

UMCU

CIMH

UMM Universitätsmedizin Mannheim, Medizinische Ethik Kommission II (UMM University Medical Mannheim, Medical Ethics Commission II)

2014-540N-MA

UCBM

Universita Campus Bio-Medica De Roma Comitato Etico (University Campus Bio-Medical Ethics Committee De Roma)

18/14 PAR

ComET CBM

KI

Centrala Etikprovningsnamnden (Central Ethical Review Board)

32-2010

Competing interests

The IMI-JU2 had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Any views expressed are those of the author(s) and not necessarily those of the IMI-JU2. Jan K. Buitelaar has been in the past 3 years a consultant to/member of advisory board of/and/or speaker for Janssen Cilag BV, Takeda/Shire, Roche, Novartis, Medice, Angelini, 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, and royalties. Tobias Banaschewski has served in an advisory or consultancy role for Actelion, Hexal Pharma, Lilly, Medice, Novartis, Oxford outcomes, Otsuka, PCM scientific, Shire and Viforpharma. He received conference support or speaker’s fee by Medice, Novartis and Shire. He is/has been involved in clinical trials conducted by Shire and Viforpharma. He received royalties from Hogrefe, Kohlhammer, CIP Medien, and Oxford University Press. The present work is unrelated to the above grants and relationships. Julian Tillmann is a consultant for F. Hoffmann-La Roche Ltd. Chris H. Chatham is an employee of F. Hoffmann–La Roche. Tony Charman has received research grant support from the Medical Research Council (UK), the National Institute for Health Research, Horizon 2020 and the Innovative Medicines Initiative (European Commission), MQ, Autistica, FP7 (European Commission), the Charles Hawkins Fund, and the Waterloo Foundation. He has served as a consultant to F. Hoffmann-La Roche Ltd and Servier. He has received royalties from Sage Publications and Guilford Publications. Declan G. Murphy has received consultancy fees from Roche and Servier, and grant support from the Medical Research Council (UK), the National Institute for Health Research, Horizon 2020 and the Innovative Medicines Initiative (European Commission). Lisa M. Berg, Tim Schaefer, Caroline Gurr, Johanna Leyhausen, Hanna Seelemeyer, Charlotte M. Pretzsch, Bethany Oakley, Eva Loth, Dorothea L. Floris, Christian F. Beckmann, Emily J.H. Jones, Thomas Bourgeron, and Christine Ecker reported no biomedical financial interests or potential conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Methods. 1. Sample Description. 2. MRI Data Quality Assessments. 3. Effects of Intellectual Disability. 4. Site effects. 5. General Least Square (GLS)-decoding accounting for autocorrelations in spatially embedded transcriptomic maps. Figure S1A. Effects of Intellectual Disability (ID) for measures of cortical thickness (CT). Figure S1B. Effects of Intellectual Disability (ID) for measures of cortical thickness (CT). Figure S1C. Effects of Intellectual Disability (ID) for measures of surface area (SA). Figure S1D. Effects of Intellectual Disability (ID) for measures of surface area (SA). Figure S2A. Effects of acquisition sites for measures of cortical thickness (CT). Figure S2B. Effects of acquisition sites for measures of cortical thickness (CT). Figure S2C. Effects of acquisition sites for measures of surface area (SA). Figure S2D. Effects of acquisition sites for measures of surface area (SA). Figure S3. Significant differences of cortical thickness (CT) and surface area (SA) for the main effects of ASD, ADHD and the ASD-by-ADHD interaction term. Figure S4. Effect size plots for individual model terms. Figure S5A. Interaction plots for all significant clusters for measures of cortical thickness (CT). Figure S5B. Interaction plots for all significant clusters for measures of surface area (SA). Figure S6. Gene enrichments as resulting from the GLS-decoding approach Figure S7. Overlap between significant decoded genes. Table S1. MRI Acquisition Parameters across sites. Table S2. Vertex-wise differences in CT and SA for the main effect of ASD. Table S3. Vertex-wise differences in CT and SA for the main effect of ADHD. Table S4. Vertex-wise differences in CT and SA for the ASD-by-ADHD interaction term. Table S5. Vertex-wise differences in CT for the main effect of ASD + ADHD.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Berg, L.M., Gurr, C., Leyhausen, J. et al. The neuroanatomical substrates of autism and ADHD and their link to putative genomic underpinnings. Molecular Autism 14, 36 (2023). https://doi.org/10.1186/s13229-023-00568-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13229-023-00568-z

Keywords