Skip to main content
  • Research article
  • Open access
  • Published:

Brain functional networks in syndromic and non-syndromic autism: a graph theoretical study of EEG connectivity

Abstract

Background

Graph theory has been recently introduced to characterize complex brain networks, making it highly suitable to investigate altered connectivity in neurologic disorders. A current model proposes autism spectrum disorder (ASD) as a developmental disconnection syndrome, supported by converging evidence in both non-syndromic and syndromic ASD. However, the effects of abnormal connectivity on network properties have not been well studied, particularly in syndromic ASD. To close this gap, brain functional networks of electroencephalographic (EEG) connectivity were studied through graph measures in patients with Tuberous Sclerosis Complex (TSC), a disorder with a high prevalence of ASD, as well as in patients with non-syndromic ASD.

Methods

EEG data were collected from TSC patients with ASD (n = 14) and without ASD (n = 29), from patients with non-syndromic ASD (n = 16), and from controls (n = 46). First, EEG connectivity was characterized by the mean coherence, the ratio of inter- over intra-hemispheric coherence and the ratio of long- over short-range coherence. Next, graph measures of the functional networks were computed and a resilience analysis was conducted. To distinguish effects related to ASD from those related to TSC, a two-way analysis of covariance (ANCOVA) was applied, using age as a covariate.

Results

Analysis of network properties revealed differences specific to TSC and ASD, and these differences were very consistent across subgroups. In TSC, both with and without a concurrent diagnosis of ASD, mean coherence, global efficiency, and clustering coefficient were decreased and the average path length was increased. These findings indicate an altered network topology. In ASD, both with and without a concurrent diagnosis of TSC, decreased long- over short-range coherence and markedly increased network resilience were found.

Conclusions

The altered network topology in TSC represents a functional correlate of structural abnormalities and may play a role in the pathogenesis of neurological deficits. The increased resilience in ASD may reflect an excessively degenerate network with local overconnection and decreased functional specialization. This joint study of TSC and ASD networks provides a unique window to common neurobiological mechanisms in autism.

Peer Review reports

Background

Tuberous Sclerosis Complex (TSC) is a genetic neurocutaneous disorder, with highly variable, unpredictable and potentially devastating neurological outcome [1], and approximately 40% of these patients develop autism spectrum disorders (ASD) [2]. No conventional magnetic resonance imaging (MRI) biomarker can reliably predict intractable epilepsy, cognitive impairment or autism in this population [3]. Research has conventionally focused on non-syndromic ASD, but now consensus is emerging that single gene disorders with high penetrance of ASD (for example, TSC, Fragile X syndrome, Rett syndrome) can be used to understand better the cellular and circuitry bases of ASD [4–6]. Moreover, to advance the understanding of common neurobiological mechanisms in ASD, these should be present in subjects with ASD regardless of an underlying neurogenetic abnormality. For example, using diffusion tensor imaging (DTI), we have recently demonstrated abnormalities in structural connectivity of the corpus callosum of children with TSC and a co-morbid diagnosis of ASD, adding to a growing body of evidence of callosal microstructural deficits in subjects with ASD alone [3, 7–10].

Although such structural data from DTI provide insight into the architecture of interregional connections, to understand how neurophysiological function is supported by this architecture, functional networks should be analyzed as well [11]. Functional networks are implicated in cognitive functioning [12] and may form the physiological basis of information processing and mental representations [11]. They are made up by brief states of coordinated activity between physiological signals from neuronal aggregates in spatially distributed and specialized brain regions [13–16]. Functional connections form the building blocks of a functional network, and can be studied with neurophysiological techniques (for example, electroencephalography, EEG) and by neuroimaging (for example, functional MRI, fMRI).

Compared to fMRI, EEG has poor spatial resolution and is subject to volume conduction. However, it has a better signal-to-noise ratio and a significantly better temporal resolution. Moreover, EEG connectivity is directly related to neural activity, whereas fMRI is derived from the cerebral hemodynamic response to an increased metabolic demand [17], with a lag of 1 to 2 seconds from the neuronal activation. Recently, intermittent motion of the head during fMRI acquisition was shown to generate an artifactual reduction in long range connectivity and an increase in short range connectivity. This artifact may mask alterations in functional connectivity associated with autism, and complicate appropriate interpretation of functional connectivity MRI studies [18]. The best way to compensate for this artifact after the acquisition is completed remains unclear and the acquisition of MRI scans of children with autism without motion is a challenging task. On the contrary, artifact assessment is part of routine EEG interpretation by the clinical neurophysiologist, and common post-processing techniques allow for motion rejection or correction.

Thus, the main advantage of EEG is the high temporal resolution, allowing for direct characterization of higher frequency coordinated activity [19]. Recently, two data-driven analyses of EEG signals allowed for robust classification of subjects with (or at high risk for) autism and controls [20, 21]. Although these EEG studies reflect functional connectivity, they do not measure the complex network properties.

To characterize these complex networks with quantitative measures, graph analysis can be applied [11]. Graph theory has been recently introduced to characterize biological systems, particularly in the brain. Graph analysis of fMRI, magnetoencephalography (MEG) and EEG signals has revealed fundamental insights into the large-scale functional organization of the human brain in health and disease. Using EEG and MEG, syndrome-specific patterns of abnormal functional networks have been described in epilepsy [22], Alzheimer's disease [14] and in adult subjects with ASD [23–25].

For ASD, graph theoretical measures of brain networks are particularly well suited as ASD is, like Alzheimer's disease, considered a disconnection syndrome [26–30]. In disconnection syndromes, functional impairment is theoretically related to the disruption or abnormal integration of spatially distributed brain regions that would normally constitute a large-scale network subserving function [11, 28]. In ASD specifically, the developmental disconnection theory proposes a decreased long-range integration accompanied by increased local connectivity [7]. To synthesize the apparent inconsistencies of various long-range deficits or local surfeits in physical (DTI) and functional (fMRI) connectivity reported in ASD, a network approach may also be used [26, 27, 31–33].

We studied connectivity in ASD and the effects of abnormal connectivity on network properties. As the study of a homogeneous group of cooperative, high-functioning (young) adults precludes generalization of findings to the entire autism spectrum, a two-way study population was chosen: patients with and without TSC, and patients with and without ASD. To include the early developmental period of accelerated brain growth, during which autism symptoms become apparent [7] and secondary, maladaptive developmental changes have not yet occurred [34], a wide age-range was included in our study.

We hypothesized that (micro-) structural deficits in connection in TSC and ASD affect functional network properties, quantifiable by neurobiologically meaningful graph measures [35] of conventional EEG coherence. In particular, we hypothesized a widespread disconnectivity in TSC based on structural imaging findings, and decreased long-range and increased short-range connectivity in ASD, in agreement with the current model of developmental disconnection.

Methods

Subjects

TSC patients were identified through the Boston Children's Hospital Multidisciplinary Tuberous Sclerosis Program and were diagnosed with definite TSC based on clinical criteria described by the Tuberous Sclerosis Consensus Conference [36]. All patients with TSC were neurologically examined, and clinical data were obtained during office visits and from review of medical records. Genetic confirmatory testing included TSC1 and TSC2 gene sequencing and micro-deletion analysis at Athena Diagnostics (Worcester, MA, USA) or Boston University School of Medicine Center for Human Genetics (Boston, MA, USA). The ASD diagnoses were based on clinical assessment by a board-certified pediatric neurologist (MS and SSJ) using the Diagnostic and Statistical Manual (DSM-IV-TR), supplemented in most with the Autism Diagnostic Observation Schedule (ADOS) [37] by clinical- or research-ADOS certified specialists.

Autistic subjects without TSC (non-syndromic ASD group) were recruited from the Early Childhood Partial Hospitalization Program (ECPHP), an intensive, multidisciplinary, and highly specialized intervention program for children, two- to five-years old, with ASD via the Center for Autism Research and Treatment, Semel Institute, University of California, Los Angeles, Los Angeles, CA. ASD diagnoses were made as described for the TSC population.

