Complexity of Multi-Channel Electroencephalogram Signal Analysis in Childhood Absence Epilepsy

Absence epilepsy is an important epileptic syndrome in children. Multiscale entropy (MSE), an entropy-based method to measure dynamic complexity at multiple temporal scales, is helpful to disclose the information of brain connectivity. This study investigated the complexity of electroencephalogram (EEG) signals using MSE in children with absence epilepsy. In this research, EEG signals from 19 channels of the entire brain in 21 children aged 5-12 years with absence epilepsy were analyzed. The EEG signals of pre-ictal (before seizure) and ictal states (during seizure) were analyzed by sample entropy (SamEn) and MSE methods. Variations of complexity index (CI), which was calculated from MSE, from the pre-ictal to the ictal states were also analyzed. The entropy values in the pre-ictal state were significantly higher than those in the ictal state. The MSE revealed more differences in analysis compared to the SamEn. The occurrence of absence seizures decreased the CI in all channels. Changes in CI were also significantly greater in the frontal and central parts of the brain, indicating fronto-central cortical involvement of “cortico-thalamo-cortical network” in the occurrence of generalized spike and wave discharges during absence seizures. Moreover, higher sampling frequency was more sensitive in detecting functional changes in the ictal state. There was significantly higher correlation in ictal states in the same patient in different seizures but there were great differences in CI among different patients, indicating that CI changes were consistent in different absence seizures in the same patient but not from patient to patient. This implies that the brain stays in a homogeneous activation state during the absence seizures. In conclusion, MSE analysis is better than SamEn analysis to analyze complexity of EEG, and CI can be used to investigate the functional brain changes during absence seizures.     

Long-Term Effects of Temporal Lobe Epilepsy on Local Neural Networks: A Graph Theoretical Analysis of Corticography Recordings

Purpose

Pharmaco-resistant temporal lobe epilepsy (TLE) is often treated with surgical intervention at some point. As epilepsy surgery is considered a last resort by most physicians, a long history of epileptic seizures prior to surgery is not uncommon. Little is known about the effects of ongoing TLE on neural functioning. A better understanding of these effects might influence the moment of surgical intervention. Functional connectivity (interaction between spatially distributed brain areas) and network structure (integration and segregation of information processing) are thought to be essential for optimal brain functioning. We report on the impact of TLE duration on temporal lobe functional connectivity and network characteristics.

Methods

Functional connectivity of the temporal lobe at the time of surgery was assessed by means of interictal electrocorticography (ECoG) recordings of 27 TLE patients by using the phase lag index (PLI). Graphs (abstract network representations) were reconstructed from the PLI matrix and characterized by the clustering coefficient C (local clustering), the path length L (overall network interconnectedness), and the “small world index” S (network configuration).

Results

Functional connectivity (average PLI), clustering coefficients, and the small world index were negatively correlated with TLE duration in the broad frequency band (0.5–48 Hz).

Discussion

Temporal lobe functional connectivity is lower in patients with longer TLE history, and longer TLE duration is correlated with more random network configuration. Our findings suggest that the neural networks of TLE patients become more pathological over time, possibly due to temporal lobe changes associated with long-standing lesional epilepsy.     

Widespread EEG Changes Precede Focal Seizures

The process by which the brain transitions into an epileptic seizure is unknown. In this study, we investigated whether the transition to seizure is associated with changes in brain dynamics detectable in the wideband EEG, and whether differences exist across underlying pathologies. Depth electrode ictal EEG recordings from 40 consecutive patients with pharmacoresistant lesional focal epilepsy were low-pass filtered at 500 Hz and sampled at 2,000 Hz. Predefined EEG sections were selected immediately before (immediate preictal), and 30 seconds before the earliest EEG sign suggestive of seizure activity (baseline). Spectral analysis, visual inspection and discrete wavelet transform were used to detect standard (delta, theta, alpha, beta and gamma) and high-frequency bands (ripples and fast ripples). At the group level, each EEG frequency band activity increased significantly from baseline to the immediate preictal section, mostly in a progressive manner and independently of any modification in the state of vigilance. Preictal increases in each frequency band activity were widespread, being observed in the seizure-onset zone and lesional tissue, as well as in remote regions. These changes occurred in all the investigated pathologies (mesial temporal atrophy/sclerosis, local/regional cortical atrophy, and malformations of cortical development), but were more pronounced in mesial temporal atrophy/sclerosis. Our findings indicate that a brain state change with distinctive features, in the form of unidirectional changes across the entire EEG bandwidth, occurs immediately prior to seizure onset. We postulate that these changes might reflect a facilitating state of the brain which enables a susceptible region to generate seizures.     

