- 1Pediatric Neurology Unit, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
- 2Dipartimento di Diagnostica per Immagini, Istituto di Radiologia, Università Cattolica del Sacro Cuore, Rome, Italy
- 3UOC Radiologia e Neuroradiologia, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
- 4Medical Physics Unit, Bambino Gesù Children’s Hospital, IRCCS, Rome, Italy
- 5Centro Clinico Nemo, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
- 6Department of Clinical and Experimental Medicine, University of Pisa, Pisa, Italy
- 7Department of Developmental Neuroscience, IRCCS Fondazione Stella Maris, Pisa, Italy
The term autism spectrum disorder (ASD) includes a wide variability of clinical presentation, and this clinical heterogeneity seems to reflect a still unclear multifactorial etiopathogenesis, encompassing different genetic risk factors and susceptibility to environmental factors. Several studies and many theories recognize as mechanisms of autism a disruption of brain development and maturation time course, suggesting the existence of common neurobiological substrates, such as defective synaptic structure and aberrant brain connectivity. Magnetic resonance imaging (MRI) plays an important role in both assessment of region-specific structural changes and quantification of specific alterations in gray or white matter, which could lead to the identification of an MRI biomarker. In this study, we performed measurement of cortical thickness in a selected well-known group of preschool ASD subjects with the aim of finding correlation between cortical metrics and clinical scores to understand the underlying mechanism of symptoms and to support early clinical diagnosis. Our results confirm that recent brain MRI techniques combined with clinical data can provide some useful information in defining the cerebral regions involved in ASD although large sample studies with homogeneous analytical and multisite approaches are needed.
Introduction
According to the Diagnostic and Statistical Manual of Mental Disorders (DSM-5), autism spectrum disorder (ASD) is a neurodevelopmental disorder, characterized by the presence of impaired social communication and unusually repetitive behaviors or restricted interests (American Psychiatric Association, 2013).
However, the term ASD includes a wide variability of clinical presentation (Lord et al., 2000), which seems to reflect a multifactorial etiopathogenesis, still unclear, encompassing different genetic risk factors and susceptibility to environmental factors (Dietert et al., 2011; Tordjman et al., 2014).
Over the years, disruption of brain development and the maturation time course has been recognized as mechanism of autism (Amaral et al., 2008), suggesting the existence of common neurobiological substrates, such as defective synaptic structure and aberrant brain connectivity (Ecker et al., 2016).
Therefore, several studies have assessed the utility of brain magnetic resonance imaging (MRI) in the evaluation of region-specific structural changes and quantification of specific alterations in gray (GM) or white matter (WM) [in terms of cortical thickness, gyrification index, GM surface and volume (Pagnozzi et al., 2018)].
The identification of an MRI biomarker could provide insight into the underlying mechanisms of symptoms and could additionally provide crucial support in early clinical diagnosis and dividing patients, sharing common features, into different groups to tailor specific interventions.
In literature, consensus exists on increased growth of total GM and WM volume (Courchesne et al., 2001; Sparks et al., 2002; Hardan et al., 2006; Hazlett et al., 2011; Nordahl et al., 2012; Mahajan and Mostofsky, 2015; Ismail et al., 2016; Libero et al., 2016; Lucibello et al., 2019), which involves specific brain regions, such as temporal lobes (Lai et al., 2013; Ecker et al., 2014; Riddle et al., 2017; Pappaianni et al., 2018), frontal lobes (Ecker et al., 2014; Zhou et al., 2014; Eilam-Stock et al., 2016; Patriquin et al., 2016), or both (Katuwal et al., 2015; Postema et al., 2019).
A comparable number of studies finds opposite results, reporting decreased volumes of the frontotemporal (Hardan et al., 2006; Ecker et al., 2010; Mak-Fan et al., 2012), and parietal cortex (Ecker et al., 2010), significantly smaller right (Haznedar et al., 1997), and bilateral (Ecker et al., 2010; Jiao et al., 2010), anterior cingulate gyrus, and decreased cortical volume in the orbitofrontal cortex bilaterally (Ecker et al., 2012).
Recently, many studies have focused their efforts on measurement of cortical thickness (Doyle-Thomas et al., 2013; Khundrakpam et al., 2017; Pereira et al., 2018; Prigge et al., 2018), and among them, the one with the largest ASD population (Van Rooij et al., 2018), showed promising results. Specifically, the authors found complex developmental trajectories involving different brain regions with significant differences in terms of increased cortical thickness in the frontal cortex and decreased thickness in the temporal cortex between ASD patients and controls during adolescence. Other authors focus their study on cortical gyrification, suggesting the involvement of both genetic and non-genetic factors in definition of cortical gyral and sulcal patterns (Lohmann et al., 1999; White et al., 2002; Kates et al., 2009; Kremen et al., 2010; Mata et al., 2010; Docherty et al., 2015; Kuhn et al., 2016; Bernardoni et al., 2018; Duan et al., 2020; Kruggel and Solodkin, 2020). Although analysis of gyrification has led to the identification of different cortical regions involved in ASD, there is high heterogeneity across studies.
However, no studies have investigated the correlation between clinical, genetic, and radiological findings in a well-selected pediatric population.
Aim of the Study
As part of an ongoing project collecting ASD preschooler data sets, we retrospectively selected MRI data sets and applied a semiautomatic brain segmentation methodology to investigate cortical thickness and gyrification. Correlations of cortical indexes and clinical features were subsequently performed.
Materials and Methods
Study Design and Clinical Assessment
A retrospective study was designed. The study group included 39 preschool-age children regularly followed at the Child Neurology Unit of the Gemelli Hospital (Rome, Italy) with a diagnosis of ASD. As part of our routine assessment since March 2016, patients referred to our unit with a suspect clinical diagnosis of ASD undergo a detailed clinical assessment, including a neurological examination, Leiter or Wechsler scales according to age and cooperation, a comprehensive neuropsychiatric assessment using parent-reported questionnaires’ (Child Behavior Checklist, CBCL), and autism-specific diagnostic tools, specifically Autism Diagnostic Interview Revised (ADI-R) and Autism Diagnostic Observation Schedule, second edition (ADOS2) (Lord et al., 1989, 1994). In addition, MRI is also routinely performed.
To have a relatively homogeneous cohort, the following exclusion criteria were used: (1) dysmorphic features with specific genetic syndrome identification, (2) presence of severe epilepsy, (3) presence of cerebral palsy or other major neurological signs, (4) malformations or other lesions at MRI.
The scores of all subscales of the ADI-R (social interaction, communication and language, restricted and repetitive behaviors) and ADOS2 (social affect score, restricted and repetitive behavior score, total score) were correlated with cerebral cortical thickness and cerebral cortical gyrification with the aim of finding a plausible neural substrate for the domains of deficit: impaired social communication and unusually repetitive behaviors or restricted interests.
Study Design and Neuroradiological Protocol
All participants underwent MRI with a 1.5T Philips Ingenia Scanner (Philips Healthcare, Eindhoven, Netherlands). The sequence used for postprocessing was a T1-weighted 3-D-TFE in a sagittal orientation (TR = 9.8 ms, TE = 4.6 ms with a delay time of 650 ms after a 180° prepulse, flip angle = 10°, FOV = 200 mm × 222 mm, 1.0 mm slice thickness with no gaps, total of 150 slices per slab, matrix size = 200 × 222, NSA = 2 with an in-plane resolution of 1.0 × 1.0 mm2). This sequence is routinely acquired in our MRI protocol for children older than 2 years of age. ASD subjects were sedated with a general anesthesia with a halogenated agent while spontaneously breathing; contrast agent injection was never required. The written informed consent from a parent or guardian of children was obtained.
The quality of the structural MRI data was rated by an experienced neuroimaging researcher (TV) on a 3-point rating scale: 0 = no motion artifacts, excellent quality; 1 = few motion artifacts, fair quality; and 2 = moderate/severe motion artifacts, poor quality. Only data sets with scores of 0 were considered of adequate quality for research purposes. In case of imaging artifacts (ghosting, aliasing, chemical shift, and distortion), patients were excluded.
The raw 3-D T1 MRI data underwent automated processing for surface-based cortex reconstruction and volumetric segmentation using Freesurfer image analysis software (version 6.0.0), which is documented and freely available for download online1, installed on an OSX El Capitan 10.11.6.
The technical details of these procedures are described in prior publications; briefly, the processing pipeline includes motion correction and averaging (Reuter et al., 2010) of volumetric T1 weighted images, removal of non-brain tissue (Ségonne et al., 2004), automated Talairach transformation, segmentation of the subcortical WM and deep GM volumetric structure (including hippocampus, amygdala, caudate, putamen, ventricles) (Fischl et al., 2002, 2004) intensity normalization (Sled et al., 1998), tessellation of the GM/WM boundary, automated topology correction (Fischl et al., 2001; Ségonne et al., 2007), and surface deformation following intensity gradients to optimally place the GM/WM and GM/CSF borders (Dale and Sereno, 1993; Dale et al., 1999; Fischl and Dale, 2000). Cortical thickness was quantified as the closest distance from the GM/WM boundary to the GM/CSF boundary at each vertex. Cortical parcelation and thickness estimations were based on the Desikan–Killiany Atlas (Desikan et al., 2006), resulting in average cortical thickness in 34 cortical parcels per hemisphere. Procedures for the measurement of cortical thickness have been validated against histological analysis (Rosas et al., 2002) and manual measurements (Kuperberg et al., 2003; Salat et al., 2004).
Local gyrification index (LGI) was computed using the approach proposed by Lyu et al. (2018)2 that quantifies cortical gyrification within sulcal and gyral regions using a spatially varying kernel shape, able to adaptively encode cortical folding patterns. The proposed LGI is then computed within the adaptive kernel as a ratio of the cortical surface area and a fixed area on the outer hull (ρ = 316 mm2) (Lyu et al., 2018; Płonka et al., 2020).
The Ethics Committee of Policlinico Gemelli, Catholic University of Sacred Hearth, examined and approved the project study (February 2018).
Statistical Analysis
We investigated differences in cortical parameter distributions (i.e., cortical thickness, CT, and local gyrification, LGI) among subjects grouped on the basis of ADOS2 and ADI-R domains scores. To this purpose, we first mapped vertex-wise CT and LGI values on a common spherical coordinate system using spherical transformation, and then we assessed differences among groups using permutation tests (1,000 permutations for all tests) based on the t-statistics, performed with the Permutation Analysis of Linear Models (PALM) FSL package (3 version 6.0). In particular, we used group age as a covariate to produce threshold-free cluster enhancement (TFCE) statistical maps, where the initial raw statistical images were enhanced using both the intensity of the data point and information from neighboring voxels (Mensen and Khatami, 2013). We detected group differences on both family-wise error (FWE) corrected and uncorrected p-value maps. Moreover, correlation analyses were evaluated vertex-wise between cortical parameters (CT and LGI) and several clinical variables, including social interaction (SI), communication (COM), and total scores (TOT, defined as SI + COM), together with their autism diagnostic interview (ADI) counterparts, i.e., ADI_SI, ADI_COM, repetition ADI (ADI_rep) and ADI_TOT as sum of the previously mentioned scores. Correlation analyses were performed testing Pearson correlation with PALM permutation test (1000 permutations).
Results
After image quality check, five patients were excluded. Data from 34 patients, 4 females (11.8%) and 30 males (88.2%), were eventually included in the analysis. Their age range was between 3.4 and 6 years. All subjects of the sample selected had normal MRI findings and normal FMR1 and CGH array analysis.
Seven patient (20.6%) had a normal intellectual quotient, 12 (35.3%) a mild intellectual disability, and three (8.8%) a moderate intellectual disability.
About ADOS valuation, 5/34 ASD patients (14.7%) had a low level of severity, 20/34 (58.8%) had a medium level of severity, and 5/34 (14.7%) had a high level of severity (Table 1).
Correlations
Autism Diagnostic Interview – Social Interaction Domain
Statistical analysis showed a correlation between the subtest of ADI that investigates social interaction aims and thickness of different left-brain regions: fusiform (p = 0.046), lingual (0.049), posterior cingulate (0.048), pre-cuneus (p = 0.046), superior parietal (p = 0.033), inferior parietal (p = 0.033), superior temporal (p = 0.043), inferior temporal (p = 0.041), middle temporal (0.037), temporal pole (p = 0.041), and lateral occipital (0.041).
A correlation was found also between the subtest of ADI that investigates social interaction aims and gyrification of different cortical regions: in left hemisphere insula (p = 0.01), rostral anterior cingulate (0.008), pars orbitalis (p = 0.009), superior frontal (p = 0.007), medial orbito-frontal (p = 0.006), lateral orbito-frontal (p = 0.008), precentral (p = 0.008), postcentral (p = 0.004), supramarginal (p = 0.007), transverse temporal (p = 0.009); in right hemisphere, superior frontal (p = 0.008), precentral (p = 0.007), paracentral (p = 0.009) (Figure 1).
Figure 1. Statistical results mapped on the surface template. (A) Family-wise error (FWE) corrected p-value maps for cortical thickness ADI-social interaction correlation. Significance was set at 0.05. (B) Uncorrected p-value maps for gyrification ADI-social interaction correlation. Significance was set at 0.01.
Autism Diagnostic Interview – Communication Domain
Statistical analysis showed a correlation between the subtest of ADI that investigates communication aims and gyrification of different left cortical regions: caudal anterior cingulate (p = 0.007), posterior cingulate (p = 0.008) (Figure 2).
Figure 2. Uncorrected p-value maps for gyrification ADI-communication correlation. Significance was set at 0.01.
Autism Diagnostic Interview – Repetitive Behavior Domain
Statistical analysis showed a correlation between the subtest of ADI that investigates social interaction aims and gyrification of different left regions: fusiform (p = 0.008), lingual (p = 0.008), para-hippocampal (p = 0.008), middle temporal (p = 0.006), inferior temporal (p = 0.005), lateral occipital (p = 0.008), but also with different right region: superior frontal (p = 0.006), frontal pole (p = 0.008), rostral middle frontal (p = 0.004), caudal middle frontal (p = 0.007), and precentral (p = 0.006) (Figure 3).
Figure 3. Uncorrected p-value maps for gyrification ADI-repetitive behavior correlation. Significance was set at 0.01.
Autism Diagnostic Observation Schedule – Social Affect Score
Statistical analysis showed a correlation between the subtest of ADOS that investigates restricted and repetitive behavior aims and thickness of different left cortical regions: pre-cuneus (p = 0.007), posterior cingulate (p = 0.006), lateral orbito-frontal (p = 0.006), medial orbito-frontal (p = 0.005), precentral (p = 0.007), paracentral (p = 0.009), postcentral (p = 0.007), and superior parietal (p = 0.007).
A correlation was found also between the subtest of ADOS that investigates restricted and repetitive behavior aims and gyrification of left cortical regions: middle temporal (p = 0.008) and superior temporal (p = 0.008).
Autism Diagnostic Observation Schedule – Restricted and Repetitive Behaviors Score
Statistical analysis showed a correlation between the subtest of ADOS that investigates social affect aims and gyrification of different left cortical regions: fusiform (p = 0.008), lingual (p = 0.005), entorhinal (p = 0.008), para-hippocampal (p = 0.008), posterior cingulate (p = 0.009), superior frontal (p = 0.009), caudal middle frontal (p = 0.004), rostral middle frontal (p = 0.004), inferior temporal (p = 0.008) and also with different right regions: lingual (p = 0.008), para-hippocampal (p = 0.009), pars triangularis (p = 0.007), pre-cuneus (p = 0.009), posterior cingulate (p = 0.008), rostral anterior cingulate (p = 0.006), caudal anterior cingulate (p = 0.009), pars orbitalis (p = 0.004), superior frontal (p = 0.009), rostral middle frontal (p = 0.006), medial orbito-frontal (p = 0.006), and caudal middle frontal (p = 0.008) (Figure 4).
Figure 4. Uncorrected p-value maps for gyrification ADOS-restricted and repetitive behavior score correlation. Significance was set at 0.01.
Autism Diagnostic Observation Schedule – Total Score
Statistical analysis showed a correlation between the total score of ADOS and gyrification of left cortical regions: middle temporal (p = 0.009) and inferior temporal (p = 0.007) (Figure 5).
Figure 5. Statistical results mapped on the surface template. (A) Uncorrected p-value maps for cortical thickness-ADOS_total score correlation. Significance was set at 0.05. (B) Uncorrected p-value maps for gyrification-ADOS-total score correlation. Significance was set at 0.01.
Discussion
Autism is a complex and heterogeneous neurodevelopmental disorder with widely varied clinical features and a multitude of possible etiological factors.
Emerging evidence supports that ASD undergoes an atypical trajectory of brain maturation (Doyle-Thomas et al., 2013; Ecker et al., 2014) or abnormal lateralization of brain networks (Conti et al., 2015). Findings of atypical brain morphology in ASD, however, are highly heterogeneous and there are still not clear neuroanatomical markers to accurately identify individuals with ASD.
Whereas most neuroradiological studies on autism have focused on finding a specific pattern of brain involvement, only a few studies aimed to establish a correlation between clinical severity and cortical thickness (Doyle-Thomas et al., 2013). Additionally, many of these studies exclusively investigate well-known cortical areas involved in specific function, whereas we have tried to assess the correlation of several cortical areas to verify some correspondence with already known brain circuits. Whereas some results have reinforced well-known concepts already widely described in literature, others are quite difficult to frame.
The correlation between cingulate and frontal regions, particularly the left medial orbito-frontal area, with a social affect score of both ADOS and ADI-R scales is not surprising as the involvement of this area in social-affective processing is already described as part of the ventral social affective processing system (Belger et al., 2011). These results are consistent with what is already known about cortical gyrification morphology in neurodevelopmental disorders (Kohli et al., 2019).
Less clear is the correlation with other areas, such as the left lateral-occipital area, even if other authors have already shown this result in general in ASD patients (Kohli et al., 2019).
Among the results obtained, it is useful to underline a partial congruence between the results of the ADOS2 and ADI-R subtests investigating social affective behaviors and cortical gyrification index. This highlights how the two diagnostic tests support each other and how it is essential to always try to perform a complete patient assessment.
In the same way, it is easy to understand the correlation found between the restricted and repetitive behavior indices of both scales with gyrification of para-hippocampal, temporal, and middle frontal area.
Many human neuroimaging studies have indeed shown increased activation in ventral striatum and ventromedial frontal cortex in response to unexpected negative feedback that implies a change of behavior strategy (Glascher et al., 2009; D’Cruz et al., 2016; Sharda et al., 2017).
Conversely, it is unclear how the cingulate cortex could be involved in communication tasks although some neuroimaging studies suggest the association with this region and indirect language and social inferential capacity, such as the comprehension of non-literal language or construction-based pragmatic information (D’Cruz et al., 2011).
In general, the assessment of cortical areas based on the ADOS and ADI scores cannot be univocal. These scores, indeed, assess other skills, such as language or different non-verbal communication ability that could be underpinned by more cortical areas. Certainly, many of these areas, such as temporal and frontal or middle temporal areas, are already described as involved in patients with ASD.
This study confirms and expands previous results while reducing the variability of cohort patients as, in other studies a systematic selection of ASD subjects was not performed, and combining imaging with more detailed clinical data. The well-selected population of our study allowed us to reduce the variability observed between different ages.
Exclusion criteria, such as dysmorphic features or epilepsy or other major neurological signs, have led to less variability related to secondary factors, and focusing on only preschool-age patients has reduced the bias linked to the noted changes in the thickness of the cerebral cortex related to age (Shibata et al., 2010).
Neuroimaging studies indeed report numerous region-specific alterations in cortical thickness in patients with ASD. However, there are many inconsistent findings, and this is probably due to atypical CT developmental trajectories in ASD, characterized by decreased cortical thinning during early adolescence and increased thinning at later stages, involving mostly frontal and parietal areas (Mensen and Khatami, 2013; Ecker et al., 2014).
Limitations of our study include a small sample size, which reduces the power of statistical analysis.
Though the study has investigated the brain structure–function correlation with the aim of catching a possible clinical-related sensitivity marker at this age, the lack of controls constitutes a limitation of the study. The need for sedation during MRI scans, indeed, limits the opportunity to collect age-matched typical subjects.
Moreover, because pediatric brains are different in size and shape from adult brains commonly used as frameworks for spatial normalization (e.g., Talairach space), specific cortical areas could show correlation of deformation with age. In this context, several studies assessed the influence of age on various spatial normalization parameters (Wilke et al., 2002; Fonov et al., 2011), and consequently, our results should be taken with caution.
In conclusion, our results confirm that recent brain MRI techniques integrated to clinical data can provide some useful information in defining the cerebral regions involved in ASD although large sample studies characterized by homogeneous analytical and multisite approaches are needed.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of Fondazione Policlinico Gemelli Hospital – Catholic University, Rome, Italy. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
RB, SL, and EM contributed to conception and design of the study. SL, TV, and GB organized the database. EP, MC, DC, and RD’A contributed to investigation. AN and ML performed the statistical analysis. SL wrote the first draft of the manuscript. TV and AN wrote sections of the manuscript. PM, CC, and EM contributed to the supervision of manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was partially supported by the Current Research (RC) of Italian Ministry of Health, IRCCS Stella Maris Foundation (FSM), entitled “Early functional and neurobiological markers of diseases associated to neurodevelopmental disorders” and 5 for thousands of IRCCS FSM (RB).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
- ^ http://surfer.nmr.mgh.harvard.edu/
- ^ https://github.com/ilwoolyu/LocalGyrificationIndex
- ^ https://fsl.fmrib.ox.ac.uk/fsl/fslwiki
References
Amaral, D. G., MillsSchumann, C., and WuNordahl, C. (2008). Neuroanatomy of autism. Trends Neurosci. 31, 137–145.
American Psychiatric Association (2013). Diagnostic and Statistical Manual of Mental Disorders, 5th Edn. Washington DC: American Psychiatric Association.
Belger, A., Carpenter, K. L. H., Yucel, G., Cleary, K. M., and Donkers, F. (2011). The neural circuitry of autism. Neurotox. Res. 20, 201–214.
Bernardoni, F., King, J. A., Geisler, D., Birkenstock, J., Tam, F. I., Weidner, K., et al. (2018). Nutritional status affects cortical folding: lessons learned from anorexia nervosa. Biol. Psychiatry 84, 692–701. doi: 10.1016/j.biopsych.2018.05.008
Conti, S., Calderoni, S., Gaglianese, A., Pannek, K., Mazzotti, S., and Rose, S. (2015). Lateralization of brain networks and clinical severity in toddlers with autism spectrum disorder: a HARDI diffusion MRI study. Autism Res. 9, 382–392. doi: 10.1002/aur.1533
Courchesne, E., Karns, C. M., Davis, H. R., Ziccardi, R., Carper, R. A., Tigue, Z. D., et al. (2001). Unusual brain growth patterns in early life in patients with autistic disorder: an MRI study. Neurology 57, 245–254. doi: 10.1212/wnl.57.2.245
Dale, A. M., Fischl, B., and Sereno, M. I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9, 179–194. doi: 10.1006/nimg.1998.0395
Dale, A. M., and Sereno, M. I. (1993). Improved localizadon of cortical activity by combining EEG and MEG with MRI cortical surface reconstruction: a linear approach. J. Cogn. Neurosci. 5, 162–176. doi: 10.1162/jocn.1993.5.2.162
D’Cruz, A., Ragozzino, M. E., Mosconi, M. W., Pavuluri, M. N., and Sweeney, J. A. (2011). Human reversal learning under conditions of certain versus uncertain outcomes. Neuroimage 56, 315–322. doi: 10.1016/j.neuroimage.2011.01.068
D’Cruz, A.-M., Mosconi, M. W., Ragozzino, M. E., Cook, E. H., and Sweeney, J. A. (2016). Alterations in the functional neural circuitry supporting flexible choice behavior in autism spectrum disorders. Transl. Psychiatry 6:e916. doi: 10.1038/tp.2016.161
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021
Dietert, R. R., Dietert, J. M., and Dewitt, J. C. (2011). Environmental risk factors for autism. Emerg. Health Threats J. 4:7111. doi: 10.3402/ehtj.v4i0.7111
Docherty, A. R., Hagler, D. J. Jr., Panizzon, M. S., Neale, M. C., Eyler, L. T., Fennema-Notestine, C., et al. (2015). Does degree of gyrification underlie the phenotypic and genetic associations between cortical surface area and cognitive ability? Neuroimage 106, 154–160. doi: 10.1016/j.neuroimage.2014.11.040
Doyle-Thomas, K. A. R., Duerden, E. G., Taylor, M. J., Lerchc, J. P., Sooryad, L. V., Wangd, A. T., et al. (2013). Effects of age and symptomatology on cortical thickness in autism spectrum disorders. Res. Autism Spectr. Disord. 7, 141–150. doi: 10.1016/j.rasd.2012.08.004
Duan, D., Xia, S., Rekik, I., Wu, Z., Wang, L., Lin, W., et al. (2020). Share individual identification and individual variability analysis based on cortical folding features in developing infant singletons and twins. Hum. Brain Mapp. 41, 1985–2003. doi: 10.1002/hbm.24924
Ecker, C., Andrews, D., Dell’Acqua, F., Daly, E., Murphy, C., Catani, M., et al. (2016). Relationship between cortical gyrification, white matter connectivity, and autism spectrum disorder. Cereb. Cortex 26, 3297–3309. doi: 10.1093/cercor/bhw098
Ecker, C., Marquand, A., Mourão-Miranda, J., Johnston, P., Daly, E. M., Brammer, M. J., et al. (2010). Describing the brain in autism in five dimensions–magnetic resonance imaging-assisted diagnosis of autism spectrum disorder using a multiparameter classification approach. J. Neurosci. 30, 10612–10623. doi: 10.1523/JNEUROSCI.5413-09.2010
Ecker, C., Shahidiani, A., Feng, Y., Daly, E., Murphy, C., D’Almeida, V., et al. (2014). The effect of age, diagnosis, and their interaction on vertex-based measures of cortical thickness and surface area in autism spectrum disorder. J. Neural Transm. 121, 1157–1170. doi: 10.1007/s00702-014-1207-1
Ecker, C., Suckling, J., Deoni, S. C., Lombardo, M. V., Bullmore, E. T., Baron-Cohen, S., et al. (2012). Brain anatomy and its relationship to behavior in adults with autism spectrum disorder: a multicenter magnetic resonance imaging study. Arch. Gen. Psychiatry 69, 195–209. doi: 10.1001/archgenpsychiatry.2011.1251
Eilam-Stock, T., Wu, T., Spagna, A., Egan, L. J., and Fan, J. (2016). Neuroanatomical alterations in high-functioning adults with autism spectrum disorder. Front. Neurosci. 10:237. doi: 10.3389/fnins.2016.00237
Fischl, B., and Dale, A. M. (2000). Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc. Natl. Acad. Sci. U.S.A. 97, 11050–11055. doi: 10.1073/pnas.200033797
Fischl, B., Liu, A., and Dale, A. M. (2001). Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. IEEE Trans. Med. Imaging 20, 70–80. doi: 10.1109/42.906426
Fischl, B., Salat, D. H., Busa, E., Albert, M., Dieterich, M., Haselgrove, C., et al. (2002). Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron 33, 341–355. doi: 10.1016/s0896-6273(02)00569-x
Fischl, B., Salat, D. H., van der Kouwe, A. J., Makris, N., Ségonne, F., Quinn, B. T., et al. (2004). Sequence-independent segmentation of magnetic resonance images. Neuroimage 23(Suppl. 1) S69–S84. doi: 10.1016/j.neuroimage.2004.07.016
Fonov, V., Evans, A. C., Botteron, K., Robert Almli, C., McKinstry, R. C., and Collins, D. L. (2011). Unbiased average age-appropriate atlases for pediatric studies. Neuroimage 54, 313–327. doi: 10.1016/j.neuroimage.2010.07.033
Glascher, J., Hampton, A. N., and O’Doherty, J. P. (2009). Determining a role for ventromedial prefrontal cortex in encoding action-based value signals during reward-related decision making. Cereb. Cortex 19, 483–495. doi: 10.1093/cercor/bhn098
Hardan, A. Y., Muddasani, S., Vemulapalli, M., Keshavan, M. S., and Minshew, N. J. (2006). An MRI study of increased cortical thickness in autism. Am. J. Psychiatry 163, 1290–1292. doi: 10.1176/ajp.2006.163.7.1290
Hazlett, H. C., Poe, M. D., Gerig, G., Styner, M., Chappell, C., Smith, R. G., et al. (2011). Early brain overgrowth in autism associated with an increase in cortical surface area before age 2 years. Arch. Gen. Psychiatry 68, 467–476. doi: 10.1001/archgenpsychiatry.2011.39
Haznedar, M. M., Buchsbaum, M. S., Metzger, M., Solimando, A., Spiegel-Cohen, J., and Hollander, E. (1997). Anterior cingulate gyrus volume and glucose metabolism in autistic disorder. Am. J. Psychiatry 154, 1047–1050. doi: 10.1176/ajp.154.8.1047
Ismail, M. M., Keynton, R. S., Mostapha, M. M., ElTanboly, A. H., Casanova, M. F., Gimel’farb, G. L., et al. (2016). Studying autism spectrum disorder with structural and diffusion magnetic resonance imaging: a survey. Front. Hum. Neurosci. 10:211. doi: 10.3389/fnhum.2016.00211
Jiao, Y., Chen, R., Ke, X., Chu, K., Lu, Z., and Herskovits, E. H. (2010). Predictive models of autism spectrum disorder based on brain regional cortical thickness. Neuroimage 50, 589–599. doi: 10.1016/j.neuroimage.2009.12.047
Kates, W. R., Ikuta, I., and Burnette, C. P. (2009). Gyrification patterns in monozygotic twin pairs varying in discordance for autism. Autism Res. 2, 267–278. doi: 10.1002/aur.98
Katuwal, G. J., Cahill, N. D., Baum, S. A., and Michael, A. M. (2015). The predictive power of structural MRI in autism diagnosis. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2015, 4270–4273. doi: 10.1109/EMBC.2015.7319338
Khundrakpam, B. S., Lewis, J. D., Kostopoulos, P., Carbonell, F., and Evans, A. C. (2017). Cortical thickness abnormalities in autism spectrum disorders through late childhood, adolescence, and adulthood: a large-scale MRI study. Cereb. Cortex 27, 1721–1731. doi: 10.1093/cercor/bhx038
Kohli, J. S., Kinnear, M. K., Fong, C. H., Fishman, I., Carper, R. A., and Müller, R.-A. (2019). Local cortical gyrification is increased in children with autism spectrum disorders, but decreases rapidly in adolescents. Cereb. Cortex 29, 2412–2423. doi: 10.1093/cercor/bhy111
Kremen, W. S., Prom-Wormley, E., Panizzon, M. S., Eyler, L. T., Fischl, B., Neale, M. C., et al. (2010). Genetic and environmental influences on the size and specific brain regions in midlife: the VETSA MRI study. Neuroimage 49, 1213–1223. doi: 10.1016/j.neuroimage.2009.09.043
Kruggel, F., and Solodkin, A. (2020). Heritability of structural patterning in the human cerebral cortex. Neuroimage 221:117169. doi: 10.1016/j.neuroimage.2020.117169
Kuhn, S., Witt, C., Banaschewski, T., Barbot, A., Barker, G. J., Buchel, C., et al. (2016). From mother to child: orbitofrontal cortex gyrification and changes of drinking behaviour during adolescence. Addict. Biol. 21, 700–708. doi: 10.1111/adb.12240
Kuperberg, G. R., Broome, M. R., McGuire, P. K., David, A. S., Eddy, M., Ozawa, F., et al. (2003). Regionally localized thinning of the cerebral cortex in schizophrenia. Arch. Gen. Psychiatry 60, 878–888. doi: 10.1001/archpsyc.60.9.878
Lai, M. C., Lombardo, M. V., Suckling, J., Ruigrok, A. N., Chakrabarti, B., Ecker, C., et al. (2013). Biological sex affects the neurobiology of autism. Brain 136(Pt 9), 2799–2815. doi: 10.1093/brain/awt216
Libero, L. E., Burge, W. K., Deshpande, H. D., Pestilli, F., and Kana, R. K. (2016). White matter diffusion of major fiber tracts implicated in autism spectrum disorder. Brain Connect. 6, 691–699. doi: 10.1089/brain.2016.0442
Lohmann, G., Von Cramon, D. Y., and Steinmetz, H. (1999). Sulcal variability in twins. Cereb. Cortex 9, 754–763. doi: 10.1093/cercor/9.7.754
Lord, C., Risi, S., Lambrecht, L., Cook, E. H. Jr., Leventhal, B. L., DiLavore, P. C., et al. (2000). The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J. Autism Dev. Disord. 30, 205–223.
Lord, C., Rutter, M., and Le Couteur, A. (1994). Autism diagnostic interview-revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J. Autism Dev. Disord. 24, 659–685. doi: 10.1007/BF02172145
Lord, C., Rutter, M., Goode, S., Heemsbergen, J., Jordan, H., Mawhood, L., et al. (1989). Autism diagnostic observation schedule: a standardized observation of communicative and social behavior. J. Autism Dev. Disord. 19, 185–212. doi: 10.1007/BF02211841
Lucibello, S., Verdolotti, T., Giordano, F. M., Lapenta, L., Infante, A., Piludu, F., et al. (2019). Brain morphometry of preschool age children affected by autism spectrum disorder: correlation with clinical findings. Clin. Anat. 32, 143–150. doi: 10.1002/ca.23252
Lyu, I., Kim, S. H., Girault, J. B., Gilmore, J. H., and Styner, M. A. (2018). A cortical shape-adaptive approach to local gyrification index. Med. Image Anal. 48, 244–258. doi: 10.1016/j.media.2018.06.009
Mahajan, R., and Mostofsky, S. H. (2015). Neuroimaging endophenotypes in autism spectrum disorder. CNS Spectr. 20, 412–426.
Mak-Fan, K. M., Taylor, M. J., Roberts, W., and Lerch, J. P. (2012). Measures of cortical grey matter structure and development in children with autism spectrum disorder. J. Autism Dev. Disord. 42, 419–427. doi: 10.1007/s10803-011-1261-6
Mata, I., Perez-Iglesias, R., Roiz-Santianez, R., Tordesillas-Gutierrez, D., Pazos, A., Gutierrez, A., et al. (2010). Gyrification brain abnormalities associated with adolescence and early-adulthood cannabis use. Brain Res. 1317, 297–304. doi: 10.1016/j.brainres.2009.12.069
Mensen, A., and Khatami, R. (2013). Advanced EEG analysis using threshold-free cluster-enhancement and non-parametric statistics. Neuroimage 67, 111–118. doi: 10.1016/j.neuroimage.2012.10.027
Nordahl, C. W., Scholz, R., Yang, X., Buonocore, M. H., Simon, T., Rogers, S., et al. (2012). Increased rate of amygdala growth in children aged 2 to 4 years with autism spectrum disorders: a longitudinal study. Arch. Gen. Psychiatry 69, 53–61. doi: 10.1001/archgenpsychiatry.2011.145
Pagnozzi, A. M., Conti, E., Calderoni, S., Fripp, J., and Rose, S. E. (2018). A systematic review of structural MRI biomarkers in autism spectrum disorder: a machine learning perspective. Int. J. Dev. Neurosci. 71, 68–82. doi: 10.1016/j.ijdevneu.2018.08.010
Pappaianni, E., Siugzdaite, R., Vettori, S., Venuti, P., Job, R., and Grecucci, A. (2018). Three shades of grey: detecting brain abnormalities in children with autism using source-, voxel- and surface based morphometry. Eur. J. Neurosci. 47, 690–700. doi: 10.1111/ejn.13704
Patriquin, M. A., DeRamus, T., Libero, L. E., Laird, A., and Kana, R. K. (2016). Neuroanatomical and neurofunctional markers of social cognition in autism spectrum disorder. Hum. Brain Mapp. 37, 3957–3978. doi: 10.1002/hbm.23288
Pereira, A. M., Campos, B. M., Coan, A. C., Pegoraro, L. F., de Rezende, T. J. R., Obeso, I., et al. (2018). Differences in cortical structure and functional MRI connectivity in high functioning autism. Front. Neurol. 9:539. doi: 10.3389/fneur.2018.00539
Płonka, O., Krześniak, A., and Adamczyk, P. (2020). Analysis of local gyrification index using a novel shape-adaptive kernel and the standard FreeSurfer spherical kernel—evidence from chronic schizophrenia outpatients. Heliyon 6:e04172. doi: 10.1016/j.heliyon.2020.e04172
Postema, M. C., Rooij, D. V., Anagnostou, E., Arango, C., Auzias, G., and Behrmann, M. (2019). Altered structural brain asymmetry in autism spectrum disorder in a study of 54 datasets. Nat. Commun. 10:4958.
Prigge, M. B. D., Bigler, E. D., Travers, B. G., Froehlich, A., Abildskov, T., Anderson, J. S., et al. (2018). Social Responsiveness Scale (SRS) in relation to longitudinal cortical thickness changes in autism spectrum disorder. J. Autism Dev. Disord. 48, 3319–3329. doi: 10.1007/s10803-018-3566-1
Reuter, M., Rosas, H. D., and Fischl, B. (2010). Highly accurate inverse consistent registration: a robust approach. Neuroimage 53, 1181–1196. doi: 10.1016/j.neuroimage.2010.07.020
Riddle, K., Cascio, C. J., and Woodward, N. D. (2017). Brain structure in autism: a voxel-based morphometry analysis of the Autism Brain Imaging Database Exchange (ABIDE). Brain Imaging Behav. 11, 541–551. doi: 10.1007/s11682-016-9534-5
Rosas, H. D., Liu, A. K., Hersch, S., Glessner, M., Ferrante, R. J., Salat, D. H., et al. (2002). Regional and progressive thinning of the cortical ribbon in Huntington’s disease. Neurology 58, 695–701. doi: 10.1212/wnl.58.5.695
Salat, D. H., Buckner, R. L., Snyder, A. Z., Greve, D. N., Desikan, R. S., Busa, E., et al. (2004). Thinning of the cerebral cortex in aging. Cereb. Cortex. 14, 721–730. doi: 10.1093/cercor/bhh032
Ségonne, F., Dale, A. M., Busa, E., Glessner, M., Salat, D., Hahn, H. K., et al. (2004). A hybrid approach to the skull stripping problem in MRI. Neuroimage 22, 1060–1075. doi: 10.1016/j.neuroimage.2004.03.032
Ségonne, F., Pacheco, J., and Fischl, B. (2007). Geometrically accurate topology-correction of cortical surfaces using nonseparating loops. IEEE Trans. Med. Imaging 26, 518–529. doi: 10.1109/TMI.2006.887364
Sharda, M., Foster, N. E. V., Tryfon, A., Doyle-Thomas, K. A. R., Ouimet, T., Anagnostou, E., et al. (2017). Language ability predicts cortical structure and covariance in boys with autism spectrum disorder. Cereb. Cortex 27, 1849–1862. doi: 10.1093/cercor/bhw024
Shibata, M., Toyomura, A., Itoh, H., and Abe, J. (2010). Neural substrates of irony comprehension: a functional MRI study. Brain Res. 1308, 114–123. doi: 10.1016/j.brainres.2009.10.030
Sled, J. G., Zijdenbos, A. P., and Evans, A. C. (1998). A nonparametric method for automatic correction of intensity nonuniformity in MRI data. IEEE Trans. Med. Imaging 17, 87–97. doi: 10.1109/42.668698
Sparks, B. F., Friedman, S. D., Shaw, D. W., Aylward, E. H., Echelard, D., Artru, A. A., et al. (2002). Brain structural abnormalities in young children with autism spectrum disorder. Neurology 59, 184–192. doi: 10.1212/wnl.59.2.184
Tordjman, S., Somogyi, E., Coulon, N., Kermarrec, S., Cohen, D., Bronsard, G., et al. (2014). Gene × Environment interactions in autism spectrum disorders: role of epigenetic mechanisms. Front. Psychiatry 4:53. doi: 10.3389/fpsyt.2014.00053
Van Rooij, D., Anagnostou, E., Arango, C., Auzias, G., Behrmann, M., Busatto, G. F., et al. (2018). 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 175, 359–369. doi: 10.1176/appi.ajp.2017.17010100
White, T., Andreasen, N. C., and Nopoulos, P. (2002). Brain volumes and surface morphology in monozygotic twins. Cereb. Cortex 12, 486–493. doi: 10.1093/cercor/12.5.486
Wilke, M., Schmithorst, V. J., and Holland, S. K. (2002). Assessment of spatial normalization of whole-brain magnetic resonance images in children. Hum. Brain Mapp. 17, 48–60. doi: 10.1002/hbm.10053
Keywords: autism spectrum disorder (ASD), cortical thickness, pre-scholar child, MRI, neuropsychological, local gyrification index
Citation: Lucibello S, Bertè G, Verdolotti T, Lucignani M, Napolitano A, D’Abronzo R, Cicala MG, Pede E, Chieffo D, Mariotti P, Colosimo C, Mercuri E and Battini R (2022) Cortical Thickness and Clinical Findings in Prescholar Children With Autism Spectrum Disorder. Front. Neurosci. 15:776860. doi: 10.3389/fnins.2021.776860
Received: 14 September 2021; Accepted: 16 December 2021;
Published: 07 February 2022.
Edited by:
Francesco Craig, University of Calabria, ItalyReviewed by:
Uicheul Yoon, Catholic University of Daegu, South KoreaFabrizio Pizzagalli, University of Turin, Italy
Copyright © 2022 Lucibello, Bertè, Verdolotti, Lucignani, Napolitano, D’Abronzo, Cicala, Pede, Chieffo, Mariotti, Colosimo, Mercuri and Battini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roberta Battini, rbattini@fsm.unipi.it