Controls were selected from the general neurology clinic at Boston Children's Hospital in 2010, and were considered when they would have an EEG prompted by a single clinical event of moderate-to-low suspicion for epilepsy (for example, syncope, tics, behavioral outbursts, headache, and prominent startle). Only those subjects were included who had normal neurological development for age, a normal physical examination, a normal EEG both during wakefulness and sleep and a clinical follow-up of at least one month to confirm the trivial nature of the EEG referral. Fifteen controls had a normal imaging study, others were not imaged. Subject recruitment, data collection, retrieval and analysis were conducted with informed consent for the participation of children in the study by the parents when appropriate (for example, waived for use of retrospective EEG data), using protocols approved by the Institutional Review Boards from Boston Children's Hospital and the Semel Institute, University of California, Los Angeles.

The study populations are represented in (Figure 1 A).

Figure 1
figure 1

Group structure and functional networks. (A) In the two-way representation of our population, group status is defined by the subject being diagnosed with Tuberous Sclerosis Complex (TSC) or not and with Autism Spectrum Disorder (ASD) or not. This structure allows an independent attribution of effects specific to TSC and ASD. (B) Electrode locations from the international 10-20 system of electrode placement are used as nodes in the network. (C) Illustrations of the functional networks of a control subject, a TSC patient, a non-syndromic autistic patient, and a TSC patient diagnosed with autism. Colors on the connection terminations correspond to the colors in (B).

Previous literature has implicated an important role of the corpus callosum in ASD [8–10]. For illustrative purposes, and to ensure our functional connectivity measure was representative of callosal integrity, the effects of an absent or severely malformed corpus callosum (ACC) on coherence were assessed. These subjects were retrieved through a search of radiological reports and confirmed by review of the images and electronic medical records.

EEG recording and artifact rejection

In Boston, through review of electronic medical records, digital EEG data were identified and retrieved from the archives. Both routine electroencephalographic data and inpatient data from long-term monitoring with video-EEG were used, utilizing the 10-20 International System of electrode placement (Figure 1 B). If multiple recordings were present, a single record was chosen based on proximity to the acquisition of imaging data, for future correlation of imaging and EEG findings. If data were of insufficient quality, the next closest EEG was chosen. EEGs were recorded on Biologic recording systems, 256 to 512 Hz sampling rate, 1 to 100 Hz bandpass, or on a Natus Neuroworks® EEG system, 200 Hz sampling rate, 0.1 to 100 Hz bandpass. Data in Los Angeles were collected using a 128 Hydrocel Geodesic Sensor Net System (EGI®, Inc). Data were collected and recorded using NetAmps Amplifiers and NetStation software, sampled at 250 Hz, and digitized with a National Instruments Board (12 bit).

All raw data were imported, pruned, notch filtered at 60 Hz and, if necessary, spatially down-sampled to the standard clinical 19 electrodes (Figure 1B). An average reference was created using the BESA® Research 3.5 software package. Next, data were imported into EEGlab for band-pass filtering (FIR filter, 1-70 Hz), rejection of artifact-ridden epochs and selection of awake task-free data, with a minimum of two minutes. Epochs with evidence of muscle artifact were, where possible, rejected. Independent Component Analysis (ICA) was used for semi-automated artifact rejection of eye blinks and lateral eye movements, according to previously described methods [38]. The average reference was used for calculation of connectivity.

After artifact removal, for each subject several segments of continuous EEG signal were available for analysis. While these segments varied in their number and length, no group difference was observed in the mean, minimum and maximum length of these segments (P > 0.15 for both TSC and ASD conditions). The total duration of EEG data analyzed was higher in TSC subjects (P < 0.01) with an average time of 646 seconds compared to 439 seconds for non-TSC subjects. No difference related to ASD was observed in the total duration of EEG data analyzed. This difference in total length was accounted for in our definition of the connectivity measure.

Connectivity measure

The connections between brain regions that make up functional networks are measures of linear or non-linear statistical interdependence between two time-series [39, 40]. Coherence is a measure of the stability of phase correlation over time and is sensitive to both changes in power and phase relationships, although the former is typically negligible [17]. High coherence values between two signals are taken as a measure of strong connectivity between the responsible brain regions [41, 42]. Advantages of this measure include ample experience across the literature, the intuitive intelligibility by clinician-scientists and the description of statistically consistent, recurrent connections over a longer period of time. Drawbacks include the neglect of shorter interactions in the time domain and of non-linear relations, the assumption of stationarity of the signal, and sensitivity to volume conductance through skull, scalp and cerebrospinal fluid [14, 17, 42].

Combining the data from the different segments of continuous EEG recording by concatenation would span transitions which translate into artifactual high frequency content of the power spectrum. Thus, coherence was calculated for each segment individually. Our connectivity measure was obtained by computing the average of these coherences weighted by segment length. This method weighs longer segments more, and gives negligible weight to short segments. Specifically, let S i (t) be the signal in the i-th segment, L i be its length, N be the number of segments, Coh(S,f) be the coherence of signal S at frequency f and Ï• be the frequency band of interest, the connectivity measure read:

C = ∫ ϕ ∑ i = 1 N L i Coh ( S i , f ) ∑ i = 1 N L i d f .

Group comparison of coherences

To illustrate the validity of pre-processing methods and coherence calculations, data were analyzed on a sample of 16 patients with an ACC. The corpus callosum is the largest inter-hemispheric white matter pathway, critical to direct long-range information transfer between homotopic cortical regions and is implicated in autism [3, 8–10]. We calculated the ratio of the mean coherence of all corresponding inter-hemispheric electrode pairs over all non-midline intrahemispheric electrode pairs. As anticipated [43, 44], decreased inter-hemispheric coherence was found in the group with an abnormal or ACC (one-tailed two-sample t-test: P < .006 in all three bands).

For comparison of long- versus short-range coherences, neighboring electrode pairs were ignored because of excessive volume conductance [42, 45]. Short-distance mean coherence was calculated from all intra- and inter-hemispheric electrode pairs not immediately adjacent, that is, with a Euclidean distance of two. Long-distance pairs were defined as a Euclidean distance of three or more on the grid (Figure 1 B), that is, 75% or more of the maximum distance between aligned electrodes.

For analysis of TSC and ASD populations, we controlled for baseline coherences and volume conduction by comparison to healthy subjects and for maturational changes by including age as a covariate into the regression model.

