Simple and difficult mathematics in children: A minimum spanning tree EEG network analysis

Comparison of data-driven thresholding methods using directed functional brain networks

Over the past two centuries, intensive empirical research has been conducted on the human brain. As an electroencephalogram (EEG) records millisecond-to-millisecond changes in the electrical potentials of the brain, it has enormous potential for identifying useful information about neuronal transactions. The EEG data can be modelled as graphs by considering the electrode sites as nodes and the linear and nonlinear statistical dependencies among them as edges (with weights). The graph theoretical modelling of EEG data results in functional brain networks (FBNs), which are fully connected (complete) weighted undirected/directed networks. Since various brain regions are interconnected via sparse anatomical connections, the weak links can be filtered out from the fully connected networks using a process called thresholding. Multiple researchers in the past decades proposed many thresholding methods to gather more insights about the influential neuronal connections of FBNs. This paper reviews various thresholding methods used in the literature for FBN analysis. The analysis showed that data-driven methods are unbiased since no arbitrary user-specified threshold is required. The efficacy of four data-driven thresholding methods, namely minimum spanning tree (MST), minimum connected component (MCC), union of shortest path trees (USPT), and orthogonal minimum spanning tree (OMST), in characterizing cognitive behavior of the normal human brain is analysed using directed FBNs constructed from EEG data of different cognitive load states. The experimental results indicate that both MCC and OMST thresholding methods can detect cognitive load-induced changes in the directed functional brain networks.

Potential long consequences from internal and external ecology: loss of gut microbiota antifragility in children from an industrialized population compared with an indigenous rural lifestyle

Human health is strongly mediated by the gut microbiota ecosystem, which, in turn, depends not only on its state but also on its dynamics and how it responds to perturbations. Healthy microbiota ecosystems tend to be in criticality and antifragile dynamics corresponding to a maximum complexity configuration, which may be assessed with information and network theory analysis. Under this complex system perspective, we used a new analysis of published data to show that a children's population with an industrialized urban lifestyle from Mexico City exhibits informational and network characteristics similar to parasitized children from a rural indigenous population in the remote mountainous region of Guerrero, México. We propose then, that in this critical age for gut microbiota maturation, the industrialized urban lifestyle could be thought of as an external perturbation to the gut microbiota ecosystem, and we show that it produces a similar loss in criticality/antifragility as the one observed by internal perturbation due to parasitosis by the helminth A. lumbricoides. Finally, several general complexity-based guidelines to prevent or restore gut ecosystem antifragility are discussed.

Spontaneous brain activity, graph metrics, and head motion related to prospective post-traumatic stress disorder trauma-focused therapy response

Introduction

Trauma-focused psychotherapy for post-traumatic stress disorder (PTSD) is effective in about half of all patients. Investigating biological systems related to prospective treatment response is important to gain insight in mechanisms predisposing patients for successful intervention. We studied if spontaneous brain activity, brain network characteristics and head motion during the resting state are associated with future treatment success.

Methods

Functional magnetic resonance imaging scans were acquired from 46 veterans with PTSD around the start of treatment. Psychotherapy consisted of trauma-focused cognitive behavioral therapy (tf-CBT), eye movement desensitization and reprocessing (EMDR), or a combination thereof. After intervention, 24 patients were classified as treatment responders and 22 as treatment resistant. Differences between groups in spontaneous brain activity were evaluated using amplitude of low-frequency fluctuations (ALFF), while global and regional brain network characteristics were assessed using a minimum spanning tree (MST) approach. In addition, in-scanner head motion was assessed.

Results

No differences in spontaneous brain activity and global network characteristics were observed between the responder and non-responder group. The right inferior parietal lobule, right putamen and left superior parietal lobule had a more central position in the network in the responder group compared to the non-responder group, while the right dorsolateral prefrontal cortex (DLPFC), right inferior frontal gyrus and left inferior temporal gyrus had a less central position. In addition, responders showed less head motion.

Discussion

These results show that areas involved in executive functioning, attentional and action processes, learning, and visual-object processing, are related to prospective PTSD treatment response in veterans. In addition, these findings suggest that involuntary micromovements may be related to future treatment success.