Characterization of synchronization in interacting groups of oscillators: application to seizures

We investigate the emergence of synchronization in two groups of oscillators; one group acts as a synchronization source, and the other as the target. Based on phase model simulations, we construct a synchrony index (SI): a combination of intra- and intergroup synchronies. The SI characterizes the extent of induced synchrony in the population. We demonstrate the usefulness of the measure in a test case of mesial temporal lobe epilepsy: the SI can be readily calculated from standard electroencephalographic measurements. We show that the synchrony index has a statistically significant increased value for the ictal periods and that the epileptic focus can be located by identifying the most synchronous pairs of electrodes during the initial part of ictal period of the seizure. We also show that it is possible in this pilot case to differentiate clinical and subclinical seizures based on the dynamical features of the synchronization. The synchronization index was found to be a useful quantity for the characterization of “pathological hypersynchronization” within a well-characterized patient with mesial temporal lobe epilepsy and thus has potential medical value in seizure detection, localizing ability, and association with later surgical outcome.     

The Role of Resting State Networks in Focal Neocortical Seizures

Objective

The role of resting state functional networks in epilepsy is incompletely understood. While some pathologic diagnoses have been shown to have maintained but altered resting state connectivity, others have implicated resting state connectivity in disease progression. However little is known about how these resting state networks influence the behavior of a focal neocortical seizure.

Methods

Using data taken from invasively monitored patients with intractable focal neocortical epilepsy, we evaluated network connectivity (as determined by oscillatory covariance of the slow cortical potential (<0.5 Hz)) as it relates to neocortical seizure foci both in the interictal and ictal states.

Results

Similar to what has been shown in the past for sleep and anesthesia, electophysiologic resting state networks that are defined by this slow cortical potential covariance maintain their topographic correlation structure throughout an ictal event. Moreover, in the context of focal epilepsy in which the seizure has a specific site of onset, seizure propagation is not chaotic or random. Rather, the seizure (reflected by an elevation of high frequency power) preferentially propagates along the network that contains the seizure onset zone.

Significance

Taken together, these findings further undergird the fundamental role of resting state networks, provide novel insights into the network-influenced behavior of seizures, and potentially identify additional targets for surgical disconnection including informing the location for the completion of multiple subpial transections (MSPTs).     

Scale Invariance Properties of Intracerebral EEG Improve Seizure Prediction in Mesial Temporal Lobe Epilepsy

Although treatment for epilepsy is available and effective for nearly 70 percent of patients, many remain in need of new therapeutic approaches. Predicting the impending seizures in these patients could significantly enhance their quality of life if the prediction performance is clinically practical. In this study, we investigate the improvement of the performance of a seizure prediction algorithm in 17 patients with mesial temporal lobe epilepsy by means of a novel measure. Scale-free dynamics of the intracerebral EEG are quantified through robust estimates of the scaling exponents—the first cumulants—derived from a wavelet leader and bootstrap based multifractal analysis. The cumulants are investigated for the discriminability between preictal and interictal epochs. The performance of our recently published patient-specific seizure prediction algorithm is then out-of-sample tested on long-lasting data using combinations of cumulants and state similarity measures previously introduced. By using the first cumulant in combination with state similarity measures, up to 13 of 17 patients had seizures predicted above chance with clinically practical levels of sensitivity (80.5%) and specificity (25.1% of total time under warning) for prediction horizons above 25 min. These results indicate that the scale-free dynamics of the preictal state are different from those of the interictal state. Quantifiers of these dynamics may carry a predictive power that can be used to improve seizure prediction performance.     

Disruption of rolandic gamma-band functional connectivity by seizures is associated with motor impairments in children with epilepsy

