Autism spectrum disorders (ASDs) are complex disorders of neurodevelopment defined by impaired social communication and restricted, repetitive behaviors. ASDs represent a significant public health concern, affecting upward of one in 110 children, with a recurrence rate among at-risk families of nearly one in five (
1,
2). Findings from prospective studies of infant siblings of children with ASDs, who are at higher than average risk for the disorder, indicate that a number of the defining behavioral features of ASDs first emerge around 12 months of age after a period of relatively typical postnatal development (
3–
5).
Although several studies have documented the early behavioral course of infants later diagnosed with ASDs, we know of no published neuroimaging data on such infants before toddlerhood. Existing studies have provided evidence from magnetic resonance imaging (MRI) of significantly larger than normal brain volume in 2- and 3-year olds with ASDs (
6–
10). These findings are consistent with reports of higher than normal brain weight in postmortem brains of individuals with autism (
11), as well as numerous reports associating ASDs with large head circumference (
12). Two large retrospective studies suggest that the onset of this larger head circumference likely occurs in the latter part of the first year of life (
7,
13). Taken together, findings from the behavioral studies and the studies of brain and head size growth in high-risk infant siblings suggest that the latter half of the first year of life is a pivotal time for both brain changes and symptom onset in infants later diagnosed with an ASD. The concurrent timing of these phenomena suggests that brain changes during this period may have an important role in the pathogenesis of autistic behavior.
Autism is increasingly considered a disorder characterized in part by aberrant neural circuitry (
14,
15). Functional neuroimaging studies have revealed patterns of disrupted connectivity in adults (
16) and children (
17,
18) with ASDs. A growing body of work has employed diffusion tensor imaging to gauge the microstructural properties of white matter circuitry (
19,
20). Studies of ASDs using diffusion tensor imaging have identified evidence of widespread abnormalities in white matter fiber tract integrity (
21–
25), with cross-sectional studies indicating less age-related change in children and adolescents with ASDs than in typical comparison subjects (
26,
27). While the preponderance of existing work suggests altered white matter development in ASDs, the extent, direction, and developmental course of these differences remain unclear (
28).
Fundamentally a disorder of development, autism emerges early in life and is generally associated with lifelong disability (
29). When development itself is intrinsic to the phenomenon in question, it is prudent to seek answers from developmental trajectories (
30–
32). This is particularly so in the earliest periods of infancy, when dramatic changes in behavior are paralleled by dramatic changes in the brain. White matter pathways, for instance, rapidly develop during the first years of life, after which change is attenuated (
33). However, to our knowledge, there have been no longitudinal studies of neural circuitry in ASDs, and but a few concerning volumetric brain development (
8,
34,
35). Although numerous cross-sectional studies of children with ASDs posit key neurological changes relevant to the disorder, few have employed a true developmental perspective. As Karmiloff-Smith (
36) aptly observed, cross-sectional neuroimaging studies of children are not tantamount to neuroimaging studies of development.
Method
Participants
This study included data from an Autism Center of Excellence funded by the National Institutes of Health. The Infant Brain Imaging Study (IBIS), the parent network, is an ongoing study of infants at risk for autism. Four clinical data collection sites are associated with the study: University of North Carolina, Chapel Hill; University of Washington, Seattle; Children's Hospital of Philadelphia; and Washington University, St. Louis. Data were coordinated through the Montreal Neurological Institute (MNI) at McGill University, and data processing was performed at the University of North Carolina and the Scientific Computing and Imaging Institute at the University of Utah. The parent study enrolled and assessed 6-month-old high-risk infants who were seen for follow-up assessments at 12 and 24 months of age. Written informed consent was obtained from parents or legal guardians before enrollment, and the study procedures were approved by institutional review boards at each site.
The exclusionary criteria were 1) a significant medical condition known to affect brain development, 2) sensory impairment, 3) low birth weight (<2,200 g) or prematurity (<36 weeks gestation), 4) perinatal brain injury secondary to maternal complications or exposure to specific medications or neurotoxins (e.g., alcohol) during gestation, 5) non-English-speaking family, 6) contraindication for MRI (e.g., metal implant), 7) adoption, and 8) first-degree relative with idiopathic intellectual disability, psychosis, schizophrenia, or bipolar disorder.
The present study group included all high-risk infant siblings who received diffusion-weighted MRI scans at 6 months and for whom behavioral assessments were completed at age 24 months as of June 2011. Symptoms of ASDs were measured at 24 months by using the Autism Diagnostic Observation Schedule (
37), which was completed by research-reliable administrators to maximize agreement across sites. On the basis of classifications from this scale, the high-risk infants were divided into two groups: ASD-negative (below the ASD cutoff) and ASD-positive (above the cutoff). At 24 months, 28 infants met the criteria for ASDs and 64 did not. To characterize differences in autism symptoms between groups, a severity score was generated for each participant. Symptom severity scores on the Autism Diagnostic Observation Schedule range from 1 (least severe) to 10 (most severe), with scores of 4 or higher consistent with the presence of an ASD (
38). The Mullen Scales of Early Learning (
39) were administered at each visit, and the early learning composite score at 24 months of age was used to characterize general cognitive ability. The composite score at 12 months was substituted in the data for four subjects who were missing complete 24-month data for the Mullen scales.
Image Acquisition
MRI brain scans were completed at the clinical sites on identical 3-T Siemens TIM Trio scanners (Siemens Medical Solutions, Malvern, Pa.) equipped with 12-channel head coils during natural sleep. The diffusion tensor imaging sequence was acquired as an ep2d_diff pulse sequence with a field of view of 190 mm (6 and 12 months) or 209 mm (24 months), 75–81 transversal slices, a slice thickness of 2 mm isotropic, 2×2×2-mm
3 voxel resolution, a TR of 12,800–13,300 ms, a TE of 102 ms, variable b values between 0 and 1,000 s/mm
2, 25 gradient directions, and a scan time of 5–6 minutes. Intra- and intersite reliabilities were initially established and regularly evaluated by scanning traveling volunteers, or “phantoms,” at all sites within the same week (
40).
Image Preprocessing
Data from diffusion-weighted imaging were screened by using DTIprep software (
41), which automatically detects artifacts, corrects for motion and eddy current deformations, excludes images with artifacts, and generates a full report. Expert raters manually removed scans with clear residual artifacts. Data sets with fewer than 18 (72%) gradient diffusion-weighted images after this quality procedure were excluded from further processing owing to a low signal-to-noise ratio.
Processing Pipeline
Group analysis of the data from diffusion-weighted imaging, processed by means of diffusion tensor estimates, employed an improved processing pipeline (
42). This processing overcomes the major challenge in implementing tract-oriented statistics in large study groups, which is finding consistent spatial parametrization within and between groups. This includes a computational anatomy approach for nonlinear coregistration of the diffusion tensor imaging data to a template reference coordinate frame, a process to parameterize fiber tracts to functions of length, and the mapping of individual tract geometries into common coordinates.
Computational Anatomy Mapping
Unbiased atlas building (
43) was used to provide one-to-one mapping between the image data and the template atlas, wherein the atlas is built from the population of data as the centered image with the smallest deformation distances. Registration proceeds in two steps (
42). The first applies linear, affine registration of diffusion-weighted imaging baseline images to a structural weighted T
2 atlas by using B-spline registration and normalized mutual information (
44). This is followed by an unbiased, deformable atlas-building procedure (
43) that applies large-deformation diffeomorphic metric mapping transformations. The procedure relates individual data sets to the study-specific atlas template space by means of nonlinear, invertible transformation. Tensor maps were calculated from the diffusion-weighted imaging data sets by using weighted least-squares estimation and were transformed into the atlas space with tensor reorientation by the finite strain approach (
45). The transformed tensor images were averaged by using the Riemannian framework (
46), resulting in a final three-dimensional average tensor atlas for tractography and tract parameterization. Longitudinal diffusion-weighted imaging data for the subjects covering ages 6 to 24 months were mapped into a common atlas space by using the preceding procedure.
Fiber Tractography
Seed label maps were created according to existing tractography methods (
47,
48) and drawn in the combined atlas for regions of interest by using 3D Slicer (
www.slicer.org). A secondary check of regions of interest was made by using an atlas for early childhood (
33). Label maps were created for the following fiber tracts: genu, body, and splenium of the corpus callosum; fornix; inferior longitudinal fasciculus; uncinate fasciculus; anterior thalamic radiation; and anterior and posterior limbs of the internal capsule. Label maps were created bilaterally for all tracts except the corpus callosum. Fiber tracts generated in 3D Slicer were processed for spurious or incomplete streamlines by means of open source software developed in-house (FiberViewer;
http://www.ia.unc.edu/dev/). Fractional anisotropy values were generated for each fiber tract. Fractional anisotropy is an index measuring the degree of anisotropy of local diffusivity, ranging from 0, for isotropic diffusion in fluid, to 1, for strongly directional diffusivity in highly structured axonal bundles (
18,
19). Axial (λ
1) and radial [(λ
2+λ
3)/2] diffusivity values, which represent diffusion parallel and transverse to axonal directions, were also produced.
Statistical Analysis
Demographic characteristics were recorded at 6, 12, and 24 months. Potential group differences between ASD-positive and ASD-negative groups were tested for age, sex, and Mullen Scales of Early Learning composite score at each time point by using either t tests (age and early learning composite score) or Fisher's exact test (sex).
Longitudinal trajectories of mean fractional anisotropy values for fiber tracts in the ASD-positive and -negative groups were compared by using random coefficient linear growth curve models. The random coefficient model fits a group developmental trend while accounting for variability in individual growth trajectories. Among the 92 high-risk subjects included in the analysis, 14 subjects had data from one visit, 40 had data from two visits, and 38 had data from three visits. The mixed-model framework accommodates different patterns of missing data and unbalanced designs.
Although white matter is known to develop more rapidly in the first year of life than in the second, a linear model was determined to best fit the change in fractional anisotropy across time given the limited number of available time points. Separate random coefficient growth curve models were fit for each fiber tract with age, group, and group-by-age interaction as fixed effects and with intercept and age slope as random effects. All growth curve models included the 24-month early learning composite score as a covariate. The primary hypothesis of different age trajectories in the two groups was tested in the interaction of group and age for each tract. With parameters from the random coefficient models, we were able to estimate separate growth slopes for the ASD-positive and ASD-negative groups. Least-squares mean fractional anisotropy values were estimated for each group and contrasted at 6, 12, and 24 months to highlight differences between groups at each individual time point. For the primary analysis of fractional anisotropy trajectories, the false discovery rate was used to correct for multiple comparisons across multiple tracts. Axial and radial diffusivity were examined in a secondary analysis to fully characterize fractional anisotropy results by means of the growth model described in the preceding.
Though the sex ratios in the two groups were not significantly different, to study for potential confounding effects of sex, we ran sensitivity analyses by adding sex and sex-by-age interaction in the preceding models. The p values for the primary hypothesis of age-by-group interactions decreased slightly with no change to the conclusions (results not presented).
All analyses were done by means of SAS software, version 9.2 (SAS Institute, Cary, N.C.).
Results
Descriptive and demographic data are presented in
Table 1. There were no differences between the ASD-negative and ASD-positive infants in age or sex ratio at the 6-, 12-, or 24-month scan (p>0.10) or age at the 24-month behavioral assessment. The final Mullen early learning composite score was approximately 12 points lower in the ASD-positive group than in the ASD-negative group, as those groups were represented at each time point (6 months: t=2.8, df=90, p=0.007; 12 months: t=2.5, df=64, p=0.01; 24 months: t=2.1, df=48, p=0.04). Severity scores on the Autism Diagnostic Observation Schedule at 24 months of age likewise differed significantly between groups as assessed at each time point (6 months: t=19.2, df=90, p<0.001; 12 months: t=18.2, df=64, p<0.001; 24 months: t=14.6, df=48, p<0.001).
Estimates of the fractional anisotropy slope parameters with standard errors are presented for the ASD-positive and -negative groups for all tracts in
Table 2. Both groups showed significant increases in fractional anisotropy from 6 to 24 months, though the rate of change for the ASD-negative group was significantly greater than that for the ASD-positive group in the bilateral limbic (fornix) and association (inferior longitudinal fasciculus and uncinate) fiber tracts. Individual and mean group trajectories for these tracts are presented in
Figure 1. The changes from 6 to 24 months in fractional anisotropy for the corpus callosum subdivisions are shown in
Figure 2; the change for the body was significantly different in the two groups. For projection tracts, the growth trajectories of the left anterior thalamic radiation and all internal capsule divisions were significantly steeper for the ASD-negative infants (
Figure 3).
Follow-up cross-sectional analyses for mean fractional anisotropy are presented in
Table 3. The mean fractional anisotropy value was significantly higher for ASD-positive than for ASD-negative infants at 6 months for the left fornix (t=2.1, df=89, p=0.04), left inferior longitudinal fasciculus (t=2.6, df=89, p=0.01), and left uncinate (t=2.2, df=89, p=0.03). The mean for the body of the corpus callosum was also significantly higher in the ASD-positive infants at 6 months (t=2.0, df=89, p=0.04). The mean value for the right posterior limb of the internal capsule was significantly higher in the ASD-positive infants at 6 months (t=2.1, df=89, p=0.04), and the mean value for the left anterior limb of the internal capsule was significantly higher in the ASD-negative infants at 24 months (t=–2.1, df=89, p=0.04). For the left anterior thalamic radiation, the value was significantly higher in the ASD-negative group at both 12 months (t=–2.0, df=89, p=0.04) and 24 months (t=–3.0, df=89, p=0.003).
Results from secondary analyses of axial and radial diffusivities are presented in
Table 4. Linear growth trajectories did not differ significantly for these indices of tract diffusion, with the exception of the posterior limbs of the internal capsule, for which radial diffusivity was significantly higher across the 6–24-month interval for the ASD-positive group.
Discussion
These preliminary findings suggest that a distinct and pervasive course of white matter fiber tract development characterizes high-risk infants who go on to develop autistic symptoms. Trajectories of the fractional anisotropy values for 12 of the 15 tracts examined in the present study differed significantly between groups. Most fiber tracts for the ASD-positive infants were characterized by higher fractional anisotropy at 6 months followed by blunted developmental trajectories such that fractional anisotropy was lower by 24 months. Around 12 months of age, the groups appeared similar across all tracts except the anterior thalamic radiation; this exception potentially indicates differential timing of development for this specific fiber bundle. Understanding why fractional anisotropy values for multiple fiber tracts are higher at 6 months but show less change over time may be critical to understanding the development of ASDs and shed light on the neural mechanisms underlying their onset.
In agreement with findings from typically developing infants, the fractional anisotropy values for the ASD-negative infants were characterized by rapid change during the 6–24-month interval (
33). The blunted trajectories seen in the ASD-positive group are consistent with the lower fractional anisotropy values and comparative absence of age-related associations observed in older children and adolescents with ASDs (
24,
26,
27), although these patterns may vary by fiber pathway throughout development (
25,
28). The differences between the high-risk groups in this study are particularly striking considering evidence that individuals with ASDs and nonaffected family members of individuals with ASDs share a neural phenotype consisting of specific structural and functional brain abnormalities (
49,
50). In a follow-up analysis, we examined axial and radial diffusivity measures to elucidate the fractional anisotropy results. With a few exceptions, the trajectories for these diffusivity measures did not significantly differ between groups, suggesting that the fractional anisotropy results stem in part from the proportional relationship of axial and radial diffusivity and not from one measure alone.
The altered trajectories of development seen here ostensibly begin in advance of the onset of clinical symptoms (
3–
5), suggesting that core behavioral features of ASDs may arise from an altered neurobiological foundation. Instantiated by these preliminary results, the organization of neural networks underlying ASDs appears to be characterized by atypical patterns of connectivity that differ across systems and time (
14–
17) and are not specific to any single brain region or behavioral domain. The relevance of developmental trajectories to understanding these dynamic processes is undeniable (
30–
32). Had the present study included only cross-sectional data centered at 12 months, for instance, it might have followed that the anterior thalamic radiation is uniquely relevant to the early development of ASDs. The wider view afforded by longitudinal data indicates that thalamic afferents are but one of many tracts implicated in the disorder, with associations at 12 months likely resulting from differential development.
Both highly experience-dependent and less environmentally mediated processes contribute to the functional and structural organization of the brain, and the dynamic interplay of these processes over time yields specialized cortical circuits designed to optimally process complex information (
51). For example, differences in structural organization prior to a period of experience-dependent development related to social cognition (
52–
54) may decrease neural plasticity through limitations on environmental input, preventing typical neural specialization (
52). These alterations could have a ripple effect through decreasing environmental responsiveness and escalating invariance, thus canalizing a specific neural trajectory that results in the behavioral phenotype that defines ASDs. In typical development, the selective refinement of neural connections through axonal pruning (
55) along with constructive processes such as myelination (
56) combine to yield efficient signal transmission among brain regions. One or both of these mechanisms may underlie the widespread differences in white matter fiber pathways observed in the current study. Although diffusion tensor imaging indirectly measures axonal organization, the present findings may reflect connective properties secondary to atypical activity-dependent neuronal modulation and development resulting from disordered molecular mechanisms (
57). It is intriguing that our data are consistent with recent evidence from postmortem stereological research (
58) and mouse models of autism (
59) suggesting that axonal plasticity may be implicated in the development of ASDs.
It is clear that the neurodevelopmental story of ASDs neither begins at 6 months of age nor ends at 24. Extending neuroimaging downward to infants younger than 6 months would help clarify the temporal origin of diverging trajectories, while extending neuroimaging to later ages would capture increasingly stable neurobehavioral outcomes. Additional data points would likewise refine the calibration of trajectories beyond linear models. Because the present results are limited to at-risk infants, the inclusion of a low-risk comparison group is necessary to frame group trajectories in the context of typical development. The addition of a multimodal approach to neuroimaging could sharpen our understanding of early brain changes in ASDs, allowing for the investigation of functional and structural covariance. Linking the structural integrity of specific fiber tracts to behaviors relevant to autism during infancy would similarly clarify mutually dependent developmental processes. To capture the mechanisms that engender and maintain the trajectories evidenced here, future researchers might also consider the role of genetic and epigenetic variables in the development of structural neural circuitry. For instance,
CNTNAP2, which has been associated with ASDs (
60), is necessary to neuron-glia interactions and may mediate axonal development by way of activity-dependent oligodendrocyte functioning (
56).
Finally, the presence of significant differences in fractional anisotropy at 6 months raises the exciting possibility of developing imaging biomarkers for risk of ASDs in advance of symptom onset. Future work might investigate the potential of predictive models for ASDs in early infancy, a process that could include refined imaging techniques or combined biobehavioral markers. Identifying infants at highest risk for ASDs before the full syndrome is manifest offers the possibility of implementing interventions that could reduce or even prevent the manifestation of the full syndrome (
61).
There are several limitations to the present study. The absence of a low-risk contrast group limits interpretation of the results outside of a familial background for ASDs. It is possible, for example, that trajectories for the ASD-negative group reflect developmental resilience to or compensation for neuropathology related to familial risk. Although fractional anisotropy corresponds to organizational properties indicative of development (
19,
33) as well as histological findings (
62), it is an imperfect index of white matter microstructure and may reflect the effects of axonal packing, crossing fibers, or partial volumes related to noise in the data from diffusion tensor imaging (
20,
63). In this study, ASD status was based on Autism Diagnostic Observation Schedule scores at 24 months. Follow-up assessment at 36–48 months of age would provide greater assurance of diagnostic outcomes.
To our knowledge, this longitudinal study has identified the earliest brain differences related to later symptoms of autism to date. The aberrant development of multiple white matter pathways seen here, along with previously reported brain and behavioral change during infancy (
3–
9), suggests a period of critical importance to the pathogenesis of ASDs. A developmental approach to brain and behavioral changes during this time, sensitive to characterizing longitudinal trajectories, is crucial to understanding the complexities inherent to the neurodevelopmental processes implicated in the emergence of ASDs.
Acknowledgments
The Infant Brain Imaging Study (IBIS) Network clinical sites are located at the University of North Carolina (J. Piven, IBIS Network primary investigator; H.C. Hazlett, C. Chappell); the University of Washington (S.R. Dager, A.M. Estes, D. Shaw); Washington University (K.N. Botteron, R.C. McKinstry, J. Constantino, J. Pruett); Children's Hospital of Philadelphia (R.T. Schultz, S.J. Paterson); and the University of Alberta (L. Zwaigenbaum). The data coordinating center is at the Montreal Neurological Institute (A.C. Evans, D.L. Collins, G.B. Pike, P. Kostopoulos, S. Das). The image processing core is at the University of Utah (G. Gerig) and the University of North Carolina (M. Styner). The statistical analysis core is at the University of North Carolina (H. Gu). The genetics analysis core is at the University of North Carolina (P. Sullivan, F. Wright).
The authors thank the IBIS children and their families for their participation; Varun Puvanesarajah and Eric Maltbie for assisting with tractography; and Rachel G. Smith, Cheryl Dietrich, and Mahshid Farzinfar for the DWI/DTI correction and quality control efforts.