Altered brain network topology during speech tracking in developmental dyslexia

Developmental dyslexia is often accompanied by altered phonological processing of speech. Underlying neural changes have typically been characterized in terms of stimulus- and/or task-related responses within individual brain regions or their functional connectivity. Less is known about potential changes in the more global functional organization of brain networks. Here we recorded electroencephalography (EEG) in typical and dyslexic readers while they listened to (a) a random sequence of syllables and (b) a series of tri-syllabic real words. The network topology of the phase synchronization of evoked cortical oscillations was investigated in four frequency bands (delta, theta, alpha and beta) using minimum spanning tree graphs. We found that, compared to syllable tracking, word tracking triggered a shift toward a more integrated network topology in the theta band in both groups. Importantly, this change was significantly stronger in the dyslexic readers, who also showed increased reliance on a right frontal cluster of electrodes for word tracking. The current findings point towards an altered effect of word-level processing on the functional brain network organization that may be associated with less efficient phonological and reading skills in dyslexia.

Aberrant dynamic minimal spanning tree parameters within default mode network in patients with autism spectrum disorder

The altered functional connectivity (FC) level and its temporal characteristics within certain cortical networks, such as the default mode network (DMN), could provide a possible explanatory framework for Autism spectrum disorder (ASD). In the current study, we hypothesized that the topographical organization along with its temporal dynamics of the autistic brain measured by temporal mean and variance of complex network measures, respectively, were significantly altered, which may further explain the autistic symptom severity in patients with ASD. To validate these hypotheses, the precise FCs between DMN regions at each time point were calculated using the resting-state functional magnetic resonance imaging (fMRI) datasets from the Autism Brain Imaging Data Exchange (ABIDE) project. Then, the minimal spanning tree (MST) technique was applied to construct a time-varying complex network of DMN. By analyzing the temporal mean and variance of MST parameters and their relationship with autistic symptom severity, we found that in persons with ASD, the information exchange efficiencies between cortical regions within DMN were significantly lower and more volatile compared with those in typical developing participants. Moreover, these alterations within DMN were closely associated with the autistic symptom severity of the ASD group.

Resting-state EEG reveals global network deficiency in prelingually deaf children with late cochlear implantation

There are individual differences in rehabilitation after cochlear implantation that can be explained by brain plasticity. However, from the perspective of brain networks, the effect of implantation age on brain plasticity is unclear. The present study investigated electroencephalography functional networks in the resting state, including eyes-closed and eyes-open conditions, in 31 children with early cochlear implantation, 24 children with late cochlear implantation, and 29 children with normal hearing. Resting-state functional connectivity was measured with phase lag index, and we investigated the connectivity between the sensory regions for each frequency band. Network topology was examined using minimum spanning tree to obtain the network backbone characteristics. The results showed stronger connectivity between auditory and visual regions but reduced global network efficiency in children with late cochlear implantation in the theta and alpha bands. Significant correlations were observed between functional backbone characteristics and speech perception scores in children with cochlear implantation. Collectively, these results reveal an important effect of implantation age on the extent of brain plasticity from a network perspective and indicate that characteristics of the brain network can reflect the extent of rehabilitation of children with cochlear implantation.

Graph Analysis of EEG Functional Connectivity Networks During a Letter-Speech Sound Binding Task in Adult Dyslexics

We performed an EEG graph analysis on data from 31 typical readers (22.27 ± 2.53 y/o) and 24 dyslexics (22.99 ± 2.29 y/o), recorded while they were engaged in an audiovisual task and during resting-state. The task simulates reading acquisition as participants learned new letter-sound mappings via feedback. EEG data was filtered for the delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), and beta (13–30 Hz) bands. We computed the Phase Lag Index (PLI) to provide an estimate of the functional connectivity between all pairs of electrodes per band. Then, networks were constructed using a Minimum Spanning Tree (MST), a unique sub-graph connecting all nodes (electrodes) without loops, aimed at minimizing bias in between groups and conditions comparisons. Both groups showed a comparable accuracy increase during task blocks, indicating that they correctly learned the new associations. The EEG results revealed lower task-specific theta connectivity, and lower theta degree correlation over both rest and task recordings, indicating less network integration in dyslexics compared to typical readers. This pattern suggests a role of theta oscillations in dyslexia and may reflect differences in task engagement between the groups, although robust correlations between MST metrics and performance indices were lacking.

