Energy landscape of resting magnetoencephalography reveals fronto-parietal network impairments in epilepsy

Reliability of energy landscape analysis of resting-state functional MRI data

Energy landscape analysis is a data-driven method to analyze multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e., within-participant reliability) than across different sets of sessions from different participants (i.e., between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Reliability of energy landscape analysis of resting-state functional MRI data

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Functional and structural networks decoupling in generalized tonic-clonic seizures and its reorganization by drugs

OBJECTIVE

To investigate potential functional and structural large-scale network disturbances in untreated patients with generalized tonic-clonic seizures (GTCS) and the effects of antiseizure drugs.

METHODS

In this study, 41 patients with GTCS, comprising 21 untreated patients and 20 patients who received antiseizure medications (ASMs), and 29 healthy controls were recruited to construct large-scale brain networks based on resting-state functional magnetic resonance imaging and diffusion tensor imaging. Structural and functional connectivity and network-level weighted correlation probability (NWCP) were further investigated to identify network features that corresponded to response to ASMs.

RESULTS

Untreated patients showed more extensive enhancement of functional and structural connections than controls. Specifically, we observed abnormally enhanced connections between the default mode network (DMN) and the frontal-parietal network. In addition, treated patients showed similar functional connection strength to that of the control group. However, all patients exhibited similar structural network alterations. Moreover, the NWCP value was lower for connections within the DMN and between the DMN and other networks in the untreated patients; receiving ASMs could reverse this pattern.

SIGNIFICANCE

Our study identified alterations in structural and functional connectivity in patients with GTCS. The influence of ASMs may be more noticeable within the functional network; moreover, abnormalities in both the functional and structural coupling state may be improved by ASM treatment. Therefore, the coupling state of structural and functional connectivity may be used as an indicator of the efficacy of ASMs.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

Contributions of Magnetoencephalography to Understanding Mechanisms of Generalized Epilepsies: Blurring the Boundary Between Focal and Generalized Epilepsies?

According to the latest operational 2017 ILAE classification of epileptic seizures, the generalized epileptic seizure is still conceptualized as "originating at some point within and rapidly engaging, bilaterally distributed networks." In contrast, the focal epileptic seizure is defined as "originating within networks limited to one hemisphere." Hence, one of the main concepts of "generalized" and "focal" epilepsy comes from EEG descriptions before the era of source localization, and a presumed simultaneous bilateral onset and bi-synchrony of epileptiform discharges remains a hallmark for generalized seizures. Current literature on the pathophysiology of generalized epilepsy supports the concept of a cortical epileptogenic focus triggering rapidly generalized epileptic discharges involving intact corticothalamic and corticocortical networks, known as the cortical focus theory. Likewise, focal epilepsy with rich connectivity can give rise to generalized spike and wave discharges resulting from widespread bilateral synchronization. Therefore, making this key distinction between generalized and focal epilepsy may be challenging in some cases, and for the first time, a combined generalized and focal epilepsy is categorized in the 2017 ILAE classification. Nevertheless, treatment options, such as the choice of antiseizure medications or surgical treatment, are the reason behind the importance of accurate epilepsy classification. Over the past several decades, plentiful scientific research on the pathophysiology of generalized epilepsy has been conducted using non-invasive neuroimaging and postprocessing of the electromagnetic neural signal by measuring the spatiotemporal and interhemispheric latency of bi-synchronous or generalized epileptiform discharges as well as network analysis to identify diagnostic and prognostic biomarkers for accurate diagnosis of the two major types of epilepsy. Among all the advanced techniques, magnetoencephalography (MEG) and multiple other methods provide excellent temporal and spatial resolution, inherently suited to analyzing and visualizing the propagation of generalized EEG activities. This article aims to provide a comprehensive literature review of recent innovations in MEG methodology using source localization and network analysis techniques that contributed to the literature of idiopathic generalized epilepsy in terms of pathophysiology and clinical prognosis, thus further blurring the boundary between focal and generalized epilepsy.

Emotion Recognition Based on Dynamic Energy Features Using a Bi-LSTM Network

Among electroencephalogram (EEG) signal emotion recognition methods based on deep learning, most methods have difficulty in using a high-quality model due to the low resolution and the small sample size of EEG images. To solve this problem, this study proposes a deep network model based on dynamic energy features. In this method, first, to reduce the noise superposition caused by feature analysis and extraction, the concept of an energy sequence is proposed. Second, to obtain the feature set reflecting the time persistence and multicomponent complexity of EEG signals, the construction method of the dynamic energy feature set is given. Finally, to make the network model suitable for small datasets, we used fully connected layers and bidirectional long short-term memory (Bi-LSTM) networks. To verify the effectiveness of the proposed method, we used leave one subject out (LOSO) and 10-fold cross validation (CV) strategies to carry out experiments on the SEED and DEAP datasets. The experimental results show that the accuracy of the proposed method can reach 89.42% (SEED) and 77.34% (DEAP).

Control energy assessment of spatial interactions among macro‐scale brain networks

Many recent studies have revealed that spatial interactions of functional brain networks derived from fMRI data can well model functional connectomes of the human brain. However, it has been rarely explored what the energy consumption characteristics are for such spatial interactions of macro‐scale functional networks, which remains crucial for the understanding of brain organization, behavior, and dynamics. To explore this unanswered question, this article presents a novel framework for quantitative assessment of energy consumptions of macro‐scale functional brain network's spatial interactions via two main effective computational methodologies. First, we designed a novel scheme combining dictionary learning and hierarchical clustering to derive macro‐scale consistent brain network templates that can be used to define a common reference space for brain network interactions and energy assessments. Second, the control energy consumption for driving the brain networks during their spatial interactions is computed from the viewpoint of the linear network control theory. Especially, the energetically favorable brain networks were identified and their energy characteristics were comprehensively analyzed. Experimental results on the Human Connectome Project (HCP) task‐based fMRI (tfMRI) data showed that the proposed methods can reveal meaningful, diverse energy consumption patterns of macro‐scale network interactions. In particular, those networks present remarkable differences in energy consumption. The energetically least favorable brain networks are stable and consistent across HCP tasks such as motor, language, social, and working memory tasks. In general, our framework provides a new perspective to characterize human brain functional connectomes by quantitative assessment for the energy consumption of spatial interactions of macro‐scale brain networks.

The ascending arousal system shapes neural dynamics to mediate awareness of cognitive states

Models of cognitive function typically focus on the cerebral cortex and hence overlook functional links to subcortical structures. This view does not consider the role of the highly-conserved ascending arousal system's role and the computational capacities it provides the brain. We test the hypothesis that the ascending arousal system modulates cortical neural gain to alter the low-dimensional energy landscape of cortical dynamics. Here we use spontaneous functional magnetic resonance imaging data to study phasic bursts in both locus coeruleus and basal forebrain, demonstrating precise time-locked relationships between brainstem activity, low-dimensional energy landscapes, network topology, and spatiotemporal travelling waves. We extend our analysis to a cohort of experienced meditators and demonstrate locus coeruleus-mediated network dynamics were associated with internal shifts in conscious awareness. Together, these results present a view of brain organization that highlights the ascending arousal system's role in shaping both the dynamics of the cerebral cortex and conscious awareness.

Resting-state functional connectivity in the idiopathic generalized epilepsies: A systematic review and meta-analysis of EEG and MEG studies

For idiopathic generalized epilepsies (IGE), brain network analysis is emerging as a biomarker for potential use in clinical care. To determine whether people with IGE show alterations in resting-state brain connectivity compared to healthy controls, and to quantify these differences, we conducted a systematic review and meta-analysis of EEG and magnetoencephalography (MEG) functional connectivity and network studies. The review was conducted according to PRISMA guidelines. Twenty-two studies were eligible for inclusion. Outcomes from individual studies supported hypotheses for interictal, resting-state brain connectivity alterations in IGE patients compared to healthy controls. In contrast, meta-analysis from six studies of common network metrics clustering coefficient, path length, mean degree and nodal strength showed no significant differences between IGE and control groups (effect sizes ranged from -0.151 -1.78). The null findings of the meta-analysis and the heterogeneity of the included studies highlights the importance of developing standardized, validated methodologies for future research. Network neuroscience has significant potential as both a diagnostic and prognostic biomarker in epilepsy, though individual variability in network dynamics needs to be better understood and accounted for.

Characterising Alzheimer's Disease with EEG-based Energy Landscape Analysis

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases, with around 50 million patients worldwide. Accessible and non-invasive methods of diagnosing and characterising AD are therefore urgently required. Electroencephalography (EEG) fulfils these criteria and is often used when studying AD. Several features derived from EEG were shown to predict AD with high accuracy, e.g. signal complexity and synchronisation. However, the dynamics of how the brain transitions between stable states have not been properly studied in the case of AD and EEG. Energy landscape analysis is a method that can be used to quantify these dynamics. This work presents the first application of this method to both AD and EEG. Energy landscape assigns energy value to each possible state, i.e. pattern of activations across brain regions. The energy is inversely proportional to the probability of occurrence. By studying the features of energy landscapes of 20 AD patients and 20 age-matched healthy counterparts (HC), significant differences are found. The dynamics of AD patients' EEG are shown to be more constrained - with more local minima, less variation in basin size, and smaller basins. We show that energy landscapes can predict AD with high accuracy, performing significantly better than baseline models. Moreover, these findings are replicated in a separate dataset including 9 AD and 10 HC above 70 years old.

Alterations of the intrinsic amygdala-hippocampal network in juvenile myoclonic epilepsy

INTRODUCTION

Several lines of evidence suggest that the amygdala-hippocampus is involved in the epileptogenic network of juvenile myoclonic epilepsy (JME). The aim of this study was to investigate the alterations in the individual nuclei of the amygdala and hippocampal subfields, and the intrinsic amygdala-hippocampal network of patients with JME compared to healthy controls.

METHODS

This retrospective study conducted at a single tertiary hospital involved 35 patients with newly diagnosed JME, and 34 healthy subjects. We calculated the individual structural volumes of 18 nuclei in the amygdala, and 38 hippocampal subfields using three-dimensional volumetric T1-weighted imaging and FreeSurfer program. We also performed an analysis of the intrinsic amygdala-hippocampal global and local network based on these volumes using a graph theory and Brain Analysis using Graph Theory (BRAPH) program. We investigated the differences in these volumes and network measures between patients with JME and healthy controls.

RESULTS

There were no significant volume differences in the nuclei of the amygdala and hippocampal subfields between patients with JME and healthy controls. However, we found significant differences in the global network between patients with JME and healthy controls. The mean clustering coefficient was significantly decreased in patients with JME compared to healthy controls (0.473 vs. 0.653, p = .047). In addition, specific regions in the hippocampal subfields showed significant differences in the local network between the two groups. The betweenness centrality of the right CA1-head, right hippocampus-amygdala-transition area, left hippocampal fissure, left fimbria, and left CA3-head, was increased in patients with JME compared to healthy controls.

CONCLUSION

The intrinsic amygdala-hippocampal global and local networks differed in patients with JME compared to healthy controls, which may be related to the pathogenesis of JME, and memory dysfunction in patients with JME.

A computational biomarker of juvenile myoclonic epilepsy from resting-state MEG

•

Computational modelling is combined with MEG to differentiate people with juvenile myoclonic epilepsy from healthy controls.

•

Brain network ictogenicity (BNI) was found higher in people with juvenile myoclonic epilepsy relative to healthy controls.

•

BNI’s classification accuracy in our cohort was 73%.

Objective

For people with idiopathic generalized epilepsy, functional networks derived from their resting-state scalp electrophysiological recordings have shown an inherent higher propensity to generate seizures than those from healthy controls when assessed using the concept of brain network ictogenicity (BNI). Herein we tested whether the BNI framework is applicable to resting-state magnetoencephalography (MEG) from people with juvenile myoclonic epilepsy (JME).

Methods

The BNI framework consists in deriving a functional network from apparently normal brain activity, placing a mathematical model of ictogenicity into the network and then computing how often such network generates seizures in silico. We considered data from 26 people with JME and 26 healthy controls.

Results

We found that resting-state MEG functional networks from people with JME are characterized by a higher propensity to generate seizures (i.e., higher BNI) than those from healthy controls. We found a classification accuracy of 73%.

Conclusions

The BNI framework is applicable to MEG and was capable of differentiating people with epilepsy from healthy controls.

Significance

The BNI framework may be applied to resting-state MEG to aid in epilepsy diagnosis.

Pairwise maximum entropy model explains the role of white matter structure in shaping emergent co-activation states

A major challenge in neuroscience is determining a quantitative relationship between the brain's white matter structural connectivity and emergent activity. We seek to uncover the intrinsic relationship among brain regions fundamental to their functional activity by constructing a pairwise maximum entropy model (MEM) of the inter-ictal activation patterns of five patients with medically refractory epilepsy over an average of ~14 hours of band-passed intracranial EEG (iEEG) recordings per patient. We find that the pairwise MEM accurately predicts iEEG electrodes' activation patterns' probability and their pairwise correlations. We demonstrate that the estimated pairwise MEM's interaction weights predict structural connectivity and its strength over several frequencies significantly beyond what is expected based solely on sampled regions' distance in most patients. Together, the pairwise MEM offers a framework for explaining iEEG functional connectivity and provides insight into how the brain's structural connectome gives rise to large-scale activation patterns by promoting co-activation between connected structures.

Recurrence quantification analysis of dynamic brain networks

Evidence suggests that brain network dynamics are a key determinant of brain function and dysfunction. Here we propose a new framework to assess the dynamics of brain networks based on recurrence analysis. Our framework uses recurrence plots and recurrence quantification analysis to characterize dynamic networks. For resting-state magnetoencephalographic dynamic functional networks (dFNs), we have found that functional networks recur more quickly in people with epilepsy than in healthy controls. This suggests that recurrence of dFNs may be used as a biomarker of epilepsy. For stereo electroencephalography data, we have found that dFNs involved in epileptic seizures emerge before seizure onset, and recurrence analysis allows us to detect seizures. We further observe distinct dFNs before and after seizures, which may inform neurostimulation strategies to prevent seizures. Our framework can also be used for understanding dFNs in healthy brain function and in other neurological disorders besides epilepsy.