Analysis of Brain Functional Connectivity Neural Circuits in Children With Autism Based on Persistent Homology
Di Liang1 
Shengxiang Xia
Weiwei Zhang- 1School of Science, Shandong Jianzhu University, Jinan, China
- 2School of Control Science and Engineering, Shandong University, Jinan, China
Autism spectrum disorder (ASD) is a complex neuropsychiatric disorder with a complex and unknown etiology. Statistics demonstrate that the number of people diagnosed with ASD is increasing in countries around the world. Currently, although many neuroimaging studies indicate that ASD is characterized by abnormal functional connectivity (FC) patterns within brain networks rather than local functional or structural abnormalities, the FC characteristics of ASD are still poorly understood. In this study, a Vietoris-Rips (VR) complex filtration model of the brain functional network was established by using resting-state functional magnetic resonance imaging (fMRI) data of children aged 6–13 years old [including 54 ASD patients and 52 typical development (TD) controls] from the Autism Brain Imaging Data Exchange (ABIDE) public database. VR complex filtration barcodes are calculated by using persistent homology to describe the changes in the FC neural circuits of brain networks. The number of FC neural circuits with different length ranges at different threshold values is calculated by using the barcodes, the different brain regions participating in FC neural circuits are discussed, and the connectivity characteristics of brain FC neural circuits in the two groups are compared and analyzed. Our results show that the number of FC neural circuits with lengths of 8–12 is significantly decreased in the ASD group compared with the TD control group at threshold values of 0.7, 0.8 and 0.9, and there is no significant difference in the number of FC neural circuits with lengths of 4–7 and 13–16 and lengths 16. When the thresholds are 0.7, 0.8, and 0.9, the number of FC neural circuits in some brain regions, such as the right orbital part of the superior frontal gyrus, the left supplementary motor area, the left hippocampus, and the right caudate nucleus, involved in the study is significantly decreased in the ASD group compared with the TD control group. The results of this study indicate that there are significant differences in the FC neural circuits of brain networks in the ASD group compared with the TD control group.
1. Introduction
Autism spectrum disorder (ASD), also known as autism, is a complex neurodevelopmental disorder. It usually develops before the age of 3 years and includes the typical clinical features of impaired language communication and social communication skills, a narrow range of interests, and repetitive and stereotyped behaviors (American-Psychiatric-Association, 2013). In 2019, the “China autism education rehabilitation III industry development status report” showed that there are more than 10 million autistic people in China, including more than 2 million children with ASD, and the number is increasing at a rate of ~2% annually. According to the Centers for Disease Control and Prevention (CDC), one out of every 59 children in the U.S. was diagnosed with ASD in 2018. Statistics demonstrate that the number of people diagnosed with ASD is increasing worldwide. The etiology of ASD is complex, the pathogenesis is still unclear, and the early symptoms are not typical. Over the years, the clinical diagnosis of ASD has been mainly based on typical clinical symptoms, the Diagnostic and Statistical Manual of Mental Disorders (Fifth Edition) and the Childhood Autism Rating Scale. Identification and diagnosis are subjective and less stable. Professor Laurent Mottron's research team (Rdgaard et al., 2019) found that the objective performance differences between individuals with and without autism are actually narrowing due to changes in and the updating of diagnostic practices. Therefore, it is becoming increasingly important to carry out quantitative assisted diagnosis through the study of the neurophysiological basis of ASD.
The emergence of functional magnetic resonance imaging (fMRI) has brought new developments to cognitive neuroscience. Functional connectivity (FC) is described as the temporal correlation of a neurophysiological index measured in different brain regions (Friston et al., 1993). FMRI realizes functional brain imaging by examining magnetic field changes caused by the blood oxygen level-dependent (BOLD) response entering brain cells (Ogawa et al., 1990). In the resting state, spontaneous neural activities in spatially separated brain regions involved in the execution of separate cognitive functions, such as default mode and executive control, are highly correlated (Salvador et al., 2005; He et al., 2009). Resting-state FC can refer to the relationship between neuronal activation patterns in brain regions and reflect the information exchange between brain regions (Cole et al., 2016; Tomasi and Volkow, 2019). Research (Lee et al., 2013; Xu et al., 2013) on resting-state fMRI can reveal the functional integration characteristics, dynamic change rules of spontaneous brain activities, and topological properties of the FC neural circuits of brain networks.
A considerable amount of neuroimaging literature (Cherkassky et al., 2006; Kennedy and Adolphs, 2012; Martino et al., 2014) indicates that ASD is characterized by abnormal overall connectivity patterns within brain networks. Just et al. (2004) believed that the deficit in verbal information-integration ability in ASD patients is related to decreased functional brain connectivity. Wicker et al. (2008) believes that abnormalities in FC among the amygdala, prefrontal cortex, and posterior association cortex are the main obstacles for integrating emotional information in ASD patients. Zhang (2015) believed that the deficits in controlling and regulating attention and alertness of ASD patients may be caused by damage to the bottom-up neural circuits with the basal nucleus of the amygdala-primary visual cortex. Analysis of brain FC networks (Kana et al., 2011; Supekar et al., 2013) based on resting-state fMRI will provide an important biological basis for the diagnosis and treatment of ASD from the perspective of FC.
Persistent homology, a newly developed method in topological data analysis (Carlsson, 2009), is an effective tool for analyzing high-dimensional nonlinear data and exploring its nonlinear structures. In fact, the concept of persistence was introduced by Robins (1999), and refined by Edelsbrunner et al. (2002) and Zomorodian and Carlsson (2005). The most important advantage of persistent homology is that it is not restricted by threshold values. Persistent homology has been applied in many kinds of data (including brain data) as described in the section 1 (Carlsson, 2020). For example, the research results of Victor et al. (2018) show that the persistent homology method is a powerful new tool for analyzing brain networks; Lee et al. (2018) used persistent homology to detect abnormal holes in brain network connectivity and provided another representation of the shape of brain networks; Pereira and de Mello (2015) developed a framework for clustering time-series and spatial data based on topological properties, which could detect similarities in recurrent behavior for time series data sets and spatial structures in spatial data sets, such as electroencephalography data from alcoholic subjects; Kanari et al. (2018) introduced the Topological Morphology Descriptor, derived from principles of persistent homology, which could reliably describe the shape of neurons and advance our understanding of the anatomy and diversity of the neuronal morphologies; Cassidy et al. (2015) have developed a new distance measure to discriminate resting state fMRI connectivity using persistent homology. In this study, with the resting-state fMRI data of some children aged 6–13 from the “Autism Brain Imaging Data Exchange (ABIDE)” public database, Vietoris-Rips (VR) complex filtration models were established for ASD patients and typical development (TD) controls based on their brain FC. The barcodes of VR complex filtration were calculated by the persistent homology method. Therefore, FC neural circuits of ASD patients and TD controls at different threshold values can be obtained from the barcodes. The change rules of FC neural circuits in ASD patients are found by analyzing the number of FC neural circuits of different length ranges and brain regions involved at different threshold values; furthermore, it is expected to provide objective biological markers for the early diagnosis of ASD.
2. Materials and Methods
2.1. Data Acquisition
The resting-state fMRI data used in this paper are from the four collection sites, ABIDEII-NYU_1, ABIDEII-NYU_2, ABIDEII-SU_2, and ABIDEII-EMC_1, in the ASD Public Dataset (ABIDE) (Martino and Milham, 2017). The imaging data conditions of the ASD group and the TD control group are as follows: (1) there (Moran et al., 2014) was a strong gender effect in ASD patients, so the subjects were all male; (2) all subjects were dextral; (3) there was an age effect on brain FC, so the age of the subjects was 6–13 years old, including 54 ASD patients [mean (9.3 ± 1.9) years old] and 52 TD controls [mean (9.6 ± 1.9) years old], and there was no significant difference in age between the two groups (p = 0.31).
2.2. Resting-State fMRI Processing
In this study, the tool DPABI (Yan et al., 2016) (http://rfmri.org/DPABI) is used for preprocessing resting-state fMRI data. The preprocessing of data mainly includes the following steps: (1) Removing the first 10 items from the time series and realigning; (2) Cleaning low frequency scanner drift and applying bandpass filtering (0.01–0.1 Hz); (3) Normalizing spatial data by DARTEL; (4) Selecting anatomical automatic labeling (AAL) (Tzourio-Mazoyer et al., 2002) to define regions of interest (ROIs); and (5) Calculating the FC matrix based on BOLD time series for each pair of brain regions. For more data preprocessing details, please refer to the Supplementary Material.
2.3. A VR Complex Filtration Model of Brain Networks
Let (X, D) be a metric space and a Vietoris-Rips (VR) complex with ε parameters be a simple complex denoted by VR(X, ε), whose vertex is the set X, and {x0, x1, ..., xk} generates a k-complex if and only if d(xi, xj) ≤ ε for all 0 ≤ i, j ≤ k.
A correlation matrix of the brain functional network is obtained by using DPABI. Each brain region is a vertex in the brain functional network. We use the Pearson correlation coefficient to define the distance between two vertices Pi and Pj,; therefore, the metric space (X, dX) is formed from the brain functional network. The correlation matrix is sorted from smallest to largest by the distance between the two vertices: ε1 < ε2 < … < εM. We can obtain a VR filtration model of the brain network:
JavaPlex software is used to calculate the VR filtration barcode (the barcode-interval set). These intervals represent the lifetime of nontrivial cycles in the growing process of complexes. The left end of the interval represents the birth of a new topological property, and the right end of the interval represents its death. From this, we can describe the change rules of the Betti number of a brain network. The 0-dimensional Betti number β0 represents the number of components of the brain functional network, and the one-dimensional Betti number β1 represents the number of one-dimensional holes in the brain functional network, that is, the number of FC neural circuits in the brain network. VR complexes constructed from a point cloud data are described in the following example to help understand persistence homology.
For a set X consisting of 5 points on the plane, we construct VR complexes at different values of ε. Let's draw five circles centered at the points with the radius ε/2. The VR complexes are shown in Figure 1. We explain each subfigure in detail to help understand Figure 1. In Figure 1A, there is no one-dimensional simplex formed at ε = 0.60. There are 5 components, this means β0 = 5. In Figure 1B, when the distance between the two points is less than ε = 1.37, the two lower right points form 1-dimensional simplex. There are 4 components, this means β0 = 4. In Figure 1C, the two bottom points are connected, and three bottom points form a new component when ε = 1.62, that is β0 = 3. In Figure 1D, in this case ε = 1.78, β0 = 2. In Figure 1E, a loop is formed, i.e., β1 = 1 and β0 = 1 when ε = 1.97. Then we calculate the persistence barcode according to persistence homology. The persistence barcodes of the data are shown in Figure 2.
Figure 1. (A–E) are VR (X, ε) complexes of the point cloud data with ε = 0.60, 1.37, 1.62, 1.78, and 1.97 respectively.
Figure 2. 0-D and 1-D persistence barcodes of the point cloud data.
We select AAL to define the ROIs, and a brain network has 90 vertices. Figure 3 presents the results obtained by calculating two original data points of an ASD patient and a TD control in ABIDE II with the above steps. A vertical line is drawn across the x axis, and the vertical coordinate of the intersection point of the vertical line and the 0-dimensional (1-dimensional) Betti curve is the 0-dimensional (1-dimensional) Betti number of the network corresponding to the threshold value. Figure 3 shows the curves of the 0-dimensional and 1-dimensional Betti numbers of VR filtration for the ASD patient and the TD control. It depicts the differences in the topological features of the ASD brain network and the TD brain network at different threshold values.
Figure 3. 0-D and 1-D Betti number curves of the ASD brain network and the TD brain network.
2.3.1. Comparison of FC in Brain Networks
The number of components of the brain network (i.e., β0) can be easily calculated from the 0-dimensional Betti curve. For example, the brain network based on AAL consists 90 brain regions. When the threshold value is 0.7, some brain regions are connected and form 83 edges, there are 7 components in the ASD patient. All brain regions are connected and form 89 edges, there is 1 component in the TD control. Figure 4 shows that there are more components in the ASD brain network than in the TD brain network. This can be interpreted as FC between brain regions of the ASD patient not being as tight as that of the TD control, and the deficit in FC between brain regions means that they cannot transmit information.
Figure 4. The ASD brain network (left) including 83 edges and the TD brain network (right) including 89 edges at a threshold value of 0.7.
2.3.2. Comparison of FC Neural Circuits in Brain Networks
The number of FC neural circuits in the brain network (i.e., β1) can be easily calculated from the one-dimensional Betti curve. For example, when the threshold value is 0.7, there are 32 FC neural circuits in the ASD patient and 51 FC neural circuits in the TD control. The number of FC neural circuits in the brain network at a threshold value of 0.7 means the total number of FC neural circuits formed in the brain network when the threshold value is no more than 0.7. Figure 5 shows part of the FC neural circuits in the ASD brain network and in the TD brain network. Figure 6 shows the average numbers of FC neural circuits in each brain region involved in the ASD group and the TD control group when the threshold value is 0.7. The brain regions and the number of brain regions involved in FC neural circuits are different. Differences in FC neural circuits may account for some of the dysfunction experienced by ASD patients.
Figure 5. Comparison of the FC neural circuits in the ASD brain network (left) and the TD brain network (right) at a threshold value of 0.7. The abbreviations for names of brain regions in the AAL atlas are listed in the Supplementary Material. The differences of the colors of edges are just for easy to observe.
Figure 6. The histogram of the average number of FC neural circuits in each brain region in AAL in the ASD group and the TD control group.
3. Results
Two-sample Mann-Whitney U-tests are performed on the data from the ASD group and the TD control group to assess the group differences in the number of FC neural circuits with different circuit length ranges. FDR adjustments are performed. The statistical results show that when the threshold value is 0.7, 0.8, or 0.9, the number of FC neural circuits with a length of 8–12 is significantly decreased in the ASD group compared with in the TD control group, and there is no significant difference in the number of short neural circuits with a length of 4–7, long neural circuits with a length of 13–16, and ultralong neural circuits with a length 16 (see Figure 7 and Table 1).
Figure 7. Violin pictures of the comparison of the numbers of FC neural circuits with different circuit length ranges of the two groups.
Table 1. Statistically significant differences in the number of FC neural circuits (8–12) between the two groups.
Two-sample Mann-Whitney U-tests are performed on the data from the ASD group and the TD control group to assess the group differences in the number of FC neural circuits in each involved brain region. FDR adjustments are performed. The statistical results show that at different threshold values, there are significant differences in the numbers of FC neural circuits in some involved brain regions. Compared with the TD control group, at different threshold values, the number of FC neural circuits in some involved brain regions is significantly decreased in the ASD group. These brain regions include the right orbital part of the superior frontal gyrus, bilateral orbital part of the middle frontal gyrus, left opercular part of the inferior frontal gyrus, left supplementary motor area, bilateral posterior cingulate gyrus, left hippocampus, right parahippocampal gyrus, left angular gyrus, right caudate nucleus, and bilateral Heschl gyrus (see Figure 8 and Table 2).
Figure 8. Violin pictures of the comparison of the numbers of FC neural circuits in significantly different brain regions involved in the two groups.
Table 2. Statistically significant differences in the numbers of FC neural circuits in the involved brain regions of the two groups.
4. Discussion
The lengths and number of FC neural circuits determine the efficiency of information transmission processing. By comparative studies, it was found that the mid-length FC neural circuits of spontaneous activities in the ASD group were significantly decreased compared with those in the TD control group, which is speculated to result in the weakened ability of information interaction between brain regions in ASD patients.
4.1. The Frontal Lobes
The frontal lobes participate in high levels of cognition and are extremely important for the development of social behavior. The orbitofrontal cortex, located on the ventral side of the frontal lobe, plays an important role in social emotional cognition and self-regulating behavior in social communication. This provides a basis for effective social communication. Some scholars (Mori et al., 2013; Jonker et al., 2015) believe that neurological functional deficits in the orbitofrontal cortex can cause ASD in patients with social function disorders. The posterior cingulate gyrus is an important part of emotional circuits (Miao et al., 2011), which participate in the processes of emotion and self-evaluation. Some scholars (Lai et al., 2010; Lee et al., 2013; Sophia et al., 2013) have found that the functions of the posterior cingulate gyrus are reduced in ASD patients. Loveland et al. (2008) found that the functional activity of brain regions (including the orbitofrontal cortex and the posterior cingulate gyrus) related to the judgment of facial expressions in emotional tasks was decreased in ASD patients. The opercular part of the inferior frontal gyrus is connected to the amygdala through the insula lobe so that when observing emotional expressions, an unconscious imitation of facial expression can promote the subjective emotional experience as well as the understanding of others' emotional intentions, which is closely related to emotional face processing (Fusar-Poli et al., 2009). Compared with the TD control group, the ASD group had abnormally fewer FC neural circuits in some involved brain regions, including the right orbital part of the superior frontal gyrus, bilateral orbital part of the middle frontal gyrus, left opercular part of the inferior frontal gyrus, and bilateral posterior cingulate gyrus. This indicates that there are FC deficits in these brain regions, which may be the cause of the cognitive dysfunction experienced by ASD patients and their inability to effectively implement social interaction behaviors. The supplementary motor area is related to the voluntary movement control of the trunk muscles. Huang et al. (2015) found that the FC deficit in the premotor cortex in ASD patients may be one of the mechanisms of habituate behaviors. The number of FC neural circuits involved in the left supplementary motor area was abnormally decreased in the ASD group. This finding indicates that the dysfunction of planning and executing movement in ASD patients may result in abnormal cognitive integration.
4.2. The Temporal Lobes
The temporal lobes process auditory information and participate in emotion and memory functions. The hippocampus is mainly responsible for the storage, conversion, and orientation of long-term memory, which is related to social episodic memory and self-awareness of memory (Yonelinas et al., 2001; Greenberg et al., 2005). In 1937, Papez found that the “parahippocampal gyrus—hippocampus—papillary body—anterior thalamic nucleus—cingulate gyrus—parahippocampal gyrus” formed the hippocampal circuit, which is related to emotion, learning, memory, and other higher functions. Rudie et al. (2012) showed that the hippocampus and the parahippocampal gyrus in ASD patients are strongly connected, and the habituate behaviors of ASD patients are related to FC of the parahippocampal gyrus. Deficts in the parahippocampal gyrus may cause the dysfunction in emotional and cognitive behaviors. Compared with the TD control group, the ASD group had abnormally fewer FC neural circuits in the left hippocampus and the right parahippocampal gyrus. This indicates that memory dysfunction in ASD patients may be a potential cause of social dysfunction. The Heschl gyrus is located on the superior surface of the temporal lobe. The researcher in Long et al. (2018) found that the Heschl gyrus is responsible for processing the complex spectrum of acoustic information, such as pitch, intensity, and timbre. This finding indicates that speech recognition disorder in ASD patients may result in cognitive dysfunctions.
4.3. The Parietal Lobes
The parietal lobes participate in memory processing, which is related to sensory and linguistic functions (Roth and Courtney, 2007; Johnson et al., 2011). The angular gyrus is the visual language center, which facilitates the understanding of the meaning of symbols and words. The left angular gyrus connects the visual, auditory, and somatosensory cortices, and various forms of information generate integrated language activities through the function of the angular gyrus (Joseph, 1993; Hofmann et al., 2008). Compared with that in the TD control group, the number of FC neural circuits in the left angular gyrus was abnormally decreased in the ASD group. This may lead to the cognitive dysfunction in ASD patients and the loss of the connection between what they view and what they hear.
4.4. The Subcortical Regions
The caudate nucleus (Sesack and Grace, 2010) is a key brain region in the reward network and plays an important role in behavior monitoring. Many scholars (Portmann et al., 2014; Staal, 2015; Qiu et al., 2016) believe that the caudate nucleus is related to the habituated behavior of ASD patients. Compared with that in the TD control group, the number of FC neural circuits in the right caudate nucleus was abnormally decreased in the ASD group. This indicates that the dysfunction of the caudate nucleus is associated with repetitive and stereotypical behaviors in ASD patients.
In addition, this study shows that most brain regions are unilateral rather than bilateral abnormalities, which indicates that there is brain laterality for FC in ASD patients.
In conclusion, the present study analyzes the FC neural circuits of ASD patients aged 6–13 years in the resting state and finds that there are significant differences in the topological properties of the FC neural circuits of ASD patients compared with those of the TD controls. These abilities of information transmission and the functional integration of brain networks are reduced. The number of FC neural circuits with different length ranges and the number of FC neural circuits in some involved brain regions are different, particularly those associated with social cognition, executive functioning, and memory. These abnormalities indicate that impairments in social cognition and social communication abilities in ASD patients are related to abnormal changes in FC neural circuits. The shortcomings of this study are the small sample size and cross-sectional data. In view of these shortcomings, the next step is to increase the sample size for research and analysis on individual longitudinal data. In addition, networks of different scales influence brain networks differently. Therefore, topographic maps of different scales should be combined for network construction to obtain more abundant local and global information. Therefore, more reliable biological markers for the diagnosis of ASD patients can be found, and more effective neuroimaging evidence for treatment will be provided.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: https://fcon_1000.projects.nitrc.org/indi/abide/.
Author Contributions
DL and SX were responsible for the study concept and design. WZ contributed to the acquisition of the data. XZ assisted with the analyses and interpretation of the findings. DL drafted the initial manuscript. DL, SX, XZ, and WZ critically reviewed the content and approved the final version of the manuscript for publication. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Natural Science Foundation of Shandong Province of China (ZR2020ZD25) and the Soft Science Research Project of Shandong Province of China (Grant No. 2020RKB01671).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2021.745671/full#supplementary-material
References
American-Psychiatric-Association (2013). Diagnostic and Statistical Manual of Mental Disorders (DSM-5Ⓡ). Arlington, VA: American Psychiatric Publishing. doi: 10.1176/appi.books.9780890425596
Carlsson, G. (2009). Topology and data. Bull. Am. Math. Soc. 46, 255–308. doi: 10.1090/S0273-0979-09-01249-X
Carlsson, G. (2020). Persistent homology and applied homotopy theory. arXiv:2004.00738v1. doi: 10.1201/9781351251624-8
Cassidy, B., Rae, C., and Solo, V. (2015). “Brain activity: conditional dissimilarity and persistent homology,” in 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI) (Brooklyn, NY), 1356–1359. doi: 10.1109/ISBI.2015.7164127
Cherkassky, V. L., Kana, R. K., Keller, T. A., and Just, M. A. (2006). Functional connectivity in a baseline resting-state network in autism. Neuroreport 17, 1687–1690. doi: 10.1097/01.wnr.0000239956.45448.4c
Cole, M. W., Ito, T., Bassett, D. S., and Schultz, D. H. (2016). Activity flow over resting-state networks shapes cognitive task activations. Nat. Neurosci. 19, 1718–1726. doi: 10.1038/nn.4406
Edelsbrunner, H., Letscher, D., and Zomorodian, A. (2002). Topological persistence and simplification. Discr. Comput. Geomet. 28, 511–533. doi: 10.1007/s00454-002-2885-2
Friston, K. J., Frith, C. D., Liddle, P. F., and Frackowiak, R. S. J. (1993). Functional connectivity: the principal-component analysis of large (pet) data sets. J. Cereb. Blood Flow Metab. 13, 5–14. doi: 10.1038/jcbfm.1993.4
Fusar-Poli, P., Placentino, A., Carletti, F., Landi, P., Allen, P., Surguladze, S., et al. (2009). Functional atlas of emotional faces processing: a voxel-based meta-analysis of 105 functional magnetic resonance imaging studies. J. Psychiatry Neurosci. 34, 418–432.
Greenberg, D. L., Rice, H. J., Cooper, J. J., Cabeza, R., Rubin, D. C., and Labar, K. S. (2005). Co-activation of the amygdala, hippocampus and inferior frontal gyrus during autobiographical memory retrieval. Neuropsychologia 43, 659–674. doi: 10.1016/j.neuropsychologia.2004.09.002
He, Y., Wang, J., Wang, L., Chen, Z. J., Yan, C., Yang, H., et al. (2009). Uncovering intrinsic modular organization of spontaneous brain activity in humans. PLoS ONE 4:e5226. doi: 10.1371/journal.pone.0005226
Hofmann, M. J., Herrmann, M. J., Dan, I., Obrig, H., Conrad, M., Kuchinke, L., et al. (2008). Differential activation of frontal and parietal regions during visual word recognition: an optical topography study. NeuroImage 40, 1340–1349. doi: 10.1016/j.neuroimage.2007.12.037
Huang, M., Zhang, Z., Gao, L., Wang, D., Chen, H., Xia, K., et al. (2015). Does premotor connectivity on resting fMRI reflect stereotyped behavior in ASD children. J. Pract. Radiol. 31, 1337–1340. doi: 10.3969/j.issn.1002-1671.2015.08.027
Johnson, R., Simon, E. J., Henkell, H., and Zhu, J. (2011). The role of episodic memory in controlled evaluative judgments about attitudes: an event-related potential study. Neuropsychologia 49, 945–960. doi: 10.1016/j.neuropsychologia.2011.01.028
Jonker, F. A., Jonker, C., Scheltens, P., and Scherder, E. J. A. (2015). The role of the orbitofrontal cortex in cognition and behavior. Rev. Neurosci. 26, 1–11. doi: 10.1515/revneuro-2014-0043
Joseph, R. (1993). The Naked Neuron: Evolution and the Languages of the Body and Brain. New York, NY: Springer. doi: 10.1007/978-1-4899-6008-5
Just, M. A., Cherkassky, V. L., Keller, T. A., and Minshew, N. J. (2004). Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity. Brain 127, 1811–1821. doi: 10.1093/brain/awh199
Kana, R. K., Libero, L. E., and Moore, M. S. (2011). Disrupted cortical connectivity theory as an explanatory model for autism spectrum disorders. Phys. Life Rev. 8, 410–437. doi: 10.1016/j.plrev.2011.10.001
Kanari, L., Dłotko, P., Scolamiero, M., Levi, R., Shillcock, J., Hess, K., and Markram, H. (2018). A topological representation of branching neuronal morphologies. Neuroinformatics 16, 3–13. doi: 10.1007/s12021-017-9341-1
Kennedy, D. P., and Adolphs, R. (2012). The social brain in psychiatric and neurological disorders. Trends Cogn. Sci. 16, 559–572. doi: 10.1016/j.tics.2012.09.006
Lai, M. C., Lombardo, M. V., Chakrabarti, B., Sadek, S. A., Pasco, G., Wheelwright, S. J., et al. (2010). A shift to randomness of brain oscillations in people with autism. Biol. Psychiatry 68, 1092–1099. doi: 10.1016/j.biopsych.2010.06.027
Lee, H., Chung, M. K., Kang, H., Choi, H., and Kim, Y. K. Lee, D. S. (2018). “Abnormal hole detection in brain connectivity by kernel density of persistence diagram and Hodge Laplacian,” in 2018 IEEE 15th International Symposium on Biomedical Imaging (ISBI 2018) (Washington, DC), 20–23. doi: 10.1109/ISBI.2018.8363514
Lee, M. H., Smyser, C. D., and Shimony, J. S. (2013). Resting-state fMRI: a review of methods and clinical applications. Am. J. Neuroradiol. 34, 1866–1872. doi: 10.3174/ajnr.A3263
Long, P., Wan, G., Roberts, M. T., and Corfas, G. (2018). Myelin development, plasticity, and pathology in the auditory system. Dev. Neurobiol. 78, 80–92. doi: 10.1002/dneu.22538
Loveland, K. A., Steinberg, J. L., Pearson, D. A., Mansour, R., and Reddoch, S. (2008). Judgments of auditory-visual affective congruence in adolescents with and without autism: a pilot study of a new task using fMRI. Percept. Motor Skills 107, 557–575. doi: 10.2466/pms.107.2.557-575
Martino, A. D., and Milham, M. (2017). Autism brain imaging data exchange II abide II. Available online at: http://fcon_1000.projects.nitrc.org/indi/abide/abide_II.html
Martino, A. D., Yan, C. G., Li, Q., Denio, E., Castellanos, F. X., Alaerts, K., et al. (2014). The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. Mol. Psychiatry 19, 659–667. doi: 10.1038/mp.2013.78
Miao, X., Wu, X., Li, R., Chen, K., and Yao, L. (2011). Altered connectivity pattern of hubs in default-mode network with Alzheimer's disease: an granger causality modeling approach. PLoS ONE 6:e25546. doi: 10.1371/journal.pone.0025546
Moran, J. S., Casey, C. G., Rutledge, T. F., Johnson, D. C., Lansdowne, C. B., Boyd, M. F., et al. (2014). Prevalence of Autism Spectrum Disorder Among Children Aged 8 Years - Autism and Developmental Disabilities Monitoring Network, 11 sites, United States, 2010. MMWR. Survllance summaries: Morbidity and mortality weekly report. Survllance summaries/CDC.
Mori, K., Toda, Y., Ito, H., Mori, T., Goji, A., Fujii, E., et al. (2013). A proton magnetic resonance spectroscopic study in autism spectrum disorders: amygdala and orbito-frontal cortex. Brain Dev. 35, 139–145. doi: 10.1016/j.braindev.2012.09.016
Ogawa, S., Lee, T. M., Kay, A. R., and Tank, D. W. (1990). Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc. Natl. Acad. Sci. U.S.A. 87, 9868–9872. doi: 10.1073/pnas.87.24.9868
Pereira, C. M. M., and de Mello, R. F. (2015). Persistent homology for time series and spatial data clustering. Expert Syst. Appl. 42, 6026–6038. doi: 10.1016/j.eswa.2015.04.010
Portmann, T., Yang, M., Mao, R., Panagiotakos, G., Ellegood, J., Dolen, G., et al. (2014). Behavioral abnormalities and circuit defects in the basal ganglia of a mouse model of 16p11.2 deletion syndrome. Cell Rep. 7, 1077–1092. doi: 10.1016/j.celrep.2014.03.036
Qiu, T., Chang, C., Li, Y., Qian, L., Xiao, C. Y., Xiao, T., et al. (2016). Two years changes in the development of caudate nucleus are involved in restricted repetitive behaviors in 2-5-year-old children with autism spectrum disorder. Dev. Cogn. Neurosci. 19, 137–143. doi: 10.1016/j.dcn.2016.02.010
Rdgaard, E. M., Jensen, K., Vergnes, J. N., Soulières, I., and Mottron, L. (2019). Temporal changes in effect sizes of studies comparing individuals with and with out autism: a meta-analysis. JAMA Psychiatry 76, 1124–1132. doi: 10.1001/jamapsychiatry.2019.1956
Roth, J. K., and Courtney, S. M. (2007). Neural system for updating object working memory from different sources: sensory stimuli or long-term memory. NeuroImage 38, 617–630. doi: 10.1016/j.neuroimage.2007.06.037
Rudie, J. D., Hernandez, L. M., Brown, J. A., Beck-Pancer, D., Colich, N. L., Gorrindo, P., et al. (2012). Autism-associated promoter variant in met impacts functional and structural brain networks. Neuron 75, 904–915. doi: 10.1016/j.neuron.2012.07.010
Salvador, R., Suckling, J., Coleman, M. R., Pickard, J. D., Menon, D., and Bullmore, E. (2005). Neurophysiological architecture of functional magnetic resonance images of human brain. Cereb. Cortex 15, 1332–1342. doi: 10.1093/cercor/bhi016
Sesack, S. R., and Grace, A. A. (2010). Cortico-basal ganglia reward network: microcircuitry. Neuropsychopharmacology 35, 27–47. doi: 10.1038/npp.2009.93
Sophia, M., Daniel, K., Samson, A. C., Valerie, K., Janusch, B., Michel, G., et al. (2013). Convergent findings of altered functional and structural brain connectivity in individuals with high functioning autism: a multimodal MRI study. PLoS ONE 8:e67329. doi: 10.1371/journal.pone.0067329
Staal, W. G. (2015). Autism, drd3 and repetitive and stereotyped behavior, an overview of the current knowledge. Eur. Neuropsychopharmacol. 25, 1421–1426. doi: 10.1016/j.euroneuro.2014.08.011
Supekar, K., Uddin, L. Q., Khouzam, A., Phillips, J., Gaillard, W. D., Kenworthy, L. E., et al. (2013). Brain hyperconnectivity in children with autism and its links to social deficits. Cell Rep. 5, 738–747. doi: 10.1016/j.celrep.2013.10.001
Tomasi, D., and Volkow, N. D. (2019). Association between brain activation and functional connectivity. Cereb. Cortex 29, 1984–1996. doi: 10.1093/cercor/bhy077
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage 15, 273–289. doi: 10.1006/nimg.2001.0978
Victor, S., Jean-Baptiste, P., Lindquist, M. A., Simpson, S. L., Dubois, B. F., Chung, M. K., et al. (2018). Connectivity in fMRI: blind spots and breakthroughs. IEEE Trans. Med. Imag. 37, 1537–1550. doi: 10.1109/TMI.2018.2831261
Wicker, B., Fonlupt, P., Hubert, B., Tardif, C., Gepner, B., and Deruelle, C. (2008). Abnormal cerebral effective connectivity during explicit emotional processing in adults with autism spectrum disorder. Soc. Cogn. Affect. Neurosci. 3, 135–143. doi: 10.1093/scan/nsn007
Xu, K., Liu, H., Li, H., Tang, Y., Womer, F., Jiang, X., et al. (2013). Amplitude of low-frequency fluctuations in bipolar disorder: a resting state fMRI study. J. Affect. Disord. 152, 237–242. doi: 10.1016/j.jad.2013.09.017
Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4
Yonelinas, A. P., Hopfinger, J. B., Buonocore, M. H., Kroll, N. E., and Baynes, K. (2001). Hippocampal, parahippocampal and occipital-temporal contributions to associative and item recognition memory: an fMRI study. Neuroreport 12, 359–363. doi: 10.1097/00001756-200102120-00035
Zhang, P. (2015). Attention bias to social threatening in autism spectrum disorders children (Master Thesis). School of Psychology, Northwest Normal University, Lanzhou, China.
Keywords: ASD, resting-state fMRI, FC neural circuits, persistent homology, VR complex
Citation: Liang D, Xia S, Zhang X and Zhang W (2021) Analysis of Brain Functional Connectivity Neural Circuits in Children With Autism Based on Persistent Homology. Front. Hum. Neurosci. 15:745671. doi: 10.3389/fnhum.2021.745671
Received: 22 July 2021; Accepted: 19 August 2021;
Published: 13 September 2021.
Edited by:
Zhen Yuan, University of Macau, ChinaReviewed by:
Jie Wu, Hebei Normal University, ChinaChengyuan Wu, Agency for Science, Technology and Research (A*STAR), Singapore
Copyright © 2021 Liang, Xia, Zhang and Zhang. 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: Shengxiang Xia, xias@sdjzu.edu.cn