Frontal Functional Network Disruption Associated with Amyotrophic Lateral Sclerosis: An fNIRS-Based Minimum Spanning Tree Analysis

Recent evidence increasingly associates network disruption in brain organization with multiple neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS), a rare terminal disease. However, the comparability of brain network characteristics across different studies remains a challenge for conventional graph theoretical methods. One suggested method to address this issue is minimum spanning tree (MST) analysis, which provides a less biased comparison. Here, we assessed the novel application of MST network analysis to hemodynamic responses recorded by functional near-infrared spectroscopy (fNIRS) neuroimaging modality, during an activity-based paradigm to investigate hypothetical disruptions in frontal functional brain network topology as a marker of the executive dysfunction, one of the most prevalent cognitive deficit reported across ALS studies. We analyzed data recorded from nine participants with ALS and ten age-matched healthy controls by first estimating functional connectivity, using phase-locking value (PLV) analysis, and then constructing the corresponding individual and group MSTs. Our results showed significant between-group differences in several MST topological properties, including leaf fraction, maximum degree, diameter, eccentricity, and degree divergence. We further observed a global shift toward more centralized frontal network organizations in the ALS group, interpreted as a more random or dysregulated network in this cohort. Moreover, the similarity analysis demonstrated marginally significantly increased overlap in the individual MSTs from the control group, implying a reference network with lower topological variation in the healthy cohort. Our nodal analysis characterized the main local hubs in healthy controls as distributed more evenly over the frontal cortex, with slightly higher occurrence in the left prefrontal cortex (PFC), while in the ALS group, the most frequent hubs were asymmetrical, observed primarily in the right prefrontal cortex. Furthermore, it was demonstrated that the global PLV (gPLV) synchronization metric is associated with disease progression, and a few topological properties, including leaf fraction and tree hierarchy, are linked to disease duration. These results suggest that dysregulation, centralization, and asymmetry of the hemodynamic-based frontal functional network during activity are potential neuro-topological markers of ALS pathogenesis. Our findings can possibly support new bedside assessments of the functional status of ALS' brain network and could hypothetically extend to applications in other neurodegenerative diseases.

EEG based emotion recognition using minimum spanning tree

Emotion is a fundamental factor that influences human cognition, motivation, decision making and social interactions. This psychological state arises spontaneously and goes with physiological changes that can be recognized by computational methods. In this study, changes in minimum spanning tree (MST) structure of brain functional connectome were used for emotion classification based on EEG data and the obtained results were employed for interpretation about the most informative frequency content of emotional states. For estimation of interaction between different brain regions, several connectivity metrics were applied and interactions were calculated in different frequency bands. Subsequently, the MST graph was extracted from the functional connectivity matrix and its features were used for emotion recognition. The results showed that the accuracy of the proposed method for separating emotions with different arousal levels was 88.28%, while for different valence levels it was 81.25%. Interestingly, the system performance for binary classification of emotions based on quadrants of arousal-valence space was also higher than 80%. The MST approach allowed us to study the change of brain complexity and dynamics in various emotional states. This capability provided us enough knowledge to claim lower-alpha and gamma bands contain the main information for discrimination of emotional states.

Androgyny eliminates sex differences in emotional reactivity: ERP and network coupling evidences

The role of gender role in interpreting sex differences in emotion is unknown. The present study examined how gender role moderates sex differences in emotional reactivity. Event-related potentials (ERP) were recorded in sixty-eight subjects with typical or androgynous gender roles when they passively observed neutral and negative pictures. Behaviorally, typical females (feminine females) reported higher emotional rating than typical males (masculine males), while androgynous males and androgynous females reported no significant differences. Electrophysiologically, we found higher late positive potential (LPP) amplitude in typical females compared to typical males, while this pattern of sex difference was absent in androgynous subjects. The network analysis of EEG data indicates that typical males showed enhanced network coupling strengths between frontal/prefrontal and parietal areas than typical females, which was again absent in androgynous subjects. These findings suggest that gender role is an important determinant in the interpretation of sex differences in emotional reactivity.