The theta band (4 to 8 Hz) and the lower- and upper alpha bands (8 to 10 and 10 to 12 Hz, respectively) were chosen on the basis of previous findings in disconnection syndromes (for example, Alzheimer's disease [14] and autism (for example, [45, 46])), the higher power and signal-to-noise ratio in these bands and the increased susceptibility of beta- and gamma-bands to contamination by muscle artifact in routine clinical EEG. This also allowed for limiting the number of statistical analyses.

In addition, graph analysis allows the avoidance of the multiple comparisons typically needed for group analyses at the connection level (19 electrodes have 171 possible connections in each subject). The actual correction for the comparisons of the few graph measures used would require knowing the correlation between these measures.

Graph analysis

Mathematically, networks are represented by graphs, which consist of nodes connected by edges. A graph based on a connection measure without directionality (for example, coherence) is referred to as undirected. Graphs can also be weighted or unweighted. In an unweighted graph, edges represent the presence or absence of a connection between two nodes regardless of its strength. By contrast, weighted graphs also encode the strength of the connections within the edges. For this study, an undirected weighted graph was built using the 19 electrodes as nodes and inter-electrode coherence values as edges (Figure 1 B-C). Edge strength can be mapped to a functional distance by applying a function f to it [35]. We choose f to be the negative logarithm as it will associate a functional distance of 0 to time series in perfect synchronization (coherence equal to 1) and an infinite functional distance to incoherent time series (coherence equal to 0). Functional distances then allow the definition of functional path lengths being the sum of the functional distances along a particular path [35].

Graphs can be characterized by various global measures. It is not yet established which measures are most appropriate for the analysis of brain networks [11]. Three important ones are the characteristic path length, clustering coefficient, and global efficiency [14, 35]. The characteristic path length is the average length of the shortest paths that must be traversed to go from one node to another. The clustering coefficient indicates the likelihood that two nodes strongly connected to a third node are also strongly connected to each other, forming a strongly connected triangular cluster [11]. As such, the clustering coefficient is a measure of the network segregation. The global efficiency is the average of the inverse path lengths. As a result, the global efficiency is primarily driven by shorter paths (stronger connections) while characteristic path length is primarily driven by longer paths (weaker connections). In particular, the characteristic path length of a disconnected network is infinite while its global efficiency is finite. Both characteristic path lengths and global efficiency are measures of network integration. A high clustering coefficient and a low average path length form a network with 'small-world' characteristics. Small-world architecture suggests a network with connections that are neither regular nor random, and is found ubiquitously in natural and technological systems [11, 35].

Additional file 1 provides an accessible introduction to graph theory and our network measures used, as well as examples of brain and airline networks.

Another interesting property of networks is its resilience to the removal of random or highly connected nodes, known as Random Failure and Targeted Attack, respectively [43, 47]. In technological networks, the resilience is typically enforced by structurally replicating the nodes, inducing a redundancy in the network. In biological systems (and in the brain in particular), nodes typically cannot be replicated and resilience may indicate that structurally different components can perform similar functions, known as functional degeneracy. Thus, while a resilient functional network may reflect the ability to preserve system function in neuropathological conditions [14, 43, 47], an excess of degeneracy indicates a decreased functional specialization [48–50]. To measure resilience, attacks and failures are simulated by removing nodes and their connections from the graph. The global efficiency is computed for the resulting damaged network and compared to its initial value. Global efficiency is chosen to investigate resilience, as suggested in [35].

Additional file 2 contains an entry-level description of network resilience, and the main methods of assessing resilience through the two modes of network attack. Again examples are provided for airline networks and brain networks.

All graph measures were computed using the NetworkX toolbox in Python [51] except for the global efficiency which was developed in-house.

Statistical analysis

To distinguish the influence of ASD from that of TSC, a two-way ANCOVA was applied. This statistical model can assess effects specific to TSC and specific to ASD, and allows the inclusion of age as a covariate. However, it cannot account for differences in the group that are not additive. For example, if some measure is larger in both the ASD and TSC group but the effects do not add up in the TSC with ASD subgroup, then the group differences may not be shown by the model. Conversely, it cannot account for situations in which a group difference is attributable to a single subgroup.

Our unique two-way study population structure is somewhat akin to a repeated analysis, since, for each hypothesis tested, we have two subgroups to study. If the hypothesis is consistent for both subgroups (that is, both ASD with and without TSC, or both TSC with and without ASD), the finding is more intrinsic to ASD (or TSC).

The subject's age was used as a covariate, given the maturational changes in both EEG coherence [52, 53] and in graph measures of brain networks [54, 55]. The corresponding generalized linear model for each of the measures y is as follows:

y = y ¯ + β A S D A S D + β T S C T S C + β a g e a g e

where y ¯ is the baseline value, ASD and TSC are binary group variables indicating the presence or absence of ASD and TSC, respectively, and age is the subject's age in years. The two-way ANCOVA then allows us to assess whether β A S D and/or β T S C is significantly different from zero, indicating the influence of ASD and/or TSC on the measured properties.

Statistical and graph analyses were done with in-house developed software and standard issue toolbox on a MatLab platform (2009a, MatLab Inc., Natick, MA).

Results

All results of statistical tests are displayed in Table 1.

Table 1 P-values of the differences associated with autism spectrum disorder (ASD) and tuberous sclerosis complex (TSC).

Demographic data

A total of 43 subjects with TSC were included (27 male, mean age 6.9 years, range 0.7 to 25.6) and 46 age-matched control subjects (19 male, mean age 7.1 years, range 0.08 to 17.4). Fourteen TSC subjects were diagnosed with ASD (9 male, mean age 9.3 years, range 1.0 to 25.6) and 29 were not (17 male, mean age 6.0 years, range 0.7 to 23.4). Sixteen subjects with non-syndromic autism were included (12 male, mean age 4.1 years, range 2.2 to 5.5). Using Fisher's exact test (binary variables) and student t-test (continuous variables), no group differences were found in gender and age between all TSC subjects and controls (P = 0.45 and 0.06, respectively). Age of TSC subjects with and without autism did not differ from controls (P = 0.48 and 0.16, respectively). No age difference was found between all ASD subjects and controls (P = 0.29), but there was a slight male predominance (P = 0.02). Non-syndromic ASD subjects were younger than controls (mean 4.1 years +/- 1.1 versus 7.9 years +/- 5.6). Age differences were controlled for all four groups in all subsequent analyses, through the incorporation of age as a covariate into the ANCOVA model.

ASD was not associated with TSC1 or TSC2 mutations (P = 1.0). In all patients, there was no association between significant cognitive impairment (clinical assessment or, if available, full scale intelligence quotient < 70) and ASD (P = 0.15). There was no difference in the prevalence of significant cognitive impairment between patients with ASD alone and patients with ASD and TSC (P = 1.0). In TSC patients, there was no association between ASD and epilepsy, or ASD and infantile spasms (P = 1.0 and 0.19, respectively), perhaps reflecting an inclusion bias of those patients who underwent EEG recordings.

Coherence measures

The age-related increase in mean coherence (data not shown) results from developmental changes in brain maturation [52, 53]. On a network level, it represents increasing integration and decreasing segregation of structural and functional network hubs found by DTI and fMRI studies [54–56].

In the TSC group, mean coherence between all electrode pairs was significantly decreased in the lower alpha band. This decrease indicates a significant global underconnectivity specific to TSC, corrected for age and regardless of the presence of ASD. In the ASD group, no difference in mean coherence was observed (Figure 2 A). However, note that the mean coherence does reflect the distribution of long- and short-range connections (Figure 2 C).

Figure 2
figure 2

Connectivity by measures of conventional coherence and network topology. (A-C) Conventional coherence measures: (A) Mean coherence over all pairs of electrodes shows a significantly global under-connectivity for the Tuberous Sclerosis Complex (TSC) group regardless of the presence of Autism Spectrum Disorder (ASD). (B) Inter- versus intra-hemispheric coherence ratio. The significantly smaller value for patients with an absent corpus callosum (ACC) illustrates the validity of coherence as a connectivity measure. (C) Long- over short-range connectivity ratio is significantly smaller in the ASD group. As the mean coherence in (A) is not altered, this indicates a short-range overconnectivity and long-range underconnectivity in patients with ASD, evident in both subgroups (in ASD related to TSC and in ASD alone). (D-F) Network topology measures: TSC is characterized by a higher average path length (D) and lower clustering coefficient (E) and global efficiency (F). This departure from the small-world network topology results in a less functionally integrated and segregated network. No ASD group effect was found for these network topology measures. Only the lower alpha band is shown, see Table 1 for details. All measures are corrected for age. * P < 0.05, ** P < 0.01, *** P < 0.001.

For TSC, despite our prior report on microstructural deficits of the corpus callosum [3], there was no difference in the ratio of inter-hemispheric over intra-hemispheric coherence (Figure 2 B).

For ASD, with the exception of a group effect in the theta band, this ratio was unchanged as well. To illustrate coherence as a measure of connectivity, this ratio indeed showed a reduction of inter-hemispheric connections in the group with an ACC [43, 44].

In the TSC group, the ratio of long- over short-distance coherence trended lower but did not reach significance (Figure 2 C).

In the ASD group, this ratio was significantly and consistently decreased over all examined frequency bands. As the mean coherence (Figure 2 A) was not altered, the decreased ratio indicates a local overconnectivity accompanied by a proportional long-range underconnectivity in patients with ASD. This pattern was evident in both subgroups, that is, in both ASD with and without TSC.

Graph measures

In patients with TSC, we found both a longer average path length (weak, long connections are weaker) and a decreased global efficiency (strong, short connections are weaker), indicating less integration through both short and long network paths (Figure 2 D-F). The clustering coefficient was decreased, indicating a decreased local connectedness in the graph. Together, the increased path length and the decreased clustering coefficient represent a network that departs from small-world topology in TSC, independent of a co-morbid diagnosis of ASD.

In patients with ASD, regardless of the presence of TSC, no significant group difference was found for the three topological measures. Since in ASD the ratio of long- over short-range connectivity is significantly lower, the absence of topological differences in this population suggests that functional networks are altered while maintaining an unaltered distribution of connection strengths (Figure 3).

Figure 3
figure 3

Distribution of all connection strengths (coherences between electrode pairs). In Tuberous Sclerosis Complex (TSC), the distribution is left-skewed from a relative increase of lower coherence values. In Autism Spectrum Disorder (ASD), the distribution of connection strengths is similar to controls. This indicates that the alteration of the functional network in ASD stems from the allocation of similar connection strengths in a different electrode pairing scheme. Only the lower alpha band is shown.

Resilience measures

For TSC, there was no group effect for either the targeted attack or the random failure in all three spectral bands. The decreased mean coherence in TSC (Figure 2 A) does not affect resilience measures, as these measures reflect a percentage change relative to the baseline global efficiency of a network.

For ASD, with the targeted attack, a significantly decreased decline of the global efficiency was found. This group effect was not present with the random failure. In ASD, regardless of the presence of TSC, this significantly increased resilience to targeted attack was found in all frequency bands (Figure 4). As increased resilience could be related to an altered organization of hubs (highly connected nodes) in the network, we calculated the connection degree of the three highest connected nodes, compared to the sum of the degree of all nodes. In ASD, this normalized degree was significantly decreased across all bands (P < 0.005 for the first hub; P < 0.001 for second hub except in the upper alpha band; P < 0.01 for the third hub except in the upper alpha band).

Figure 4
figure 4

Resilience to targeted attacks. Autism Spectrum Disorder (ASD) is characterized by significantly higher network resilience to targeted attacks, indicating a higher degree of functional degeneracy (see discussion for details). This difference is consistent in both ASD subgroups (with and without TSC). Only the lower alpha band is shown. A sample bar plot of the global efficiency of the network after removal of three nodes is inserted (top right).

Subgroup analysis

Results in Figures 2, 3, 4 and 5 and Table 1 suggest that the observed differences are specific to a condition (TSC or ASD) rather than a subgroup. To validate this observation, we performed post-hoc t-tests on the differences between the subgroups. The contributions of TSC with and without ASD to the findings related to TSC (mean coherence and graph measures) were not significantly different. Similarly, the contributions of ASD with and without TSC to the findings related to ASD (long- over short-range connectivity and resilience to targeted attack) were not significantly different across all bands (except in the theta band for 1 and 2 nodes removed, P = 0.02 and 0.04, respectively). The difference in sub-group age between idiopathic ASD and controls did not in turn lead to the identification of a subgroup difference on these measures. There were no differences in the findings between the groups and the subgroups on these measures.

Figure 5
figure 5

Summary of network analysis. In summary, functional networks in Autism Spectrum Disorder (ASD) are characterized by a high resilience to targeted attacks, and this is consistent for both ASD alone and ASD related to Tuberous Sclerosis Complex (TSC). Functional networks in the TSC group have an altered topology with a decreased clustering coefficient and increased average path length, regardless of a co-morbid diagnosis of ASD. The consistency of our findings is highlighted by the alignment of the data-cloud on the Y-axis (increased resilience for ASD) and on the X-axis (altered network topology for TSC). Groups are displayed with a span of one standard deviation.

For clarity, the findings characterizing TSC by an altered network topology and ASD by an increased resilience are summarized in Figure 5.

Discussion

We analyzed functional connectivity through EEG coherence in a large sample of children with TSC with and without ASD. Incorporating subjects with ASD alone allowed us to study connectivity abnormalities common to autism, regardless of etiology. In TSC, a pattern of global underconnectivity and altered network topology was found. ASD was characterized by a decreased long- over short-range connectivity and a markedly increased resilience to targeted attack.

Coherence measures

In TSC, mean coherence was significantly lower, suggesting a global connection deficit. On a structural level, diffuse deficits in connectivity have also been described in TSC. First, animal models of TSC have demonstrated aberrant structural connectivity on a neural level, stemming from abnormalities in myelination, guidance and specification of the axon [57–59]. Second, in human subjects, DTI studies have revealed widespread decreased white matter microstructural integrity (for a brief overview, see [3]). This study is the first to demonstrate a functional correlate of structurally aberrant connectivity in TSC.

In ASD as well as in TSC, the ratio of inter- over intra-hemispheric connectivity did not reflect abnormalities of the corpus callosum, evident from volumetric, microstructural and functional imaging studies [3, 7–9, 60]. This either suggests a more subtle decrease in interhemispheric functional connectivity or more widespread distribution of altered connectivity in these disorders.

In ASD, a consistent and significantly decreased ratio of long- over short-range coherence, in the setting of an unaltered mean coherence, supports the current model of autism as a developmental disconnection syndrome. A comprehensive synthesis of prior findings in EEG and MEG coherence studies of ASD is challenging due to methodological differences [21, 27]. Nonetheless, patterns of regional under-connectivity and local over-connectivity were found in several studies. Eighteen autistic adults were found to have locally increased frontal and temporal resting state theta coherence, and decreased coherence between frontal lobe and all other regions in the lower alpha band [45]. In 20 children with ASD, decreased intra-hemispheric and inter-hemispheric delta and theta coherences were reported [46]. Similar to our findings, Mathewson and colleagues found no significant difference of coherence at rest in the alpha band in adults with autism compared to controls [61]. Barttfeld et al. studied EEGs of 10 autistic adults with a measure of synchronization and reported a prominent deficit in long-range and an increase in short-range connectivity [25]. Our data for the first time demonstrate similar findings in both ASD associated with TSC, and in ASD alone, in support of a common mechanism.

Graph measures

In TSC, despite extensive neurological involvement, functional connectivity has not been studied before. In ASD, the classic autistic cognitive profile of superior simple information processing and impaired higher order information processing stresses the importance of examining functional network as a whole and not only specific connections between specific regions [10]. To investigate properties of the entire functional network both in TSC and in ASD, we applied graph theoretical analysis.

The widespread deficits in local and regional connectivity in TSC are reflected both in conventional and graph measures of coherence. The aberrant network topology results in a decreased efficiency of information processing. The miswiring of axonal connections may contribute to the pathogenesis of neurological deficits in TSC [3, 62, 63], and our EEG study demonstrates the functional implications of these structural abnormalities on a network level.

The comparison of TSC to other disorders is complicated by the developmental rather than neurodegenerative nature of the disconnection. Using MEG, a recent study on connectivity and demyelination from multiple sclerosis reported an increased path length and clustering coefficient, suggesting a more regular network topology [64]. A similar MEG study of connectivity in Alzheimer's disease, considered a disconnection syndrome, found a decreased path length and clustering coefficient, indicating a more random network [14]. Both studies found associations between neuropsychological performance measures and graph connectivity measures, underscoring the neurobiological relevance of network analysis.

In patients with ASD, no significant differences in network topology measures were found. In the model of autism as a developmental disconnection syndrome, decreased small-worldness could be anticipated, as demonstrated in a MEG study using the synchronization likelihood in young adults with high-functioning autism [25]. Decreased small-worldness may then represent a decrease in local specialization (lower clustering coefficient) and regional integration (longer average path length) [65]. Our data may not have shown this due to the young age of the study population, where such refinement has not taken place yet [56]. Also, the absence of a higher clustering coefficient and of a longer path length for ASD despite the decrease in long-over-short range coherence demonstrates that nodes are spatially more clustered (conventional analysis) but not functionally more clustered (network analysis). This discrepancy is based on the conceptual difference between 'physical distance' and 'network distance'.

Resilience

In the ASD group, a key finding of significantly increased resilience to targeted attack was found. Different explanations to this observation can be posited, each reflecting different aspects of proposed neurobiological mechanisms of ASD.

First, increased resilience in the networks of autistic subjects could be related to redundant connectivity patterns. An abundance of connections makes a network highly resilient to attacks. However, connections in the brain are formed at a high physical cost [11, 66] and the brain constantly negotiates the trade-off between wiring costs and topological efficiency [67]. In particular, normal early developmental overconnectivity is followed by a pruning of connections in the maturing brain [55], suggesting network refinement [56].

Physiologically, in autism, the redundancy of connections could be explained by an impaired pruning of connections in the aforementioned dynamic process. The remaining overconnected network operates at different scales, from the neuronal level to the system level, and is consistent with studies of early cerebral overgrowth in autism (for a summary, see [7, 26]).

Cognitively, overconnectivity can result in a poor signal-to-noise ratio where the system is flooded with noise that swamps the signal [34]. With a poor signal-to-noise ratio, the output of a network may not be sufficiently distinct to achieve the necessary information processing [30]. Thus, overconnectivity can create an abnormally undifferentiated response to any stimulus. This excess of information gets equal rather than selective attention, creating an overstimulated, inefficient and delayed processing bottleneck [34].

Second, increased resilience could imply decreased functional specialization of brain regions. In technological systems, resilience refers to redundancy, as the same function is performed by identical elements. In biological networks, however, it refers to degeneracy, as structurally different elements can perform the same function [48, 50]. A degenerate system implies less specialization as the same output can be generated by different elements. Therefore, the increased resilience found in ASD could indicate an excessively degenerate system, where the removal of targeted nodes does not much affect the global properties of the network. Their presence is apparently less critical to the network, providing evidence of decreased functional specialization of these nodes. Our finding of a decreased level of connectivity of the main three hubs in the ASD population adds support to this interpretation.

In summary, the integration of primary order perceptions into higher order concepts is altered in ASD, but whether this is a top-down deficit (developmental disconnection syndrome) or due to heightened primary processing remains unclear [28]. A study using local coherence measures has argued that decreased responsiveness of autistic subjects to external stimuli may stem from a signal reduction through excess dampening [21]. Our network approach, however, suggests it rather comes from an excess of information processed in an overconnected, less specialized network [7, 31].

Future directions

Nodes should best represent brain regions with coherent patterns of extrinsic anatomical or functional connections which is problematic with EEG [35]. Specifically, node definition with only 19 electrodes is problematic as the locations do not match well-defined functional regions. Higher density EEG with 128 or 256 channels in an experimental setting can overcome this problem in part, although, inherent to the technique, only superficial aspects of the brain network can be modeled. In addition, electrode positions relative to underlying anatomical structures and functional areas are subject to variability, to changes related to growth and maturation and to methods of electrode placement used (for example, high density electrode cap or net). Thus, the level of anatomic accuracy of the current study does not allow for examining in detail the relation between functional and structural connectivity. Volume conduction will result in lattice-like graphs with highly clustered connections between neighboring electrodes, potentially confounding analysis of network properties [66]. Other measures of connectivity are less sensitive to this problem such as the phase lag index, synchronization likelihood [14, 68, 69], and partial directed coherence [17]. Nevertheless, it is encouraging that many 'headline' results seem to be robust to methodological details at several steps of network generation [70].

Frequency bands clearly have different associations with different aspects of cognitive activity, with different roles in pathology and with different biophysical mechanisms, reviewed in [71, 72]. Also, connectivity levels between regions are different for each frequency band [73]. Some authors have proposed that long distance communication may be mainly reflected by synchronization in low frequency bands (alpha and theta range) while shorter distance local communication is supported by synchronization in beta and gamma frequency bands [74]. As a result, it is imperative that any study should try to be as complete as possible in investigating connectivity in the different bands. However, to limit the number of comparisons, but more importantly because of potential muscle artifact (interference with beta and gamma bands) and potential residual motion artifact (interference with the delta band) in the EEG data of this challenging population, our study was restricted to the theta, upper and lower alpha bands.

Bands not studied include the delta band and the beta and gamma bands. Slower brain oscillations in the delta and sub-delta range appear to have a physiological role in sensory processing and cognition, even in the absence of environmental stimulation (for example, the default mode network in resting EEG and fMRI studies) [75]. In autistic patients, Coben and colleagues found the most significant coherence changes in autistic patients in the delta and theta bands compared to controls [46]. Gamma frequency oscillatory activity has been implicated in local cognitive processes and in the development of distributed cortical networks through both resting-state and task-related neural synchrony in this band [76]. Unfortunately, while analysis of these bands would be no additional challenge to execute, with the current limitations of the clinically acquired data, the findings would not be meaningful.

Finally, the heterogeneity of the study population should be emphasized. First, the use of clinical EEG data as controls could have introduced a bias from subtle EEG abnormalities that have escaped routine interpretation by the clinical neurophysiologist, and controls were recruited from a population with neurological complaints. Second, autism is a spectrum disorder with a wide range of severity, not well reflected by a binary variable (presence or absence of ASD). In future studies continuous variables could be used, such as the calibrated severity score of the ADOS [77], or the Social Responsiveness Score (SRS), which was recently used in another cross-disorder approach of autism [78]. Third, patients with TSC have wide variability in their phenotypical presentation. We did not incorporate anti-epileptic or psychoactive medication, and epilepsy-severity variables into the model, while in TSC up to 90% experience seizures in their lifetime [1]. With this high co-occurrence of epilepsy, cognitive impairment and autism in TSC [1], differences found may represent more global neurocognitive and behavioral dysfunction in TSC [78], although in our patients we did not find an association between severe cognitive impairment and autism. While EEG segments with epileptic discharges were excluded from analysis, it remains possible that epilepsy impacted the network analysis, in particular of the TSC population. However, interictal functional networks of patients with epilepsy are characterized by increased connectivity (especially in the theta band) and topological changes including increased regularity and hub-like organization [79–82] which we did not find.

In summary, several limitations of our study including its retrospective nature, the sources of EEG data and the possibility of different cognitive states of the subjects [27] can largely be addressed by a prospective study design. The burden associated with an EEG procedure is especially prominent in the young, low-functioning autistic population and may only be justified by clinical indication [21] warranting the exploratory use of already collected data. Still, a large, prospective, multicenter endeavor for determination of advanced neuroimaging and EEG correlates of autism in TSC has been launched. In addition, graph theoretical analysis of resting state fMRI connectivity in our population could validate our findings with much higher spatial resolution.

Conclusions

Connectivity analysis can provide fundamental insights into temporal functional coupling of spatially separate, specialized brain regions. Our EEG coherence study demonstrates decreased functional connectivity related to TSC in a global manner and to ASD in a more complex pattern.

In TSC, this study is the first to demonstrate altered functional connectivity, both through direct measurement of EEG coherence and on a network level. These results may represent a functional correlate of structural connectivity abnormalities in TSC and contribute to the neurological pathogenesis in TSC.

In ASD, a decreased long-over short-range coherence and markedly increased resilience to targeted attack renders an excessively degenerate network with local overconnection and decreased functional specialization.

Authors' information

JMP is a pediatric neurologist, epileptologist and clinical neurophysiologist with an interest in EEG signal processing. He aims to increase objectivity of EEG interpretation in the clinical setting, and improve availability to the non-neurophysiologist by translating the signal to colors, numbers or graphs. In the Computational Radiology Laboratory, his research focuses on establishing imaging and neurophysiological correlates of the neurological phenotype in Tuberous Sclerosis Complex (TSC), a neurocutaneous disorder strongly but unpredictably associated with epilepsy, cognitive disability and autism. MT is an electrical engineer and graduate student in neuroimaging. His work focuses on the development of novel image processing and statistical techniques to carry out population studies of the brain structure and function. In particular, he aims to better understand and analyze the brain network and its resilience. His contributions have been applied to the analysis of neurological differences in tuberous sclerosis complex and autism.

Abbreviations

ACC:

absent corpus callosum

ADOS:

Autism Diagnostic Observation Schedule

ANCOVA:

analysis of covariance

ASD:

Autism Spectrum Disorder

C:

control

DSM:

Diagnostic and Statistical Manual

DTI:

diffusion tensor imaging

EEG:

electroencephalography

(f)MRI:

(functional) magnetic resonance imaging

MEG:

magnetoencephalography

TSC:

Tuberous Sclerosis Complex.

References

  1. Crino PB, Nathanson KL, Henske EP: The tuberous sclerosis complex. N Engl J Med. 2006, 355: 1345-1356. 10.1056/NEJMra055323.

    Article  CAS  PubMed  Google Scholar 

  2. Jeste SS, Sahin M, Bolton P, Ploubidis GB, Humphrey A: Characterization of autism in young children with tuberous sclerosis complex. J Child Neurol. 2008, 23: 520-525. 10.1177/0883073807309788.

    Article  PubMed  Google Scholar 

  3. Peters JM, Sahin M, Vogel-Farley VK, Jeste SS, Nelson CA, Gregas MC, Prabhu SP, Scherrer B, Warfield SK: Loss of white matter microstructural integrity is associated with adverse neurological outcome in tuberous sclerosis complex. Acad Radiol. 2012, 19: 17-25. 10.1016/j.acra.2011.08.016.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Khwaja OS, Sahin M: Translational research: Rett syndrome and tuberous sclerosis complex. Curr Opin Pediatr. 2011, 23: 633-639. 10.1097/MOP.0b013e32834c9251.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Sahin M: Targeted treatment trials for tuberous sclerosis and autism: no longer a dream. Curr Opin Neurobiol. 2012, 22: 895-901. 10.1016/j.conb.2012.04.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Hampson DR, Gholizadeh S, Pacey LK: Pathways to drug development for autism spectrum disorders. Clin Pharmacol Ther. 2012, 91: 189-200. 10.1038/clpt.2011.245.

    Article  CAS  PubMed  Google Scholar 

  7. Courchesne E, Pierce K, Schumann CM, Redcay E, Buckwalter JA, Kennedy DP, Morgan J: Mapping early brain development in autism. Neuron. 2007, 56: 399-413. 10.1016/j.neuron.2007.10.016.

    Article  CAS  PubMed  Google Scholar 

  8. Keller TA, Kana RK, Just MA: A developmental study of the structural integrity of white matter in autism. Neuroreport. 2007, 18: 23-27. 10.1097/01.wnr.0000239965.21685.99.

    Article  PubMed  Google Scholar 

  9. Alexander AL, Lee JE, Lazar M, Boudos R, DuBray MB, Oakes TR, Miller JN, Lu J, Jeong EK, McMahon WM, Bigler ED, Lainhart JE: Diffusion tensor imaging of the corpus callosum in Autism. Neuroimage. 2007, 34: 61-73. 10.1016/j.neuroimage.2006.08.032.

    Article  PubMed  Google Scholar 

  10. Travers BG, Adluru N, Ennis C, Tromp DP, Destiche D, Doran S, Bigler ED, Lange N, Lainhart JE, Alexander AL: Diffusion tensor imaging in autism spectrum disorder: a review. Autism Res. 2012, 5: 289-313. 10.1002/aur.1243.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Bullmore E, Sporns O: Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009, 10: 186-198. 10.1038/nrn2575.

    Article  CAS  PubMed  Google Scholar 

  12. Varela F, Lachaux JP, Rodriguez E, Martinerie J: The brainweb: phase synchronization and large-scale integration. Nat Rev Neurosci. 2001, 2: 229-239.

    Article  CAS  PubMed  Google Scholar 

  13. Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL: A default mode of brain function. Proc Natl Acad Sci USA. 2001, 98: 676-682. 10.1073/pnas.98.2.676.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Stam CJ, de Haan W, Daffertshofer A, Jones BF, Manshanden I, van Cappellen van Walsum AM, Montez T, Verbunt JP, de Munck JC, van Dijk BW, Berendse HW, Scheltens P: Graph theoretical analysis of magnetoencephalographic functional connectivity in Alzheimer's disease. Brain. 2009, 132: 213-224.

    Article  CAS  PubMed  Google Scholar 

  15. Truccolo W, Donoghue JA, Hochberg LR, Eskandar EN, Madsen JR, Anderson WS, Brown EN, Halgren E, Cash SS: Single-neuron dynamics in human focal epilepsy. Nat Neurosci. 2011, 14: 635-641. 10.1038/nn.2782.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. van Putten MJ: Proposed link rates in the human brain. J Neurosci Methods. 2003, 127: 1-10. 10.1016/S0165-0270(03)00090-6.

    Article  PubMed  Google Scholar 

  17. Sakkalis V: Review of advanced techniques for the estimation of brain connectivity measured with EEG/MEG. Comput Biol Med. 2011, 41: 1110-1117. 10.1016/j.compbiomed.2011.06.020.

    Article  CAS  PubMed  Google Scholar 

  18. Deen B, Pelphrey K: Perspective: brain scans need a rethink. Nature. 2012, 491: S20-10.1038/491S20a.

    Article  CAS  PubMed  Google Scholar 

  19. Van de Ville D, Britz J, Michel CM: EEG microstate sequences in healthy humans at rest reveal scale-free dynamics. Proc Natl Acad Sci USA. 2010, 107: 18179-18184. 10.1073/pnas.1007841107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bosl W, Tierney A, Tager-Flusberg H, Nelson C: EEG complexity as a biomarker for autism spectrum disorder risk. BMC Med. 2011, 9: 18-10.1186/1741-7015-9-18.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Duffy FH, Als H: A stable pattern of EEG spectral coherence distinguishes children with autism from neuro-typical controls-a large case control study. BMC Med. 2012, 10: 64-10.1186/1741-7015-10-64.

    Article  PubMed  PubMed Central  Google Scholar 

  22. van Dellen E, Douw L, Baayen JC, Heimans JJ, Ponten SC, Vandertop WP, Velis DN, Stam CJ, Reijneveld JC: Long-term effects of temporal lobe epilepsy on local neural networks: a graph theoretical analysis of corticography recordings. PLoS One. 2009, 4: e8081-10.1371/journal.pone.0008081.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Tsiaras V, Simos PG, Rezaie R, Sheth BR, Garyfallidis E, Castillo EM, Papanicolaou AC: Extracting biomarkers of autism from MEG resting-state functional connectivity networks. Comput Biol Med. 2011, 41: 1166-1177. 10.1016/j.compbiomed.2011.04.004.

    Article  PubMed  Google Scholar 

  24. Pollonini L, Patidar U, Situ N, Rezaie R, Papanicolaou AC, Zouridakis G: Functional connectivity networks in the autistic and healthy brain assessed using Granger causality. Conf Proc IEEE Eng Med Biol Soc. 2010, 2010: 1730-1733.

    PubMed  Google Scholar 

  25. Barttfeld P, Wicker B, Cukier S, Navarta S, Lew S, Sigman M: A 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-263. 10.1016/j.neuropsychologia.2010.11.024.

    Article  PubMed  Google Scholar 

  26. Belmonte MK, Allen G, Beckel-Mitchener A, Boulanger LM, Carper RA, Webb SJ: Autism and abnormal development of brain connectivity. J Neurosci. 2004, 24: 9228-9231. 10.1523/JNEUROSCI.3340-04.2004.

    Article  CAS  PubMed  Google Scholar 

  27. Vissers ME, Cohen MX, Geurts HM: Brain connectivity and high functioning autism: a promising path of research that needs refined models, methodological convergence, and stronger behavioral links. Neurosci Biobehav Rev. 2012, 36: 604-625. 10.1016/j.neubiorev.2011.09.003.

    Article  PubMed  Google Scholar 

  28. Geschwind DH, Levitt P: Autism spectrum disorders: developmental disconnection syndromes. Curr Opin Neurobiol. 2007, 17: 103-111. 10.1016/j.conb.2007.01.009.

    Article  CAS  PubMed  Google Scholar 

  29. Minshew NJ, Williams DL: The new neurobiology of autism: cortex, connectivity, and neuronal organization. Arch Neurol. 2007, 64: 945-950. 10.1001/archneur.64.7.945.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Rippon G, Brock J, Brown C, Boucher J: Disordered connectivity in the autistic brain: challenges for the "new psychophysiology". Int J Psychophysiol. 2007, 63: 164-172. 10.1016/j.ijpsycho.2006.03.012.

    Article  PubMed  Google Scholar 

  31. Rubenstein JL, Merzenich MM: Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes Brain Behav. 2003, 2: 255-267. 10.1034/j.1601-183X.2003.00037.x.

    Article  CAS  PubMed  Google Scholar 

  32. Anderson JS, Druzgal TJ, Froehlich A, DuBray MB, Lange N, Alexander AL, Abildskov T, Nielsen JA, Cariello AN, Cooperrider JR, Bigler ED, Lainhart JE: Decreased interhemispheric functional connectivity in autism. Cereb Cortex. 2011, 21: 1134-1146. 10.1093/cercor/bhq190.

    Article  PubMed  Google Scholar 

  33. Dinstein I, Pierce K, Eyler L, Solso S, Malach R, Behrmann M, Courchesne E: Disrupted neural synchronization in toddlers with autism. Neuron. 2011, 70: 1218-1225. 10.1016/j.neuron.2011.04.018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Belmonte MK, Cook EH, Anderson GM, Rubenstein JL, Greenough WT, Beckel-Mitchener A, Courchesne E, Boulanger LM, Powell SB, Levitt PR, Perry EK, Jiang YH, DeLorey TM, Tierney E: Autism as a disorder of neural information processing: directions for research and targets for therapy. Mol Psychiatry. 2004, 9: 646-663.

    CAS  PubMed  Google Scholar 

  35. Rubinov M, Sporns O: Complex network measures of brain connectivity: uses and interpretations. Neuroimage. 2010, 52: 1059-1069. 10.1016/j.neuroimage.2009.10.003.

    Article  PubMed  Google Scholar 

  36. Roach ES, Gomez MR, Northrup H: Tuberous sclerosis complex consensus conference: revised clinical diagnostic criteria. J Child Neurol. 1998, 13: 624-628. 10.1177/088307389801301206.

    Article  CAS  PubMed  Google Scholar 

  37. Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, Pickles A, Rutter M: The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000, 30: 205-223. 10.1023/A:1005592401947.

    Article  CAS  PubMed  Google Scholar 

  38. Jung TP, Makeig S, Humphries C, Lee TW, McKeown MJ, Iragui V, Sejnowski TJ: Removing electroencephalographic artifacts by blind source separation. Psychophysiology. 2000, 37: 163-178. 10.1016/S0167-8760(00)00088-X.

    Article  CAS  PubMed  Google Scholar 

  39. Friston KJ, Tononi G, Reeke GN, Sporns O, Edelman GM: Value-dependent selection in the brain: simulation in a synthetic neural model. Neuroscience. 1994, 59: 229-243. 10.1016/0306-4522(94)90592-4.

    Article  CAS  PubMed  Google Scholar 

  40. Fingelkurts AA, Kahkonen S: Functional connectivity in the brain--is it an elusive concept?. Neurosci Biobehav Rev. 2005, 28: 827-836. 10.1016/j.neubiorev.2004.10.009.

    Article  PubMed  Google Scholar 

  41. van Drongelen W: Signal Processing for Neuroscientists: An Introduction to the Analysis of Physiological Signals. 2006, Burlington, MA: Academic Press

    Google Scholar 

  42. Srinivasan R, Winter WR, Ding J, Nunez PL: EEG and MEG coherence: measures of functional connectivity at distinct spatial scales of neocortical dynamics. J Neurosci Methods. 2007, 166: 41-52. 10.1016/j.jneumeth.2007.06.026.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Alstott J, Breakspear M, Hagmann P, Cammoun L, Sporns O: Modeling the impact of lesions in the human brain. PLoS Comput Biol. 2009, 5: e1000408-10.1371/journal.pcbi.1000408.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Nielsen T, Montplaisir J, Lassonde M: Decreased interhemispheric EEG coherence during sleep in agenesis of the corpus callosum. Eur Neurol. 1993, 33: 173-176. 10.1159/000116928.

    Article  CAS  PubMed  Google Scholar 

  45. Murias M, Webb SJ, Greenson J, Dawson G: Resting state cortical connectivity reflected in EEG coherence in individuals with autism. Biol Psychiatry. 2007, 62: 270-273. 10.1016/j.biopsych.2006.11.012.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Coben R, Clarke AR, Hudspeth W, Barry RJ: EEG power and coherence in autistic spectrum disorder. Clin Neurophysiol. 2008, 119: 1002-1009. 10.1016/j.clinph.2008.01.013.

    Article  PubMed  Google Scholar 

  47. Achard S, Salvador R, Whitcher B, Suckling J, Bullmore E: A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J Neurosci. 2006, 26: 63-72. 10.1523/JNEUROSCI.3874-05.2006.

    Article  CAS  PubMed  Google Scholar 

  48. Price CJ, Friston KJ: Degeneracy and cognitive anatomy. Trends Cogn Sci. 2002, 6: 416-421. 10.1016/S1364-6613(02)01976-9.

    Article  PubMed  Google Scholar 

  49. Sole RV, Ferrer Cancho R, Montoya JM, Valverde S: Selection, tinkering, and emergence in complex networks. Complexity. 2003, 8: 20-33.

    Article  Google Scholar 

  50. Tononi G, Sporns O, Edelman GM: Measures of degeneracy and redundancy in biological networks. Proc Natl Acad Sci USA. 1999, 96: 3257-3262. 10.1073/pnas.96.6.3257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Hagberg AA, Schult DA, Swart PJ: Exploring network structure, dynamics, and function using NetworkX. 7th Python in Science Conference (SciPy 2008): Aug 2008; Pasadena, CA. 2008, 11-15.

    Google Scholar 

  52. Tarokh L, Carskadon MA, Achermann P: Developmental changes in brain connectivity assessed using the sleep EEG. Neuroscience. 2010, 171: 622-634. 10.1016/j.neuroscience.2010.08.071.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Thatcher RW, Walker RA, Giudice S: Human cerebral hemispheres develop at different rates and ages. Science. 1987, 236: 1110-1113. 10.1126/science.3576224.

    Article  CAS  PubMed  Google Scholar 

  54. Micheloyannis S, Vourkas M, Tsirka V, Karakonstantaki E, Kanatsouli K, Stam CJ: The influence of ageing on complex brain networks: a graph theoretical analysis. Hum Brain Mapp. 2009, 30: 200-208. 10.1002/hbm.20492.

    Article  PubMed  Google Scholar 

  55. Supekar K, Musen M, Menon V: Development of large-scale functional brain networks in children. PLoS Biol. 2009, 7: e1000157-10.1371/journal.pbio.1000157.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Hagmann P, Sporns O, Madan N, Cammoun L, Pienaar R, Wedeen VJ, Meuli R, Thiran JP, Grant PE: White matter maturation reshapes structural connectivity in the late developing human brain. Proc Natl Acad Sci USA. 2010, 107: 19067-19072. 10.1073/pnas.1009073107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Choi YJ, Di Nardo A, Kramvis I, Meikle L, Kwiatkowski DJ, Sahin M, He X: Tuberous sclerosis complex proteins control axon formation. Genes Dev. 2008, 22: 2485-2495. 10.1101/gad.1685008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Meikle L, Talos DM, Onda H, Pollizzi K, Rotenberg A, Sahin M, Jensen FE, Kwiatkowski DJ: A mouse model of tuberous sclerosis: neuronal loss of Tsc1 causes dysplastic and ectopic neurons, reduced myelination, seizure activity, and limited survival. J Neurosci. 2007, 27: 5546-5558. 10.1523/JNEUROSCI.5540-06.2007.

    Article  CAS  PubMed  Google Scholar 

  59. Nie D, Di Nardo A, Han JM, Baharanyi H, Kramvis I, Huynh T, Dabora S, Codeluppi S, Pandolfi PP, Pasquale EB, Sahin M: Tsc2-Rheb signaling regulates EphA-mediated axon guidance. Nat Neurosci. 2010, 13: 163-172. 10.1038/nn.2477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Just MA, Cherkassky VL, Keller TA, Kana RK, Minshew NJ: Functional and anatomical cortical underconnectivity in autism: evidence from an FMRI study of an executive function task and corpus callosum morphometry. Cereb Cortex. 2007, 17: 951-961.

    Article  PubMed  Google Scholar 

  61. Mathewson KJ, Jetha MK, Drmic IE, Bryson SE, Goldberg JO, Schmidt LA: Regional EEG alpha power, coherence, and behavioral symptomatology in autism spectrum disorder. Clin Neurophysiol. 123: 1798-1809.

  62. Tsai P, Sahin M: Mechanisms of neurocognitive dysfunction and therapeutic considerations in tuberous sclerosis complex. Curr Opin Neurol. 2011, 24: 106-113. 10.1097/WCO.0b013e32834451c4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Lewis WW, Sahin M, Scherrer B, Peters JM, Suarez RO, Vogel-Farley VK, Jeste SS, Gregas MC, Prabhu SP, Nelson CA, Warfield SK: Impaired language pathways in tuberous sclerosis complex patients with autism spectrum disorders. Cereb Cortex.

  64. Schoonheim MM, Geurts JJ, Landi D, Douw L, van der Meer ML, Vrenken H, Polman CH, Barkhof F, Stam CJ: Functional connectivity changes in multiple sclerosis patients: a graph analytical study of MEG resting state data. Hum Brain Mapp. 2013, 34: 52-61. 10.1002/hbm.21424.

    Article  PubMed  Google Scholar 

  65. Sporns O, Zwi JD: The small world of the cerebral cortex. Neuroinformatics. 2004, 2: 145-162. 10.1385/NI:2:2:145.

    Article  PubMed  Google Scholar 

  66. Bullmore E, Bassett DS: Brain graphs: graphical models of the human brain connectome. Annu Rev Clin Psychol. 2011, 7: 113-140. 10.1146/annurev-clinpsy-040510-143934.

    Article  PubMed  Google Scholar 

  67. Bullmore E, Sporns O: The economy of brain network organization. Nat Rev Neurosci. 2012, 13: 336-349.

    CAS  PubMed  Google Scholar 

  68. Stam CJ, Nolte G, Daffertshofer A: Phase lag index: assessment of functional connectivity from multi channel EEG and MEG with diminished bias from common sources. Hum Brain Mapp. 2007, 28: 1178-1193. 10.1002/hbm.20346.

    Article  PubMed  Google Scholar 

  69. Stam CJ, van Straaten EC: Go with the flow: use of a directed phase lag index (dPLI) to characterize patterns of phase relations in a large-scale model of brain dynamics. Neuroimage. 2012, 62: 1415-1428. 10.1016/j.neuroimage.2012.05.050.

    Article  CAS  PubMed  Google Scholar 

  70. Bassett DS, Bullmore ET: Human brain networks in health and disease. Curr Opin Neurol. 2009, 22: 340-347. 10.1097/WCO.0b013e32832d93dd.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Basar E, Basar-Eroglu C, Karakas S, Schurmann M: Gamma, alpha, delta, and theta oscillations govern cognitive processes. Int J Psychophysiol. 2001, 39: 241-248. 10.1016/S0167-8760(00)00145-8.

    Article  CAS  PubMed  Google Scholar 

  72. Knyazev GG: Motivation, emotion, and their inhibitory control mirrored in brain oscillations. Neurosci Biobehav Rev. 2007, 31: 377-395. 10.1016/j.neubiorev.2006.10.004.

    Article  PubMed  Google Scholar 

  73. Stam CJ, van Straaten EC: The organization of physiological brain networks. Clin Neurophysiol. 2012, 123: 1067-1087. 10.1016/j.clinph.2012.01.011.

    Article  CAS  PubMed  Google Scholar 

  74. von Stein A, Sarnthein J: Different frequencies for different scales of cortical integration: from local gamma to long range alpha/theta synchronization. Int J Psychophysiol. 2000, 38: 301-313. 10.1016/S0167-8760(00)00172-0.

    Article  CAS  PubMed  Google Scholar 

  75. Van Someren EJ, Van Der Werf YD, Roelfsema PR, Mansvelder HD, da Silva FH: Slow brain oscillations of sleep, resting state, and vigilance. Prog Brain Res. 2011, 193: 3-15.

    Article  CAS  PubMed  Google Scholar 

  76. Uhlhaas PJ, Roux F, Rodriguez E, Rotarska-Jagiela A, Singer W: Neural synchrony and the development of cortical networks. Trends Cogn Sci. 14: 72-80.

  77. Gotham K, Pickles A, Lord C: Standardizing ADOS scores for a measure of severity in autism spectrum disorders. J Autism Dev Disord. 2009, 39: 693-705. 10.1007/s10803-008-0674-3.

    Article  PubMed  Google Scholar 

  78. van Eeghen AM, Pulsifer MB, Merker VL, Neumeyer AM, van Eeghen EE, Thibert RL, Cole AJ, Leigh FA, Plotkin SR, Thiele EA: Understanding relationships between autism, intelligence, and epilepsy: a cross-disorder approach. Dev Med Child Neurol. 2013, 55: 146-153. 10.1111/dmcn.12044.

    Article  PubMed  Google Scholar 

  79. Chavez M, Valencia M, Navarro V, Latora V, Martinerie J: Functional modularity of background activities in normal and epileptic brain networks. Phys Rev Lett. 2010, 104: 118701.

    Article  CAS  PubMed  Google Scholar 

  80. Douw L, van Dellen E, de Groot M, Heimans JJ, Klein M, Stam CJ, Reijneveld JC: Epilepsy is related to theta band brain connectivity and network topology in brain tumor patients. BMC Neurosci. 2010, 11: 103-10.1186/1471-2202-11-103.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Horstmann MT, Bialonski S, Noennig N, Mai H, Prusseit J, Wellmer J, Hinrichs H, Lehnertz K: State dependent properties of epileptic brain networks: comparative graph-theoretical analyses of simultaneously recorded EEG and MEG. Clin Neurophysiol. 2010, 121: 172-185. 10.1016/j.clinph.2009.10.013.

    Article  PubMed  Google Scholar 

  82. Wilke C, Worrell G, He B: Graph analysis of epileptogenic networks in human partial epilepsy. Epilepsia. 2010, 52: 84-93.

    Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This work was supported in part by the National Institutes of Health (grant numbers R01 RR021885, R01 LM010033, R03 EB008680, UL1 RR025758 to S.K.W., P20 RFA-NS-12-006, 1U01NS082320-01 to M.S. and J.M.P.), the National Institute of Mental Health (grant number K23MH094517 to S.S.J.), the National Institute on Deafness and Other Communication Disorders (grant number DC 10290 to C.A.N.), and the Department of Defense (grant number W81XWH-11-1-0365 to C.A.N.) In addition, JMP is supported by a Faculty Development Fellowship from the 'Eleanor and Miles Shore 50th Anniversary Fellowship Program for Scholars in Medicine', Boston Children's Hospital, Department of Neurology, 2012-2013 and performs video-EEG long-term monitoring, EEGs, and other electrophysiological studies at Boston Children's Hospital and bills for these procedures. MT is supported by the Fonds de la Recherche Scientifique-FNRS, and the Belgian American Educational Foundation. CAN is supported by the Simons Foundation. MS is supported by the John Merck Fund and a Junior Investigator Award from the Boston Children's Hospital Translational Research Program. The authors thank the EEG technologists of BCH for their high quality data collection in this challenging population, Tina Shimizu and Amanda Noroña for their assistance with the Los Angeles data, and Caterina Stamoulis, PhD, for her help with questions related to coherence. We are most indebted to the patients and their families for their support and time contributed.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jurriaan M Peters or Maxime Taquet.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JMP and MT conceived and designed the experimental approach, performed data collection, data analysis and prepared the manuscript. CV was critical for analysis and interpretation of phenotypical data. ISF and JT assisted with data collection and pre-processing. CN and SSJ assisted with data collection and significant edits of the manuscript. MS and SKW supervised the project, and contributed to experimental design and analysis. All authors read and approved the final version of the manuscript.

Jurriaan M Peters, Maxime Taquet contributed equally to this work.

Electronic supplementary material

12916_2012_739_MOESM1_ESM.PDF

Additional file 1: Graph Analysis: An introduction. Description: An easily accessible introduction to graph theory and network analysis. Each measure is explained with both text and illustration, and for each an example is given of airline networks and brain networks. (PDF 144 KB)

12916_2012_739_MOESM2_ESM.PDF

Additional file 2: Graph Analysis: Resilience. Description: An entry-level description of the concept of network resilience and the main methods of assessing resilience through two modes of attack of the network. Examples are provided for airline networks and brain networks. (PDF 70 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Peters, J.M., Taquet, M., Vega, C. et al. Brain functional networks in syndromic and non-syndromic autism: a graph theoretical study of EEG connectivity. BMC Med 11, 54 (2013). https://doi.org/10.1186/1741-7015-11-54

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1741-7015-11-54

Keywords