Abstract
Understanding the brain differences present at the earliest possible diagnostic age for autism spectrum disorder (ASD) is crucial for delineating the underlying neuropathology of the disorder. However, knowledge of brain structural network changes in the early important developmental period between 2 and 7 years of age is limited in children with ASD. In this study, we aimed to fill the knowledge gap by characterizing age-related brain structural network changes in ASD from 2 to 7 years of age, and identify sensitive network-based imaging biomarkers that are significantly correlated with the symptom severity. Diffusion MRI was acquired in 30 children with ASD and 21 typically developmental (TD) children. With diffusion MRI and quantified clinical assessment, we conducted network-based analysis and correlation between graph-theory-based measurements and symptom severity. Significant age-by-group interaction was found in global network measures and nodal efficiencies during the developmental period of 2–7 years old. Compared with significant age-related growth of the structural network in TD, relatively flattened maturational trends were observed in ASD. Hyper-connectivity in the structural network with higher global efficiency, global network strength, and nodal efficiency were observed in children with ASD. Network edge strength in ASD also demonstrated hyper-connectivity in widespread anatomical connections, including those in default-mode, frontoparietal, and sensorimotor networks. Importantly, identified higher nodal efficiencies and higher network edge strengths were significantly correlated with symptom severity in ASD. Collectively, structural networks in ASD during this early developmental period of 2–7 years of age are characterized by hyper-connectivity and slower maturation, with aberrant hyper-connectivity significantly correlated with symptom severity. These aberrant network measures may serve as imaging biomarkers for ASD from 2 to 7 years of age.
Introduction
Autism spectrum disorder (ASD) is characterized by deficits in social communication and interactions, repetitive patterns of behaviors, and restricted interests (). The age range of 2–7 years is a critical developmental period for understanding ASD as it is a time frame immediately after the earliest possible ASD diagnosis around 2–3 years of age based on current clinical diagnostic methods (). However, nearly all the current literature on brain structure and function in ASD focuses on late childhood periods after 7 years of age (e.g., ; Supekar et al., 2013; ), or adolescent and adult groups (e.g., Rudie et al., 2013; ; ). Recently, a few retrospective neuroimaging studies included very young children before 2–3 years who went on to be diagnosed with ASD (, ; Shiohama et al., 2021). Diffusion MRI studies on the important early developmental period of 2–7 years in ASD are relatively scarce and were conducted to investigate white matter microstructural changes (e.g., ; Walker et al., 2012; ). Since ASD has numerous implications for patient function and integration in society, early interventions for ASD immediately after its diagnosis at around 2 years of age are more likely to reduce symptoms and to positively affect long-term neurodevelopment (Zwaigenbaum et al., 2015). Current lack of understanding about ASD in the age range of 2–7 years limits potential approaches for early intervention.
Connectomics, a framework to comprehensively map brain organization into a network (), has been applied to studying a variety of brain disorders (van den Heuvel and Sporns, 2019). Topological properties of complex brain networks can be quantified with graph theoretical analysis. The network is represented as nodes of a graph with connecting edges based on measures of connectivity. A growing body of literature suggests that ASD is associated with alterations in multiple interconnected brain systems rather than isolated regions (e.g., ; Rudie et al., 2013; Uddin et al., 2013; ). Meanwhile, large-scale brain systems affected in ASD may underlie patients’ complex phenotype of behavioral impairments. For instance, prior studies reported that the aberrant brain connectivity of the default mode network (DMN) is likely linked to social deficits in ASD (e.g., ; ). Using connectome analysis to systematically delineate alterations in brain organization of ASD may contribute to understanding its underlying neuropathology and identifying sensitive network-based imaging biomarkers that are correlated with its symptom severity.
Heterogenous patterns of aberrant brain connectivity, characterized by both hyper- and hypo-connectivity, have been reported in prior studies on ASD (e.g., ; Rudie et al., 2013; Supekar et al., 2013; Uddin et al., 2013; ; ; ). Specifically, brain hyper-connectivity is more common in children with ASD, while hypo-connectivity is prevalent in adolescents and adults with ASD compared to typically developing (TD) individuals (; Supekar et al., 2013; Uddin et al., 2013). The discrepancies between findings of hyper- and hypo-connectivity are likely dependent on altered age-related trajectories associated with ASD in different developmental stages (e.g., Uddin et al., 2013; ). With connectomic analysis as an effective tool to probe brain organizational changes, characterizing age-related brain network and connectivity changes in ASD in different developmental stages may elucidate the non-uniform changes associated with ASD. Structural networks, usually constructed by white matter (WM) bundles traced with diffusion MRI (dMRI) tractography, serve as the structural substrate of brain functional organization. So far, the majority of the network and connectivity research in ASD has adopted functional magnetic resonance imaging (fMRI) approaches (e.g., ; Supekar et al., 2013; Uddin et al., 2013; ; ). Relatively few studies have investigated structural network changes in children with ASD (e.g., Rudie et al., 2013; ). To our knowledge, none of the structural network studies were focused on early childhood of 2–7 years. Understanding of brain structural networks in ASD in this early developmental period of 2–7 years may shed light on initial structural connectional configurations that can evolve into both hyper- and hypo-connectivity organizations across different developmental stages and provide a structural basis for functional network alterations.
In this study, we aimed to characterize age-related brain structural network changes in young children with ASD from 2 to 7 years of age, and to identify sensitive network-based imaging biomarkers that are significantly correlated with the symptom severity during this early developmental period. Based on our previous findings of flattened WM microstructural changes quantified by dMRI microstructural measures in children with ASD from the same age range (), we hypothesized aberrant hyper-connectivity and slower maturation of structural networks in children with ASD during this same age period. Diffusion MRI was acquired to map the brain structural connectome in 30 children with ASD and 21 TD children. All brain structural networks were constructed with dMRI tractography. We compared graph-theory-based network measurements between TD and ASD individuals and delineated age-related characteristics and alterations of these network properties. We also conducted correlations between graph-theory-based network measurements and symptom severity in children with ASD.
Materials and Methods
Children With Autism Spectrum Disorder and Children With Typical Development
All participants were children recruited at Beijing’s Children’s Hospital. The study was approved by the institute Research Ethics Committee, and informed/parental consent was obtained. Thirty male children with ASD aged 2.33–7.00 years (4.15 ± 1.42 years) and 21 male TD children aged 1.99–5.96 years (3.90 ± 1.11 years) participated in this study. TD children at the time of MR imaging were referred for seizures with fever (n = 9), intermittent headache (n = 10), and strabismus (n = 2). All TD children had normal neurological examinations documented in their medical records. The exclusion criteria for TD consisted of known nervous system disease or a history of psychiatric, neurodevelopmental, or systemic illness. Children with ASD were not receiving any CNS-active medications before the MRI studies. The diagnosis of ASD was established using the Autism Diagnostic Interview-Revised (ADI-R) (), Childhood Autism Rating Scale (CARS), Clancy Autism Behavior Scale (CABS) (), and Autism Behavior Checklist (ABC) (), and confirmed based on expert opinion according to Diagnostic and Statistical Manual of Mental Disorder criteria (). All diagnoses were performed by two experienced pediatricians with 11 and 12 years of experience in clinical neuropsychology, respectively. These clinical assessments were not performed in TD. Detailed demographics and clinical characteristics of participants are provided in Table 1. Thirty children with ASD have been reported in a prior article () focused on atypical WM microstructure in ASD, whereas in this study we examined completely different brain properties (i.e., the brain macrostructural network) using completely different analysis (i.e., graph-theory analysis and network-based statistics) in a mostly overlapped ASD cohort and a larger TD cohort.
TABLE 1
| Parameter | Children with ASD (n = 30) | Children with TD (n = 21) |
| Age (years) mean ± SD | 4.15 ± 1.42 | 3.90 ± 1.11 |
| Median (min-max) | 3.50 (2.33–7.00) | 3.84 (1.99–5.96) |
| Gender (male/female) | 30/0 | 21/0 |
| Clinical assessment score | ||
| Autism behavior checklist (ABC) | ||
| Total score | 94.39 ± 7.41 | – |
| Sensory | 8.03 ± 3.03 | – |
| Relating | 30.97 ± 4.56 | – |
| Stereotypes and object use | 9.81 ± 4.23 | – |
| Language | 27.29 ± 2.51 | – |
| Self-help and social | 18.23 ± 2.90 | – |
| Autism diagnostic interview (ADI-R) | 52.77 ± 6.92 | – |
| Childhood autism rating scale (CARS) | 41.03 ± 3.79 | – |
| Clancy autism behavior scale (CABS) | 18.03 ± 1.96 | – |
Demographics and clinical assessments of participants.
Data are presented as mean ± standard deviation. All of the scores are raw values. ASD, autism spectrum disorder. TD, typically developing.
MRI Data Acquisition
All children were scanned using a 3T Philips Achieva MR System with sedation. Axial diffusion MRI data was acquired with the anterior-posterior commissure (AC-PC) line parallel to the phase-encoding direction. Single-shot, echo-planar imaging (EPI) sequence was used with Sensitivity Encoding parallel imaging scheme (SENSE, reduction factor = 2.5). Eight-channel SENSE head coil was used. Other imaging parameters were as follows: repetition time (TR) = 7.96 s, echo time (TE) = 83 ms, field of view (FOV) = 256 × 256 mm2, imaging matrix = 128 × 128, voxel size = 2 × 2 × 2 mm3, slice number = 70 covering the entire brain without a slice gap. Diffusion weighting was encoded along 30 independent directions, and b-value was 1,000 s/mm2. The acquisition was repeated twice to improve signal-to-noise ratio (SNR), resulting in a scan time of 11.5 min. T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) image was also acquired. Imaging parameters for MPRAGE were as follows: TR/TE = 8.34/3.83 ms, flip angle = 12°, FOV = 256 × 256 × 160 mm3, voxel size = 1 × 1 × 1 mm3. T1-weighted images (T1w) have superior gray and white matter contrast for cortical parcellation. Visual inspection was conducted for all MRI data by the pediatric radiologists (Y.P. and H.C.), and no apparent motion artifacts were found.
Data Preprocessing
Diffusion MRI data preprocessing including eddy current and motion correction, tensor fitting, and estimation of diffusion MRI-derived measures was conducted with DTIStudio (). Specifically, eddy current distortion and head motion were corrected by registering all raw diffusion weighted images (DWIs) to a b0 image using a 12-parameter (affine) automated image registration (AIR) algorithm (Woods et al., 1998; ). Few motion artifacts were observed in the dMRI datasets. Head motions in dMRI data were quantified for all subjects using the methods described in the literature () and in Supplementary Material. As shown in Supplementary Figure 1, few motion artifacts were observed in both TD and ASD subjects with DWI volume-by-volume translation less than 0.5 mm and rotation around 0.12 degrees for all subjects. There were no significant differences in translation (p = 0.719) and rotation (p = 0.227) between ASD and TD groups. There were no significant correlations between age and translation (p = 0.672) or rotation (p = 0.373) in this cohort, either. Standard tensor fitting was conducted to generate dMRI-derived metrics, including fractional anisotropy (FA).
Network Construction
The two fundamental elements of a network, nodes and edges, were defined using the following procedures to construct an individual structural network.
Network Node Definition
Network nodes of each subject in the native dMRI space were obtained by transferring their parcellated cortical regions from their native T1w space (Figures 1A–F). Briefly, contrasts of individual subject’s T1w image (Figure 1D) and b0 image (Figure 1A) in their native dMRI space were used to drive linear registration with a transformation matrix T. T1w image of each subject (Figure 1D) was parcellated into 68 cortical gyri (Figure 1E and Supplementary Table 1) using Freesurfer software () with Desikan–Kiliany atlas (). Visual inspection was conducted for Freesurfer cortical parcellation by the pediatric radiologists (Y.P. and H.C.) to ensure the parcellation quality as well as to make sure there were no clearly misregistered gyral labels across all subjects in automatic cortical parcellation. The inverse transformation (T–1) was used to warp the individual’s cortical parcellations (Figure 1E) into individual native dMRI space. The 68 cortical regions of interest (ROIs) with 34 ROIs in each hemisphere (see Supplementary Table 1 for detail) represented 68 nodes of subject’s structural network (Figure 1). Of note, parcellated cortical ribbon in the dMRI space was dilated by 8 mm toward the deep WM direction (Figure 1F), using an in-house code to penetrate superficial WM and reach deep WM to initiate fiber tracking (). All registrations were performed using the SPM8 package1.
FIGURE 1
Network Edge Definition
Network edge was defined as the number of fiber streamlines connecting two regions (Figure 1F). A brain structural network (Figure 1H) was constructed for each participant, represented by a symmetric 68 × 68 connectivity matrix (Figure 1G). Network edges were defined with reconstructed whole brain white matter fibers. Whole-brain tractography was performed within a brain mask derived from the subject’s b0 image, using diffusion MRI Brute-force deterministic fiber tractography (; ) in Diffusion Toolkit.2 In this approach, a streamline was initiated from each voxel of the brain mask. Due to relatively low FA in brains of young children, FA threshold was set to 0.15 and angle threshold was 40° for tractography (e.g., Zhao et al., 2019). Only reconstructed fibers with two end points located in the pair of dilated cortical regions, network nodes (Figure 1F), were kept to define network edges (Figure 1C). Number of fiber streamlines (FN) connecting two regions was defined as edge weight (wij). To remove spurious connections, we used a minimum threshold of fiber streamlines number (wij = 5) between two regions. As a result, we constructed a weighted structural network (Figure 1H) for each subject, represented by a symmetric 68 × 68 connectivity matrix (Figure 1G).
Network Analysis
A brain structural network graph G (Figure 1H) is composed by N nodes and K edges. Both global and regional network measures were calculated to provide a summarized scalar for topological characterization of individual brain structural network. Both global and regional network measures were quantified with graph-theory approaches. These network measures characterize the topological organization of brain structural network. Specifically, the following global graph measures were calculated: network strength, global and local efficiency (Eg and Eloc), and shortest path length (Lp) (). For regional properties, we calculated nodal efficiency of each node (). Based on nodal efficiency, we identified network hubs in both groups, as a node with efficiency that was at least 1 standard deviation larger than averaged nodal efficiency across all nodes. All network analysis was performed using GRETNA software (Wang et al., 2015), and results were visualized using BrainNet Viewer software (Xia et al., 2013). Detailed definitions of these network properties, including network strength, global and local efficiency (Eg and Eloc), shortest path length (Lp) (), nodal efficiency Enodal, and hub, are described in Supplementary Material.
Statistical Analysis
Network-Based Statistic
To localize specific altered structural network edges strength in ASD, we used a network-based statistic (NBS)3 approach (Zalesky et al., 2012) with age as covariate to identify between-group differences in pairwise edge (or connection). The NBS analysis was performed in three steps. First, a threshold of p < 0.05 (before correction) was used to yield t statistic (two-sample t-tests) matrix of suprathreshold connections, among which any connected components and their size (the number of connections) were determined. Second, a nonparametric permutation approach (10,000 permutations) was used to estimate statistical significance of observed component sizes in the un-corrected connection matrix, controlling family wise error. Briefly, in each permutation, all participants were randomly shuffled into two groups, and two-sample t-tests were recomputed to examine group differences in network edges after controlling for age. The same primary threshold (p < 0.05) was used to produce suprathreshold connections among which the size of the maximal connected component was recorded. This permutation approach derived the empirical null distribution of connected component size for estimating the significance of observed component sizes. Finally, for a connected component of size N found in the real grouping of control and ASD, its corrected p-value was determined by finding the proportion of the 10,000 permutations for which the maximal connected component was larger than N. Interconnected subnetwork components with a corrected p < 0.05 were considered statistically significant.
Between-Group Differences
To examine between-group difference in global and nodal network properties, a general linear model (GLM) was performed with age as covariate.
Age-by-Group Interaction
To assess age-related alteration in global network properties and nodal efficiency, the age-by-group interaction term was added into GLM as the main effect with age and group as covariates. If interaction effect was significant, rates of age-dependent trendline in network measures would be significantly different between groups.
Correlation With Clinical Assessments
For the network nodes or edges with significant group differences, relationships between network properties and clinical assessments in ASD were explored using GLM with network measures as dependent variables and clinical assessments (Table 1) as independent variables. Subject age was treated as confounding covariate in the GLM models.
All GLM analyses were performed in R statistic-software (version 3.5.1).4 For global network measures, a p < 0.05 was considered significant. For regional network metrics, multiple comparisons were corrected with false discovery rate correction within each hemisphere (p < 0.05). To capture more subtle differences in group comparisons of regional network metrics, we used a less strict false-discovery rate (FDR) correction for multiple comparisons within a hemisphere instead of the whole brain.
Results
The demographic and clinical assessment data for children with ASD and TD are shown in Table 1. There was no significant difference in median age between ASD and TD (Wilcoxon rank-sum test, p = 0.8). Assessed clinical scores of 30 children with ASD ranged from 79 to 107 (mean ± sd: 94.39 ± 7.41) for total score of ABC scale, 41–67 (52.77 ± 6.92) for ADI-R, 36–50 (41.03 ± 3.79) for CARS, and 14–22 (18.03 ± 1.96) for CABS (Table 1).
Flattened Age-Related Characteristics of Structural Network in Autism Spectrum Disorder at Whole-Brain, Global Level
Both TD and ASD groups showed significantly age-related increases in global efficiency (TD: r = 0.75, p < 0.0001; ASD: r = 0.39, p = 0.03), local efficiency (TD: r = 0.73, p = 0.0002; ASD: r = 0.40, p = 0.03), network strength (TD: r = 0.75, p < 0.0001; ASD: r = 0.44, p = 0.01), and decrease in network shortest path length (TD: r = −0.73, p = 0.0002; ASD: r = −0.36, p = 0.05) from 2 to 7 years old (Figure 2A). However, the growth rates of global network measures were significantly slower in ASD compared with TD, revealed by age-by-group interaction analyses (Figure 2A, global efficiency p = 0.006, network strength p = 0.009, local efficiency p = 0.01, and shortest path length p = 0.0005). Overall, from 2 to 7 years of age, children with ASD exhibited a significant hyper-connectivity pattern across multiple global network measures, including increased network strength (t = 2.597, p = 0.012), increased global efficiency (t = 3.161, p = 0.003), and decreased shortest path length (Lp, t = −3.34, p = 0.002) (Figure 2B). We also investigated motion effects on age-by-group interaction and group comparison results at the global level and found that the above findings were not affected after statistically controlling for motion estimates in dMRI scans in the statistical analyses, as shown in Supplementary Tables 2, 3.
FIGURE 2
Flattened Age-Related Trends of Nodal Efficiency in Structural Network of Autism Spectrum Disorder at Regional Level
Network hub regions, serving pivotal roles for communication between any pairs of network nodes, were identified for both ASD and TD groups (Figure 3A, left and central panels). Similar hub distributions, with core regions mainly in the frontal and parietal cortices, were found for both groups. Specifically, the bilateral caudal middle frontal gyrus (cMFG), Isthmus cingulate cortex (IsC), precuneus cortex (PCUN), superior frontal gyrus (SFG), and superior parietal cortex (SPC) are common hubs for both groups. The left superior temporal gyrus (STG) was identified as hub only in TD children.
FIGURE 3
Similar to global network properties, significant age-by-group interactions were found in nodal efficiency across seven brain regions (Figure 3A, right panel, all p < 0.006), including bilateral SPC, left IsC, left SFG, left PCAL, right supramarginal gyrus (SMG), and right temporal pole (Tpole). Four regions were network hubs (i.e., bilateral SPC, left IsC, and left SFG) and displayed as yellow nodes in Figure 3A (right panel). Scatter plots clearly show initial high nodal efficiencies in children with ASD, with subsequent flattened and slower development. Nodal efficiencies of structural networks in TD children significantly increased with age from 2 to 7 years old (Figure 3B, all p < 0.008). However, relatively flattened age-dependent trendlines of nodal efficiency were observed in ASD, indicating an atypical slower maturation of the structural network based on nodal efficiency. The correlation r and p-values from linear regression between age and nodal efficiency for both groups are provided in Figure 3B. Similarly, Supplementary Table 3 demonstrates that motion did not contaminate the significant age-by-group interaction in nodal efficiency found in this study.
Hyper-Connectivity of Structural Network of Autism Spectrum Disorder at Regional Level
Children with ASD also exhibit a hyper-connective pattern in structural networks at the regional level. This hyper-connective pattern is reflected by higher nodal efficiencies in ASD and involves multiple brain subnetworks across several network hub regions (Figure 4A). Significantly higher nodal efficiencies were found in 7 hub regions (i.e., bilateral cMFG, bilateral SPC, bilateral PCUN, and right IsC) and 5 non-hub regions [i.e., left caudal anterior cingulate cortex (cACC), left posterior cingulate cortex (PCC), left parahippocampal gyrus (PHG), left Tpole, and right precentral gyrus (PrCG)] (Figure 4A and Table 2, corrected p < 0.05, and Supplementary Table 2 after motion adjustment). Notably, among these hyper-connective regions, bilateral cMFG, bilateral PCUN, left PCC, and right IsC are key regions of functional default-mode network. Left cACC and right PrCG are critical regions of frontoparietal and sensorimotor networks, respectively.
FIGURE 4
TABLE 2
| Regions | Category | Enodal (mean ± SD) | t-value | p-value (FDR corrected) | |
| Children with TD | Children with ASD | ||||
| L.cMFG | Hub | 55.81 ± 9.93 | 65.84 ± 8.58 | −3.838 | 0.0003 |
| L.cACC | Non-hub | 39.85 ± 8.02 | 47.51 ± 6.68 | −3.710 | 0.0005 |
| L.PCC | Non-hub | 49.90 ± 8.37 | 57.98 ± 8.00 | −3.533 | 0.0009 |
| L.SPC | Hub | 77.90 ± 15.89 | 90.28 ± 10.56 | −3.400 | 0.0013 |
| L.PCUN | Hub | 57.20 ± 9.57 | 66.10 ± 9.27 | −3.264 | 0.0020 |
| L.Tpole | Non-hub | 46.68 ± 8.91 | 54.98 ± 8.85 | −3.196 | 0.0025 |
| L.PHG | Non-hub | 41.61 ± 6.27 | 47.01 ± 6.38 | −2.949 | 0.0049 |
| R.PrCG | Non-hub | 37.29 ± 7.70 | 45.14 ± 6.04 | −3.957 | 0.0002 |
| R.SPC | Hub | 76.26 ± 15.84 | 89.44 ± 10.48 | −3.651 | 0.0006 |
| R.PCUN | Hub | 56.49 ± 10.00 | 65.54 ± 7.88 | −3.505 | 0.0010 |
| R.IsC | Hub | 51.92 ± 9.54 | 59.65 ± 8.44 | −2.938 | 0.0051 |
| R.cMFG | Hub | 53.75 ± 10.72 | 61.19 ± 7.15 | −2.898 | 0.0056 |
Brain regions with significant group difference between TD and ASD in nodal efficiency.
The regions with significant increased nodal efficiency (p < 0.05, corrected) are listed in ascending order by absolute t scores in each hemisphere. L, left hemisphere; R, right hemisphere. Cortical regions were classified into hub and non-hub regions. FDR, false discovery rate. See Supplementary Table 1 for abbreviations of brain regions.
Significantly higher strength of network edges between widespread brain regions also indicates hyper-connectivity in structural network in ASD (Figure 5). When compared with TD, children with ASD showed significantly increased connection strengths in 61 network edges connecting 51 nodes, widely distributed across the whole brain (Figure 5A, p < 0.05 corrected). Higher edge strengths were found in many connections between network hubs demonstrated as yellow nodes in the circle view of Figure 5A, including but not limited to the connections between left cMFG and left PCUN (L.cMFG-L.PCUN), between right IsC and right PCUN (R.IsC-R.PCUN), and between bilateral SPC. Supplementary Figure 2 demonstrates that motion did not contaminate the significantly higher strength of network edges between widespread brain regions in ASD found in this study.
FIGURE 5
Hyper-Connectivity in Autism Spectrum Disorder Correlated With Symptom Severity
We further examined the relationship between widely observed hyper-connectivity of structural networks in ASD and symptom severity. Among the 12 nodes with significantly higher nodal efficiency in ASD, two were positively correlated with symptom severity (Figure 4B, p < 0.05). Specifically, nodal efficiencies from right IsC of functional default-mode network (r = 0.47, p = 0.008) and right PrCG of sensorimotor network (r = 0.44, p = 0.015) in ASD were significantly correlated with the ABC and CABS scores, respectively. In addition, five connections with higher edge strength in ASD were positively correlated with their symptom severity (Figure 5B, p < 0.05). Particularly, L.cMFG-L.PCUN (r = 0.39, p = 0.03) and R.IsC-R.PCUN (r = 0.57, p = 0.001), connections within functional default-mode network, significantly correlated with ABC and CARS scores, respectively. Connection between left and right rostral middle frontal gyrus (L.rMFG-R.rMFG), as part of frontoparietal network, was significantly correlated with CARS score (r = 0.39, p = 0.03). Other sensorimotor network related connections between right PrCG and right SPC (R.PrCG-R.SPC; r = 0.40, p = 0.03), and between left banks superior temporal sulcus (BSTS) and the left Tpole (L.BSTS-L.Tpole; r = 0.54, p = 0.002) were significantly correlated with ABC and CABS scores, respectively. These findings between clinical scores and nodal efficiency or edge strength demonstrate the significant association between hyper-connectivity from aberrant structural network of ASD and severity of symptoms.
Discussion
The present study found flattened age-related structural network changes associated with ASD in children from 2 to 7 years of age and revealed a significant correlation between structural network alterations and symptom severity. It sheds light on the brain structural organizational alterations in children with ASD in this critical and rarely studied developmental age of 2–7 years, immediately after the earliest possible ASD diagnosis based on current clinical diagnostic methods. Structural hyper-connectivity in the pivotal default-mode, frontoparietal, and sensorimotor networks underlies the flattened and slower maturation of these structural networks in ASD. Consistent with functions of the default-mode, frontoparietal, and sensorimotor networks, hyper-connectivity in these networks was significantly correlated with quantified symptom severity based on behavioral assessment, complementing the existing structure-function-behavior framework for ASD in early childhood. These network-based hyperconnectivity in ASD over the studied age range may serve as biomarkers potentially predicting the disorder outcomes and setting the stage for possible early intervention.
Flattened Age-Related Structural Network Changes and Hyper-Connectivity in Children With Autism Spectrum Disorder
Non-uniform brain connectivity findings across different age-ranges support the need to delineate age-related characteristics of structural networks in ASD. The pattern of widespread hyper-connectivity in networks, with relatively slow age-related progression, corroborates the previous microstructural observation of higher WM integrity at an earlier stage and subsequent slower WM maturation in ASD (). During typical development from infancy to childhood, WM maturation is characterized by increased WM FA (; Yu et al., 2020) and structural network reconfiguration toward stronger and more efficient connectivity (). This pattern of structural development facilitates neural information integration and transformation across brain regions (). Structural hyper-connectivity found in 2–7-year-old children with ASD is consistent with a previous structural network study of ASD (). Enhanced WM microstructural properties with higher FA values across WM bundles in young children with ASD (e.g., ; ) may contribute to the observed hyper-connectivity in their structural networks as demonstrated by stronger network edges and higher network efficiencies. The pattern of local or short-range hyperconnectivity has been frequently suggested in the brains of individuals with ASD (e.g., ; ), especially in younger age groups (e.g., Rudie and Dapretto, 2013). Such pattern was even observed in neonates with high risk of ASD (e.g., ). However, the progression from hyper-connectivity observed in young childhood ASD to the hypo-connectivity observed in other studies from adolescents and adults with ASD implies a pattern of overall connectivity decrease with development, relative to TD children (e.g., Rudie et al., 2013; ; ).
Importantly, multiple network hubs exhibited significant hyper-connectivity and atypical network property changes, further indicating the crucial role of network hubs in information transfer and their relative vulnerability in brain disorders (). Given the evolution from hyper- to hypo-connectivity, it is plausible that these atypical network hubs revert into non-hubs with age. As such, the connectivity state of these hubs may be helpful in anticipating the progression of ASD pathology and symptoms. For instance, the left STG, or primary auditory cortex, was a network hub for TD children, whereas children with ASD did not have a network hub in this location. This finding is supported by prior magnetic encephalography (MEG) studies, which also revealed atypical development of the primary auditory cortex in children with ASD, as evidenced by a characteristic delay in their auditory processing (). Aberrant age-related trends of topological network properties, indicating an important developmental discontinuity, may serve as a neural signature in young childhood ASD.
The human brain network is reconfigured toward an optimal global balance between information segregation and integration during development (; ), resulting in the reorganization of the brain network architecture from a relatively randomized configuration to a well-organized one. When examining global network measures of structural networks, we found that individuals with ASD had higher levels of global efficiency and lower levels of shortest path length. Global efficiency reflects the ability to integrate distributed information between distant brain regions and is highly related to long-range connections. Enhanced structural integrity of long-range WM bundles found in ASD (e.g., ) may contribute to making the brain networks more globally integrated (). Since a random network usually has a short characteristic path length (Sporns, 2011), the altered global network measures in ASD may suggest a less organized or more random distribution of network connections. This is consistent with previous network findings in children with ASD (e.g., Rudie et al., 2013) as well as the increased randomness of brain oscillations findings in resting-state fMRI of individuals with ASD (). Furthermore, structural networks displayed similar levels of local efficiency across both groups in our study. The local efficiency indicates the network segregation capacity, which facilitates functional specialization. Our findings suggest that the network integration and segregation in children with ASD are not appropriately balanced during development.
Correlation of Hyper-Connectivity in Autism Spectrum Disorder With Symptom Severity
We established the relationship between severe ASD symptoms and greater structural connectivity. This is consistent with prior functional network research that discovered associations between aberrant functional hyper-connectivity and symptom severity in children with ASD (; ; Supekar et al., 2013). Notably, structural hyper-connectivity, which was significantly associated with symptom severity in the present study, was found in regions or connections in functional default-mode, frontoparietal, and sensorimotor networks. Specifically, IsC, cMFG, and PCUN in Figures 4B, 5B are critical hubs of the default-mode network. The default-mode network is considered to play a critical role in self-referential and high-order social cognitive processes, such as Theory of Mind (). Impairments in this type of mental process are part of the social-cognitive symptoms that characterize ASD (). Existing research indicates that disruptions in the default-mode network might significantly contribute to social deficits in ASD across the age span (e.g., ; ; ; Uddin et al., 2013; Yerys et al., 2015; ). For instance, functional hyper-connectivity in IsC and PCC predicted social communication deficits in children with ASD aged 7–12 years (), and abnormal PCUN and cMFG functional connectivity significantly correlated with the severity of social and communication deficits in adolescents with ASD aged 11–20 years (). Abnormal gray matter volumes potentially related to their symptom severity (Uddin et al., 2011) in default-mode network areas were also reported in children and adolescents with ASD. The association in Figure 5B between symptom severity and the banks of superior temporal sulcus (STS) connection might be attributed to social deficits in ASD as well. STS, a multimodal association region involved in cortical integration of both sensory and limbic information, has been recognized as a key cortical area of the social brain (). Anatomical and functional abnormalities in the STS are highly implicated in ASD with social-interaction impairment (Zilbovicius et al., 2006). The observed correlation between structural hyper-connectivity in the frontoparietal network with symptom severity in Figure 5B may be related to executive dysfunction in the ASD population. The frontoparietal network, also known as the central executive network, is involved in highly adaptive cognitive control processes and is critical to the completion of executive functions (). A meta-analysis of fMRI studies on executive functions in ASD aged 7–52 years revealed that abnormal functional connectivity in the frontoparietal network may underline some executive dysfunctions such as rigidity, preservation, and repeated behaviors, which are commonly seen in ASD (). Lastly, identified regions and connections from the sensorimotor network in Figures 4B, 5B are likely due to the motor abnormality in ASD. Previous studies have demonstrated that motor impairments linked with ASD can be observed as early as infancy () and are prevalent in children with ASD (). Disruption of functional organization within PrCG, a key component of sensorimotor network, was related to ASD diagnosis and to the severity of ASD traits in children (). Taken together, our findings further elucidate the structural connectional basis for these established functional alterations and revealed the relationship between altered structural connectivity and the phenotype of behavioral impairments in ASD.
Limitations and Future Directions
We are aware of a few limitations for our study. First, dMRI-based deterministic tractography was adopted to construct network edges. There is potential loss of tracing existing fibers in areas with fiber crossings. We have considered different tractography algorithms in the construction of structural networks in the aspects of connectome specificity and sensitivity. Although deterministic methods yield sparse connectomes with false negatives, more sophisticated probabilistic methods (e.g., ) yield dense connectomes with low specificity due to false positives (Zalesky et al., 2016). Given that only single shell dMRI data with a b-value of 1,000 s/mm2 along 30 directions was acquired, we adopted deterministic tractography in the current study. Since the same deterministic tractography was applied to all subjects, false negatives in connectomes across subjects were offset as relative connectome metric changes instead of absolute connectome measures were focused on in this study. Second, Freesurfer parcellation of individual T1w images was used to define network nodes. Topological organization of brain networks could be affected by different parcellation strategies when defining network nodes (Zalesky et al., 2010). Despite that choice of network nodes varying across studies, Freesurfer parcellation has been widely used in previous brain network studies (e.g., ) and is able to reliably define brain anatomic regions. Future studies with high-spatial-resolution parcellation strategies in network analysis may further evaluate the reproducibility of our findings. Third, we recognize that atypical development of structural networks in ASD comes from cross-sectional datasets. Future longitudinal studies are warranted to confirm these findings and to fully characterize developmental trajectories of structural network in ASD. Fourth, we only examined the cortico-cortical structural connectivity but not the cortical-subcortical connectivity, such as thalamocortical connectivity. Given the critical role of the thalamus in information processing and cortical functioning of the brain (Sherman, 2016), more research focusing on systematically examining thalamocortical connections is necessary for a more comprehensive understanding of brain connectivity patterns in young children with ASD (e.g., ). In addition, considering ASD is more prevalent in males within a general population (), current study only involved male subjects. Although statistical significance was detected in this study, we acknowledge that larger sample sizes with more homogeneous age distribution and balanced male-to-female ratio are needed for better characterizing such a heterogeneous neurodevelopmental disorder. This work focused on young children aged 2–7 years can also be extended in the future to infants at risk for ASD under 2 years of age for identifying sensitive imaging markers that can predict ASD diagnosis and allowing early intervention at a time prior to the age of diagnosis (; ). Specifically, changes in the brain have been observed in infants at risk for ASD under 2 years of age (e.g., Wolff et al., 2012; ), preceding clinical manifestations. The changes in brain were also reported for infants under 2 years of age later developing ASD and identified through retrospective study (e.g., , ; Shiohama et al., 2021). A machine learning model that leverages sensitive imaging biomarkers may be capable of predicting the clinical diagnostic outcome of individual infants before the development of the full syndrome. The prediction makes early detection and intervention possible and has a significant potential impact on improving outcomes.
Conclusion
In conclusion, our results suggest that hyper-connectivity with slower subsequent maturation of structural networks is a key component of the underlying neurobiology of young childhood ASD. These networks with hyper-connectivity are significantly correlated with ASD symptom severity in the 2–7-year-old range. Measurements of structural connectivity as indicated by network efficiency and edge strength may serve as early predictors of abnormal developmental trajectories, thereby designating patients with a need for more targeted early intervention.
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.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. The MRI data is available from the website http://brainmrimap.org, a public website maintained by Huang lab, or upon reasonable request to the corresponding author.
Ethics statement
This study involving human participants was reviewed and approved by the Beijing Children’s Hospital, Capital Medical University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
HH and YP designed the study. MO, YP, SS, DH, TZ, HC, and HH performed the research and wrote the manuscript. MO analyzed the data. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Institute of Health to HH (R01MH092535, R01MH125333, R01EB031284, and P50HD105354) and MO (R21MH123930).
Acknowledgments
We would like to thank the lab member, Jessica L. Hyland, for contribution to English editing.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.757838/full#supplementary-material
References
1
AchardS.BullmoreE. (2007). Efficiency and cost of economical brain functional networks.PLoS Comput. Biol.3:e17. 10.1371/journal.pcbi.0030017
2
AdolphsR. (2003). Cognitive neuroscience of human social behaviour.Nat. Rev. Neurosci.4165–178. 10.1038/nrn1056
3
American Psychiatric Association [APA] (2013). Diagnostic and Statistical Manual of Mental Disorders, 5th Edn. Arlington, VA: American Psychiatric Association.
4
AssafM.JagannathanK.CalhounV. D.MillerL.StevensM. C.SahlR.et al (2010). Abnormal functional connectivity of default mode sub-networks in autism spectrum disorder patients.Neuroimage53247–256. 10.1016/j.neuroimage.2010.05.067
5
BehrensT. E.BergH. J.JbabdiS.RushworthM. F.WoolrichM. W. (2007). Probabilistic diffusion tractography with multiple fibre orientations: what can we gain?Neuroimage34144–155. 10.1016/j.neuroimage.2006.09.018
6
Ben BashatD.Kronfeld-DueniasV.ZachorD. A.EksteinP. M.HendlerT.TarraschR.et al (2007). Accelerated maturation of white matter in young children with autism: a high b value DWI study.Neuroimage3740–47. 10.1016/j.neuroimage.2007.04.060
7
BucknerR. L.Andrews-HannaJ. R.SchacterD. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. NY Acad. Sci.1124, 1–38. 10.1196/annals.1440.011
8
BullmoreE.SpornsO. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems.Nat. Rev. Neurosci.10186–198. 10.1038/nrn2575
9
CaoM.HuangH.HeY. (2017). Developmental connectomics from infancy through early childhood.Trends Neurosci.40494–506. 10.1016/j.tins.2017.06.003
10
CastelliF.FrithC.HappéF.FrithU. (2002). Autism, Asperger syndrome and brain mechanisms for the attribution of mental states to animated shapes.Brain1251839–1849. 10.1093/brain/awf189
11
CiarrustaJ.DimitrovaR.BatalleD.O’MuircheartaighJ.Cordero-GrandeL.PriceA.et al (2020). Emerging functional connectivity differences in newborn infants vulnerable to autism spectrum disorders.Transl. Psychiatry10:131. 10.1038/s41398-020-0805-y
12
ClancyH.DugdaleA.Rendle-ShortJ. (1969). The diagnosis of infantile autism.Dev. Med. Child Neurol.11432–442. 10.1111/j.1469-8749.1969.tb01461.x
13
ColeM. W.RepovšG.AnticevicA. (2014). The frontoparietal control system: a central role in mental health.Neuroscientist20652–664. 10.1177/1073858414525995
14
ConstantinoJ. N.CharmanT. (2016). Diagnosis of autism spectrum disorder: reconciling the syndrome, its diverse origins, and variation in expression.Lancet Neurol.15279–291. 10.1016/S1474-4422(15)00151-9
15
CourchesneE.PierceK. (2005). Why the frontal cortex in autism might be talking only to itself: local over-connectivity but long-distance disconnection.Curr. Opin. Neurobiol.15225–230. 10.1016/j.conb.2005.03.001
16
CrossleyN. A.MechelliA.ScottJ.CarlettiF.FoxP. T.McGuireP.et al (2014). The hubs of the human connectome are generally implicated in the anatomy of brain disorders.Brain1372382–2395. 10.1093/brain/awu132
17
DennisE. L.JahanshadN.McMahonK. L.de ZubicarayG. I.MartinN. G.HickieI. B.et al (2013). Development of brain structural connectivity between ages 12 and 30: a 4-Tesla diffusion imaging study in 439 adolescents and adults.Neuroimage64671–684. 10.1016/j.neuroimage.2012.09.004
18
DesikanR. S.SégonneF.FischlB.QuinnB. T.DickersonB. C.BlackerD.et al (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest.Neuroimage31968–980. 10.1016/j.neuroimage.2006.01.021
19
Di MartinoA.YanC. G.LiQ.DenioE.CastellanosF. X.AlaertsK.et al (2014). The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism.Mol. Psychiatry19659–667. 10.1038/mp.2013.78
20
FischlB.Van Der KouweA.DestrieuxC.HalgrenE.SégonneF.SalatD. H.et al (2004). Automatically parcellating the human cerebral cortex.Cereb. Cortex1411–22. 10.1093/cercor/bhg087
21
HazlettH. C.GuH.MunsellB. C.KimS. H.StynerM.WolffJ. J.et al (2017). Early brain development in infants at high risk for autism spectrum disorder.Nature542348–351. 10.1038/nature21369
22
HuangH.ShuN.MishraV.JeonT.ChalakL.WangZ. J.et al (2015). Development of human brain structural networks through infancy and childhood.Cereb. Cortex251389–1404. 10.1093/cercor/bht335
23
HuangH.ZhangJ.Van ZijlP. C.MoriS. (2004). Analysis of noise effects on DTI-based tractography using the brute-force and multi-ROI approach.Magn. Reson. Med.52559–565. 10.1002/mrm.20147
24
JansiewiczE. M.GoldbergM. C.NewschafferC. J.DencklaM. B.LandaR.MostofskyS. H. (2006). Motor signs distinguish children with high functioning autism and Asperger’s syndrome from controls.J. Autism Dev. Disord.36613–621. 10.1007/s10803-006-0109-y
25
JeonT.MishraV.OuyangM.ChenM.HuangH. (2015). Synchronous changes of cortical thickness and corresponding white matter microstructure during brain development accessed by diffusion MRI tractography from parcellated cortex.Front. Neuroanat.9:158. 10.3389/fnana.2015.00158
26
JiangH.van ZijlP. C.KimJ.PearlsonG. D.MoriS. (2006). DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking.Comput. Methods Programs Biomed.81106–116. 10.1016/j.cmpb.2005.08.004
27
KeownC. L.ShihP.NairA.PetersonN.MulveyM. E.MüllerR. A. (2013). Local functional overconnectivity in posterior brain regions is associated with symptom severity in autism spectrum disorders.Cell Rep.5567–572. 10.1016/j.celrep.2013.10.003
28
KrugD. A.ArickJ.AlmondP. (1980). Behavior checklist for identifying severely handicapped individuals with high levels autistic behavior.J. Child Psychol. Psychiatry21221–229. 10.1111/j.1469-7610.1980.tb01797.x
29
LaiM. C.LombardoM. V.ChakrabartiB.SadekS. A.PascoG.WheelwrightS. J.et al (2010). A shift to randomness of brain oscillations in people with autism.Biol. Psychiatry681092–1099. 10.1016/j.biopsych.2010.06.027
30
LandaR.Garrett-MayerE. (2006). Development in infants with autism spectrum disorders: a prospective study.J. Child Psychol. Psychiatry47629–638. 10.1111/j.1469-7610.2006.01531.x
31
LevmanJ.MacDonaldP.RowleyS.StewartN.LimA.EwensonB.et al (2019). Structural magnetic resonance imaging demonstrates abnormal regionally-differential cortical thickness variability in autism: from newborns to adults.Front. Hum. Neurosci.13:75. 10.3389/fnhum.2019.00075
32
LevmanJ.VasungL.MacDonaldP.RowleyS.StewartN.LimA.et al (2018). Regional volumetric abnormalities in pediatric autism revealed by structural magnetic resonance imaging.Int. J. Dev. Neurosci.7134–45. 10.1016/j.ijdevneu.2018.08.001
33
LiH.XueZ.EllmoreT. M.FryeR. E.WongS. T. (2014). Network-based analysis reveals stronger local diffusion-based connectivity and different correlations with oral language skills in brains of children with high functioning autism spectrum disorders.Hum. Brain Mapp.35396–413. 10.1002/hbm.22185
34
LiS. J.WangY.QianL.LiuG.LiuS. F.ZouL. P.et al (2018). Alterations of white matter connectivity in preschool children with autism spectrum disorder.Radiology288209–217. 10.1148/radiol.2018170059
35
LoomesR.HullL.MandyW. P. L. (2017). What is the male-to-female ratio in autism spectrum disorder? A systematic review and meta-analysis.J. Am. Acad. Child Adolesc. Psychiatry56466–474. 10.1016/j.jaac.2017.03.013
36
LordC.RutterM.Le CouteurA. (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.24659–685. 10.1007/BF02172145
37
LynchC. J.UddinL. Q.SupekarK.KhouzamA.PhillipsJ.MenonV. (2013). Default mode network in childhood autism: posteromedial cortex heterogeneity and relationship with social deficits.Biol. Psychiatry74212–219. 10.1016/j.biopsych.2012.12.013
38
MayK. E.KanaR. K. (2020). Frontoparietal network in executive functioning in autism spectrum disorder.Autism Res.131762–1777. 10.1002/aur.2403
39
MinshewN. J.WilliamsD. L. (2007). The new neurobiology of autism: cortex, connectivity, and neuronal organization.Arch. Neurol.64945–950. 10.1001/archneur.64.7.945
40
MonkC. S.PeltierS. J.WigginsJ. L.WengS. J.CarrascoM.RisiS.et al (2009). Abnormalities of intrinsic functional connectivity in autism spectrum disorders.Neuroimage47764–772. 10.1016/j.neuroimage.2009.04.069
41
MoriS.TournierJ. D. (2014). Introduction to Diffusion Tensor Imaging: And Higher Order Models.Amsterdam: Elsevier Science, 10.1016/C2011-0-07607-X
42
MoriS.CrainB. J.ChackoV. P.van ZijlP. C. (1999). Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging.Ann. Neurol.45265–269. 10.1002/1531-8249(199902)45:2<265:AID-ANA21<3.0.CO;2-3
43
NairA.JalalR.LiuJ.TsangT.McDonaldN. M.JacksonL.et al (2021). Altered thalamocortical connectivity in 6-week-old infants at high familial risk for autism spectrum disorder. Cereb. Cortex31, 4191–4205. 10.1093/cercor/bhab078
44
NebelM. B.EloyanA.BarberA. D.MostofskyS. H. (2014). Precentral gyrus functional connectivity signatures of autism.Front. Syst. Neurosci.8:80. 10.3389/fnsys.2014.00080
45
NomiJ. S.UddinL. Q. (2015). Developmental changes in large-scale network connectivity in autism.NeuroImage Clin.7732–741. 10.1016/j.nicl.2015.02.024
46
OuyangM.ChengH.MishraV.GongG.MosconiM. W.SweeneyJ.et al (2016). Atypical age-dependent effects of autism on white matter microstructure in children of 2–7 years.Hum. Brain Mapp.37819–832. 10.1002/hbm.23073
47
OuyangM.DuboisJ.YuQ.MukherjeeP.HuangH. (2019). Delineation of early brain development from fetuses to infants with diffusion MRI and beyond.Neuroimage185836–850. 10.1016/j.neuroimage.2018.04.017
48
OuyangM.KangH.DetreJ. A.RobertsT. P.HuangH. (2017). Short-range connections in the developmental connectome during typical and atypical brain maturation.Neurosci. Biobehav. Rev.83109–122. 10.1016/j.neubiorev.2017.10.007
49
OuyangM.PengQ.JeonT.HeyneR.ChalakL.HuangH. (2020). Diffusion-MRI-based regional cortical microstructure at birth for predicting neurodevelopmental outcomes of 2-year-olds.Elife9:e58116. 10.7554/eLife.58116
50
OzonoffS.IosifA. M.BaguioF.CookI. C.HillM. M.HutmanT.et al (2010). A prospective study of the emergence of early behavioral signs of autism.J. Am. Acad. Child Adolesc. Psychiatry49256–266. 10.1016/j.jaac.2009.11.009
51
PadmanabhanA.LynchC. J.SchaerM.MenonV. (2017). The default mode network in autism.Biol. Psychiatry Cogn. Neurosci. Neuroimaging2476–486. 10.1016/j.bpsc.2017.04.004
52
PayabvashS.PalaciosE. M.OwenJ. P.WangM. B.TavassoliT.GerdesM.et al (2019). White matter connectome edge density in children with autism spectrum disorders: potential imaging biomarkers using machine-learning models.Brain Connect.9209–220. 10.1089/brain.2018.0658
53
RobertsT. P.KhanS. Y.ReyM.MonroeJ. F.CannonK.BlaskeyL.et al (2010). MEG detection of delayed auditory evoked responses in autism spectrum disorders: towards an imaging biomarker for autism.Autism Res.38–18. 10.1002/aur.111
54
RubinovM.SpornsO. (2010). Complex network measures of brain connectivity: uses and interpretations.Neuroimage521059–1069. 10.1016/j.neuroimage.2009.10.003
55
RudieJ. D.DaprettoM. (2013). Convergent evidence of brain overconnectivity in children with autism?Cell Rep.5565–566. 10.1016/j.celrep.2013.10.043
56
RudieJ. D.BrownJ. A.Beck-PancerD.HernandezL. M.DennisE. L.ThompsonP. M.et al (2013). Altered functional and structural brain network organization in autism.NeuroImage Clin.279–94. 10.1016/j.nicl.2012.11.006
57
ShermanS. M. (2016). Thalamus plays a central role in ongoing cortical functioning.Nature neuroscience19533–541. 10.1038/nn.4269
58
ShiohamaT.OrtugA.WarrenJ. L. A.ValliB.LevmanJ.FajaS. K.et al (2021). Small nucleus accumbens and large cerebral ventricles in infants and toddlers prior to receiving diagnoses of autism spectrum disorder.Cereb. Cortexbhab283. 10.1093/cercor/bhab283[Epub ahead of print].
59
SpornsO. (2011). The non-random brain: efficiency, economy, and complex dynamics.Front. Comput. Neurosci.5:5. 10.3389/fncom.2011.00005
60
SupekarK.UddinL. Q.KhouzamA.PhillipsJ.GaillardW. D.KenworthyL. E.et al (2013). Brain hyperconnectivity in children with autism and its links to social deficits.Cell Rep.5738–747. 10.1016/j.celrep.2013.10.001
61
UddinL. Q.MenonV.YoungC. B.RyaliS.ChenT.KhouzamA.et al (2011). Multivariatete searchlight classification of structural magnetic resonance imaging in children and adolescents with autism.Biol. Psychiatry70833–841. 10.1016/j.biopsych.2011.07.014
62
UddinL. Q.SupekarK.LynchC. J.KhouzamA.PhillipsJ.FeinsteinC.et al (2013). Salience network–based classification and prediction of symptom severity in children with autism.JAMA Psychiatry70869–879. 10.1001/jamapsychiatry.2013.104
63
van den HeuvelM. P.SpornsO. (2019). A cross-disorder connectome landscape of brain dysconnectivity.Nat. Rev. Neurosci.20435–446. 10.1038/s41583-019-0177-6
64
WalkerL.GozziM.LenrootR.ThurmA.BehsetaB.SwedoS.et al (2012). Diffusion tensor imaging in young children with autism: biological effects and potential confounds.Biol. Psychiatry721043–1051. 10.1016/j.biopsych.2012.08.001
65
WangJ.WangX.XiaM.LiaoX.EvansA.HeY. (2015). GRETNA: a graph theoretical network analysis toolbox for imaging connectomics.Front. Hum. Neurosci.9:386. 10.3389/fnhum.2015.00386
66
WolffJ. J.GuH.GerigG.ElisonJ. T.StynerM.GouttardS.et al (2012). Differences in white matter fiber tract development present from 6 to 24 months in infants with autism.Am. J. Psychiatry169589–600. 10.1176/appi.ajp.2011.11091447
67
WoodsR. P.GraftonS. T.HolmesC. J.CherryS. R.MazziottaJ. C. (1998). Automated image registration: I. General methods and intrasubject, intramodality validation.J. Comput. Assist. Tomogr.22139–152. 10.1097/00004728-199801000-00027
68
XiaM.WangJ.HeY. (2013). BrainNet viewer: a network visualization tool for human brain connectomics.PLoS One8:e68910. 10.1371/journal.pone.0068910
69
YerysB. E.GordonE. M.AbramsD. N.SatterthwaiteT. D.WeinblattR.JankowskiK. F.et al (2015). Default mode network segregation and social deficits in autism spectrum disorder: evidence from non-medicated children.NeuroImage: Clin.9223–232. 10.1016/j.nicl.2015.07.018
70
YuQ.PengY.KangH.PengQ.OuyangM.SlingerM.et al (2020). Differential white matter maturation from birth to 8 years of age.Cereb. Cortex.302674–2690. 10.1093/cercor/bhz268
71
ZaleskyA.CocchiL.FornitoA.MurrayM. M.BullmoreE. D. (2012). Connectivity differences in brain networks.Neuroimage601055–1062. 10.1016/j.neuroimage.2012.01.068
72
ZaleskyA.FornitoA.CocchiL.GolloL. L.van den HeuvelM. P.BreakspearM. (2016). Connectome sensitivity or specificity: which is more important?Neuroimage142407–420. 10.1016/j.neuroimage.2016.06.035
73
ZaleskyA.FornitoA.HardingI. H.CocchiL.YücelM.PantelisC.et al (2010). Whole-brain anatomical networks: does the choice of nodes matter?Neuroimage50970–983. 10.1016/j.neuroimage.2009.12.027
74
ZhaoT.MishraV.JeonT.OuyangM.PengQ.ChalakL.et al (2019). Structural network maturation of the preterm human brain.NeuroImage185699–710. 10.1016/j.neuroimage.2018.06.047
75
ZilboviciusM.MeresseI.ChabaneN.BrunelleF.SamsonY.BoddaertN. (2006). Autism, the superior temporal sulcus and social perception.Trends Neurosci.29359–366. 10.1016/j.tins.2006.06.004
76
ZwaigenbaumL.BaumanM. L.ChoueiriR.KasariC.CarterA.GranpeeshehD.et al (2015). Early intervention for children with autism spectrum disorder under 3 years of age: recommendations for practice and research.Pediatrics136S60–S81. 10.1542/peds.2014-3667E
Summary
Keywords
autism spectrum disorder, brain development, structural network, hyper-connectivity, early childhood, symptom severity
Citation
Ouyang M, Peng Y, Sotardi S, Hu D, Zhu T, Cheng H and Huang H (2022) Flattened Structural Network Changes and Association of Hyperconnectivity With Symptom Severity in 2–7-Year-Old Children With Autism. Front. Neurosci. 15:757838. doi: 10.3389/fnins.2021.757838
Received
12 August 2021
Accepted
21 December 2021
Published
14 February 2022
Volume
15 - 2021
Edited by
Jacob Levman, St. Francis Xavier University, Canada
Reviewed by
Ian Max Andolina, Institute of Neuroscience, Shanghai Institute for Biological Sciences, Chinese Academy of Sciences (CAS), China; Janice Hau, San Diego State University, United States
Updates
Copyright
© 2022 Ouyang, Peng, Sotardi, Hu, Zhu, Cheng and Huang.
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: Yun Peng, ppengyun@hotmail.comHao Huang, huangh6@email.chop.edu
This article was submitted to Neurodevelopment, a section of the journal Frontiers in Neuroscience
Disclaimer
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.