Skip to main content

Organization of brain networks governed by long-range connections index autistic traits in the general population

Abstract

Background

The dimensional approach to autism spectrum disorder (ASD) considers ASD as the extreme of a dimension traversing through the entire population. We explored the potential utility of electroencephalography (EEG) functional connectivity as a biomarker. We hypothesized that individual differences in autistic traits of typical subjects would involve a long-range connectivity diminution within the delta band.

Methods

Resting-state EEG functional connectivity was measured for 74 neurotypical subjects. All participants also provided a questionnaire (Social Responsiveness Scale, SRS) that was completed by an informant who knows the participant in social settings. We conducted multivariate regression between the SRS score and functional connectivity in all EEG frequency bands. We explored modulations of network graph metrics characterizing the optimality of a network using the SRS score.

Results

Our results show a decay in functional connectivity mainly within the delta and theta bands (the lower part of the EEG spectrum) associated with an increasing number of autistic traits. When inspecting the impact of autistic traits on the global organization of the functional network, we found that the optimal properties of the network are inversely related to the number of autistic traits, suggesting that the autistic dimension, throughout the entire population, modulates the efficiency of functional brain networks.

Conclusions

EEG functional connectivity at low frequencies and its associated network properties may be associated with some autistic traits in the general population.

Background

In the seminal paper in which Kanner first characterized autism, he noticed that some of the symptoms observed in children were shared at sub-threshold levels by their parents [1]. This original observation has been confirmed by several studies [2, 3]. Today, the predominant view is that autism spectrum disorder (ASD) is not an all-or-nothing condition; instead, its severity is graded and can be quantified through several diagnostic assessments [46]. The graded nature of ASD and the progression of autistic traits in the general population suggest that ASD constitutes the extreme of a dimension traversing through the entire population [7, 8]. This has led to the development of behavioural measures seeking to evaluate related traits in the general population, such as the Social Responsiveness Scale (SRS) [7], and to use this dimensional approach to study ASD [79].

There is increasing evidence that ASD could be a condition of altered brain connectivity [1013], including a reduced corpus callosum [14] and diminished long-range functional connectivity [12, 15, 16], producing a system that is ineffective for integrating information [1719]. Functional brain networks of ASD compared to typical subjects in the resting state (that is, during free thought) showed qualitatively different organizations at the group level, which broadly reflects a deficit in long-range connectivity, especially along the long-distance fronto-posterior axis [16, 18]; see [19] for a review. Most of the knowledge about ASD brain connectivity has been acquired through the use of functional magnetic resonance imaging (fMRI); however, changes in connectivity based on stationary electroencephalography (EEG) measures can also reliably discriminate between ASD and control populations [16]. Compared to fMRI, EEG has remarkable economical and practical advantages [20], such as the much less stressful set-up and a short application time, making it better for large-scale screening, especially with a resting-state paradigm. There are theoretical and practical motivations for using the resting state to assess clinical populations [20, 21]. Of these there are the relatively short data-acquisition times, a simple set-up, the subject does not need to be stimulated, behavioural responses do not need to be collected (allowing a broader sampling of patient populations) and a better signal-to-noise ratio compared to task-related protocols [21].

To assist in the diagnosis of ASD, it is essential to find robust brain biomarkers that characterize ASD as the upper extreme of a dimension across the entire population. To date only one piece of research has addressed this issue: it showed that fMRI connectivity in a single link connecting the anterior cingulate cortex and the mid-insula diminishes in strength as the number of autistic traits increases [22]. Here, instead, we study how global brain organization, measured through network properties derived from EEG [23, 24], relates to autistic traits in the general population.

In a previous study we showed a diminution of long-range connectivity within the (low-frequency) delta band in an ASD population compared to control groups, leading to a ‘big-world’ organization of brain connectivity in ASD [16]. Here we investigate whether this progression is gradually modulated by the autistic traits of typical subjects by studying resting-state EEG functional networks. Specifically, we hypothesize that: (1) markers for autistic traits in the general population are indexed by the strength of long-range connections (compared to short-range connections), predominantly fronto-occipital connections, (2) changes in connectivity with autistic traits are most prominent in low-frequency bands, which are decoupled in the autistic population [16] and (3) the small-world index significantly decreases with autistic traits in the general population.

Methods

Participants and assessment

In this study, there were 74 subjects of similar educational and cultural backgrounds (37 male, 37 female; mean age = 27.33, SD = 5.10; educational level = 19.2 years; SD = 2.94 years). None of the volunteers had a history of neurological or psychiatric conditions as determined by a semi-structured interview (Schedules of Clinical Assessment in Neuropsychiatry) [25]. The interview was conducted by a trained physician in order to exclude any subjects with psychiatric, neurological or sensory impairment, addictions (such as alcohol and drug abuse) or general cognitive impairment. After the subjects were given a complete description of the study, written informed consent was obtained in agreement with the Declaration of Helsinki and the Institution’s ethical committee, which approved this work.