Although children with epilepsy exhibit numerous neurological and cognitive deficits, the mechanisms underlying these impairments remain unclear. Synchronization of oscillatory neural activity in the gamma frequency range (>30 Hz) is purported to be a mechanism mediating functional integration within neuronal networks supporting cognition, perception and action. Here, we tested the hypothesis that seizure-induced alterations in gamma synchronization are associated with functional deficits. By calculating synchrony among electrodes and performing graph theoretical analysis, we assessed functional connectivity and local network structure of the hand motor area of children with focal epilepsy from intracranial electroencephalographic recordings. A local decrease in inter-electrode phase synchrony in the gamma bands during ictal periods, relative to interictal periods, within the motor cortex was strongly associated with clinical motor weakness. Gamma-band ictal desychronization was a stronger predictor of deficits than the presence of the seizure-onset zone or lesion within the motor cortex. There was a positive correlation between the magnitude of ictal desychronization and impairment of motor dexterity in the contralateral, but not ipsilateral hand. There was no association between ictal desynchronization within the hand motor area and non-motor deficits. This study uniquely demonstrates that seizure-induced disturbances in cortical functional connectivity are associated with network-specific neurological deficits.     

Increased Intrinsic Connectivity of the Default Mode Network in Temporal Lobe Epilepsy: Evidence from Resting-State MEG Recordings

The electrophysiological signature of resting state oscillatory functional connectivity within the default mode network (DMN) during spike-free periods in temporal lobe epilepsy (TLE) remains unclear. Using magnetoencephalographic (MEG) recordings, this study investigated how the connectivity within the DMN was altered in TLE, and we examined the effect of lateralized TLE on functional connectivity. Sixteen medically intractable TLE patients and 22 controls participated in this study. Whole-scalp 306-channel MEG epochs without interictal spikes generated from both MEG and EEG data were analyzed using a minimum norm estimate (MNE) and source-based imaginary coherence analysis. With this processing, we obtained the cortical activation and functional connectivity within the DMN. The functional connectivity was increased between DMN and the right medial temporal (MT) region at the delta band and between DMN and the bilateral anterior cingulate cortex (ACC) regions at the theta band. The functional change was associated with the lateralization of TLE. The right TLE showed enhanced DMN connectivity with the right MT while the left TLE demonstrated increased DMN connectivity with the bilateral MT. There was no lateralization effect of TLE upon the DMN connectivity with ACC. These findings suggest that the resting-state functional connectivity within the DMN is reinforced in temporal lobe epilepsy during spike-free periods. Future studies are needed to examine if the altered functional connectivity can be used as a biomarker for treatment responses, cognitive dysfunction and prognosis in patients with TLE.     

颞叶癫痫脑“默认模式”的fMRI研究

功能磁共振成像(functional magnetic resonance imaging,fMRI)被用于检测静息时脑功能神经网络.作者运用静息fMRI检测海马硬化颞叶癫痫(temporal lobe epilepsy,TLE)脑"默认模式",采用感兴趣区域功能连接分析检测16例TLE患者和16名正常对照静息时脑的"默认模式",并进行组内和组间分析.研究发现,与正常对照相比,TLE静息时海马、颞极、额叶、颞叶、壳核及楔前叶等脑区与后扣带回的功能连接增强.研究结果表明TLE患者的固有脑功能组织模式有可能出现紊乱.这一研究将有助于从脑功能的角度了解癫痫患者某些临床症状的发病机理,为今后癫痫诊治的发展提供一定的帮助.     

Disruptions in Brain Networks of Older Fallers Are Associated with Subsequent Cognitive Decline: A 12-Month Prospective Exploratory Study

Cognitive impairment and impaired mobility are major public health concerns. There is growing recognition that impaired mobility is an early biomarker of cognitive impairment and dementia. The neural basis for this association is currently unclear. We propose disrupted functional connectivity as a potential mechanism. In this 12-month prospective exploratory study, we compared functional connectivity of four brain networks– the default mode network (DMN), fronto-executive network (FEN), fronto-parietal network (FPN), and the primary motor sensory network (SMN) – between community-dwelling older adults with ≥ two falls in the last 12 months and their non-falling counterparts (≤ one fall in the last 12 months). Functional connectivity was examined both at rest and during a simple motor tapping task. Compared with non-fallers, fallers showed more connectivity between the DMN and FPN during right finger tapping (p = 0.04), and significantly less functional connectivity between the SMN and FPN during rest (p≤0.05). Less connectivity between the SMN and FPN during rest was significantly associated with greater decline in both cognitive function and mobility over the12-month period (r = −0.32 and 0.33 respectively; p≤0.04). Thus, a recent history of multiple falls among older adults without a diagnosis of dementia may indicate sub-clinical changes in brain function and increased risk for subsequent decline.     

