Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder characterized by impairments in social interaction and communication that begin at a young age, which commonly exhibits repetitive behaviors and restricted interests. Currently, ASD occurs in as much as 0.6 % of the global population (Elsabbagh et al. 2012
), and is more common in males (Baron-Cohen et al. 2011
; Baird et al. 2006
). There is increasing evidence that the heterogeneity and complexity of core symptoms in ASD can be explained as a fundamental impairment in both grey and white matter structure (Minshew and Williams 2007
). In this perspective, a recent noninvasive magnetic resonance-based method of combined structure imaging and diffusion tensor imaging (DTI) specifically sensitive to the presence of brain structure and WM integrity alterations, has allowed for a better understanding of the neuroanatomical abnormalities in ASD.
Although the etiology of ASD remains unclear, a growing number of efforts have been made on studies of the neural bases in ASD. There was variability in the literature on brain imaging findings in ASD. Early brain overgrowth was probably the most replicated finding in this population (Bryńska 2012
; Courchesne et al. 2003
). There were some of the available evidence of both grey matter (GM) and WM volume abnormalities in ASD group compare to normal people (Courchesne et al. 2001
; Ulay and Ertuğrul 2009
). Additionally some specific brain regions were particularly implicated, including the frontal lobe, temporal lobe, limbic structure, amygdale, basal ganglia and cerebellar regions (Carper et al. 2002
; Hazlett et al. 2005
; Nordahl et al. 2012
; Rojas et al. 2006
; Schumann et al. 2010
; Ulay and Ertuğrul 2009
). Moreover, few papers based on DTI have shed light on microstructural WM integrity in certain brain regions in ASD (Aoki et al. 2013
; Barnea et al. 2004
).Despite these existing findings, it’s still difficult to compare those reports with each other, depending on substantial differences not only as to the clinical and demographic characteristics of the samples, but also as to the methods of analysis. To date, as far as we know, few studies have combined neuroimaging data from DTI and structural imaging with voxel-based morphometry (VBM) to investigate the neuroanatomical changes in 2–3 years old toddlers with ASD. The goal of current research is to investigate the available brain imaging data and examine their implication for understanding the neuroanatomical bases of ASD. In order to assess these changes, the whole brain volumes, regional anatomy of GM and WM, and WM integrity in 2–3 years old toddlers with ASD and other DD controls were measured using magnetic resonance imaging (MRI) with VBM and DTI. In particular, the purpose of this research is twofold: first, to examine possible brain imaging abnormalities in the structure and WM integrity between ASD children and DD controls; second, to analyze voxel-based brain volume differences between ASD and DD group.
Participants and Assessments
A total of 50 toddlers with ASD (mean age 29.92 ± 5.54 months) and 28 age, gender, and developmental quotient (DQ) matched toddlers with DD (mean age 28.25 ± 4.38 months) were included in a case–control study conducted from May 2010 to November 2011 at the Child Mental Health Research Center of the Nanjing Brain Hospital affiliated of Nanjing Medical University. The study protocol was fully disclosed to all participants and their guardians, and written informed consent was obtained from each participant’s guardian according the provision of the Declaration of Helsinki. The study protocol was approved by the Institutional Review Board of the Nanjing Brain Hospital affiliated of Nanjing Medical University.
The diagnoses of ASD were based on the DSM-IV-TR criteria of pervasive developmental disorders (Rett’s Syndrome and Heller’s Syndrome was excluded) by two licensed child psychiatrists. All patients were assessed by childhood autism rating scale (CARS) (Schopler et al. 1980
) and autism diagnostic inventory-revised (ADI-R) (Lord et al. 1997
). The inclusion criteria of DD was mental retardation or communication disorders on DSM-IV-TR. Mental retardation caused by chromosome abnormalities or gene defect s such as: Down syndrome, Klinefelter’s syndrome, Fragile X syndrome, congenital hypothyroidism, Phenylketonuria and Prader–Willi syndrome was excluded from DD group by genetic tests.
We also excluded patients who had a history of head injury, neurological disorders, or other major medical problems from ASD and DD group. DQ scores of all participants were obtained using the Bayley Scales of Infant Development—Chinese Version (BSID-C) (Shou et al. 1993
Imaging Data Acquisition
Structural MRI images were acquired using a standard quadrature head coil on a 3.0 T Verio MRI system (Siemens Medical Systems, Germany). Before MRI scanning, each subject was sedated by using chloral hydrate with parental consent. All the participants took chloral hydrate (0.5 g in 10 ml) for oral administration around 12 ml each time (50 mg/kg). The maximum dosage was usually <20 ml. About 85 % of the subjects successfully took chloral hydrate for oral administration at first time, and the rest took enema instead at next time. So each participant was sedated before MRI scanning.
During scanning, the patient’s head was gently restrained by foam cushions. High-resolution MRI images for volumetric analysis were obtained with a T1-weighted (T1W) three-dimensional (3D) spoiled gradient (SPGR) sequence with the following parameters: TR = 2,530 ms, TE = 3.34 ms, flip angle = 7°, field of view (FOV) = 256 mm × 256 mm2, in-plane resolution = 256 × 192, inversion time = 1,100 ms, and slice thickness = 1.33 mm. Image orientation parallel to the anterior commissure–posterior commissure (AC–PC) plane was used.
Diffusion tensor imaging (DTI) were obtained with single shot echo planar (SE-EPI) sequences with diffusion gradients applied in 30 non-collinear directions and b = 1,000 s/mm2. The thickness of each slice was 2.5 mm without gap. The sequence parameters for DTI were: TE = 104 ms; repetition time (TR) = 9,000 ms; FOV = 230 × 230 mm2, and acquisition matrix = 128 × 128. Total DTI scanning time was 5.1 min. Routine clinical MRI scans (T1W, T2W, and fluid-attenuated inversion recovery) were conducted to further reveal incidental pathological abnormalities.
Analysis of Structural Images
Customized template was created based on all participants during the optimized VBM8 protocol, which was then used as a prior image to normalize, segment and modulate all images. Briefly, T1W MRI images were bias field-corrected and GM, WM, and cerebrospinal fluid (CSF) were identified. Global GM, WM, CSF, and total intracranial volume (TIV) were then calculated in native space for each subject and used for between-group comparisons. Individual native space GM and WM segments were normalized with an affine registration. From these affine-registered GM and WM segments, an average diffeomorphic anatomical registration through exponentiated lie algebra (DARTEL) template of all subjects in Montreal Neurological Institute (MNI) space was constructed. Next, affine-registered GM segments were warped against the average template using high-dimensional DARTEL means and modulated. Resultant GM maps code the relative local volume of GM in MNI space (Chapman et al. 2006
). To distinguish possible artifacts and failed segmentation or normalization, a quality check was performed using the SPM “check reg” function and MRIcron (http://www.cabiatl.com/mricro/mricron/index.html
). When no unexpected abnormalities were reported, GM maps were smoothed with an 8 × 8 × 8 mm3
Analysis of Diffusion Tensor Images
SPM8, DTI studio 2.4.01 (https://www.mristudio.org/
), and a custom-developed in-house software system based on MatLab 2009 (http://www.mathworks.cn/
) were used for preprocessing and analysis of DTI data. In short, data was denoised, and images were used to calculate FA and mean diffusion (MD) in DTI studio. Because stock templates may result in mistakes in FA and MD difference localization, a custom template for participants was created, and data was spatially normalized into this customized template. Resultant normalized maps were smoothed with 8 mm full-width-at-half-maximum isotropic Gaussian kernels for SPM8 analysis. Group analyses were carried out using random effects voxel-wise multiple regression analyses with a threshold at t
> 3.20 (p
< 0.001, uncorrected) and a spatial extent threshold of k = 50 voxels. These parameters were designed to remove isolated small regions of differences, allowing significantly different regions of interest (ROIs) in FA and MD to be clearly delineated. All coordinates were provided in MNI space.
Between-group differences in age, gender, DQ, CARS score, and ADI-R scores were calculated using t
tests or χ
analyses. Global differences in whole brain structure were assessed by independent-sample t
test in SPSS version 19.0 (IBM, USA). GM and WM volumes as well as total intracranial volume (TIV) were included as test variables, and diagnoses as ASD and DD were the primary grouping variables. Differences in regional GM and WM volumes as well as differences in FA and MD were thresholded at p
value of <0.001 (two-tailed tests) uncorrected for multiple comparisons, with an extent threshold of 50 voxels to eliminate isolated small areas of volume differences using two-sample t
tests in SPM8 (http://www.fil.ion.ucl.ac.uk/spm
In this study, brain structure anatomy and WM integrity in ASD toddlers were investigated following a combined approach of structure imaging and DTI. The neurophysiological results provided three important findings: (1) These young ASD children possessed larger global GM and WM volumes than DD group of the same ages, genders and DQ scores; (2) these ASD children exhibited elevated GM volumes most commonly observed in the right superior temporal gyrus of the temporal lobe and elevated WM volumes most commonly observed in the right superior temporal gyrus, left middle temporal gyrus, right insular cortex, and right Heschl’s gyrus; (3) ASD also exhibited higher FA and lower MD in the CC, posterior cingulate cortex, and limbic lobe as well as lower MD in the insular cortex.
The first result not only agreed with the initial evidence of brain and cerebral enlargement in toddlers with AS reported by Courchesne and his colleagues in 2001, but also extends the stage into Chinese population at which total brain volumes were increased in ASD patients at early age. Courchesne et al. (2001
) indicated that by age 2–3 years over 90 % of autistic toddlers exhibited abnormally larger total brain, and the autistic 2- to 3-year-olds had more cerebral (18 %) and cerebellar (39 %) WM, and more cerebral cortical GM (12 %) than normal boys. Furthermore, similar findings were described in a recent study concerning brain WM volumes in autistic boys of ASD which also suggested that there was an overall increase in brain volumes compared with controls (Herbert et al. 2003
). The current study demonstrated results consistent with a large body of literature supporting that possible brain imaging abnormalities in the structure indeed existed in ASD children at the early age (Courchesne et al. 2001
; Cosgrove et al. 2007
; Hazlett et al. 2005
; Schumann et al. 2010
; Sparks et al. 2002
Secondly, greater regional GM volumes in the right superior temporal gyrus and higher WM volumes in the right superior temporal gyrus, left middle temporal gyrus, right insular cortex, and right Heschl’s gyrus were observed in ASD toddlers in the current study by using optimized VBM. Several previous studies performed on ASD also reported that autistic children exhibited extreme enlargement of GM and WM volumes in the frontal and temporal lobes, a slightly variant result from many other studies (Carper et al. 2002
; Cosgrove et al. 2007
; Courchesne et al. 2004
; Hazlett et al. 2005
; Rojas et al. 2006
; Schumann et al. 2010
). Notably, the structure found in our report were predominantly localized in the temporal regions, which were implicated in social perception, language, and “theory-of-mind”, an ability to put oneself into someone else’s shoes, to imagine their thoughts and feelings, all of which were impaired in ASD. In accordance with above findings, we confirmed the GM and WM volumes increase in the temporal regions of 2- to 3-year-old toddlers with ASD.
Diffusion tensor imaging (DTI) data from the current study also demonstrated that changes in microstructural WM organization with increased FA and reduced MD were apparent in the CC and cingulum of the ASD group. The CC is the largest inter-hemispheric tract in the brain and it is considered to be involved in emotion and social functioning, as well as in higher cognitive processes, such as decoding non-literal meaning, affective prosody, and understanding humour (Badaruddin et al. 2007
; Brown et al. 2005a
). The cingulate cortex participates in tasks related to higher-level cognitive processes, like empathic cognition, social behavior and pain perception in which is more likely to be hypoactivated in ASD (Di Martino et al. 2009
; Thakkar et al. 2008
). Most DTI researches on adolescents or adults with ASD showed reduced FA in several brain regions, including the CC and cingulum (Alexander et al. 2007
; Keller et al. 2007
; Weinstein et al. 2011
However, in a study (Cheng et al. 2010
) carried out on adolescents with ASD, an imbalance of FA was found, as some areas showed an increased FA, while others a decreased FA. Moreover, Billeci et al. (2012
) found an increased FA in about 5 years old ASD in numerous WM tracts like CC and cingulum. These contradictory results could reflect a change of the trend of FA in young children compared with adolescents and adults. Furthermore, it may be linked to the alteration in ASD brain volume, which has a macrostructural expansion in the first years of life and then shifts to an abnormally slow growth (Jones and Cercignani 2010
; Liu et al. 2010
; Muratori et al. 2012
; Vos et al. 2011
; Westerhausen et al. 2011
). In fact, the weak concordance among DTI investigations largely depends on the wide heterogeneity both of ASD clinical phenotype and of methods of analysis, which makes it more difficult for comparison of results. Furthermore, experimental artifacts and confounding factors should be taken into consideration in order to enhance the reliability and comparability of DTI results (Walker et al. 2012
). Caution is needed in interpreting high FA levels, because this change could represent an index of various WM microstructural alterations, like increased myelination, axon size and density, path geometry, and the presence of crossing fiber pathway. To date, as our knowledge, the voxel-based analysis has been applied to limited number of ASD studies on young autistic children, especially in this early stage of 2- to 3-year-old toddlers.
The results of this study should be considered in lieu of the matching strategy, which involved the use of DD group as controls, a group which may not accurately reflect characteristics of the general population. ASD commonly co-occurs with other neurodevelopmental disorders. Many ASD individuals have language deficits, ranging from complete lack of speech through language delays, poor comprehension of speech, or stilted and overly literal language. So it is impossible to exclude the language impairment in our ASD group. Specific learning disabilities are also common in individuals with an ASD. But it is difficult to identify the specific learning disorder during 2–3 years old. Thus, further studies with larger and more diverse cohorts will be required to verify these findings. Furthermore, when comparing the results of this relatively small study with much larger studies without adjustment, it’s impossible to avoid errors. Although customized template was created based on all participants, the use of VBM was still unable to overcome the problem of technical difficulties for segmentation at this early age. Lastly, not including the autism diagnostic observation schedule (ADOS) as part of the diagnostic process should be included as a limitation of our study. Therefore, further substantiation of brain regions, especially for microstructural organization, remains necessary to confirm these findings.