We requested participants to select someone who knew them well, preferably a close relative or partner, to complete the adult version of the SRS questionnaire (SRS-A, [7]). SRS is a 65-item questionnaire completed by an informant who knows the evaluated subject’s preferences and personality. It measures autistic traits across the entire range of severity observed in nature. SRS outputs a score that indexes the severity of social deficits. Higher scores on the SRS indicate greater severity of social impairment. Scores between 60 and 80 are associated with mild forms of autism [7]. Within our population, the scores varied between 11 and 69. Since no tool is currently available in Spanish for the quantitative assessment of autistic impairment across a wide range of severity including the identification of sub-threshold levels of autistic symptomatology, we created an Argentine version of the SRS and ran a pilot study as a first step to full validation. The translation of the adult scale (SRS-A) into Spanish required two simultaneous translations by qualified professionals (author SC and his clinical team at Buenos Aires [26]) with education, training and work experience in the diagnosis of autism spectrum conditions and in the use of psychological tests and assessments in research and clinical settings. After the translation, two back-translations were made and two consensus meetings were held. The final translation was reviewed for clarity by a panel of experts in ASD diagnostics not involved in the translation. After the translation, a pilot study was ran in which the SRS-A test was completed by relatives or close friends of ten adults with an ASD diagnosis, as assessed by the Autistic Diagnostic Observation Schedule (ADOS), and ten typical adults, all with an IQ higher than 85. Interpretation of SRS scores was centred on the total scores. The SRS mean total T-value was highest in the ASD sample (mean SRS = 121.40, SD = 23.71) and lowest in the typically developing sample (mean SRS = 36.9, SD = 9.91; T-value = 10.39; P = 4.90 × 10-9). Also, we observed a significant positive correlation between ADOS and SRS scores within the ASD group (correlation = 0.67, P = 0.03, R-squared = 0.42). All ASD subjects scored in the higher range of the mild to moderate group (mild or high functioning autism spectrum conditions) or in the severe range group (autistic disorder or severe cases of pervasive developmental disorder not otherwise specified (PDD-NOS) or Asperger’s syndrome) of the SRS profile sheet. Findings in this pilot study provide adequate support for the application of SRS-A in the assessment of autism traits of the subjects in the present study.

EEG results

EEG measurements were taken in a Faraday cage with a Biosemi Active Two 128-channel 24-bit resolution system, with active electrodes (the first amplifying stage on the electrode improves the signal-to-noise ratio), digitalized at 512 Hz and low-passed DC-1/5th of the sample rate (-3 dB) by a fifth-order digital sync anti-aliasing filter. There were no additional hardware filters during acquisition. Temporal signals between 5 and 10 minutes were recorded during an eyes-closed rest while subjects sat on a reclining chair in a sound-attenuated room with a dim light. During the experiment, participants and EEG recordings were monitored to ensure that they maintained vigilance and did not fall asleep. After the acquisition, signals were re-referenced to the average of all electrodes. Segments containing movement artefacts were manually deleted (mean length of the remaining time series = 7.90 min, SD = 2.21), followed by an ICA-based rejection of residual artefact-laden ICA-components. After this pre-processing, we filtered the EEG signals on specific frequency bands: delta (0.5 Hz to 4 Hz), theta (4 Hz to 8 Hz), alpha (8 Hz to 12 Hz), sigma (12 Hz to 15 Hz), beta (15 Hz to 25 Hz) and gamma (25 Hz to 35 Hz).

Data analysis

MATLAB (MathWorks Inc, Natick, MA) was used for the analyses. We first explored the relation between SRS score and gender. We calculated the mean SRS score for men and women, and assessed their difference statistically using a t-test for two independent samples. To study the relation between SRS score and age, we conducted both a t-test comparing groups (low SRS score and high SRS score groups, obtained after a median split on the SRS score) and a regression between SRS score and age (Figure 1a, b).

Figure 1
figure 1