Dynamic Network Drivers of Seizure Generation,Propagation and Termination in Human Neocortical Epilepsy

The epileptic network is characterized by pathologic, seizure-generating ‘foci’ embedded in a web of structural and functional connections. Clinically, seizure foci are considered optimal targets for surgery. However, poor surgical outcome suggests a complex relationship between foci and the surrounding network that drives seizure dynamics. We developed a novel technique to objectively track seizure states from dynamic functional networks constructed from intracranial recordings. Each dynamical state captures unique patterns of network connections that indicate synchronized and desynchronized hubs of neural populations. Our approach suggests that seizures are generated when synchronous relationships near foci work in tandem with rapidly changing desynchronous relationships from the surrounding epileptic network. As seizures progress, topographical and geometrical changes in network connectivity strengthen and tighten synchronous connectivity near foci—a mechanism that may aid seizure termination. Collectively, our observations implicate distributed cortical structures in seizure generation, propagation and termination, and may have practical significance in determining which circuits to modulate with implantable devices.     

Functional Connectivity Estimated from Intracranial EEG Predicts Surgical Outcome in Intractable Temporal Lobe Epilepsy

This project aimed to determine if a correlation-based measure of functional connectivity can identify epileptogenic zones from intracranial EEG signals, as well as to investigate the prognostic significance of such a measure on seizure outcome following temporal lobe lobectomy. To this end, we retrospectively analyzed 23 adult patients with intractable temporal lobe epilepsy (TLE) who underwent an invasive stereo-EEG (SEEG) evaluation between January 2009 year and January 2012. A follow-up of at least one year was required. The primary outcome measure was complete seizure-freedom at last follow-up. Functional connectivity between two areas in the temporal lobe that were sampled by two SEEG electrode contacts was defined as Pearson’s correlation coefficient of interictal activity between those areas. SEEG signals were filtered between 5 and 50 Hz prior to computing this correlation. The mean and standard deviation of the off diagonal elements in the connectivity matrix were also calculated. Analysis of the mean and standard deviation of the functional connections for each patient reveals that 90% of the patients who had weak and homogenous connections were seizure free one year after temporal lobectomy, whereas 85% of the patients who had stronger and more heterogeneous connections within the temporal lobe had recurrence of seizures. This suggests that temporal lobectomy is ineffective in preventing seizure recurrence for patients in whom the temporal lobe is characterized by weakly connected, homogenous networks. This pilot study shows promising potential of a simple measure of functional brain connectivity to identify epileptogenicity and predict the outcome of epilepsy surgery.     

Analysis of Epileptic Discharges from Implanted Subdural Electrodes in Patients with Sturge-Weber Syndrome

ObjectiveAlmost two-thirds of patients with Sturge-Weber syndrome (SWS) have epilepsy, and half of them require surgery for it. However, it is well known that scalp electroencephalography (EEG) does not demonstrate unequivocal epileptic discharges in patients with SWS. Therefore, we analyzed interictal and ictal discharges from intracranial subdural EEG recordings in patients treated surgically for SWS to elucidate epileptogenicity in this disorder.MethodsFive intractable epileptic patients with SWS who were implanted with subdural electrodes for presurgical evaluation were enrolled in this study. We examined the following seizure parameters: seizure onset zone (SOZ), propagation speed of seizure discharges, and seizure duration by visual inspection. Additionally, power spectrogram analysis on some frequency bands at SOZ was performed from 60 s before the visually detected seizure onset using the EEG Complex Demodulation Method (CDM).ResultsWe obtained 21 seizures from five patients for evaluation, and all seizures initiated from the cortex under the leptomeningeal angioma. Most of the patients presented with motionless staring and respiratory distress as seizure symptoms. The average seizure propagation speed and duration were 3.1 ± 3.6 cm/min and 19.4 ± 33.6 min, respectively. Significant power spectrogram changes at the SOZ were detected at 10–30 Hz from 15 s before seizure onset, and at 30–80 Hz from 5 s before seizure onset.SignificanceIn patients with SWS, seizures initiate from the cortex under the leptomeningeal angioma, and seizure propagation is slow and persists for a longer period. CDM indicated beta to low gamma-ranged seizure discharges starting from shortly before the visually detected seizure onset. Our ECoG findings indicate that ischemia is a principal mechanism underlying ictogenesis and epileptogenesis in SWS.     

