- 1Department of Psychiatry and Psychotherapy, Jena University Hospital, Jena, Germany
- 2Center for Intervention and Research on Adaptive and Maladaptive Brain Circuits Underlying Mental Health (C-I-R-C), Jena, Germany
- 3Department of Psychosomatic Medicine and Psychotherapy, Jena University Hospital, Jena, Germany
- 4Clinical Affective Neuroimaging Laboratory (CANLAB), Magdeburg, Germany
- 5Department of Psychiatry and Psychotherapy, University of Tübingen, Tübingen, Germany
- 6Leibniz Institute for Neurobiology, Magdeburg, Germany
- 7Center for Behavioral Brain Sciences, Magdeburg, Germany
- 8German Center for Mental Health (DZPG), Jena, Germany
- 9Department of Psychiatry and Psychotherapy, Hannover Medical School, Hannover, Germany
- 10Department of Psychiatry and Psychotherapy, Jüdisches Krankenhaus Berlin—Berlin Jewish Hospital, Academic Teaching Hospital of the Charité—Universitätsmedizin Berlin, Berlin, Germany
Introduction: Autism spectrum disorder (ASD) is a neurodevelopmental disorder that persists into adulthood with both social and cognitive disturbances. Asperger's syndrome (AS) was a distinguished subcategory of autism in the DSM-IV-TR defined by specific symptoms including difficulties in social interactions, inflexible thinking patterns, and repetitive behaviour without any delay in language or cognitive development. Studying the functional brain organization of individuals with these specific symptoms may help to better understand Autism spectrum symptoms.
Methods: The aim of this study is therefore to investigate functional connectivity as well as functional network organization characteristics using graph-theory measures of the whole brain in male adults with AS compared to healthy controls (HC) (AS: n = 15, age range 21–55 (mean ± sd: 39.5 ± 11.6), HC: n = 15, age range 22–57 [mean ± sd: 33.5 ± 8.5]).
Results: No significant differences were found when comparing the region-by-region connectivity at the whole-brain level between the AS group and HC. However, measures of “transitivity,” which reflect local information processing and functional segregation, and “assortativity,” indicating network resilience, were reduced in the AS group compared to HC. On the other hand, global efficiency, which represents the overall effectiveness and speed of information transfer across the entire brain network, was increased in the AS group.
Discussion: Our findings suggest that individuals with AS may have alterations in the organization and functioning of brain networks, which could contribute to the distinctive cognitive and behavioural features associated with this condition. We suggest further research to explore the association between these altered functional patterns in brain networks and specific behavioral traits observed in individuals with AS, which could provide valuable insights into the underlying mechanisms of its symptomatology.
Introduction
Autism spectrum disorder (ASD) is a neurodevelopmental condition, which has been recognized as a spectrum disorder since its re-evaluation in the fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) in 2013 (1). ASD is characterized by a combination of (1) alterations in social interaction, (2) communicative behavior, and (3) repetitive patterns of behaviors and narrowly circumscribed interests. It additionally takes into account the subdivision of autistic individuals in terms of the severity of their impairments and the need for support in everyday life. Referring to the preceding classifications in the text-revised, fourth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV-TR) and the 10th edition of the International Classification of Diseases (ICD-10), autism was not only identified with all three core symptoms mentioned above but also included several subcategories including Asperger's syndrome (AS), Rett's syndrome, Kanner's syndrome, childhood disintegrative disorder, and atypical autism (2, 3). Following the DSM-IV-TR, AS comprised all the typical traits of the autistic triad but was typically not associated with delays in speech and cognitive development (2). It additionally included atypia of communicative behavior and special linguistic elements such as differing prosody or rhythm of speech (4). As a member of the group of pervasive development disorders, AS generally persists into adulthood.
Understanding the importance of AS lies in recognizing the unique strengths and challenges associated with it. Individuals with AS often exhibit remarkable focus and attention to detail impacting their fluid and general intelligence (5, 6). These strengths can contribute to their success in various domains, including academic pursuits as well as detailed and highly original problem-solving strategies (6). However, individuals with AS may also face difficulties in social interactions, especially in noisy conditions, which can be attributed to abnormal sensitivity to sensory stimuli, atypical eye contact, and altered speech perception (7–9). Recognizing the significance of Asperger's syndrome is essential for several reasons. First, it promotes a better understanding and acceptance of neurodiversity (10), acknowledging that individuals with AS have unique strengths and abilities that can contribute positively to society. Second, it allows for tailored support and accommodations to address the challenges individuals with AS may face in social interactions and sensory processing. Asperger's syndrome, as a subtype of autism spectrum disorder, holds significance in understanding the diverse range of strengths and challenges experienced by individuals on the autism spectrum. By recognizing the importance of AS and exploring the neural basis of this condition, we can base further research line, which is more adaptive to their unique brain functionality (11, 12).
Non-invasive neuroimaging studies have broadened our knowledge of both structural and functional neural alterations in ASD. The use of resting-state functional magnetic resonance imaging (rs-fMRI) enabled us to study the correspondence between autistic traits and brain functional connectome without any explicit task. In this context, age has been discussed to play an important role in ASD-related disturbances as significant changes in brain connectivity have already been hypothesized in late adolescence and early adulthood (13, 14). In addition, ASD is sometimes recognized very late or not at all in affected individuals, especially in high-functioning individuals, making research into the neural basis of this condition in adulthood particularly important (15). However, multiple studies have shown that altered functional connectivity in ASD is heterogeneous, making it difficult to reach a consensus on the nature and clinical relevance of these changes (16, 17). Despite these inconsistencies, studying functional connectivity at the resting state in individuals with Asperger's syndrome was shown to provide valuable insights into the neural networks associated with this specific subtype of ASD (14). Therefore, further studies are needed to fully elucidate the role of functional connectivity in autism and its subcategories like AS.
Graph theory has been extensively employed to investigate alterations in brain functional networks in both healthy individuals and psychiatric populations (18). It has provided valuable insights into the underlying mechanisms and potential biomarkers of various psychiatric conditions, such as schizophrenia (19), conduct disorder (20), and dissociative experiences (21). In neuroscience, graph theory can provide a framework to analyze the connectivity patterns between different brain regions, treating them as nodes and the connections between them as edges. These nodes and edges are organized into multilevel networks based on their specific functions, such as sensory or higher order cognitive functions. Ultimately, these networks come together to form the global brain network. The interactions between different regions within this network can be measured using graph theory properties, providing a deeper understanding of brain function and dysfunction. These properties include the clustering tendency of nodes within specific local brain networks and the neuronal path lengths between those networks (22). One way to quantify the organization of a complex network is to form the ratio of the “clustering coefficient” of networks to the average shortest path length within the vertices of a network, termed small-worldness. The small-worldness reflects the fact that nodes within or between networks can be reached with a few steps via their immediate neighbors, even if located in separate clusters. This feature indicates if the brain network can efficiently integrate simultaneous information of different quality across neural pathways or segregate them into modality-specific modules within a network (23). The average length of neural pathways within a network negatively correlates with the “global efficiency” of a network (18). Another metric for efficient brain processing involves the functioning of specialized and segregated clusters. In addition to global efficiency, networks must also consist of local centers with specialized functions, such as processing specific sensory stimuli. This property of networks can be assessed by measuring the “segregation.” As a component of segregation, “transitivity” refers to the clustering coefficient of brain networks and therefore reflects the interconnectivity of nodes within their direct neighborhood (18). “Assortativity” as a further metric of network segregation reflects the vulnerability of brain networks to harmful events or neuropathological processes, such as stroke, neurodegeneration, or traumatic microlesions (24, 25). High assortativity originates from robust connectivity between nodes in given networks corresponding with narrowly distributed links between major nodal hubs in those networks and is therefore interpreted as a marker for network resilience.
Previous studies have reported abnormality in the small-world configuration of the brain network in several psychiatric diseases (26–28), which resulted in the broad application of graph theory-based methods to study brain connectome in ASD (13, 26, 29–37). Because of the complex structure of its symptomatology including alterations in social interaction and communication, changes in ASD were discussed to originate from multiple parts of the brain, such as areas related to multimodal language processing or the Theory of Mind (9, 38). Furthermore, it has been suggested that because of this complex structure, changes in brain organization must predominantly occur at a global level rather than in terms of isolated functional connectivity abnormalities between or within particular brain regions (29). Age was again postulated to constitute another important covariate in the expression of ASD symptomatology, especially at critical developmental stages of the brain such as adolescence and early adulthood (13, 32, 39–42). In this context, research by Henry et al. described an initial decrease in global efficiency in childhood and adolescence with ASD, followed by an increase in early adulthood, while segregation of brain networks was found to increase in childhood and then decrease afterward (13). Confusingly, another study found that network segregation and global efficiency in ASD individuals were both diminished compared with healthy controls (39). Following the findings by Henry et al. (13), ASD can be characterized by a more extreme progression of the typically inverted U-shaped increasing curve of integrated connectome topology from childhood/adolescence into adulthood, which appears to be typical for the developing brain at this stage of life (40). From this perspective, ASD in adulthood would be mainly characterized by the presence of higher global integration (lower characteristic path length) as well as lower clustering and network resilience compared with healthy cohorts (41–43). Interestingly, Itahashi et al. (42) also stressed the need for further studies concentrating on changes in network organization from adolescence into adulthood.
Here, we hypothesized that complex network measures in adults with AS might be altered when compared to the control group. Following the existing literature on adult ASD, we expected functional connectivity alteration within sensory-motor regions as well as between frontal and parietal regions (44). We assumed that assortativity and transitivity of brain networks might be reduced in individuals with AS, whereas metrics quantifying the functional integration of the global connectome might be increased. We aimed for the whole-brain approach and its promising recent applications to psychiatric and autistic cohorts (13, 29–37). Measuring functional connectivity between each pair of regions throughout the whole brain as well as investigating graph theory measures were already subject to scientific research on autism and postulated to provide higher accuracy in diagnosing ASD in male adults (45, 46). For quantifying functional network organization, we investigated group differences in graph theory metrics, including functional segregation (transitivity), functional integration (global efficiency), and network resilience (assortativity coefficient) using the Brain Connectivity Toolbox (18).
Methods
Participants
Fifteen male participants diagnosed with Asperger's syndrome (AS) based on DSM-IV-TR (47), and fifteen healthy gender- and age-matched (only male) individuals without any psychiatric or neurologic disorders were recruited in this study as the healthy comparison group (HC). All participants were native German speakers. The final sample consisted of 15 male participants with AS with an age range of 21 to 55 (mean ± SD: 39.5 ± 11.6) and 15 male age-matched healthy controls from 22 to 57 (mean ± SD: 33.5 ± 8.5), and the group comparison using a two-sample t-test did not show any significantly different regarding age (t-value = 1.61, p-value = 0.12). The mean verbal IQ as assessed by the MWT-B—“Mehrfachwahl-Wortschatz-Intelligenztest” (48) was 30.20 for participants with AS and 31.43 for healthy controls and was not significantly different between the groups in the two-sample t-test (t-value = 0.73, p-value = 0.47). The mean of the autism spectrum quotient (AQ) for the AS group was 40.33 ± 5.45 (49), indicating a high expression of autistic traits. We explored DSM-IV-TR criteria for AS through a self-developed semi-structured interview (“Diagnostic interview: AS in adulthood”). This interview contained a general section focusing on medical anamnesis (somatic, psychiatric, and social histories, including childhood development) and continued with a special section involving AS that included the following items with respect to childhood and adulthood: social interaction and communication (e.g., friendships with/relationship to/interest in peers, and being a loner and suffering from loneliness); special interests (e.g., spending leisure time and interest in specific objects/topics); stereotypic behavior (e.g., rituals and reaction toward disturbances of rituals); and other characteristics (e.g., clumsiness and sensitivity toward noises/smells/tactile stimuli). It also included items and descriptions of all relevant criteria for the diagnosis of AS as defined in DSM-IV-TR. We confirmed the result of the interview for every AS subject by verifying the threshold value of the AQ. Additionally, eye contact, facial expression, prosody, and “mirroring” of affections and clumsiness during the interview were included in the assessment. The interview was always conducted by the same experienced psychiatrist and had a duration of ~90 min. For the whole duration of the interview, the investigator was blind to the research questions. The diagnosis was completed with information from personal interviews, gained by telephone or in written form, from observers in childhood and/or adulthood (e.g., partners, friends, and parents or siblings), and sometimes from incorporated school reports. Moreover, retrospective data on the development of speech were assessed. All DSM-IV-TR criteria had to be clearly fulfilled to confirm the diagnosis. An additional examination for axis-I co-morbidity was undertaken by using the German version of the Structured Clinical Interview for DSM-IV Axis I Disorders (SCID-I) (50). All probands of the AS group were diagnosed in adulthood and have not received any psychopharmacological treatment.
The present study was approved by the local Ethics Committee of the Hannover Medical School and has been performed following the latest version of the Declaration of Helsinki (51). The participants gave written informed consent before their participation and took part in the study for a small financial compensation for their travel expenses.
Scanner information
MR images were acquired on a 3-T Siemens Skyra Scanner (Siemens, Erlangen, Germany) equipped with a standard head coil. A total of 640 T2*-weighted volumes of the whole brain were collected in the AC-PC orientation with the following parameters: repetition time (TR) = 2,400 ms, echo time (TE) = 30 ms, flip angle = 80°, and a field of view (FOV) of 192 mm, matrix size = 64 × 64, 30 transversal slices, and voxel resolution of 3 × 3 × 3 mm3 with an interslice gap 0.33 mm. After the fMRI scan, a 3D high-resolution T1-weighted anatomical scan (MPRAGE-sequence, 192 sagittal slices, FOV = 256 mm, voxel resolution: 1 × 1 × 1 mm3, TR = 2.4 s, TE = 4.37 ms, flip angle = 7°) was recorded. The subject's head was fixed during the entire measurement to avoid head movements.
Image preprocessing
In the first step, DICOM data were converted to NIFTI files and organized in BIDS (Brain Imaging Data Structure) using dcm2bids (version: 2.1.4). Then, fMRIPrep (version:1.5.5) was used for preprocessing of the structural and functional MR data (52, 53), based on Nipype 1.4.0 (54, 55). The complete preprocessing pipeline could be found in the Supplementary material.
In brief, the functional images were slice-time corrected using 3dTshift from AFNI (56) and corrected for head motion. The T1 image was corrected for intensity non-uniformity, and the non-brain tissue was removed. Subsequently, the functional images were then co-registered to the T1w reference using bbregister (FreeSurfer), which implements boundary-based registration (57). Then, T1 images were normalized to the MNI space by applying non-linear registration using Advanced Normalization Tools (ANTs, version: 2.2.0). The Gaussian kernel of 6 mm FWHM (full-width half-maximum) was used for spatial smoothing with mcflirt (FSL 5.0.9) (58).
Head-motion parameters, i.e., transformation matrices, and six corresponding rotation and translation parameters were used to calculate the framewise displacement (FD) as implemented in Nipype (59) (https://fmriprep.readthedocs.io/en/latest/workflows.html).
Brain parcellation and identifying brain networks
To test our first hypothesis, AS and healthy control groups were compared for differences in functional connectivity at the network level. For this, we used the parcellation method by Schaefer et al. (61), which consists of 400 regions, to extract the time series of the cortical areas. These 400 regions were categorized into seven brain networks based on the classification by Thomas Yeo et al. (60), namely frontoparietal network (FPN), the default mode network (DMN), salience network (SN), limbic network (LN), dorsal attention network (DAN), sensorimotor network (SMN), and visual network (VN). Several studies have replicated the functional networks of this parcellation scheme based on extended resting-state data.
This homogeneous parcellation is composed of all cortical regions based on fMRI scans of 1,489 participants (61). The Fisher-Z transformation was applied for individual correlation coefficient matrices based on 400 parcels. Each value of these matrices is representative of functional connectivity between two nodes. Additionally, we also compared FC at the nodal level from all 400 ROIs (79,800 FCs) between AS patients and healthy controls.
Statistical models
To compare AS and control groups regarding the FC metrics at the network and nodal level, we used the two-sample t-test. To have comparable normalized variables, the dependent variables were inverse normal transformed before the two-sample t-test comparisons (62). To correct for multiple comparisons using the Bonferroni method, the alpha (0.05) was divided by 7, the number of comparisons at the network level (p-value < 0.007). However, due to a large number of comparisons at the pairwise FC level (79,800 FCs), we used the false discovery rate (FDR) correction to balance 214 between the likelihood of making a type I error and the probability of making a type II error (63).
Graph theoretical analysis
We computed both the clustering coefficient and transitivity to investigate the network segregation for different network densities. The clustering coefficient quantifies how much nodes in a network tend to cluster together. It considers the connectivity of both high- and low-degree nodes in the network. However, because the clustering coefficient is normalized at the nodal level, it gives low-degree nodes the same weight as high-degree nodes. This means that low-degree nodes are not ignored when computing the clustering coefficient, and their contribution to the overall network structure is considered. Therefore, networks with higher clustering coefficients might reflect locally efficient networks and suffer from fewer connected nodes. Meanwhile, transitivity is normalized collectively and less influenced by nodes with fewer connections (18). To measure the integration of the brain organization, we have chosen global efficiency as this graph feature is less influenced by long paths and it might be more reliable to show the information flow of a network (64). We also calculated the assortativity coefficient to examine the resilience of a network against “attacks” (65). The assortativity coefficient is a correlation coefficient between the degrees of all nodes on two ends of a link (66).
The correlation matrices of time series from 400 nodes based on the Schaefer parcellation scheme were calculated for each subject (61). These adjacency matrices were used to compute the graph theory features.
Graph theoretical measures were investigated at different network density thresholds, ensuring that all constructed networks had an equal number of edges for a given density (67). We used network densities ranging from 10 to 34% at intervals of 1%. Brain functional networks of this range show small-world properties, as described previously by Zhang et al. (68).
To perform group comparisons of graph theoretical measures, we employed a bootstrapping approach. First, we generated 1,000 bootstrap samples (with replacement) for each group and computed an average connectivity matrix for each sample. Second, we computed the clustering coefficient, transitivity, global efficiency, and assortativity from the average connectivity matrix of each bootstrap sample after thresholding and binarizing over the range of density thresholds mentioned above. Then, we calculated the area under the curve (AUC) for each graph feature. Finally, the comparisons between groups of AUC were made using a two-sample t-test and corrected for multiple comparisons using Bonferroni correction.
Results
Network-level and nodal-level functional connectivity
To compare the overall functional connectivity within each of the seven brain networks, we averaged the connections within regions of each defined functional network (DMN, FPN, SN, DAN, SMN, VN, and LN). There were no significant differences between the AS and control groups by comparing averaged FCs within seven main cortical networks.
We also compared all 79,800 FC between 400 regions. However, our analysis did not indicate any significant differences between the AS and control groups at the FDR-corrected p-level.
Graph theory properties
As presented in Figure 1, significant differences between groups were found in both small-world metrics and the assortativity across several network densities. Regarding transitivity as an indicator of brain network segregation, AS group showed significantly lower segregation than the control group (Figure 2A: AS < controls, t-value = −21.20, p-value < 0.0001, Bonferroni-corrected). Significant differences in the clustering coefficient confirmed the finding of lower segregation in the AS group (AS < controls, t-value = −29.85, p-value < 0.0001, Bonferroni-corrected). The AS group showed significantly lower network assortativity compared with healthy controls (Figure 2B: AS < controls, t-value = −15.19, p-value < 0.0001, Bonferroni-corrected).
Figure 1. Schaefer parcellation scheme with 400-area categorized in seven networks by Thomas Yeo et al. (60).
Figure 2. Comparisons of graph theory metrics between individuals with Asperger's syndrome (AS) and healthy controls. In the first row, the X-axis of each plot represents the density threshold ranging from 10 to 34%, and the Y-axis displays the values of graph theory measure, including transitivity (A), assortativity (B), and global efficiency (C). The Asperger's syndrome group (AS) is denoted by green, while the control group is indicated by orange. The second row shows boxplots of the area under the curve (AUC) for each graph measure [transitivity in (A), assortativity in (B), and global efficiency in (C)], which are compared between the two groups of AS (green) and controls (orange). The reported p-values were corrected for multiple comparisons using the Bonferroni approach (0.05/3).
Furthermore, the AS group exhibited higher network integration compared with the control group as revealed by higher global efficiency (Figure 2C: AS > controls, t-value = 2.33, p-value = 0.02, Bonferroni-corrected).
Discussion
In the present study, we investigated differences in pairwise FC at the whole-brain level as well as in the whole-brain functional organization in individuals with Asperger's syndrome (AS) and HCs. Our study revealed that individuals with AS show decreased network assortativity and transitivity compared with the HC group pointing out reduced activity in local information processing within specified brain regions and lower network resilience in AS. The global efficiency of whole-brain networks was significantly higher in participants with AS compared with HCs, which on the other hand reveals smaller shortest path lengths between distant nodes in AS. Comparing the averaged functional connectivity (FC) within the seven main functional brain networks, we did not find any significant differences between AS and HCs.
Previous studies on brain organization in AS have shown several inconsistencies and are partially competing with the current findings. Our finding of no difference in whole-brain FC in AS individuals compared with the control group is in line with a previous study by Tyszka et al. (46), which found no significant differences in whole-brain FC in adults with high-functioning autism (HFA) in comparison with HCs. However, the aforementioned study considered that mainly individuals in adolescence and with low functional abilities may exhibit FC alterations. Further studies testing this hypothesis demonstrated hyper- and hypoconnectivity in whole-brain analyses of ASD (14, 69, 70). In the case of AS, hyperconnectivity of subcortical and primary sensory networks was found to be associated with disturbed sensory perception in adults with AS (71). Social impairments, one of the core symptoms of AS, were also correlated with hypoconnectivity in insular regions in adolescents with AS (72). Our study provides evidence to support the hypothesis that in adulthood, the brain functional connectivity of individuals with AS might be similar to the functional connectivity of typically developing individuals. This convergence of functional connectivity has already been observed in other brain functions such as sensory processing, which is essential for understanding autistic symptomatology (73). Inconsistencies in previous findings might be also related to high sample heterogeneity regarding age, sex, and severity of symptoms (14, 74). Here, despite the relatively small sample size, our study can contribute important evidence, as our sample contained narrowly defined subjects with AS, exclusively male sex, and high cognitive functioning.
On the other hand, inconsistencies in studies of FC in AS samples could be also related to the complexity of AS symptomatology, which might affect the network organization of the brain. Keeping in mind, that the already mentioned core symptomatology of ASD most likely depends on the interplay of various, also sensory-related, brain networks (9, 38), we argue that the complexity of AS symptomatology could be better explained by relating them to alterations in the whole-brain connectome level analyzed by means of graph theory. Consequently, our results revealed that brain network transitivity and assortativity are globally reduced in the AS cohort, while global integration is increased. This could be interpreted as a shift toward a less efficient and therefore more randomized network configuration in adult AS, which has already been observed in previous connectome studies with other psychiatric cohorts (28). This hypothesis is furthermore in line with studies proving that lower path lengths typically occur in more randomized networks (75, 76). In further consistency with our findings, increased global efficiency was found for children with ASD in comparison with matched subjects with developmental delay (30). Moreover, increases in global integration in ASD were reported in studies with ASD samples in childhood and adolescence, showing that global efficiency in male subjects with ASD increased from adolescence into early adulthood, whereas transitivity of brain networks diminished in opposition to this (13).
Partially contradicting our findings are results reporting a decrease in clustering coefficients at the local network level in a sample with ASD individuals in childhood and adolescence (e.g., DMN) (41, 75), at least one of which was based on a highly heterogeneous sample in terms of subjects' verbal IQ, symptom severity, and global functioning (41). Interestingly, segregation deficits in children with ASD have been described for particular brain regions involved in social interaction (77) or emotional face processing (78). Here, we speculate that our finding of no significant alterations at the local network level in adult AS can be explained by the fact that clinically better-performing adults with AS or HFA might only dispose of significant deficits at the global level, whereas deficits in local network connectivity might be associated with earlier stages of brain development in ASD (13, 79). Moreover, our study is in partial contradiction with results that have postulated that significant parameters in ASD can be expected at the local level for the introduction of screening measures (34). However, this study included a sample with a relatively young subject age [age: 20.49 ± 6.16 (14–42)] and may therefore not adequately depict the neural constellation of the adults with AS or HFA.
The findings of our study, therefore, provide further evidence for a characteristic brain network configuration in adults with ASD, which has already been characterized in previous studies by increased global integration and lower clustering coefficients (41–43, 75). However, it is not possible to specify the exact symptoms that can be explained by the present findings as the study did not directly investigate the relationship between network configurations and specific symptoms. The identified network configuration of increased global integration and lower clustering coefficients in AS group may nevertheless provide insights into the underlying neural mechanisms that contribute to this specific subgroup of ASD, characterized by alterations in social communication, repetitive behaviors, and restricted interests without any delay in verbal or cognitive development (2, 3). Further research is needed to establish the relationship between network configurations and specific symptoms of ASD.
Limitations
In the first line, the sample size is small and diagnostics for Asperger's syndrome are already obsolete. The results of studies only deliver information on a subset of the current ASD population (high-functioning autism/Asperger's syndrome). Another important consideration is the potential impact of gender on brain network organization in ASD and to control the effect of gender, our study only included male participants, which has also been discussed as an independent factor in the context of altered brain functionality in autism (13, 14, 80–82). Additionally, we did not correlate the severity of the AS symptomatology to other clinical scores than the AQ. Therefore, the limitations of our investigation may potentially be the cause of the unexpected results in our study, such as the non-finding in the FC analyses.
Conclusion
In conclusion, this study showed that global FC in adult AS was significantly different from adult HC. As adults exhibited increased global integration and decreased network transitivity and assortativity. Current findings indicated more randomized brain organization in adult AS compared to HCs. We recommend two areas for future research. First, comparing brain organization measures in different groups of individuals stratified by clusters of ASD symptoms, gender, age, and severity of functional impairments could deepen our understanding of the symptomology of ASD and related disturbances in brain functional organization. Second, longitudinal studies on high-functioning individuals with ASD could help map changes in brain functional organization alongside changes in symptoms from childhood to adulthood.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The present study was approved by the Local Ethics Committee of the Hannover Medical School. The participants gave written informed consent before their participation and took part in the study.
Author contributions
F-AT and GS participated in the coordination of the study and performed the measurement. NJ, GW, and F-AT conceptualized the topic, scope of the research, and wrote the original draft of the manuscript. NJ and FC performed the data analysis. All authors reviewed the original draft and read and agreed to the published version of the manuscript.
Funding
We acknowledge the generous support provided by the German Research Foundation (Projekt-Nr. 512648189) and the Open Access Publication Fund of the Thüringer Universitäts- und Landesbibliothek Jena.
Conflict of interest
MW is a member of the following advisory boards and gave presentations to the following companies: Bayer AG, Germany; Boehringer Ingelheim, Germany; and Biologische HeilmittelHeel GmbH, Germany, conducted studies with institutional research support from HEEL and Janssen Pharmaceutical Research for a clinical trial (IIT) on ketamine in patients with MDD, unrelated to this investigation, and did not receive any financial compensation from the companies mentioned above.
The remaining 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.
References
1. American Psychiatric Association, American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders: DSM-5. 5th ed. Washington, DC: American Psychiatric Association (2013). p. 947.
2. American Psychiatric Association, American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders: DSM-IV-TR. 4th ed. Washington, DC: American Psychiatric Association (2000). p. 943.
3. World Health Organization. The ICD-10 Classification of Mental and Behavioural Disorders: Clinical Descriptions and Diagnostic Guidelines. Geneva: World Health Organization (1992). p. 362.
4. Shriberg LD, Paul R, McSweeny JL, Klin A, Cohen DJ, Volkmar FR. Speech and prosody characteristics of adolescents and adults with high-functioning Autism and Asperger syndrome. J Speech Lang Hear Res. (2001) 44:1097–115. doi: 10.1044/1092-4388(2001/087)
5. Soulières I, Dawson M, Gernsbacher MA, Mottron L. The level and nature of autistic intelligence II: what about Asperger syndrome? PLoS ONE. (2011) 6:e25372. doi: 10.1371/journal.pone.0025372
6. Pennisi P, Giallongo L, Milintenda G, Cannarozzo M. Autism, autistic traits and creativity: a systematic review and meta-analysis. Cogn Process. (2021) 22:1–36. doi: 10.1007/s10339-020-00992-6
7. Saalasti S, Tiippana K, Kätsyri J, Sams M. The effect of visual spatial attention on audiovisual speech perception in adults with Asperger syndrome. Exp Brain Res. (2011) 213:283–90. doi: 10.1007/s00221-011-2751-7
8. Pelphrey KA, Morris JP, McCarthy G. Neural basis of eye gaze processing deficits in autism. Brain. (2005) 128:1038–48. doi: 10.1093/brain/awh404
9. Tietze FA, Hundertmark L, Roy M, Zerr M, Sinke C, Wiswede D, et al. Auditory deficits in audiovisual speech perception in adult Asperger's syndrome: fMRI study. Front Psychol. (2019) 10:2286. doi: 10.3389/fpsyg.2019.02286
10. Smith O, Jones SC. ‘Coming out' with autism: identity in people with an Asperger's diagnosis after DSM-5. J Autism Dev Disord. (2020) 50:592–602. doi: 10.1007/s10803-019-04294-5
11. Hughes JA. Does the heterogeneity of autism undermine the neurodiversity paradigm? Bioethics. (2021) 35:47–60. doi: 10.1111/bioe.12780
12. Chown N, Robinson J, Beardon L, Downing J, Hughes L, Leatherland J, et al. Improving research about us, with us: a draft framework for inclusive autism research. Disabil Soc. (2017) 32:720–34. doi: 10.1080/09687599.2017.1320273
13. Henry TR, Dichter GS, Gates K. Age and gender effects on intrinsic connectivity in autism using functional integration and segregation. Biol Psychiatry. (2018) 3:414–22. doi: 10.1016/j.bpsc.2017.10.006
14. Hull JV, Dokovna LB, Jacokes ZJ, Torgerson CM, Irimia A, Van Horn JD. Resting-state functional connectivity in autism spectrum disorders: a review. Front Psychiatry. (2017) 7:205. doi: 10.3389/fpsyt.2016.00205
15. Lai MC, Baron-Cohen S. Identifying the lost generation of adults with autism spectrum conditions. Lancet Psychiatry. (2015) 2:1013–27. doi: 10.1016/S2215-0366(15)00277-1
16. Benkarim O, Paquola C, Park B, Hong SJ, Royer J, Vos de Wael R, et al. Connectivity alterations in autism reflect functional idiosyncrasy. Commun Biol. (2021) 4:1078. doi: 10.1038/s42003-021-02572-6
17. Kumar V, Garg R. Resting state functional connectivity alterations in individuals with autism spectrum disorders: a systematic review. Psychiatry Clin Psychol. (2021) 27:21261165. doi: 10.1101/2021.07.27.21261165
18. Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. (2010) 52:1059–69. doi: 10.1016/j.neuroimage.2009.10.003
19. Yuan L, Ma X, Li D, Ouyang L, Fan L, Li C, et al. Alteration of a brain network with stable and strong functional connections in subjects with schizophrenia. Schizophr. (2022) 8:91. doi: 10.1038/s41537-022-00305-0
20. Pu W, Luo Q, Jiang Y, Gao Y, Ming Q, Yao S. Alterations of brain functional architecture associated with psychopathic traits in male adolescents with conduct disorder. Sci Rep. (2017) 7:11349. doi: 10.1038/s41598-017-11775-z
21. Lotfinia S, Soorgi Z, Mertens Y, Daniels J. Structural and functional brain alterations in psychiatric patients with dissociative experiences: a systematic review of magnetic resonance imaging studies. J Psychiatr Res. (2020) 128:5–15. doi: 10.1016/j.jpsychires.2020.05.006
22. Stanley ML, Simpson SL, Dagenbach D, Lyday RG, Burdette JH, Laurienti PJ. Changes in brain network efficiency and working memory performance in aging. PLoS ONE. (2015) 10:e0123950. doi: 10.1371/journal.pone.0123950
23. Lord LD, Stevner AB, Deco G, Kringelbach ML. Understanding principles of integration and segregation using whole-brain computational connectomics: implications for neuropsychiatric disorders. Phil Trans R Soc A. (2017) 375:20160283. doi: 10.1098/rsta.2016.0283
24. Agosta F, Sala S, Valsasina P, Meani A, Canu E, Magnani G, et al. Brain network connectivity assessed using graph theory in frontotemporal dementia. Neurology. (2013) 81:134–43. doi: 10.1212/WNL.0b013e31829a33f8
25. Moreira da Silva N, Cowie CJA, Blamire AM, Forsyth R, Taylor PN. Investigating brain network changes and their association with cognitive recovery after traumatic brain injury: a longitudinal analysis. Front Neurol. (2020) 11:369. doi: 10.3389/fneur.2020.00369
26. Barttfeld P, Wicker B, Cukier S, Navarta S, Lew S, Sigman M, et al. Big-world network in ASD: dynamical connectivity analysis reflects a deficit in long-range connections and an excess of short-range connections. Neuropsychologia. (2011) 49:254–63. doi: 10.1016/j.neuropsychologia.2010.11.024
27. Brier MR, Thomas JB, Fagan AM, Hassenstab J, Holtzman DM, Benzinger TL, et al. Functional connectivity and graph theory in preclinical Alzheimer's disease. Neurobiol Aging. (2014) 35:757–68. doi: 10.1016/j.neurobiolaging.2013.10.081
28. Wagner G, de la Cruz F, Köhler S, Pereira F, Richard-Devantoy S, Turecki G, et al. Connectomics-based functional network alterations in both depressed patients with suicidal behavior and healthy relatives of suicide victims. Sci Rep. (2019) 9:14330. doi: 10.1038/s41598-019-50881-y
29. Chen L, Chen Y, Zheng H, Zhang B, Wang F, Fang J, et al. Changes in the topological organization of the default mode network in autism spectrum disorder. Brain Imaging Behav. (2021) 15:1058–67. doi: 10.1007/s11682-020-00312-8
30. Fang H, Wu Q, Li Y, Ren Y, Li C, Xiao X, et al. Structural networks in children with autism spectrum disorder with regression: a graph theory study. Behav Brain Res. (2020) 378:112262. doi: 10.1016/j.bbr.2019.112262
31. Kazeminejad A, Sotero RC. The importance of anti-correlations in graph theory based classification of autism spectrum disorder. Front Neurosci. (2020) 14:676. doi: 10.3389/fnins.2020.00676
32. Keown CL, Datko MC, Chen CP, Maximo JO, Jahedi A, Müller RA. Network organization is globally atypical in autism: a graph theory study of intrinsic functional connectivity. Biol Psychiatry. (2017) 2:66–75. doi: 10.1016/j.bpsc.2016.07.008
33. Roine U, Roine T, Salmi J, Nieminen-von Wendt T, Tani P, Leppämäki S, et al. Abnormal wiring of the connectome in adults with high-functioning autism spectrum disorder. Mol Autism. (2015) 6:65. doi: 10.1186/s13229-015-0058-4
34. Sadeghi M, Khosrowabadi R, Bakouie F, Mahdavi H, Eslahchi C, Pouretemad H. Screening of autism based on task-free fMRI using graph theoretical approach. Psychiatry Res. (2017) 263:48–56. doi: 10.1016/j.pscychresns.2017.02.004
35. Soma D, Hirosawa T, Hasegawa C, An K, Kameya M, Hino S, et al. Atypical resting state functional neural network in children with autism spectrum disorder: graph theory approach. Front Psychiatry. (2021) 12:790234. doi: 10.3389/fpsyt.2021.790234
36. Tolan E, Isik Z. Graph theory based classification of brain connectivity network for autism spectrum disorder. In:I Rojas, F Ortuño, , editors, Bioinformatics and Biomedical Engineering. Cham: Springer International Publishing (2018). p. 520–30.
37. Vagnetti R, Pino MC, Masedu F, Peretti S, Le Donne I, Rossi R, et al. Exploring the social cognition network in young adults with autism spectrum disorder using graph analysis. Brain Behav. (2020) 10:1524. doi: 10.1002/brb3.1524
38. Bolis D, Schilbach L. Observing and participating in social interactions: action perception and action control across the autistic spectrum. Dev Cogn Neurosci. (2018) 29:168–75. doi: 10.1016/j.dcn.2017.01.009
39. Bathelt J, Koolschijn PC, Geurts HM. Age-variant and age-invariant features of functional brain organization in middle-aged and older autistic adults. Mol Autism. (2020) 11:9. doi: 10.1186/s13229-020-0316-y
40. Collin G, van den Heuvel MP. The ontogeny of the human connectome: development and dynamic changes of brain connectivity across the life span. Neuroscientist. (2013) 19:616–28. doi: 10.1177/1073858413503712
41. Harlalka V, Naik S, Bapi RS, Vinod PK, Roy D. Integration and segregation in Autism Spectrum Disorders modulated by age, disease, and interaction: a graph theoretic study of intrinsic functional connectivity. Neuroscience. (2018) 2018:278846. doi: 10.1101/278846
42. Itahashi T, Yamada T, Watanabe H, Nakamura M, Jimbo D, Shioda S, et al. Altered network topologies and hub organization in adults with autism: a resting-state fMRI study. PLoS ONE. (2014) 9:e94115. doi: 10.1371/journal.pone.0094115
43. Harlalka V, Bapi RS, Vinod PK, Roy D. Age, disease, and their interaction effects on intrinsic connectivity of children and adolescents in autism spectrum disorder using functional connectomics. Brain Connect. (2018) 8:407–19. doi: 10.1089/brain.2018.0616
44. Holiga Š, Hipp JF, Chatham CH, Garces P, Spooren W, D'Ardhuy XL, et al. Patients with autism spectrum disorders display reproducible functional connectivity alterations. Sci Transl Med. (2019) 11:eaat9223. doi: 10.1126/scitranslmed.aat9223
45. Chaitra N, Vijaya PA, Deshpande G. Diagnostic prediction of autism spectrum disorder using complex network measures in a machine learning framework. Biomed Signal Process Control. (2020) 62:102099. doi: 10.1016/j.bspc.2020.102099
46. Tyszka JM, Kennedy DP, Paul LK, Adolphs R. Largely typical patterns of resting-state functional connectivity in high-functioning adults with autism. Cerebral Cortex. (2014) 24:1894–905. doi: 10.1093/cercor/bht040
47. American Psychiatric Association AP. Diagnostic and Statistical Manual of Mental Disorders (DSM-5®). American Psychiatric Pub (2013). p. 1520.
48. Lehrl S, Triebig G, Fischer B. Multiple choice vocabulary test MWT as a valid and short test to estimate premorbid intelligence. Acta Neurol Scand. (1995) 91:335–45. doi: 10.1111/j.1600-0404.1995.tb07018.x
49. Baron-Cohen S, Wheelwright S, Skinner R, Martin J, Clubley E. The autism-spectrum quotient (AQ): evidence from asperger syndrome/high-functioning autism, males and females, scientists and mathematicians. J Autism Dev Disord. (2001) 31:5–17. doi: 10.1023/A:1005653411471
50. Münster RD, Wittchen HU, Zaudig M, Fydrich T. SKID Strukturiertes Klinisches Interview für DSM-IV Achse I und II. Göttingen: Hogrefe (1997).
51. World Medical Association. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA. (2013) 20:2191–4. doi: 10.1001/jama.2013.281053
52. Esteban O, Markiewicz CJ, Blair RW, Moodie CA, Isik AI, Erramuzpe A, et al. fMRIPrep: a robust preprocessing pipeline for functional MRI. Nat Methods. (2019) 16:111–6. doi: 10.1038/s41592-018-0235-4
53. Esteban, O, Markiewicz, CJ, DuPre, E, Kent, JD, Ciric, R, Goncalves, M,, et al. fMRIPrep: A Robust Preprocessing Pipeline for Functional MRI. Zenodo (2020). Available online at: https://zenodo.org/record/3611399#.XpBTwVMzZ24 (accessed April 10, 2020).
54. Esteban O, Markiewicz, CJ, Burns, C, Johnson, H, Ziegler, E, Manhães-Savio, A,, et al. nipy/nipype: 1.4.0. Zenodo (2019). Available online at: https://zenodo.org/record/3588470#.XpBVEFMzZ24 (accessed April 10, 2020).
55. Gorgolewski K, Burns CD, Madison C, Clark D, Halchenko YO, Waskom ML, et al. Nipype: a flexible, lightweight and extensible neuroimaging data processing framework in python. Front Neuroinform. (2011) 5:13. doi: 10.3389/fninf.2011.00013
56. Cox RW, Hyde JS. Software tools for analysis and visualization of fMRI data. NMR Biomed. (1997) 10:171–8. doi: 10.1002/(SICI)1099-1492(199706/08)10:4/5<171::AID-NBM453>3.0.CO;2-L
57. Greve DN, Fischl B. Accurate and robust brain image alignment using boundary-based registration. Neuroimage. (2009) 48:63–72. doi: 10.1016/j.neuroimage.2009.06.060
58. Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage. (2002) 17:825–41. doi: 10.1006/nimg.2002.1132
59. Power JD, Mitra A, Laumann TO, Snyder AZ, Schlaggar BL, Petersen SE. Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage. (2014) 84:320–41. doi: 10.1016/j.neuroimage.2013.08.048
60. Thomas Yeo BT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol. (2011) 106:1125–65. doi: 10.1152/jn.00338.2011
61. Schaefer A, Kong R, Gordon EM, Laumann TO, Zuo XN, Holmes AJ, et al. Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI. Cereb Cortex. (2018) 28:3095–114. doi: 10.1093/cercor/bhx179
62. Bishara AJ, Hittner JB. Testing the significance of a correlation with nonnormal data: comparison of Pearson, Spearman, transformation, and resampling approaches. Psychol Methods. (2012) 17:399–417. doi: 10.1037/a0028087
63. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
64. Achard S, Bullmore E. Efficiency and cost of economical brain functional networks. PLoS Comput Biol. (2007) 3:e17. doi: 10.1371/journal.pcbi.0030017
65. Noldus R, Van Mieghem P. Assortativity in complex networks. Journal of Complex Networks. (2015) 3:507–42. doi: 10.1093/comnet/cnv005
66. Newman MEJ. Assortative mixing in networks. Phys Rev Lett. (2002) 89:208701. doi: 10.1103/PhysRevLett.89.208701
67. van den Heuvel MP, de Lange SC, Zalesky A, Seguin C, Yeo BTT, Schmidt R. Proportional thresholding in resting-state fMRI functional connectivity networks and consequences for patient-control connectome studies: issues and recommendations. Neuroimage. (2017) 152:437–49. doi: 10.1016/j.neuroimage.2017.02.005
68. Zhang J, Wang J, Wu Q, Kuang W, Huang X, He Y, et al. Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biol Psychiatry. (2011) 70:334–42. doi: 10.1016/j.biopsych.2011.05.018
69. Hahamy A, Behrmann M, Malach R. The idiosyncratic brain: distortion of spontaneous connectivity patterns in autism spectrum disorder. Nat Neurosci. (2015) 18:302–9. doi: 10.1038/nn.3919
70. Lee Y, Park B, James O, Kim SG, Park H. Autism spectrum disorder related functional connectivity changes in the language network in children, adolescents and adults. Front Hum Neurosci. (2017) 11:418. doi: 10.3389/fnhum.2017.00418
71. Cerliani L, Mennes M, Thomas RM, Di Martino A, Thioux M, Keysers C. Increased functional connectivity between subcortical and cortical resting-state networks in autism spectrum disorder. J Am Med Assoc Psychiatry. (2015) 72:767. doi: 10.1001/jamapsychiatry.2015.0101
72. Ebisch SJH, Gallese V, Willems RM, Mantini D, Groen WB, Romani GL, et al. Altered intrinsic functional connectivity of anterior and posterior insula regions in high-functioning participants with autism spectrum disorder. Hum Brain Mapp. (2011) 32:1013–28. doi: 10.1002/hbm.21085
73. Foxe JJ, Molholm S, Del Bene VA, Frey HP, Russo NN, Blanco D, et al. Severe multisensory speech integration deficits in high-functioning school-aged children with autism spectrum disorder (ASD) and their resolution during early adolescence. Cerebral Cortex. (2015) 25:298–312. doi: 10.1093/cercor/bht213
74. Feldman JI, Dunham K, Cassidy M, Wallace MT, Liu Y, Woynaroski TG. Audiovisual multisensory integration in individuals with autism spectrum disorder: a systematic review and meta-analysis. Neurosci Biobehav Rev. (2018) 95:220–34. doi: 10.1016/j.neubiorev.2018.09.020
75. Rudie JD, Brown JA, Beck-Pancer D, Hernandez LM, Dennis EL, Thompson PM, et al. Altered functional and structural brain network organization in autism. NeuroImage. (2013) 2:79–94. doi: 10.1016/j.nicl.2012.11.006
76. Sporns O. The non-random brain: efficiency, economy, and complex dynamics. Front Comput Neurosci. (2011) 2011:5. doi: 10.3389/fncom.2011.00005
77. Yerys BE, Gordon EM, Abrams DN, Satterthwaite TD, Weinblatt R, Jankowski KF, et al. Default mode network segregation and social deficits in autism spectrum disorder: evidence from non-medicated children. NeuroImage. (2015) 9:223–32. doi: 10.1016/j.nicl.2015.07.018
78. Rudie JD, Shehzad Z, Hernandez LM, Colich NL, Bookheimer SY, Iacoboni M, et al. Reduced functional integration and segregation of distributed neural systems underlying social and emotional information processing in autism spectrum disorders. Cerebral Cortex. (2012) 22:1025–37. doi: 10.1093/cercor/bhr171
79. Magiati I, Tay XW, Howlin P. Cognitive, language, social and behavioural outcomes in adults with autism spectrum disorders: a systematic review of longitudinal follow-up studies in adulthood. Clin Psychol Rev. (2014) 34:73–86. doi: 10.1016/j.cpr.2013.11.002
80. Floris DL, Lai MC, Nath T, Milham MP, Di Martino A. Network-specific sex differentiation of intrinsic brain function in males with autism. Mol Autism. (2018) 9:17. doi: 10.1186/s13229-018-0192-x
81. Lai MC, Lerch JP, Floris DL, Ruigrok ANV, Pohl A, Lombardo MV, et al. Imaging sex/gender and autism in the brain: etiological implications: imaging sex/gender and autism in the brain. J Neurosci Res. (2017) 95:380–97. doi: 10.1002/jnr.23948
Keywords: autism spectrum disorder (ASD), Asperger's syndrome (AS), graph theory, resting-state functional magnetic resonance imaging (rs-fMRI), autism spectrum traits
Citation: Javaheripour N, Wagner G, de la Cruz F, Walter M, Szycik GR and Tietze F-A (2023) Altered brain network organization in adults with Asperger's syndrome: decreased connectome transitivity and assortativity with increased global efficiency. Front. Psychiatry 14:1223147. doi: 10.3389/fpsyt.2023.1223147
Received: 16 May 2023; Accepted: 26 July 2023;
Published: 17 August 2023.
Edited by:
Xiaona Wang, Children's Hospital Affiliated to Zhengzhou University, ChinaReviewed by:
Zhengqin Zhao, Affiliated Hospital of Guangdong Medical University, ChinaBenedictor Alexander Nguchu, University of Science and Technology of China, China
Copyright © 2023 Javaheripour, Wagner, de la Cruz, Walter, Szycik and Tietze. 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: Nooshin Javaheripour, nooshin.javaheripour@med.uni-jena.de; Fabian-Alexander Tietze, fabian.tietze@jkb-online.de
†These authors have contributed equally to this work