Autistic traits and functional EEG connectivity. (a) SRS scores. The SRS score is higher for men than women. Error bars show the SEM. (b) Scatter plot showing that age does not modulate SRS within our population. (c) Cumulative β-value distributions, showing that the entire distributions are shifted to negative values. The dotted red line is a cumulative value of 0.95. (d) Mean and SEM of β-value distributions, obtained from independent multivariate regression between SRS score and SL connectivity (including gender and age as regressors of no interest). See Table 1 for mean and standard deviations of the distributions and the P-values. All frequency bands except alpha show a bias to negative values. Matrices show thresholded β-values; a blue entry in the matrix is a negative β-value, whereas a red entry is a positive β-value. (e) Scalp plot of β-values. A link is traced between two electrodes if that connection is significant (P < 0.05, uncorrected) in all frequency bands whose β-values distribution was proved to be significantly shifted to negative values (that is, delta, theta, beta and gamma).

To test whether synchronisation likelihood (SL) connectivity in spontaneous activity between electrodes covaries with SRS score, we conducted a functional connectivity analysis. The synchronisation between all pair-wise combinations of EEG channels was computed for all subjects with the SL method [27]. SL quantifies the probability that a pair of channels is synchronized. For each participant p and frequency band f we measured a 128 × 128 connectivity matrix SL f,p . A matrix entry SL f,p (i,j) indicates the temporal synchronization of the signal measured by electrodes i and j, for subject p at the frequency band f, which henceforth is referred to as the functional connectivity. All subsequent analysis and statistics were performed on these SL f,p matrices. To investigate connectivity changes associated with SRS score, we conducted an across-subjects multivariate linear regression, using least squares (as implemented in MATLAB function regstats()) between each entry of the matrix SL f,p and the SRS score for each subject, including gender and age of the subject as regressors of no interest. This lead to six matrices of beta (β) values, B f (i,j), one per frequency band f. For example, a positive value for B Delta (i,j) indicates that connectivity between electrode i and electrode j increases with SRS score for networks measured in the delta band. For visualization we projected all B f (i,j) values exceeding a threshold of P = 0.05, uncorrected, into a scalp plot (Figure 1e). This threshold is arbitrary: it was used only for visualization and played no role in any statistical analysis.

To assess the B f matrices statistically we performed a bootstrap analysis [28]. We obtained the mean value of each B f matrix, and called it the observed mean B f per frequency band f, since this was the β-value obtained experimentally. Then, we calculated the null distribution of β-values for each frequency by band shuffling the SRS scores across participants (thus breaking any possible dependence between functional connectivity and autistic traits across individuals), and repeated the whole regression analysis, to obtain a random mean β-value. We repeated this procedure 5,000 times, obtaining for each frequency band a distribution of 5,000 random mean β-values that approaches a Gaussian distribution. This distribution of random β-values is called a null distribution, or the distribution of expected β-values under the hypothesis of no relation between SRS score and functional connectivity. If any of our observed β-values truly reflects a covariation between functional connectivity and SRS score, its value should be located on the tails of the null distributions. We fitted a Gaussian to each distribution of random β-values to obtain a Z-score, by subtracting from the observed β-value the mean value of the fitted Gaussian and dividing it by the standard deviation of the fitted Gaussian. This Z-score reflects the distance between the mean of the random distribution of β-values and the observed β-value. We obtained the P-values corresponding to the Z-scores, and set the P-value threshold for significance at 0.05, Bonferroni corrected for multiple comparisons.

To further characterize the β-value distributions, we calculated the standard error of the mean (SEM) for each distribution B f through a jackknife procedure [29], repeating the regression N – 1 times (where N is the number of subjects), each time excluding a different subject from the analysis. The SEM was then calculated as std B N * N - 1 where B N is the standard deviation of the β-values over the N regressions.

To estimate the discriminative power of SL at characterizing autistic traits, we calculated a receiver operating characteristic (ROC) curve [30], categorizing as a ‘hit’ each time a high SRS subject (after a median split) was assigned to the high SRS group, and as a ‘false alarm’ each time a low SRS subject was assigned to the high SRS group. The area under the curve (AROC) quantifies how separable the two groups are: AROC = 0.50 means that the two groups completely overlap (along the variable considered), while AROC = 1 indicates that the two groups are perfectly separable by their respective SRS scores.

To address the issue of length of connections and their relation with SRS score, we defined four different regions grouping electrodes: frontal, occipital, lateral right and lateral left (Figure 2), covering the contiguous frontal, occipital and temporal electrodes. We then measured connectivity between and within these regions averaging the SL value across electrodes, and we performed a linear multivariate regression (using the least squares method) between SL and SRS score, including age and gender as covariables of no interest. We also conducted a ROC analysis as described above, to explore how well SL values separate between subjects with low and high SRS scores (after grouping subjects into low and high SRS score groups by means of a median split on the SRS score).

Figure 2
figure 2