The power of language: Functional brain network topology of deaf and hearing in relation to sign language experience

Prolonged auditory sensory deprivation leads to brain reorganization. This is indicated by functional enhancement in remaining sensory systems and known as cross-modal plasticity. In this study we investigated differences in functional brain network topology between deaf and hearing individuals. We also studied altered functional network responses between deaf and hearing individuals with a recording paradigm containing an eyes-closed and eyes-open condition. Electroencephalography activity was recorded in a group of sign language-trained deaf (N = 71) and hearing people (N = 122) living in rural Africa. Functional brain networks were constructed from the functional connectivity between fourteen electrodes distributed over the scalp. Functional connectivity was quantified with the phase lag index based on bandpass filtered epochs of brain signal. We studied the functional connectivity between the auditory, somatosensory and visual cortex and performed whole-brain minimum spanning tree analysis to capture network backbone characteristics. Functional connectivity between different regions involved in sensory information processing tended to be stronger in deaf people during the eyes-closed condition in both the alpha and beta frequency band. Furthermore, we found differences in functional backbone topology between deaf and hearing individuals. The backbone topology altered during transition from the eyes-closed to eyes-open condition irrespective of deafness, but was more pronounced in deaf individuals. The transition of backbone strength was different between individuals with congenital, pre-lingual or post-lingual deafness. Functional backbone characteristics correlated with the experience of sign language. Overall, our study revealed more insights in functional network reorganization caused by auditory deprivation and cross-modal plasticity. It further supports the idea of a brain plasticity potential in deaf and hearing people. The association between network organization and acquired sign language experience reflects the ability of ongoing brain adaptation in people with hearing disabilities.

EEG Resting State Functional Connectivity in Adult Dyslexics Using Phase Lag Index and Graph Analysis

Developmental dyslexia may involve deficits in functional connectivity across widespread brain networks that enable fluent reading. We investigated the large-scale organization of electroencephalography (EEG) functional networks at rest in 28 dyslexics and 36 typically reading adults. For each frequency band (delta, theta alpha and beta), we assessed functional connectivity strength with the phase lag index (PLI). Network topology was examined using minimum spanning tree (MST) graphs derived from the functional connectivity matrices. We found significant group differences in the alpha band (8-13 Hz). The graph analysis indicated more interconnected nodes, in dyslexics compared to typical readers. The graph metrics were significantly correlated with age in dyslexics but not in typical readers, which may indicate more heterogeneity in maturation of brain networks in dyslexics. The present findings support the involvement of alpha oscillations in higher cognition and the sensitivity of graph metrics to characterize functional networks in adult dyslexia. Finally, the current results extend our previous findings on children.

A comparison between scalp- and source-reconstructed EEG networks

EEG can be used to characterise functional networks using a variety of connectivity (FC) metrics. Unlike EEG source reconstruction, scalp analysis does not allow to make inferences about interacting regions, yet this latter approach has not been abandoned. Although the two approaches use different assumptions, conclusions drawn regarding the topology of the underlying networks should, ideally, not depend on the approach. The aim of the present work was to find an answer to the following questions: does scalp analysis provide a correct estimate of the network topology? how big are the distortions when using various pipelines in different experimental conditions? EEG recordings were analysed with amplitude- and phase-based metrics, founding a strong correlation for the global connectivity between scalp- and source-level. In contrast, network topology was only weakly correlated. The strongest correlations were obtained for MST leaf fraction, but only for FC metrics that limit the effects of volume conduction/signal leakage. These findings suggest that these effects alter the estimated EEG network organization, limiting the interpretation of results of scalp analysis. Finally, this study also suggests that the use of metrics that address the problem of zero lag correlations may give more reliable estimates of the underlying network topology.

Arithmetic in the developing brain: A review of brain imaging studies

