Probing to Observe Neural Dynamics Investigated with Networked Kuramoto Oscillators

A model integrating multiple processes of synchronization and coherence for information instantiation within a cortical area

What is the form of dynamic, e.g., sensory, information in the mammalian cortex? Information in the cortex is modeled as a coherence map of a mixed chimera state of synchronous, phasic, and disordered minicolumns. The theoretical model is built on neurophysiological evidence. Complex spatiotemporal information is instantiated through a system of interacting biological processes that generate a synchronized cortical area, a coherent aperture. Minicolumn elements are grouped in macrocolumns in an array analogous to a phased-array radar, modeled as an aperture, a "hole through which radiant energy flows." Coherence maps in a cortical area transform inputs from multiple sources into outputs to multiple targets, while reducing complexity and entropy. Coherent apertures can assume extremely large numbers of different information states as coherence maps, which can be communicated among apertures with corresponding very large bandwidths. The coherent aperture model incorporates considerable reported research, integrating five conceptually and mathematically independent processes: 1) a damped Kuramoto network model, 2) a pumped area field potential, 3) the gating of nearly coincident spikes, 4) the coherence of activity across cortical lamina, and 5) complex information formed through functions in macrocolumns. Biological processes and their interactions are described in equations and a functional circuit such that the mathematical pieces can be assembled the same way the neurophysiological ones are. The model can be conceptually convolved over the specifics of local cortical areas within and across species. A coherent aperture becomes a node in a graph of cortical areas with a corresponding distribution of information.

Active probing to highlight approaching transitions to ictal states in coupled neural mass models

The extraction of electrophysiological features that reliably forecast the occurrence of seizures is one of the most challenging goals in epilepsy research. Among possible approaches to tackle this problem is the use of active probing paradigms in which responses to stimuli are used to detect underlying system changes leading up to seizures. This work evaluates the theoretical and mechanistic underpinnings of this strategy using two coupled populations of the well-studied Wendling neural mass model. Different model settings are evaluated, shifting parameters (excitability, slow inhibition, or inter-population coupling gains) from normal towards ictal states while probing stimuli are applied every 2 seconds to the input of either one or both populations. The correlation between the extracted features and the ictogenic parameter shifting indicates if the impending transition to the ictal state may be identified in advance. Results show that not only can the response to the probing stimuli forecast seizures but this is true regardless of the altered ictogenic parameter. That is, similar feature changes are highlighted by probing stimuli responses in advance of the seizure including: increased response variance and lag-1 autocorrelation, decreased skewness, and increased mutual information between the outputs of both model subsets. These changes were mostly restricted to the stimulated population, showing a local effect of this perturbational approach. The transition latencies from normal activity to sustained discharges of spikes were not affected, suggesting that stimuli had no pro-ictal effects. However, stimuli were found to elicit interictal-like spikes just before the transition to the ictal state. Furthermore, the observed feature changes highlighted by probing the neuronal populations may reflect the phenomenon of critical slowing down, where increased recovery times from perturbations may signal the loss of a systems' resilience and are common hallmarks of an impending critical transition. These results provide more evidence that active probing approaches highlight information about underlying system changes involved in ictogenesis and may be able to play a role in assisting seizure forecasting methods which can be incorporated into early-warning systems that ultimately enable closing the loop for targeted seizure-controlling interventions.

Controllability of Networks of Multiple Coupled Neural Populations: An Analytical Method for Neuromodulation's Feasibility

Neuromodulation plays a vital role in the prevention and treatment of neurological and psychiatric disorders. Neuromodulation's feasibility is a long-standing issue because it provides the necessity for neuromodulation to realize the desired purpose. A controllability analysis of neural dynamics is necessary to ensure neuromodulation's feasibility. Here, we present such a theoretical method by using the concept of controllability from the control theory that neuromodulation's feasibility can be studied smoothly. Firstly, networks of multiple coupled neural populations with different topologies are established to mathematically model complicated neural dynamics. Secondly, an analytical method composed of a linearization method, the Kalman controllable rank condition and a controllability index is applied to analyze the controllability of the established network models. Finally, the relationship between network dynamics or topological characteristic parameters and controllability is studied by using the analytical method. The proposed method provides a new idea for the study of neuromodulation's feasibility, and the results are expected to guide us to better modulate neurodynamics by optimizing network dynamics and network topology.

Multiplex Limited Penetrable Horizontal Visibility Graph from EEG Signals for Driver Fatigue Detection