SL connectivity between and within ROIs. (a) ROC curve value for all frequency bands, between and within frontal and occipital ROIs. Fronto-occipital SL connectivity produced ROC curves showing significant P-values for the delta, theta, beta and gamma bands. See Additional file 2: Table S2 for all AROC values. (b) ROC curve value for all frequency bands, between and within right and left ROIs. Left-right SL connectivity produced ROC curves showing significant P-values for the delta and theta bands. See Additional file 4: Table S3 for all AROC values. (c) ROC curves after performing a median-split group separation, classifying subjects as high or low SRS score according to their gamma band SL value. Only the fronto-occipital SL shows significant classification (AROC fronto-occipital = 0.86; AROC occipital = 0.51; AROC frontal = 0.59). (d) Multivariate regressions obtained by averaging all gamma band SL values within (frontal and occipital) and between (fronto-occipital) ROIs, against SRS score (and gender and age as other regressors). Frontal and occipital SLs do not show significant regression with SRS score (occipital: P = 0.23, R-square = 0.01; frontal: P = 0.31, R-square = 0.001). The fronto-occipital SL, in contrast, does show good regression between SL and SRS score (P < 0.05, R-square = 0.09).

Graph theory metrics

We used graph theory metrics to summarize topological information. The connectivity matrix SL f,p defines a weighted graph where each electrode corresponds to a node and the weight of each link is determined by the SL of the electrode pair. To calculate network measures, SL f,p matrices were converted to binary undirected matrices by applying a threshold T. The arbitrary parameter T was chosen so that in all cases the resulting networks had a link density of 0.10, that is, 10% of the total number of possible links in the networks were actually present, to ensure that only the strongest links were present and that the network was not disaggregated into subcomponents [31], and to normalize networks of different subjects by size, in order to avoid spurious effects on the metrics due to network size. After transforming the SL f,p matrix to a binary undirected graph, we calculated the clustering coefficient C and the characteristic path length L using the Brain Connectivity Toolbox [32]. Combining the metrics C and L, we calculated the small-world index, C / L, as an estimate of the small-world properties of the networks. Small world refers to a ubiquitous topological network that has a relatively short (compared to random networks) L and high C[33]. The small-world index quantifies the optimality of a network in terms of information processing and storage.

To quantify the impact of autistic traits on network properties, we performed two analyses. First we conducted a median-split analysis, grouping subjects as for the SL quantification analysis (conducting a t-test and a ROC analysis). Also we conducted a multivariate linear regression analysis between the small-world index and SRS score (including gender and age as regressors of no interest). To assess the regression analysis statistically, we conducted a bootstrap analysis repeating the methods used to assess the significance of the B f matrices. We obtained the P-value corresponding to the observed small-world index, and set the P-value threshold for significance at 0.05, Bonferroni corrected for multiple comparisons.

Results

First we measured the dependence of SRS score on the demographical covariates age and gender. As expected [4], the women had on average a lower SRS score than the men (women = 34.29 ± 1.66; men = 41.00 ± 1.72; T-value = -1.99; P < 0.05), indicating that men are more likely to have autistic traits than women (Figure 1a). In contrast, age did not show an effect with SRS score, either through a median-split (younger = 38.65 ± 1.99; older = 37.77 ± 1.46; T-value = 0.24; P = 0.80) or through a regression between SRS score and age (β-value = 0.09; P = 0.78) (Figure 1b).

Next we measured the relation between SRS score and connectivity. For each participant in this study, we calculated the synchronization likelihood across all pairs of channels. The element (i,j) of the SL matrix provides an estimate of the probability that the time series of electrodes i and j are related during eyes-closed stationary EEG, which we refer to as functional connectivity. To see SL changes along the autistic dimension we conducted a multivariate regression between each entry of the SL matrix C(s,p) ij and the values for total SRS, with gender and age as regressors of no interest. A positive β-value βij indicates that the SL between electrodes i and j increased as the SRS value increased. Conversely, a negative β-value indicates that SL is greater when the SRS score diminished (Figure 1d).

The regression analysis showed an overall decrease of the mean connectivity (averaged across all electrode pairs) as the SRS score increased. This global decrease in connectivity was observed for all frequency bands, except in the alpha band and was significant for the delta, theta, beta and sigma bands (see Table 1). The effect of the decrease in connectivity with increasing SRS over the entire distribution is clearly depicted by plotting the cumulative value of the entire distribution of β-values (Figure 1d).

Table 1 Mean, standard deviation and significance of β-values from the regression between SL and SRS score