Brain imaging studies on academic achievement offer an exciting window on experience-dependent cortical plasticity, as they allow us to understand how developing brains change when children acquire culturally transmitted skills. This contribution focuses on the learning of arithmetic, which is quintessential to mathematical development. The nascent body of brain imaging studies reveals that arithmetic recruits a large set of interconnected areas, including prefrontal, posterior parietal, occipito-temporal and hippocampal areas. This network undergoes developmental changes in its function, connectivity and structure, which are not yet fully understood. This network only partially overlaps with what has been found in adults, and clear differences are observed in the recruitment of the hippocampus, which are related to the development of arithmetic fact retrieval. Despite these emerging trends, the literature remains scattered, particularly in the context of atypical development. Acknowledging the distributed nature of the arithmetic network, future studies should focus on connectivity and analytic approaches that investigate patterns of brain activity, coupled with a careful design of the arithmetic tasks and assessments of arithmetic strategies. Such studies will produce a more comprehensive understanding of how the arithmetical brain unfolds, how it changes over time, and how it is impaired in atypical development.

Topological Filtering of Dynamic Functional Brain Networks Unfolds Informative Chronnectomics: A Novel Data-Driven Thresholding Scheme Based on Orthogonal Minimal Spanning Trees (OMSTs)

The human brain is a large-scale system of functionally connected brain regions. This system can be modeled as a network, or graph, by dividing the brain into a set of regions, or "nodes," and quantifying the strength of the connections between nodes, or "edges," as the temporal correlation in their patterns of activity. Network analysis, a part of graph theory, provides a set of summary statistics that can be used to describe complex brain networks in a meaningful way. The large-scale organization of the brain has features of complex networks that can be quantified using network measures from graph theory. The adaptation of both bivariate (mutual information) and multivariate (Granger causality) connectivity estimators to quantify the synchronization between multichannel recordings yields a fully connected, weighted, (a)symmetric functional connectivity graph (FCG), representing the associations among all brain areas. The aforementioned procedure leads to an extremely dense network of tens up to a few hundreds of weights. Therefore, this FCG must be filtered out so that the "true" connectivity pattern can emerge. Here, we compared a large number of well-known topological thresholding techniques with the novel proposed data-driven scheme based on orthogonal minimal spanning trees (OMSTs). OMSTs filter brain connectivity networks based on the optimization between the global efficiency of the network and the cost preserving its wiring. We demonstrated the proposed method in a large EEG database (N = 101 subjects) with eyes-open (EO) and eyes-closed (EC) tasks by adopting a time-varying approach with the main goal to extract features that can totally distinguish each subject from the rest of the set. Additionally, the reliability of the proposed scheme was estimated in a second case study of fMRI resting-state activity with multiple scans. Our results demonstrated clearly that the proposed thresholding scheme outperformed a large list of thresholding schemes based on the recognition accuracy of each subject compared to the rest of the cohort (EEG). Additionally, the reliability of the network metrics based on the fMRI static networks was improved based on the proposed topological filtering scheme. Overall, the proposed algorithm could be used across neuroimaging and multimodal studies as a common computationally efficient standardized tool for a great number of neuroscientists and physicists working on numerous of projects.

Comparison of network analysis approaches on EEG connectivity in beta during Visual Short-term Memory binding tasks

We analyse the electroencephalogram signals in the beta band of working memory representation recorded from young healthy volunteers performing several different Visual Short-Term Memory (VSTM) tasks which have proven useful in the assessment of clinical and preclinical Alzheimer's disease. We compare network analysis using Maximum Spanning Trees (MSTs) with network analysis obtained using 20% and 25% connection thresholds on the VSTM data. MSTs are a promising method of network analysis negating the more classical use of thresholds which are so far chosen arbitrarily. However, we find that the threshold analyses outperforms MSTs for detection of functional network differences. Particularly, MSTs fail to find any significant differences. Further, the thresholds detect significant differences between shape and shape-colour binding tasks when these are tested in the left side of the display screen, but no such differences are detected when these tasks are tested for in the right side of the display screen. This provides evidence that contralateral activity is a significant factor in sensitivity for detection of cognitive task differences.