Changes in Task-Related Functional Connectivity across Multiple Spatial Scales Are Related to Reading Performance

Reading requires the interaction of a distributed set of cortical areas whose distinct patterns give rise to a wide range of individual skill. However, the nature of these neural interactions and their relation to reading performance are still poorly understood. Functional connectivity analyses of fMRI data can be used to characterize the nature of interactivity of distributed brain networks, yet most previous studies have focused on connectivity during task-free (i.e., “resting state”) conditions. Here, we report new methods for assessing task-related functional connectivity using data-driven graph theoretical methods and describe how large-scale patterns of connectivity relate to individual variability in reading performance among children. We found that connectivity patterns of subjects performing a reading task could be decomposed hierarchically into multiple sub-networks, and we observed stronger long-range interaction between sub-networks in subjects with higher task accuracy. Additionally, we found a network of hub regions known to be critical to reading that displays increased short-range synchronization in higher accuracy subjects. These individual differences in task-related functional connectivity reveal that increased interaction between distant regions, coupled with selective local integration within key regions, is associated with better reading performance. Importantly, we show that task-related neuroimaging data contains far more information than usually extracted via standard univariate analyses – information that can meaningfully relate neural connectivity patterns to cognition and task.     

Patient-Specific Epileptic Seizure Prediction in Long-Term Scalp EEG Signal Using Multivariate Statistical Process Control

An accurate epileptic seizure prediction algorithm can alleviate the problem and reduce risks in the life of a patient suffering from epilepsy. The main motive of this work is to propose a model which can predict seizures well in advance of its occurrence. Multivariate statistical process control (MSPC) has been used for seizure predictions in long-term scalp EEG signal. It has been observed that excessive neuronal activity in the preictal period of seizure changes the electrical characteristic from chaotic to rhythmic behavior. These changes have been utilized for prediction. Eight temporal based features are used for predicting the seizures by using multivariate statistical process control, which is widely known as an anomaly monitoring method. 90 seizures from the CHB-MIT EEG data of ten patients are analyzed.ResultThe results of the proposed method demonstrated that 80 seizures out of 90 in preictal period were correctly predicted prior to the seizure onset, thereby giving a sensitivity of 88.89%. The false positive rate is observed to 0.39 per hour.ConclusionThis study proposed a temporal based patient-specific epileptic seizure prediction method using MSPC in long-term scalp EEG signals. It also provides the possibility of realizing an EEG-based epileptic seizure prediction system which requires less computational power.SignificanceThe proposed method does not require preictal data for modeling. The extracted features are computationally easy. The tested result shows good accuracy on the CHB-MIT data base.     

Using network dynamic fMRI for detection of epileptogenic foci

Background

Epilepsy is one of the most prevalent neurological disorders. It remains medically intractable for about one-third of patients with focal epilepsy, for whom precise localization of the epileptogenic zone responsible for seizure initiation may be critical for successful surgery. Existing fMRI literature points to widespread network disturbances in functional connectivity. Per previous scalp and intracranial EEG studies and consistent with excessive local synchronization during interictal discharges, we hypothesized that, relative to same regions in healthy controls, epileptogenic foci would exhibit less chaotic dynamics, identifiable via entropic analyses of resting state fMRI time series.

Methods

In order to first validate this hypothesis on a cohort of patients with known ground truth, here we test individuals with well-defined epileptogenic foci (left mesial temporal lobe epilepsy). We analyzed voxel-wise resting-state fMRI time-series using the autocorrelation function (ACF), an entropic measure of regulation and feedback, and performed follow-up seed-to-voxel functional connectivity analysis. Disruptions in connectivity of the region exhibiting abnormal dynamics were examined in relation to duration of epilepsy and patients’ cognitive performance using a delayed verbal memory recall task.

Results

ACF analysis revealed constrained (less chaotic) functional dynamics in left temporal lobe epilepsy patients, primarily localized to ipsilateral temporal pole, proximal to presumed focal points. Autocorrelation decay rates differentiated, with 100 % accuracy, between patients and healthy controls on a subject-by-subject basis within a leave-one-subject out classification framework. Regions identified via ACF analysis formed a less efficient network in patients, as compared to controls. Constrained dynamics were linked with locally increased and long-range decreased connectivity that, in turn, correlated significantly with impaired memory (local left temporal connectivity) and epilepsy duration (left temporal – posterior cingulate cortex connectivity).