Next, we investigated the hypothesis that long-range connections are more informative about individual autistic traits than short connections. We averaged SL values for all pairs of electrodes within and between frontal, occipital and lateral clusters of electrodes (Figure 2). This electrode clustering broadly defines cortical regions and is not intended (due to the low resolution of the EEG) to define precise boundaries of fine cortical structures. After a median split of subjects according to their SRS score (grouping by high SRS and low SRS scores), we estimated the discriminative power of within and between cluster connectivity to characterize autistic traits, by calculating a ROC curve [30]; we categorize as a ‘hit’ each time a high SRS subject is assigned to the high SRS group, and as a ‘false alarm’ each time a low SRS subject is assigned to the high SRS group. We observed that the SL measured from electrode pairs connecting the frontal and occipital clusters classified subjects better than the SL obtained from electrode pairs within each cluster (Figure 2a). This difference was significant for the SL measured in the delta, theta, beta and gamma bands (P < 0.01, Bonferroni corrected; all AROC are listed in Additional file 1: Table S1). The SL measured for the alpha and sigma bands, as expected from our previous analysis, did not produce a significant classification of subjects. The negative relation between SRS score and between-clusters SL was also observed from a multivariate regression, using SL as the dependent variable and the SRS score as a regressor (along with age and gender as covariates of no interest). The regression between SRS score and within-clusters connectivity (Figure 2d depicts results for the gamma band; other bands gave similar results) was not significant (occipital: P = 0.23, R-square = 0.01; frontal: P = 0.31, R-square = 0.001). In contrast, the regression between SRS score and between-regions SL was significant (P < 0.05, R-square = 0.09) (Figure 2a), showing that fronto-occipital SL decreases monotonically as the SRS score increases.

The SL between electrode pairs between the left and right clusters also classified subjects better than the SL obtained from electrode pairs within left or right clusters (Figure 2b), although significant differences were only observed in the lower (delta and theta) frequencies (all AROC values are listed in Additional file 2: Table S2). This suggests that the observed relation between SRS score and between-clusters SL reflects a general effect of distance on SL (see Additional file 3: Figure S1).

Finally we examined the hypothesis that changes in connectivity result in a different network topology for participants with high or low autistic traits. Specifically we hypothesized that the small worldness of the network, a property which weights the compactness and clustering of a network using optimal architectures for information storage and propagation [17, 23], decreased with increasing prevalence of autistic traits. We calculated the small-world index, which quantifies the similarity of the network to a ubiquitous topological network usually referred to as small world [23], and conducted a median-split analysis, grouping subjects as previously, based on their SRS score. High and low SRS score groups showed different mean small-world indices only for the delta band (low SRS group = 0.38, high SRS group = 0.22; T-value = 4.17; P = 0.0001; AROC = 0.75; Figure 3a,b; see Additional file 4: Table S3 for the size effect of all frequency bands). To see whether the differences quantified with the median-split analysis were also observed in a more rigorous analysis not dependent on an explicit group definition along the autistic dimension, we conducted a multivariate regression between the small-world index and SRS score (Figure 3c). For delta band SL matrices, the regression showed a covariation between the small-world index and SRS score (P = 0.008, R-square = 0.09), showing that SRS score is negatively related to the optimality of the brain network. Relations between SRS score and the small-world index obtained from all other bands were not significant (see Additional file 5: Table S4 for all effect sizes).

Figure 3
figure 3

Topology of networks obtained from delta band SL f,p matrices. (a) Small-world indices, after a median-split of subjects into high and low SRS groups. Low SRS subjects have a higher small-world index (low SRS group = 0.38, high SRS group = 0.22; T-value = 4.17; P = 0.0001; see Additional file 4: Table S3 for the size effect of all frequency bands). (b) ROC curve for the discrimination between low and high SRS score based on the small-world index. AROC is highly significant (AROC = 0.74, P < 0.001). (c) Regression between small-world index and SRS score (and age and gender), showing the negative relation between small-world index and SRS score (β-value = -0.08; P = 0.008, R-square = 0.09).

Discussion

The main purpose of this study was to characterize and compare resting-state functional brain networks in typical subjects along the autistic dimension. In agreement with our three working hypothesis we observed that: (1) the SRS score in the general population is indexed by the strength of connections between EEG electrodes – long-range connections are more predictive than short-range connections of an individual’s SRS score, (2) the lower part of the EEG spectrum is the most informative for individual autistic traits and (3) the small worldness of the network (and hence its optimality for storage and transfer of information) increases as the SRS score diminishes.

