- Open Access
Inter-individual heterogeneity of functional brain networks in children with autism spectrum disorder
Molecular Autism volume 13, Article number: 52 (2022)
Autism spectrum disorder (ASD) is a neurodevelopmental disorder with considerable clinical heterogeneity. This study aimed to explore the heterogeneity of ASD based on inter-individual heterogeneity of functional brain networks.
Resting-state functional magnetic resonance imaging data from the Autism Brain Imaging Data Exchange database were used in this study for 105 children with ASD and 102 demographically matched typical controls (TC) children. Functional connectivity (FC) networks were first obtained for ASD and TC groups, and inter-individual deviation of functional connectivity (IDFC) from the TC group was then calculated for each individual with ASD. A k-means clustering algorithm was used to obtain ASD subtypes based on IDFC patterns. The FC patterns were further compared between ASD subtypes and the TC group from the brain region, network, and whole-brain levels. The relationship between IDFC and the severity of clinical symptoms of ASD for ASD subtypes was also analyzed using a support vector regression model.
Two ASD subtypes were identified based on the IDFC patterns. Compared with the TC group, the ASD subtype 1 group exhibited a hypoconnectivity pattern and the ASD subtype 2 group exhibited a hyperconnectivity pattern. IDFC for ASD subtype 1 and subtype 2 was found to predict the severity of social communication impairments and the severity of restricted and repetitive behaviors in ASD, respectively.
Only male children were selected for this study, which limits the ability to study the effects of gender and development on ASD heterogeneity.
These results suggest the existence of subtypes with different FC patterns in ASD and provide insight into the complex pathophysiological mechanism of clinical manifestations of ASD.
Autism spectrum disorder (ASD) is an early-onset neurodevelopmental disorder. Its core symptoms are social interaction and communication impairments, repetitive, stereotyped behaviors, and narrow interests . The etiology of ASD is unclear, and current research generally indicates that the clinical presentation of individuals with ASD is highly heterogeneous [2,3,4]. Due to the high clinical heterogeneity of ASD, the diagnostic criteria have evolved over the years and now include a significantly broader spectrum . Therefore, the interpretation of the heterogeneity pattern of ASD is a pressing issue.
There is a general consensus that ASD is related to an atypical pattern of brain functional networks [6, 7]. Abnormal resting-state brain functional connectivity (FC) patterns in ASD have been successively identified in various brain regions and networks, such as the amygdala, visual network, and default-mode network [8,9,10]. Both reduced and increased patterns of FC in the brain of individuals with ASD have been reported, which are associated with developmental stage [7, 11,12,13,14]. Although studies continue to report atypical FC in the brain of individuals with ASD, the results are often inconsistent , which hinders the study of neurobiological mechanisms of ASD and the identification of effective biomarkers to guide clinical diagnosis and treatment of ASD.
Studying inter-individual differences in brain networks is essential for understanding complex psychiatric disorders . For example, studies of inter-individual differences in functional brain networks have found that it can be used to predict schizophrenia positive symptoms . Highly heterogeneous brain structures are reported in ASD, and there is a link between heterogeneity of brain structures and clinical symptoms, which together influence the ASD classification . In addition, subtyping of ASD is a direction to study the high heterogeneity of ASD. Performance of supervised learning for prediction of Autism Diagnostic Observation Schedule (ADOS) scores can be improved by differentiating multidimensional neuroanatomical subtypes . ASD subtype’s unique brain–behavior relationships were reported in the FC-based ASD subtype study . ASD subtyping studies have been widely applied to analyze the functional and structural networks of the ASD brain and have explained to some extent the heterogeneity of ASD . However, to the best of our knowledge, existing studies on ASD subtypes in the field of ASD functional brain networks still distinguish subtypes solely rely on ASD brain network characteristics, ignoring the brain network information provided by the typical controls (TC) group. Previous studies have suggested that the heterogeneity of atypical FC patterns in ASD stems from idiosyncratic distortions of the FC pattern . Therefore, subtyping studies based on inter-individual deviation could better reveal the heterogeneous pattern of individuals with ASD.
This study aimed to explore the heterogeneity pattern of brain functional network among individuals with ASD based on the inter-individual deviation between ASD and TC. Based on the FC network for each individual, the inter-individual deviation of functional connectivity (IDFC) of each individual with ASD from the TC group is calculated. Individuals with ASD were then clustered based on IDFC to obtain ASD subtypes. The relationship between IDFC and clinical symptoms of ASD was also assessed. We predict that individuals with ASD have different functional network patterns of IDFC and can be subdivided into subtypes based on the IDFC.
Materials and methods
Resting-state functional magnetic resonance imaging (fMRI) and phenotype data from the open-access Autism Brain Imaging Data Exchange database (ABIDE, fcon_1000.projects.nitrc.org/indi/abide/) were used . The participant selection principles were the same as in our previous study . Exclusion criteria include: (a) subjects younger than 7 years of age or older than 12 years of age; (b) female subjects (female subjects were less than 10% of the total data set); (c) subjects with excessive head motion during resting-state scanning (i.e., motion greater than 2 mm translation and 2 degrees of rotation, and greater than 50% of frames with large frame-wise displacement [FD]); (d) subjects with missing information on full intelligence quotient (FIQ), handedness or eye status values; and (e) subjects with low-quality scanned structural images. The remaining subjects maximized p values for group differences in age, handedness, FIQ, eye status, and mean FD by using a data-driven algorithm to create well-matched datasets of ASD and TC groups within each site. Data from centers with fewer than 10 subjects per group were finally removed. The final remaining 105 children with ASD and 102 TC from 6 centers were applied to the study. Participant demographic details are summarized in Table 1.
The advanced edition of Data Processing Assistant for Resting-State fMRI (DPARSF A, http://rfmri.org/DPARSF) toolbox  was used to preprocess the resting-state fMRI data of the subjects. The same preprocessing steps as in our previous study were used . Main steps include: the first 10 volumes of each subject being removed, slicing time corrected, spatial realigned (participants with translational or rotational motion higher than 2 mm or 2° were excluded), normalized to standard Montreal Neurological Institute (MNI) stereotaxic space, and resampling to 3 × 3 × 3 mm3, spatial smoothing with an isotropic Gaussian kernel (full width at half maximum = 6 mm), linear trends were removed, potential motion artifacts were resolved by using the 3dDespike algorithm in Analysis of Functional NeuroImaging (https://afni.nimh.nih.gov/afni/), potential nuisance signals were regressed (Friston-24 motion parameters, white matter, and cerebrospinal fluid signals) [25,26,27], and bandpass filtering was carried out at 0.01–0.1 Hz. In addition, to reduce the effect of head movement on the FC analysis , we examined the percentage of high head movement time points for each subject. Time points with a mean head movement displacement greater than 0.5 mm and their preceding and following two time points were marked as high head movement time points . Subjects with more than 50% high head movement time points were not included in the following analysis [30, 31].
Inter-individual deviation of functional connectivity
The time series data were extracted according to a 264-region of interest (ROI) partitioning scheme, and 10 functional networks were delineated [32, 33]. The functional networks include the motor and somatosensory network (SMN), cingulo-opercular network (CO), auditory network (AN), default-mode network (DMN), visual network (VN), fronto-parietal network (FPN), salience network (SAN), subcortical network (SUB), ventral attention network (VAN), and dorsal attention network (DAN). ROIs that are not part of the 10 functional networks are included in the uncertain network (UND). A 264 × 264 Pearson correlation matrix was constructed for each subject as a measure of inter-regional FC, and we set the negative correlation to zero due to the current ambiguity about the meaning of negative correlation [34, 35]. To prevent skewing of the data distribution due to zeroing of negative correlation, we tested the reliability of the results by retesting the results using the correlation matrix without removing negative correlation (see Additional file 1). The IDFC of subject i from the ASD group at ROI k is defined as:
Fik denotes the overall FC profile of brain region k of subject i with ASD defined by the functional architecture, i.e., Fi(k,:). The vector Fi(k,:) represents the k-th column of the correlation matrix of ASD subject i, which stores the correlation coefficients of brain region k with other brain regions. Similarly, Fjk denotes the overall FC profile of brain region k of subject j in the TC group. To characterize the functional deviations of ASD from the TC group, the cosine distances of the overall FC profile of region k between subject i from the ASD group and all subjects in the TC group were calculated and the mean value was used to describe the IDFC of ASD subject i at region k. At the region level, the IDFC matrix V_roi between the ASD and TC group was obtained by calculating the IDFC for all regions within the ASD group for each subject. The IDFC matrix V_net at the network level was obtained by averaging the IDFC of regions within the same network (Fig. 1A).
Clustering analysis based on IDFC
We hypothesized the existence of subtypes of ASD with high heterogeneity of functional brain network characteristics and clustering analysis was performed on 105 individuals with ASD, as illustrated in Fig. 1B. The analysis was performed using the k-means clustering method, selecting network-level IDFC as the clustering feature, because analyzing the number of features from the region level greater than the number of subjects would lead to unreliable clustering results. The IDFC matrix V_net (a 105 × 11 matrix; 105 represents the number of ASD subjects, and 11 represents the number of networks) was calculated from the network level (including 10 functional networks and an uncertain network). A k-means clustering analysis was performed after regressing out the covariates (age, FIQ, mean FD, handedness, eye status, and sites) from the V_net matrix. For the optimal k value of clustering, it is determined by the mean silhouette value; that is, the silhouette of each point after clustering is calculated and the mean silhouette is calculated. The larger the mean silhouette is, the better the clustering effect is. The k-means clustering was repeated 19 times using k ranging from 2 to 20, and the mean silhouette value was calculated for each clustering analysis.
To determine whether clustering results were influenced by participant demographic characteristics, clinical symptom severity of ASD, and data sites, two-sample t tests were used to analyze whether there were significant differences in age, FIQ, mean FD, and ADOS scores between subtypes of ASD, χ2 tests were used to analyze whether there were significant differences in handedness and eye status between subtypes of ASD, and the distribution of sites was also examined for each subtype of ASD. To exclude the effect of potential heterogeneity among TC on the results, we randomly grouped individuals in the TC group to obtain two well-matched subgroups of TC (51 subjects in subgroup 1 and 51 subjects in subgroup 2). The cluster analysis was then replicated using the two TC subgroups to test for the reliability of this study (see Additional file 1).
Atypical FC patterns of ASD subtypes
The obtained ASD subtypes and TC were compared for FC differences to determine whether there is an atypical FC pattern in the ASD subtypes. The FC differences between each ASD subtype group and TC group and between subtype groups were analyzed at three levels: region, network, and whole-brain. The FC matrix for each group was obtained by Fisher z-transformation of the Pearson’s correlation matrix for all subjects, and two-sample t tests were then performed to analyze the FC differences. Region and network-level results were false discovery rate (FDR) corrected, and statistical significance was set at p < 0.05. Results at the whole-brain level were corrected for Bonferroni (q < 0.05, p = 0.05/3 = 0.017). Age, FIQ, mean FD, handedness, eye status, and sites were used as covariates. At the network level, we examined the t-values of the two-sample t tests in the FC analysis of differences. And we calculated the percentage of connectivity edges with significant FC differences in ROI between networks to the total connectivity edges between networks based on the results of the FC difference analysis at the region level. These analyses were used to determine which networks had anomalies. The same methods were also used to analyze whether there were differences in FC between the whole ASD and TC groups at the region, network, and whole-brain levels.
Relationship between IDFC of ASD subtypes and severity of clinical symptoms of ASD
For each ASD subtype, the relationship between the network-level IDFC (regressing out age, FIQ, mean FD, handedness, eye status, and sites as covariates), used to represent the degree of heterogeneity between individuals with ASD and TC, and the severity of ASD symptoms as assessed by ADOS subscores (i.e., communication, social, restricted and repetitive behaviors) was explored . Only the 82 ASD subjects with complete ADOS subscale score information were used for the follow-up analysis. A multivariate support vector regression method was applied for the analysis, where the ADOS subscores were used as dependent variables and the IDFC of each network was used as the independent variable . Library for Large Linear Classification (LIBLINEAR) (https://www.csie.ntu.edu.tw/~cjlin/liblinear/) toolbox was used in this analysis . The regression analysis was performed using the L2-regularized support vector regression model. The performance of the regression algorithm was evaluated using leave-one-out cross-validation (LOOCV) . Briefly, a model is constructed based on n-1 subjects and then a prediction is made for the remaining subject. After all subjects were predicted by the model, the correlation coefficient R between the predicted and observed values was calculated. Statistical significance was determined by a non-parametric permutation test , by disordering the labels to obtain a new correlation coefficient Rp based on the disordered dataset. The p value is determined by the ratio of the number of times the Rp-value is greater than the R-value in 1000 permutations to the total number of permutations (i.e., 1000). To reflect the contribution of each network, we average the feature weights of each LOOCV as the feature weights of each network.
Inter-individual FC deviation within the group
The FC heterogeneity within the group was also compared between ASD and TC groups. An algorithm similar to IDFC was used to represent the inter-individual FC deviation within the group. At the brain region level, the inter-individual FC deviation matrix Vaa_roi for the ASD group can be obtained by using the FC of the ASD group as Fik and Fjk in the IDFC calculation equation (i ≠ j). Similarly, the inter-individual FC deviation matrix Vtt_roi for the TC group can be calculated. Two-sample t tests were performed on Vaa_roi and Vtt_roi to explore the patterns of differences. Age, FIQ, mean FD, handedness, eye status, and sites were used as covariates. The results were FDR corrected, and the statistical significance was set at p < 0.05.
Subtyping ASD based on IDFC
The optimal value of k for the k-means clustering was determined as 2, since it can be seen from Fig. 2A that the maximum value of the mean silhouette is 0.61 at k = 2. The 105 ASDs were classified into two ASD subtypes (47 subjects in subtype 1 and 58 subjects in subtype 2) by k-means clustering based on IDFC at the network level, as shown in Fig. 2B. The IDFC patterns of the two ASD subtypes at the network level and at the brain region level are shown in Fig. 2C, D, and the IDFC patterns are also presented in table form as detailed in Additional file 1: Tables S1 and S2.
No significant differences were found between ASD subtype 1 and ASD subtype 2 in terms of demographics, including: age, FIQ, handedness, eye status, and mean FD. And no significant differences were found between the two ASD subtype groups in terms of clinical symptom severity by comparing ADOS scores either. Considering that we used multicenter data of subjects, we also examined the data center distribution of both ASD subtypes (Additional file 1: Fig. S1; see Additional file 1 for details).
Atypical FC patterns of ASD subtypes
FC comparison analysis at the region, network, and whole-brain levels was performed for subtype 1, subtype 2, and TC groups, using two-sample t tests, as shown in Fig. 3. Of the 34,716 possible connectivity edges at the region level, the aberrant FC pattern was found to be significantly lower for the subtype 1 group than for the TC group on 37% of the edges (i.e., 12,700 edges), significantly higher for the subtype 2 group than for the TC group on 16% of the edges (i.e., 5623 edges), and significantly higher for the subtype 2 group than for the subtype 1 group on 90% of the edges (i.e., 31,388 edges) (p < 0.05, FDR corrected) (Fig. 3A). In the comparison of the whole ASD and TC groups, no significantly different connectivity edges were found, as detailed in Additional file 1: Fig. S2.
Significant differences were also found in the network-level FC analysis (including intra-network connectivity and inter-network connectivity). The subtype 1 group showed a significant decrease in FC on all 66 connectivity edges compared to the TC group, the subtype 2 group showed a significant increase in FC on 64 connectivity edges compared to the TC group, and the subtype 2 group showed a significant increase in FC on all 66 connectivity edges compared to the subtype 1 group (p < 0.05, FDR corrected) (Fig. 3B). Analysis of the t-values of the two-sample t tests for network-level FC difference analysis revealed that the absolute value of t-values was larger between ASD subtype 1 and TC in SMN-FPN, SMN-SUB, and VAN-VAN, and between ASD subtype 2 and TC in SUB-DMN, SUB-VN, SUB-FPN, and SUB-DAN, as detailed in Additional file 1: Fig. S4A. The analysis of the percentage of connectivity edges with significant FC differences in ROI between networks showed that the percentage of connectivity edges with significant FC differences in ROI between ASD subtype 1 and TC was larger at SMN-FPN, SMN-SUB, CO-DAN, AN-FPN, AN-SUB, and VAN-VAN, and the percentage of connectivity edges with significant FC differences in ROI between ASD subtype 2 and TC was larger at SUB-DMN, SUB-VN, SUB-FPN, and SUB-DAN, as detailed in Additional file 1: Fig. S4B. In the comparison of the whole ASD and TC groups, no significantly different network connectivity edges were found, as detailed in Additional file 1: Fig. S5.
Finally, the FC at the whole-brain level was analyzed and the subtype 1 group showed a significant decrease in FC compared to the TC group, the subtype 2 group showed a significant increase in FC compared to the TC group, and the subtype 2 group showed a significant increase in FC compared to the subtype 1 group (p < 0.05, Bonferroni corrected) (Fig. 3C). In contrast, no significant differences were found comparing the whole ASD group with the TC group, as detailed in Additional file 1: Fig. S6.
Relationship between IDFC of ASD subtypes and severity of ASD clinical symptoms
For each ASD subtype, the multivariate support vector regression model was used to investigate the relationship between IDFC and ASD symptom severity. Model performance was further assessed using LOOCV, and statistical significance was determined by a non-parametric permutation test. The network-level IDFC of ASD subtype 1 predicted the ADOS communication subscore (r = 0.38, p = 0.004; Fig. 4A), and no significant relationship was observed with other subscores of ADOS. The network-level IDFC of ASD subtype 2 predicted the ADOS stereotypic behavior subscore (r = 0.48, p = 0.002; Fig. 4B), and no significant relationship was observed with other subscores of ADOS. The feature weight analysis resulted in ASD subtype 1 exhibiting higher weights for VN, VAN, and DAN and lower weights for DMN, FPN, and SUB, and ASD subtype 2 exhibiting higher weights for SMN, VN, and DAN and lower weights for CO and SUB, as shown in Additional file 1: Fig. S7.
Within-group inter-individual FC deviation differences between the ASD and TC groups
The differences in within-group inter-individual FC deviation between the ASD and TC groups were analyzed from the ROI level using two-sample t tests. Compared with the TC group, the ASD group showed higher within-group inter-individual FC deviation at 20 ROIs and lower at 5 ROIs (p < 0.05, FDR corrected), as shown in Fig. 5. Among them, the 20 ROIs with higher within-group inter-individual FC deviation include: one in SMN, five in DMN, eight in VN, one in FPN, two in SAN, one in SUB, and two in UND. The 5 ROIs with lower within-group inter-individual FC deviation included one in SMN, two in DMN, one in VAN, and one in UND. This result is also presented in table form as detailed in Additional file 1: Table S3.
In this study, IDFC constructed based on FC networks was used to analyze the heterogeneity of functional brain deviations in children with ASD. Two ASD subtypes were obtained by IDFC-based ASD clustering analysis, and aberrant FC patterns of ASD subtypes were identified. The IDFC of ASD subtypes was found to predict the severity of different ASD symptoms separately by a multivariate support vector regression model. By comparing the within-group inter-individual FC deviation of the ASD group with that of the TC group, a complex pattern of differences in brain functional network heterogeneity was found. These results further illustrate the heterogeneity of ASD from the perspective of functional brain networks and provide a new way to study brain heterogeneity in individuals with ASD.
Subtypes of individuals with ASD
ASD can be divided into different subtypes in terms of brain structure  and brain function , as reported in previous studies. Subtyping is widely used in many aspects of ASD research. While 25–40% of children with ASD showed normal brain development in early life and were found to have symptoms of ASD after 18 months of life, there were also children with ASD who were found to have developmental delays within 18 months, and this possibly suggests that there were different subtypes of ASD at the level of brain development . In a study of subtypes of ASD, two subtypes of ASD with differential FC patterns in different resting-state networks were identified, and the two subtypes differed in clinical symptom severity . Another subtype study based on FC showed that each ASD subtype had unique behavioral brain relationships, suggesting the possibility of clinical applications for differentiating ASD subtypes based on FC . Consistent with previous studies, we also support that ASD can be subdivided into subtypes, and the present study reveals two ASD subtypes with different IDFC patterns, further revealing a heterogeneous pattern of functional brain networks in ASD from the perspective of inter-individual deviation.
Although heterogeneity in ASD has been explained by subtypes from different perspectives in previous studies, studies have also attributed heterogeneity in ASD to factors such as age [43,44,45], intelligence , and gender . To exclude the effects of these factors, only male subjects were used in this study and age and FIQ were regressed out prior to subtype identification, and two highly heterogeneous ASD subtypes were finally obtained. That is, the difference in IDFC between the two ASD subtypes was not influenced by age, intelligence, and gender. This means that the results of the present study are not due to these factors, but more likely to the inherent heterogeneity of the functional network of the ASD brain, which reinforces the need for subtyping in the future studies of ASD.
Atypical FC pattern of ASD subtypes
Recent ASD subtype studies reported that different ASD subtypes have different atypical FC patterns [29, 48], and we found similar results for ASD subtype studies based on IDFC. We obtained ASD subtype 1 exhibiting a lower functional connectivity pattern, especially among SMN-FPN, SMN-SUB, and VAN-VAN networks, and ASD subtype 2 exhibiting a higher functional connectivity pattern, especially among SUB-DMN, SUB-VN, SUB-FPN, and SUB-DAN networks, which is similar to the previous findings. Previous studies on ASD brain function reported anomalous FC patterns between SUB network and other networks , and anomalous FC between SMN and FPN network was also reported . The IDFC of the two ASD subtypes predicted different ASD symptom severities possibly attributed to the different atypical FC patterns of the two subtypes, which further illustrates the important role of the two ASD subtypes obtained in this study.
In this study, a subtyping approach was used to study FC abnormality patterns, and IDFC was used as a subtyping feature to emphasize the influence of brain heterogeneity on ASD more than traditional features. The two ASD subtype groups exhibited hypoconnectivity pattern and hyperconnectivity pattern relative to the TC group, and in line with our conjecture, no significant differences were found when comparing the FC of the whole ASD group with that of the TC group. Previous studies have identified complex patterns of FC in resting-state functional brain network studies in ASD, and inconsistent results are common . Studying the entire ASD group led to some potential atypical FC patterns being overlooked or misinterpreted, which is an important reason why the results of previous studies are often inconsistent.
Relationships between IDFC and severity of ASD clinical symptoms
Previous studies have generally reported that abnormal functional networks potentially related to clinical symptoms in individuals with ASD [23, 50]. Clinical correlation analysis based on ASD subtypes to obtain a better prediction of ADOS scores has also been reported , which suggests a link between the heterogeneity of brain networks in individuals with ASD and the heterogeneity of their clinical symptom severity. This view is also supported by our study, in which both ASD subtypes of IDFC were correlated with ASD symptoms but did not behave in the same way, with subtype 1 showing a potential association with the severity of social communication impairments in ASD and subtype 2 showing a relationship with ASD symptoms in terms of severity of restricted and repetitive behaviors. Previous studies also found different potential relationships between different ASD subtypes and ASD symptoms. A previous study classified ASD into three subtypes: Asperger’s, Pervasive Developmental Disorder-Not otherwise Specified, and Autism, and each subtype-specific brain pattern is correlated with different ADOS subdomains . The correlation of one ASD subtype with core ASD symptoms was reported in a subtype study of ASD by resting-state FC and was not found in other subtypes . ASD subtyping studies hold promise as a way to reflect the clinical heterogeneity of ASD.
The two ASD subtypes showed different network feature weights in the brain–behavioral analysis, indicating that the contribution of different functional networks to the prediction of symptoms differed between the two subtypes, and this difference could potentially account for the difference between the relationships between the two subtypes and ASD symptoms; this difference in relationship could also be attributed to the hypoconnectivity and hyperconnectivity patterns of the two subtypes. Our study on the potential relationship between deviant features of brain networks and clinical symptoms may provide a promising way to explain the high heterogeneity of clinical symptoms of ASD and boost the development of clinical diagnostic markers.
Inter-individual FC deviation within the group differences between the ASD and TC groups
Brain heterogeneity in ASD is different from that in TC in the present study. The ASD group showed higher inter-individual FC deviation within the group at 20 ROIs and lower levels at 5 ROIs, i.e., greater individual differences in FC at more ROIs. The results of the restudy using two TC subgroups also found greater individual differences in FC at more ROIs for ASD compared to the TC subgroups. This suggests broader inter-individual variation among ASD compared to TC, consistent with previous findings on individuals brain heterogeneity in ASD [16, 21, 29]. Moreover, significantly different ROIs were located in multiple networks (including 7 functional networks and an uncertain network), and significantly different ROIs were also found to be located in multiple networks in a complementary analysis using two TC subgroups, which further corroborates the high complexity of ASD and the need to study it on a large scale [53, 54].
First, due to the relatively limited number of female participants in the ABIDE database, only male children were selected, which prevented analysis of the effects of gender. Because gender is a factor that influences ASD heterogeneity , there is a need to focus on collecting data on female individuals with ASD in future studies to be able to explore gender heterogeneity in ASD. Second, previous studies have generally concluded that ASD is a developmental disorder [56, 57] and the heterogeneity of functional brain networks of individuals with ASD at other developmental stages, such as adolescence and adulthood, remains unknown. In order to figure out the development of heterogeneity of functional brain networks in ASD, longitudinal data including subjects at different developmental stages are needed for future longitudinal studies. Third, recent studies have shown that resting-state FC is a dynamic process [31, 58, 59], and it is necessary to develop dynamic subtyping methods in the future, which may better reveal ASD heterogeneity. Finally, the present study was a unimodal study based on inter-individual heterogeneity of functional brain networks and the influence of structural brain networks needs to be considered in future studies. A multimodal study that introduces structural brain networks could provide a more comprehensive analysis of heterogeneity in ASD.
This study focused on the heterogeneity among individuals with ASD based on functional brain networks deviations of ASD from TC, revealing two ASD subtypes with distinct FC patterns at multiple spatial scales. Furthermore, IDFC of ASD subtypes was found to predict the severity of different ASD symptoms. These results reveal individual heterogeneity of functional brain networks in ASD and highlight the important role of individual heterogeneity in explaining the complex brain patterns of ASD.
Availability of data and materials
The data that support the findings of this study are available in Autism Brain Imaging Data Exchange (ABIDE, fcon_1000.projects.nitrc.org/indi/abide/) database.
Autism spectrum disorder
Inter-individual deviation of functional connectivity
Autism Diagnostic Observation Schedule
Communication subscore of the ADOS
Social subscore of the ADOS
Stereotyped behaviors and restricted interests subscore of the ADOS
Functional magnetic resonance imaging
Autism Brain Imaging Data Exchange database
Full intelligence quotient
Data Processing Assistant for Resting-State fMRI
Montreal Neurological Institute
Region of interest
Ventral attention network
Dorsal attention network
False discovery rate
Library for Large Linear Classification
American Psychiatric Association. Diagnostic and statistical manual of mental disorders (DSM-5®). Arlington: American Psychiatric Publishing. 2013.
Huerta M, Lord C. Diagnostic evaluation of autism spectrum disorders. Pediatr Clin. 2012;59(1):103–11.
Lombardo MV, Lai MC, Baron-Cohen S. Big data approaches to decomposing heterogeneity across the autism spectrum. Mol Psychiatry. 2019;24(10):1435–50.
Lombardo MV, Pramparo T, Gazestani V, Warrier V, Bethlehem RAI, Carter Barnes C, et al. Large-scale associations between the leukocyte transcriptome and BOLD responses to speech differ in autism early language outcome subtypes. Nat Neurosci. 2018;21(12):1680–8.
Wing L, Potter D. The epidemiology of autistic spectrum disorders: is the prevalence rising? Ment Retard Dev Disabil Res Rev. 2002;8(3):151–61.
Rane P, Cochran D, Hodge SM, Haselgrove C, Kennedy DN, Frazier JA. Connectivity in autism: a review of MRI connectivity studies. Harv Rev Psychiatry. 2015;23(4):223–44.
Muller RA, Shih P, Keehn B, Deyoe JR, Leyden KM, Shukla DK. Underconnected, but how? a survey of functional connectivity MRI studies in autism spectrum disorders. Cereb Cortex. 2011;21(10):2233–43.
Guo X, Duan X, Long Z, Chen H, Wang Y, Zheng J, et al. Decreased amygdala functional connectivity in adolescents with autism: a resting-state fMRI study. Psychiatry Res Neuroimaging. 2016;257:47–56.
McKinnon CJ, Eggebrecht AT, Todorov A, Wolff JJ, Elison JT, Adams CM, et al. Restricted and repetitive behavior and brain functional connectivity in infants at risk for developing autism spectrum disorder. Biol Psychiatry-Cognit Neurosci Neuroimaging. 2019;4(1):50–61.
Lynch CJ, Uddin LQ, Supekar K, Khouzam A, Phillips J, Menon V. Default mode network in childhood autism: posteromedial cortex heterogeneity and relationship with social deficits. Biol Psychiatry. 2013;74(3):212–9.
Just MA, Cherkassky VL, Keller TA, Minshew NJ. Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity. Brain. 2004;127(8):1811–21.
Belmonte MK, Allen G, Beckel-Mitchener A, Boulanger LM, Carper RA, Webb SJ. Autism and abnormal development of brain connectivity. J Neurosci. 2004;24(42):9228–31.
Supekar K, Uddin Lucina Q, Khouzam A, Phillips J, Gaillard William D, Kenworthy Lauren E, et al. Brain hyperconnectivity in children with autism and its links to social deficits. Cell Rep. 2013;5(3):738–47.
Padmanabhan A, Lynn A, Foran W, Luna B, O’Hearn K. Age related changes in striatal resting state functional connectivity in autism. Front Hum Neurosci. 2013;7:814.
Hull JV, Dokovna LB, Jacokes ZJ, Torgerson CM, Irimia A, Van Horn JD. Resting-state functional connectivity in autism spectrum disorders: a review. Front Psych. 2017;7:205.
Agelink van Rentergem JA, Deserno MK, Geurts HM. Validation strategies for subtypes in psychiatry: a systematic review of research on autism spectrum disorder. Clin Psychol Rev. 2021;87:102033.
Fan Y-S, Li L, Peng Y, Li H, Guo J, Li M, et al. Individual-specific functional connectome biomarkers predict schizophrenia positive symptoms during adolescent brain maturation. Hum Brain Mapp. 2020. https://doi.org/10.1002/hbm.25307.
Katuwal GJ, Baum SA, Cahill ND, Michael AM. Divide and conquer: sub-grouping of ASD improves ASD detection based on brain morphometry. PLoS ONE. 2016;11(4):e0153331.
Hong S-J, Valk SL, Di Martino A, Milham MP, Bernhardt BC. Multidimensional neuroanatomical subtyping of autism spectrum disorder. Cereb Cortex. 2018;28(10):3578–88.
Easson AK, Fatima Z, McIntosh AR. Functional connectivity-based subtypes of individuals with and without autism spectrum disorder. Netw Neurosci. 2019;3(2):344–62.
Hahamy A, Behrmann M, Malach R. The idiosyncratic brain: distortion of spontaneous connectivity patterns in autism spectrum disorder. Nat Neurosci. 2015;18(2):302–9.
Di Martino A, Yan CG, Li Q, Denio E, Castellanos FX, Alaerts K, et al. The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. Mol Psychiatry. 2014;19(6):659–67.
Guo X, Duan X, Chen H, He C, Xiao J, Han S, et al. Altered inter- and intrahemispheric functional connectivity dynamics in autistic children. Hum Brain Mapp. 2020;41(2):419–28.
Yan C, Zang Y. DPARSF: a MATLAB toolbox for" pipeline" data analysis of resting-state fMRI. Front Syst Neurosci. 2010. https://doi.org/10.3389/fnsys.2010.00013/full.
Friston KJ, Williams S, Howard R, Frackowiak RS, Turner R. Movement-related effects in fMRI time-series. Magn Reson Med. 1996;35(3):346–55.
Satterthwaite TD, Wolf DH, Loughead J, Ruparel K, Elliott MA, Hakonarson H, et al. Impact of in-scanner head motion on multiple measures of functional connectivity: relevance for studies of neurodevelopment in youth. Neuroimage. 2012;60(1):623–32.
Yan C-G, Cheung B, Kelly C, Colcombe S, Craddock RC, Di Martino A, et al. A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. Neuroimage. 2013;76:183–201.
Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage. 2012;59(3):2142–54.
Chen H, Uddin LQ, Guo X, Wang J, Wang R, Wang X, et al. Parsing brain structural heterogeneity in males with autism spectrum disorder reveals distinct clinical subtypes. Hum Brain Mapp. 2019;40(2):628–37.
Chen H, Nomi JS, Uddin LQ, Duan X, Chen H. Intrinsic functional connectivity variance and state-specific under-connectivity in autism. Hum Brain Mapp. 2017;38(11):5740–55.
Guo X, Duan X, Suckling J, Chen H, Liao W, Cui Q, et al. Partially impaired functional connectivity states between right anterior insula and default mode network in autism spectrum disorder. Hum Brain Mapp. 2019;40(4):1264–75.
Power JD, Cohen AL, Nelson SM, Wig GS, Barnes KA, Church JA, et al. Functional network organization of the human brain. Neuron. 2011;72(4):665–78.
Cole MW, Reynolds JR, Power JD, Repovs G, Anticevic A, Braver TS. Multi-task connectivity reveals flexible hubs for adaptive task control. Nat Neurosci. 2013;16(9):1348–55.
Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage. 2009;44(3):893–905.
Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci. 2005;102(27):9673–8.
Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, et al. The autism diagnostic observation schedule—generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000;30(3):205–23.
Drucker H, Burges CJ, Kaufman L, Smola A, Vapnik V. Support vector regression machines. Adv Neural Inf Process Syst. 1997;9:155-61.
Fan R-E, Chang K-W, Hsieh C-J, Wang X-R, Lin C-J. LIBLINEAR: A library for large linear classification. J Mach Learn Res. 2008;9:1871–4.
Schölkopf B, Smola AJ, Bach F. Learning with kernels: support vector machines, regularization, optimization, and beyond. Cambridge: MIT Press; 2002.
Golland P, Fischl B, editors. Permutation tests for classification: towards statistical significance in image-based studies. In: Biennial international conference on information processing in medical imaging. Springer; 2003.
Werner E, Dawson G. Validation of the phenomenon of autistic regression using home videotapes. Arch Gen Psychiatry. 2005;62(8):889–95.
Chen C, Bailey B, Muller R, editors. Towards autism subtypes? Unsupervised machine learning using fcMRI features. In: Oral presentation at the organization for human brain mapping annual meeting, Honolulu, HI, USA; 2015.
Aylward EH, Minshew NJ, Field K, Sparks B, Singh N. Effects of age on brain volume and head circumference in autism. Neurology. 2002;59(2):175–83.
Lin HY, Ni HC, Lai MC, Tseng WYI, Gau SSF. Regional brain volume differences between males with and without autism spectrum disorder are highly age-dependent. Mol Autism. 2015;6(1):1–18.
He C, Chen H, Uddin LQ, Erramuzpe A, Bonifazi P, Guo X, et al. Structure-function connectomics reveals aberrant developmental trajectory occurring at preadolescence in the autistic brain. Cereb Cortex. 2020;30(9):5028–37.
Salmond CH, Vargha-Khadem F, Gadian DG, de Haan M, Baldeweg T. Heterogeneity in the patterns of neural abnormality in autistic spectrum disorders: evidence from ERP and MRI. Cortex. 2007;43(6):686–99.
Lai M-C, Lombardo MV, Suckling J, Ruigrok AN, Chakrabarti B, Ecker C, et al. Biological sex affects the neurobiology of autism. Brain. 2013;136(9):2799–815.
Reardon AM, Li K, Langley J, Hu XP. Subtyping autism spectrum disorder via joint modeling of clinical and connectomic profiles. Brain Connect. 2022;12(2):193–205.
Maximo JO, Kana RK. Aberrant “deep connectivity” in autism: a cortico–subcortical functional connectivity magnetic resonance imaging study. Autism Res. 2019;12(3):384–400.
Gotts SJ, Simmons WK, Milbury LA, Wallace GL, Cox RW, Martin A. Fractionation of social brain circuits in autism spectrum disorders. Brain J Neurol. 2012;135(Pt 9):2711–25.
Qi S, Morris R, Turner JA, Fu Z, Jiang R, Deramus TP, et al. Common and unique multimodal covarying patterns in autism spectrum disorder subtypes. Mol Autism. 2020;11(1):1–15.
Tang S, Sun N, Floris DL, Zhang X, Di Martino A, Yeo BT. Reconciling dimensional and categorical models of autism heterogeneity: a brain connectomics and behavioral study. Biol Psychiat. 2020;87(12):1071–82.
Goch CJ, Stieltjes B, Henze R, Hering J, Poustka L, Meinzer HP, et al. Quantification of changes in language-related brain areas in autism spectrum disorders using large-scale network analysis. Int J Comput Assist Radiol Surg. 2014;9(3):357–65.
Uddin LQ, Supekar K, Lynch CJ, Khouzam A, Phillips J, Feinstein C, et al. Salience network-based classification and prediction of symptom severity in children with autism. JAMA Psychiat. 2013;70(8):869–79.
Halladay AK, Bishop S, Constantino JN, Daniels AM, Koenig K, Palmer K, et al. Sex and gender differences in autism spectrum disorder: summarizing evidence gaps and identifying emerging areas of priority. Mol Autism. 2015;6(1):1–5.
Guo X, Chen H, Long Z, Duan X, Zhang Y, Chen H. Atypical developmental trajectory of local spontaneous brain activity in autism spectrum disorder. Sci Rep. 2017;7:39822.
Henry TR, Dichter GS, Gates K. Age and gender effects on intrinsic connectivity in autism using functional integration and segregation. Biol Psychiatry Cogn Neurosci Neuroimaging. 2018;3(5):414–22.
Li Y, Zhu YY, Nguchu BA, Wang YM, Wang HJ, Qiu BS, et al. Dynamic functional connectivity reveals abnormal variability and hyper-connected pattern in autism spectrum disorder. Autism Res. 2020;13(2):230–43.
Gao Y, Sun J, Cheng L, Yang Q, Li J, Hao Z, et al. Altered resting state dynamic functional connectivity of amygdala subregions in patients with autism spectrum disorder: a multi-site fMRI study. J Affect Disord. 2022;312:69–77.
Funding sources for the datasets comprising the 1000 Functional Connectome Project are listed at fcon_1000.projects.nitrc.org/fcpClassic/FcpTable.html. Funding sources for the ABIDE dataset are listed at fcon_1000.projects.nitrc.org/indi/abide/.
This work was supported by the Natural Science Foundation of Hebei Province (No. H2021203002), Funding Project for the Returned Overseas Chinese Scholars of Hebei Province of China (No. C20220334), and Hebei Key Laboratory Project (No. 202250701010046). Junfeng Liu was supported by the National Natural Science Foundation of China (No. 81901199). Dong Cui was supported by the National Natural Science Foundation of China (No. 62173291).
Ethics approval and consent to participate
All experimental protocols were approved by the local Institutional Review Boards. Written informed consent was obtained from a parent/guardian, and assent was obtained from each child. Detailed information on scanning information, diagnostic protocols, and ethical statements can be found at fcon_1000.projects.nitrc.org/indi/abide/.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Guo, X., Zhai, G., Liu, J. et al. Inter-individual heterogeneity of functional brain networks in children with autism spectrum disorder. Molecular Autism 13, 52 (2022). https://doi.org/10.1186/s13229-022-00535-0
- Autism spectrum disorder
- Functional magnetic resonance imaging
- Functional connectivity
- k-means clustering