Conclusions

Our current results suggest that data driven functional MRI methods that target network dynamics hold promise in providing clinically valuable tools for identification of epileptic regions.

     

Transition dynamics of EEG-based network microstates during mental arithmetic and resting wakefulness reflects task-related modulations and developmental changes

We studied how maturation influences the organization of functional brain networks engaged during mental calculations and in resting state. Surface EEG measurements from 20 children (8–12 years) and 25 students (21–26 years) were analyzed. Interregional synchronization of brain activity was quantified by means of Phase Lag Index and for various frequency bands. Based on these pairwise estimates of functional connectivity, we formed graphs which were then characterized in terms of local structure [local efficiency (LE)] and overall integration (global efficiency). The overall data analytic scheme was applied twice, in a static and time-varying mode. Our results showed a characteristic trend: functional segregation dominates the network organization of younger brains. Moreover, in childhood, the overall functional network possesses more prominent small-world network characteristics than in early acorrect in xmldulthood in accordance with the Neural Efficiency Hypothesis. The above trends were intensified by the time-varying approach and identified for the whole set of tested frequency bands (from δ to low γ). By mapping the time-indexed connectivity patterns to multivariate timeseries of nodal LE measurements, we carried out an elaborate study of the functional segregation dynamics and demonstrated that the underlying network undergoes transitions between a restricted number of stable states, that can be thought of as “network-level microstates”. The rate of these transitions provided a robust marker of developmental and task-induced alterations, that was found to be insensitive to reference montage and independent component analysis denoising.

Electronic supplementary material

The online version of this article (doi:10.1007/s11571-015-9330-8) contains supplementary material, which is available to authorized users.     

Altered Resting State Brain Dynamics in Temporal Lobe Epilepsy Can Be Observed in Spectral Power,Functional Connectivity and Graph Theory Metrics

Despite a wealth of EEG epilepsy data that accumulated for over half a century, our ability to understand brain dynamics associated with epilepsy remains limited. Using EEG data from 15 controls and 9 left temporal lobe epilepsy (LTLE) patients, in this study we characterize how the dynamics of the healthy brain differ from the “dynamically balanced” state of the brain of epilepsy patients treated with anti-epileptic drugs in the context of resting state. We show that such differences can be observed in band power, synchronization and network measures, as well as deviations from the small world network (SWN) architecture of the healthy brain. The θ (4–7 Hz) and high α (10–13 Hz) bands showed the biggest deviations from healthy controls across various measures. In particular, patients demonstrated significantly higher power and synchronization than controls in the θ band, but lower synchronization and power in the high α band. Furthermore, differences between controls and patients in graph theory metrics revealed deviations from a SWN architecture. In the θ band epilepsy patients showed deviations toward an orderly network, while in the high α band they deviated toward a random network. These findings show that, despite the focal nature of LTLE, the epileptic brain differs in its global network characteristics from the healthy brain. To our knowledge, this is the only study to encompass power, connectivity and graph theory metrics to investigate the reorganization of resting state functional networks in LTLE patients.     

Development of Large-Scale Functional Brain Networks in Children

The ontogeny of large-scale functional organization of the human brain is not well understood. Here we use network analysis of intrinsic functional connectivity to characterize the organization of brain networks in 23 children (ages 7–9 y) and 22 young-adults (ages 19–22 y). Comparison of network properties, including path-length, clustering-coefficient, hierarchy, and regional connectivity, revealed that although children and young-adults' brains have similar “small-world” organization at the global level, they differ significantly in hierarchical organization and interregional connectivity. We found that subcortical areas were more strongly connected with primary sensory, association, and paralimbic areas in children, whereas young-adults showed stronger cortico-cortical connectivity between paralimbic, limbic, and association areas. Further, combined analysis of functional connectivity with wiring distance measures derived from white-matter fiber tracking revealed that the development of large-scale brain networks is characterized by weakening of short-range functional connectivity and strengthening of long-range functional connectivity. Importantly, our findings show that the dynamic process of over-connectivity followed by pruning, which rewires connectivity at the neuronal level, also operates at the systems level, helping to reconfigure and rebalance subcortical and paralimbic connectivity in the developing brain. Our study demonstrates the usefulness of network analysis of brain connectivity to elucidate key principles underlying functional brain maturation, paving the way for novel studies of disrupted brain connectivity in neurodevelopmental disorders such as autism.