Long-range intra-cortical and feedback cortico-cortical connections, which are thought to be altered in ASD, are revealed by the slow cortical potentials of the EEG [34]. Cortico-cortical connections can be roughly classified in two main groups [23, 35]: local connections linking neurons in the same cortical area (which are critical in generating functional specificity, that is, information) and long-distance connections between neurons of different cortical regions (which ensure that distant cortical sites can interact rapidly to generate dynamical patterns of temporal correlations, allowing the integration of different sources of information into coherent behavioural and cognitive states) [23, 30]. Long-range connectivity provided a good correlate of the individual level of autistic traits, suggesting that functional brain connectivity across distant cortical regions is modulated by the SRS score, diminishing as the SRS score increases and becoming indicative of ASD [16]. The reduced long-range connections may provide a physiological measure for the lack of proper integration of information observed in ASD [36].

Changes in connectivity patterns have an impact on the global organization of a network, which in turn determines the efficiency of information transfer and storage [23, 37]. Small-world networks have attracted significant attention during recent decades [33] because they are ubiquitously present in a broad range of natural phenomenon and also because they establish an optimal balance between local specialization and global integration [23]. Our results suggest that the functional networks in the general population are related to the number of individual autistic traits.

In order to become useful tools and assist ASD diagnosis, brain-imaging techniques must find robust biomarkers to characterize ASD and its relation with sub-threshold traits in the general population. Only one piece of research has addressed this issue; Di Martino et al. [22] showed that connectivity within two nodes of the saliency network diminishes in strength as the number of autistic traits increases in neurotypical adults. An estimation of a single connection, however, might not constitute a robust biomarker for characterizing such a complex and diverse condition as ASD, and its performance in actual subject classification remains to be tested. On the other hand, a network approach involving global measures of connectivity and network quality might be better for a robust biomarker [24]. In addition fMRI is not practical for large-scale clinical screening and EEG is a much more suitable, economical and practical tool [20]. Our finding that the optimality of an individual’s EEG network is markedly related to the individual’s SRS score show that it is a good candidate for a biomarker characterizing ASD with practical clinical relevance for large-scale fast screening. Moreover, most brain-imaging evidence suggests that ASD is associated with a diminished connectivity between the frontal lobe and occipitoparietal regions, typically involving default mode network (DMN) nodes such as the ventromedial prefrontal cortex and the precuneus/posterior cingulate [1013]. It is possible then to have a composite biomarker combining global metrics such as the small-world index and the underlying well-known changes in network topography, a combination that might outperform single-measure biomarkers.

One limitation of the present study is the lack of a measure of specificity regarding possible comorbidities in the autistic traits we measured. The SRS score might capture, along with autistic traits, traits characterizing other psychiatric conditions [38]. Future work should test this network approach through proper classification studies that assess its performance in realistic diagnostic situations and measure its capacity not only to characterize the autistic dimension but to detect it specifically.

Conclusions

The present study demonstrates that a resting-state EEG can identify robust and monotonic changes associated with SRS score, a graded measure of autistic traits in the general population. Our results show a decrease in functional connectivity, mainly for the delta and theta bands, is associated with an increased number of autistic traits. When inspecting the impact of autistic traits on the global organization of the functional network we found that the optimal properties of the network are inversely related to the number of autistic traits, suggesting that the autistic dimension, throughout the entire population, modulates the efficiency of functional brain networks.

Abbreviations

ADOS:

Autism Diagnostic Observation Schedule

ASD:

autism spectrum disorder

EEG:

electroencephalography

fMRI:

functional magnetic resonance imaging

ROC:

receiver operating characteristic

SEM:

standard error of the mean

SL:

synchronisation likelihood

SRS:

Social Responsiveness Scale

ICA:

independent component analysis

ROI:

region of interest.