Driver fatigue is an important contributor to road accidents, and driver fatigue detection has attracted a great deal of attention on account of its significant importance. Numerous methods have been proposed to fulfill this challenging task, though, the characterization of the fatigue mechanism still, to a large extent, remains to be investigated. To address this problem, we, in this work, develop a novel Multiplex Limited Penetrable Horizontal Visibility Graph (Multiplex LPHVG) method, which allows in not only detecting fatigue driving but also probing into the brain fatigue behavior. Importantly, we use the method to construct brain networks from EEG signals recorded from different subjects performing simulated driving tasks under alert and fatigue driving states. We then employ clustering coefficient, global efficiency and characteristic path length to characterize the topological structure of the networks generated from different brain states. In addition, we combine average edge overlap with the network measures to distinguish alert and mental fatigue states. The high-accurate classification results clearly demonstrate and validate the efficacy of our multiplex LPHVG method for the fatigue detection from EEG signals. Furthermore, our findings show a significant increase of the clustering coefficient as the brain evolves from alert state to mental fatigue state, which yields novel insights into the brain behavior associated with fatigue driving.

Strong-Weak Pruning for Brain Network Identification in Connectome-Wide Neuroimaging: Application to Amyotrophic Lateral Sclerosis Disease Stage Characterization

Magnetic resonance imaging allows acquiring functional and structural connectivity data from which high-density whole-brain networks can be derived to carry out connectome-wide analyses in normal and clinical populations. Graph theory has been widely applied to investigate the modular structure of brain connections by using centrality measures to identify the "hub" of human connectomes, and community detection methods to delineate subnetworks associated with diverse cognitive and sensorimotor functions. These analyses typically rely on a preprocessing step (pruning) to reduce computational complexity and remove the weakest edges that are most likely affected by experimental noise. However, weak links may contain relevant information about brain connectivity, therefore, the identification of the optimal trade-off between retained and discarded edges is a subject of active research. We introduce a pruning algorithm to identify edges that carry the highest information content. The algorithm selects both strong edges (i.e. edges belonging to shortest paths) and weak edges that are topologically relevant in weakly connected subnetworks. The newly developed "strong-weak" pruning (SWP) algorithm was validated on simulated networks that mimic the structure of human brain networks. It was then applied for the analysis of a real dataset of subjects affected by amyotrophic lateral sclerosis (ALS), both at the early (ALS2) and late (ALS3) stage of the disease, and of healthy control subjects. SWP preprocessing allowed identifying statistically significant differences in the path length of networks between patients and healthy subjects. ALS patients showed a decrease of connectivity between frontal cortex to temporal cortex and parietal cortex and between temporal and occipital cortex. Moreover, degree of centrality measures revealed significantly different hub and centrality scores between patient subgroups. These findings suggest a widespread alteration of network topology in ALS associated with disease progression.

Empirical Functional PCA for 3D Image Feature Extraction Through Fractal Sampling

Medical image classification is currently a challenging task that can be used to aid the diagnosis of different brain diseases. Thus, exploratory and discriminative analysis techniques aiming to obtain representative features from the images play a decisive role in the design of effective Computer Aided Diagnosis (CAD) systems, which is especially important in the early diagnosis of dementia. In this work, we present a technique that allows using specific time series analysis techniques with 3D images. This is achieved by sampling the image using a fractal-based method which preserves the spatial relationship among voxels. In addition, a method called Empirical functional PCA (EfPCA) is presented, which combines Empirical Mode Decomposition (EMD) with functional PCA to express an image in the space spanned by a basis of empirical functions, instead of using components computed by a predefined basis as in Fourier or Wavelet analysis. The devised technique has been used to classify images from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) and the Parkinson Progression Markers Initiative (PPMI), achieving accuracies up to 93% and 92% differential diagnosis tasks (AD versus controls and PD versus Controls, respectively). The results obtained validate the method, proving that the information retrieved by our methodology is significantly linked to the diseases.

Differential Effects of Simulated Cortical Network Lesions on Synchrony and EEG Complexity

Brain function has been proposed to arise as a result of the coordinated activity between distributed brain areas. An important issue in the study of brain activity is the characterization of the synchrony among these areas and the resulting complexity of the system. However, the variety of ways to define and, hence, measure brain synchrony and complexity has sometimes led to inconsistent results. Here, we study the relationship between synchrony and commonly used complexity estimators of electroencephalogram (EEG) activity and we explore how simulated lesions in anatomically based cortical networks would affect key functional measures of activity. We explored this question using different types of neural network lesions while the brain dynamics was modeled with a time-delayed set of 66 Kuramoto oscillators. Each oscillator modeled a region of the cortex (node), and the connectivity and spatial location between different areas informed the creation of a network structure (edges). Each type of lesion consisted on successive lesions of nodes or edges during the simulation of the neural dynamics. For each type of lesion, we measured the synchrony among oscillators and three complexity estimators (Higuchi's Fractal Dimension, Sample Entropy and Lempel-Ziv Complexity) of the simulated EEGs. We found a general negative correlation between EEG complexity metrics and synchrony but Sample Entropy and Lempel-Ziv showed a positive correlation with synchrony when the edges of the network were deleted. This suggests an intricate relationship between synchrony of the system and its estimated complexity. Hence, complexity seems to depend on the multiple states of interaction between the oscillators of the system. Our results can contribute to the interpretation of the functional meaning of EEG complexity.