References

  1. Kanner L: Autistic disturbances of affective contact. Acta Paedopsychiatr. 1968, 35: 100-136.

    CAS  PubMed  Google Scholar 

  2. Hoekstra RA, Bartels M, Verweij CJ, Boomsma DI: Heritability of autistic traits in the general population. Arch Pediatr Adolesc Med. 2007, 161: 372-377. 10.1001/archpedi.161.4.372.

    Article  PubMed  Google Scholar 

  3. Spiker D, Lotspeich LJ, Dimiceli S, Myers RM, Risch N: Behavioral phenotypic variation in autism multiplex families: evidence for a continuous severity gradient. Am J Med Genet. 2002, 114: 129-136. 10.1002/ajmg.10188.

    Article  PubMed  Google Scholar 

  4. Baron-Cohen S, Knickmeyer RC, Belmonte MK: Sex differences in the brain: implications for explaining autism. Science. 2005, 310: 819-823. 10.1126/science.1115455.

    Article  CAS  PubMed  Google Scholar 

  5. 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 

  6. Lord C, Rutter M, Le Couteur A: Autism diagnostic interview-revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994, 24: 659-685. 10.1007/BF02172145.

    Article  CAS  PubMed  Google Scholar 

  7. Constantino JN, Davis SA, Todd RD, Schindler MK, Gross MM, Brophy SL, Metzger LM, Shoushtari CS, Splinter R, Reich W: Validation of a brief quantitative measure of autistic traits: comparison of the social responsiveness scale with the autism diagnostic interview-revised. J Autism Dev Disord. 2003, 33: 427-433. 10.1023/A:1025014929212.

    Article  PubMed  Google Scholar 

  8. Constantino JN, Todd RD: Intergenerational transmission of subthreshold autistic traits in the general population. Biol Psychiatry. 2005, 57: 655-660. 10.1016/j.biopsych.2004.12.014.

    Article  PubMed  Google Scholar 

  9. Lau YC, Hinkley LB, Bukshpun P, Strominger ZA, Wakahiro ML, Baron-Cohen S, Allison C, Auyeung B, Jeremy RJ, Nagarajan SS: Autism traits in individuals with agenesis of the corpus callosum. J Autism Dev Disord. 2013, 43: 1106-1118. 10.1007/s10803-012-1653-2.

    Article  PubMed Central  PubMed  Google Scholar 

  10. 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 

  11. Courchesne E, Pierce K: Why the frontal cortex in autism might be talking only to itself: local over-connectivity but long-distance disconnection. Curr Opin Neurobiol. 2005, 15: 225-230. 10.1016/j.conb.2005.03.001.

    Article  CAS  PubMed  Google Scholar 

  12. 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 Central  PubMed  Google Scholar 

  13. Just MA, Cherkassky VL, Keller TA, Minshew NJ: Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity. Brain. 2004, 127: 1811-1821. 10.1093/brain/awh199.

    Article  PubMed  Google Scholar 

  14. Alexander AL, Lee JE, Lazar M, Boudos R, DuBray MB, Oakes TR, Miller JN, Lu J, Jeong EK, McMahon WM: 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 

  15. Barttfeld P, Wicker B, Cukier S, Navarta S, Lew S, Leiguarda R, Sigman M: State-dependent changes of connectivity patterns and functional brain network topology in autism spectrum disorder. Neuropsychologia. 2012, 50: 3653-3662. 10.1016/j.neuropsychologia.2012.09.047.

    Article  PubMed  Google Scholar 

  16. 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 

  17. Gallos LK, Sigman M, Makse HA: The conundrum of functional brain networks: small-world efficiency or fractal modularity. Front Physiol. 2012, 3: 123-

    Article  PubMed Central  PubMed  Google Scholar 

  18. Schipul SE, Keller TA, Just MA: Inter-regional brain communication and its disturbance in autism. Front Syst Neurosci. 2011, 5: 10-

    Article  PubMed Central  PubMed  Google Scholar 

  19. Wass S: Distortions and disconnections: disrupted brain connectivity in autism. Brain Cogn. 2011, 75: 18-28. 10.1016/j.bandc.2010.10.005.

    Article  PubMed  Google Scholar 

  20. Cruse D, Chennu S, Chatelle C, Bekinschtein TA, Fernandez-Espejo D, Pickard JD, Laureys S, Owen AM: Bedside detection of awareness in the vegetative state: a cohort study. Lancet. 2011, 378: 2088-2094. 10.1016/S0140-6736(11)61224-5.

    Article  PubMed  Google Scholar 

  21. Fox MD, Greicius M: Clinical applications of resting state functional connectivity. Front Syst Neurosci. 2010, 4: 19-

    PubMed Central  PubMed  Google Scholar 

  22. Di Martino A, Shehzad Z, Kelly C, Roy AK, Gee DG, Uddin LQ, Gotimer K, Klein DF, Castellanos FX, Milham MP: Relationship between cingulo-insular functional connectivity and autistic traits in neurotypical adults. Am J Psychiatry. 2009, 166: 891-899. 10.1176/appi.ajp.2009.08121894.

    Article  PubMed Central  PubMed  Google Scholar 

  23. 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 

  24. Stam CJ, De Haan W, Daffertshofer A, Jones BF, Manshanden I, van Walsum AM VC, Montez T, Verbunt JP, de Munck JC, van Dijk BW: Graph theoretical analysis of magnetoencephalographic functional connectivity in Alzheimer’s disease. Brain. 2009, 132: 213-224.

    Article  CAS  PubMed  Google Scholar 

  25. World Health Organization: Schedules for Clinical Assessment in Neuropsychiatry (SCAN). 1992, Geneva

    Google Scholar 

  26. Programa Argentino para Niños, Adolescentes y Adultos con Condiciones del Espectro Autista.http://www.panaacea.org,

  27. Montez T, Linkenkaer-Hansen K, van Dijk BW, Stam CJ: Synchronization likelihood with explicit time-frequency priors. NeuroImage. 2006, 33: 1117-1125. 10.1016/j.neuroimage.2006.06.066.

    Article  CAS  PubMed  Google Scholar 

  28. Efron B, Tibshirani RJ: An Introduction to the Bootstrap. 1994, New York: Chapman and Hall

    Google Scholar 

  29. Miller RG: The jackknife – a review. Biometrika. 1974, 61: 1-15.

    Google Scholar 

  30. Zweig MH, Campbell G: Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem. 1993, 39: 561-577.

    CAS  PubMed  Google Scholar 

  31. Tagliazucchi E, Von Wegner F, Morzelewski A, Brodbeck1 V, Laufs H: Dynamic BOLD functional connectivity in humans and its electrophysiological correlates. Front Hum Neurosci. 2012, 6: 339-

    Article  PubMed Central  PubMed  Google Scholar 

  32. 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 

  33. Watts DJ, Strogatz SH: Collective dynamics of ‘small-world’ networks. Nature. 1998, 393: 440-442. 10.1038/30918.

    Article  CAS  PubMed  Google Scholar 

  34. He BJ, Raichle ME: The fMRI signal, slow cortical potential and consciousness. Trends Cogn Sci. 2009, 13: 302-309. 10.1016/j.tics.2009.04.004.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Schroeder CE, Lakatos P: Low-frequency neuronal oscillations as instruments of sensory selection. Trends Neurosci. 2009, 32: 9-18. 10.1016/j.tins.2008.09.012.

    Article  CAS  PubMed  Google Scholar 

  36. Frith U: Autism: Explaining the Enigma. 1989, Oxford: Basil Blackwell

    Google Scholar 

  37. Barabasi AL: Scale-free networks: a decade and beyond. Science. 2009, 325: 412-413. 10.1126/science.1173299.

    Article  CAS  PubMed  Google Scholar 

  38. Torralva T, Gleichgerrcht E, Torrente F, Roca M, Strejilevich SA, Cetkovich M, Lischinsky A, Manes F: Neuropsychological functioning in adult bipolar disorder and ADHD patients: a comparative study. Psychiatry Res. 2011, 186: 261-266. 10.1016/j.psychres.2010.08.007.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

This work is funded by CONICET and UBACYT. MS is sponsored by the James McDonnell Foundation 21st Century Science Initiative in Understanding Human Cognition - Scholar Award. MS and PB are supported by the Human Frontiers Science Program and a post-doctoral grant (PB). AI is supported by CONICET, CONICYT/FONDECYT Regular (1130920), INECO Foundation grant, and grants FONDECYT (1130920) and PICT 2012–1309. FM is supported by a grant, PICT 2012–0412. LB is supported by a post-doctoral grant from CONICET, Argentina.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Agustín Ibanez.

Additional information

Competing interests

All authors report no financial relationships with commercial interests.

Authors’ contributions

PB and MS conceived the experiment. PB, SC, LA, JA, LB and AT collected the data. PB and MS analysed the data. PB, SC, MS, AI and FM wrote the paper. All authors read and approved the final manuscript.

Electronic supplementary material

11689_2013_41_MOESM1_ESM.doc

Additional file 1: Table S1: AROC for all frequency bands, for the combination of ROIs frontal, occipital and fronto-occipital. (DOC 28 KB)

11689_2013_41_MOESM2_ESM.doc

Additional file 2: Table S2: AROC for all frequency bands, for the combination of ROIs right, left and right-left. (DOC 27 KB)

11689_2013_41_MOESM3_ESM.pdf

Additional file 3: Figure S1: SL connectivity between long and short distances. (a) ROC curve value for all frequency bands, between and within frontal and occipital ROIs. (b) Long-distance SL connectivity produced ROC curves showing significant P-values for the delta, theta, beta and gamma bands, while short-distance SL connectivity produced ROC curves showing significant P-values only for the delta band. (PDF 415 KB)

11689_2013_41_MOESM4_ESM.doc

Additional file 4: Table S3: Significance and size effects of the t-test between low and high SRS groups for all frequency bands. (DOC 26 KB)

11689_2013_41_MOESM5_ESM.doc

Additional file 5: Table S4: Significance and size effects of the regression analysis between the small-world index and SRS score for all frequency bands. (DOC 20 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://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

Barttfeld, P., Amoruso, L., Ais, J. et al. Organization of brain networks governed by long-range connections index autistic traits in the general population. J Neurodevelop Disord 5, 16 (2013). https://doi.org/10.1186/1866-1955-5-16

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1866-1955-5-16

Keywords