Nonlinear dynamical analysis of EEG and MEG: review of an emerging field

Evaluating and Comparing Measures of Aperiodic Neural Activity

Neuro-electrophysiological recordings contain prominent aperiodic activity - meaning irregular activity, with no characteristic frequency - which has variously been referred to as 1/f (or 1/f-like activity), fractal, or 'scale-free' activity. Previous work has established that aperiodic features of neural activity is dynamic and variable, relating (between subjects) to healthy aging and to clinical diagnoses, and also (within subjects) tracking conscious states and behavioral performance. There are, however, a wide variety of conceptual frameworks and associated methods for the analyses and interpretation of aperiodic activity - for example, time domain measures such as the autocorrelation, fractal measures, and/or various complexity and entropy measures, as well as measures of the aperiodic exponent in the frequency domain. There is a lack of clear understanding of how these different measures relate to each other and to what extent they reflect the same or different properties of the data, which makes it difficult to synthesize results across approaches and complicates our overall understanding of the properties, biological significance, and demographic, clinical, and behavioral correlates of aperiodic neural activity. To address this problem, in this project we systematically survey the different approaches for measuring aperiodic neural activity, starting with an automated literature analysis to curate a collection of the most common methods. We then evaluate and compare these methods, using statistically representative time series simulations. In doing so, we establish consistent relationships between the measures, showing that much of what they capture reflects shared variance - though with some notable idiosyncrasies. Broadly, frequency domain methods are more specific to aperiodic features of the data, whereas time domain measures are more impacted by oscillatory activity. We extend this analysis by applying the measures to a series of empirical EEG and iEEG datasets, replicating the simulation results. We conclude by summarizing the relationships between the multiple methods, emphasizing opportunities for re-examining previous findings and for future work.

Quantifying instability in neurological disorders EEG based on phase space DTM function

Classifying individuals with neurological disorders and healthy subjects using EEG is a crucial area of research. The current feature extraction approach focuses on the frequency domain features in each of the EEG frequency bands and functional brain networks. In recent years, researchers have discovered and extensively studied stability differences in the electroencephalograms (EEG) of patients with neurological disorders. Based on this, this paper proposes a feature descriptor to characterize EEG instability. The proposed method starts by forming a signal point cloud through Phase Space Reconstruction (PSR). Subsequently, a pseudo-metric space is constructed, and pseudo-distances are calculated based on the consistent measure of the point cloud. Finally, Distance to Measure (DTM) Function are generated to replace the distance function in the original metric space. We calculated the relative distances in the point cloud by measuring signal similarity and, based on this, summarized the point cloud structures formed by EEG with different stabilities after PSR. This process demonstrated that Multivariate Kernel Density Estimation (MKDE) based on a Gaussian kernel can effectively separate the mappings of different stable components within the signal in the phase space. The two average DTM values are then proposed as feature descriptors for EEG instability.In the validation phase, the proposed feature descriptor is tested on three typical neurological disorders: epilepsy, Alzheimer's disease, and Parkinson's disease, using the Bonn dataset, CHB-MIT, the Florida State University dataset, and the Iowa State University dataset. DTM values are used as feature inputs for four different machine learning classifiers, and The results show that the best classification accuracy of the proposed method reaches 98.00 %, 96.25 %, 96.71 % and 95.34 % respectively, outperforming commonly used nonlinear descriptors. Finally, the proposed method is tested and analyzed using noisy signals, demonstrating its robustness compared to other methods.

Application of HFO and scaling analysis of neuronal oscillations in the presurgical evaluation of focal epilepsy

PURPOSE

To explore the utility of high frequency oscillations (HFO) and long-range temporal correlations (LRTCs) in preoperative assessment of epilepsy.

METHODS

MEG ripples were detected in 59 drug-resistant epilepsy patients, comprising 5 with parietal lobe epilepsy (PLE), 21 with frontal lobe epilepsy (FLE), 14 with lateral temporal lobe epilepsy (LTLE), and 19 with mesial temporal lobe epilepsy (MTLE) to identify the epileptogenic zone (EZ). The results were compared with clinical MEG reports and resection area. Subsequently, LRTCs were quantified at the source-level by detrended fluctuation analysis (DFA) and life/waiting -time at 5 bands for 90 cerebral cortex regions. The brain regions with larger DFA exponents and standardized life-waiting biomarkers were compared with the resection results.

RESULTS

Compared to MEG sensor-level data, ripple sources were more frequently localized within the resection area. Moreover, source-level analysis revealed a higher proportion of DFA exponents and life-waiting biomarkers with relatively higher rankings, primarily distributed within the resection area (p<0.01). Moreover, these two LRCT indices across five distinct frequency bands correlated with EZ.

CONCLUSION

HFO and source-level LRTCs are correlated with EZ. Integrating HFO and LRTCs may be an effective approach for presurgical evaluation of epilepsy.

S3T-Net: A novel electroencephalogram signals-oriented emotion recognition model

In this paper, a novel skipping spatial-spectral-temporal network (S3T-Net) is developed to handle intra-individual differences in electroencephalogram (EEG) signals for accurate, robust, and generalized emotion recognition. In particular, aiming at the 4D features extracted from the raw EEG signals, a multi-branch architecture is proposed to learn spatial-spectral cross-domain representations, which benefits enhancing the model generalization ability. Time dependency among different spatial-spectral features is further captured via a bi-directional long-short term memory module, which employs an attention mechanism to integrate context information. Moreover, a skip-change unit is designed to add another auxiliary pathway for updating model parameters, which alleviates the vanishing gradient problem in complex spatial-temporal network. Evaluation results show that the proposed S3T-Net outperforms other advanced models in terms of the emotion recognition accuracy, which yields an performance improvement of 0.23% , 0.13%, and 0.43% as compared to the sub-optimal model in three test scenes, respectively. In addition, the effectiveness and superiority of the key components of S3T-Net are demonstrated from various experiments. As a reliable and competent emotion recognition model, the proposed S3T-Net contributes to the development of intelligent sentiment analysis in human-computer interaction (HCI) realm.

Assessment of Fractal Synchronization during an Epileptic Seizure

In this paper, we define fractal synchronization (FS) based on the idea of stochastic synchronization and propose a mathematical apparatus for estimating FS. One major advantage of our proposed approach is that fractal synchronization makes it possible to estimate the aggregate strength of the connection on multiple time scales between two projections of the attractor, which are time series with a fractal structure. We believe that one of the promising uses of FS is the assessment of the interdependence of encephalograms. To demonstrate this approach in evaluating the cross-dependence between channels in a network of electroencephalograms, we evaluated the FS of encephalograms during an epileptic seizure. Fractal synchronization demonstrates the presence of desynchronization during an epileptic seizure.

Estimation of Carleman operator from a univariate time series

Reconstructing a nonlinear dynamical system from empirical time series is a fundamental task in data-driven analysis. One of the main challenges is the existence of hidden variables; we only have records for some variables, and those for hidden variables are unavailable. In this work, the techniques for Carleman linearization, phase-space embedding, and dynamic mode decomposition are integrated to rebuild an optimal dynamical system from time series for one specific variable. Using the Takens theorem, the embedding dimension is determined, which is adopted as the dynamical system's dimension. The Carleman linearization is then used to transform this finite nonlinear system into an infinite linear system, which is further truncated into a finite linear system using the dynamic mode decomposition technique. We illustrate the performance of this integrated technique using data generated by the well-known Lorenz model, the Duffing oscillator, and empirical records of electrocardiogram, electroencephalogram, and measles outbreaks. The results show that this solution accurately estimates the operators of the nonlinear dynamical systems. This work provides a new data-driven method to estimate the Carleman operator of nonlinear dynamical systems.

Complexity measures reveal age-dependent changes in electroencephalogram during working memory task

The alterations in electroencephalogram (EEG) signals are the complex outputs of functional factors, such as normal physiological aging, pathological process, which results in further cognitive decline. It is not clear that when brain aging initiates, but elderly people are vulnerable to be incipient of neurodegenerative diseases such as Alzheimer's disease. The EEG signals were recorded from 20 healthy middle age and 20 healthy elderly subjects while performing a working memory task. Higuchi's fractal dimension (HFD), Katz's fractal dimension (KFD), sample entropy and three Hjorth parameters were extracted to analyse the complexity of EEG signals. Four machine learning classifiers, multilayer perceptron (MLP), support vector machine (SVM), K-nearest neighbour (KNN), and logistic model tree (LMT) were employed to distinguish the EEG signals of middle age and elderly age groups. HFD, KFD and Hjorth complexity were found significantly correlated with age. MLP achieved the highest overall accuracy of 93.75%. For posterior region, the maximum accuracy of 92.50% was achieved using MLP. Since fractal dimension associated with the complexity of EEG signals, HFD, KFD and Hjorth complexity demonstrated the decreased complexity from middle age to elderly groups. The complexity features appear to be more appropriate indicators of monitoring EEG signal complexity in healthy aging.

Bibliometric analysis of electroencephalogram research in Parkinson's disease from 2004 to 2023

Background

Parkinson's disease (PD) is a prevalent neurodegenerative disorder affecting millions globally. It encompasses both motor and non-motor symptoms, with a notable impact on patients' quality of life. Electroencephalogram (EEG) is a non-invasive tool that is increasingly utilized to investigate neural mechanisms in PD, identify early diagnostic markers, and assess therapeutic responses.

Methods

The data were sourced from the Science Citation Index Expanded within the Web of Science Core Collection database, focusing on publications related to EEG research in PD from 2004 to 2023. A comprehensive bibliometric analysis was conducted using CiteSpace and VOSviewer software. The analysis began with an evaluation of the selected publications, identifying leading countries, institutions, authors, and journals, as well as co-cited references, to summarize the current state of EEG research in PD. Keywords are employed to identify research topics that are currently of interest in this field through the analysis of high-frequency keyword co-occurrence and cluster analysis. Finally, burst keywords were identified to uncover emerging trends and research frontiers in the field, highlighting shifts in interest and identifying future research directions.

Results

A total of 1,559 publications on EEG research in PD were identified. The United States, Germany, and England have made notable contributions to the field. The University of London is the leading institution in terms of publication output, with the University of California closely following. The most prolific authors are Brown P, Fuhr P, and Stam C In terms of total citations and per-article citations, Stam C has the highest number of citations, while Brown P has the highest H-index. In terms of the total number of publications, Clinical Neurophysiology is the leading journal, while Brain is the most highly cited. The most frequently cited articles pertain to software toolboxes for EEG analysis, neural oscillations, and PD pathophysiology. Through analyzing the keywords, four research hotspots were identified: research on the neural oscillations and connectivity, research on the innovations in EEG Analysis, impact of therapies on EEG, and research on cognitive and emotional assessments.

Conclusion

This bibliometric analysis demonstrates a growing global interest in EEG research in PD. The investigation of neural oscillations and connectivity remains a primary focus of research. The application of machine learning, deep learning, and task analysis techniques offers promising avenues for future research in EEG and PD, suggesting the potential for advancements in this field. This study offers valuable insights into the major research trends, influential contributors, and evolving themes in this field, providing a roadmap for future exploration.

Exploring the complexity of EEG patterns in Parkinson's disease

Parkinson's disease (PD) is a progressive neurodegenerative disorder primarily associated with motor dysfunctions. By the time of definitive diagnosis, about 60% of dopaminergic neurons have already been lost; moreover, even if dopaminergic drugs are highly effective in symptoms control, they only help maintaining a near-healthy condition when started as soon as possible. Therefore, interest in identifying early biomarkers of PD has grown in recent years, especially using neurophysiological techniques such as electroencephalography (EEG). This study aims to investigate brain complexity differences in PD patients compared to healthy controls, focusing on the beta band using approximate entropy (ApEn) analysis of resting-state EEG recordings. Sixty participants were recruited, including 25 PD patients and 35 healthy elderly subjects, matched for age and gender. EEG were recorded for each participant and ApEn values were computed in the beta 1 (13-20 Hz) and beta 2 (20-30 Hz) frequency bands for each EEG-channel and for ROIs. PD patients showed statistically lower ApEn values compared to controls in both beta 1 and beta 2 bands. Regarding electrodes analysis, beta 1 band alterations were found in frontocentral areas, while beta 2 band alterations were observed in centroparietal and frontocentral areas. Considering ROIs, statistically lower ApEn values for PD patients has been reported in central and parietal ROIs in the beta 2 band. Complexity reduction in these areas may underlie beta oscillatory activity dysfunction, reflecting impaired cortical mechanisms associated with motor dysfunction in PD. The results suggest that ApEn analysis of resting EEG activity may serve as a potential tool for early PD detection. Further studies are necessary to validate this approach in PD diagnosis and rehabilitation planning.

Beta-to-Theta Entropy Ratio of EEG in Aging, Frontotemporal Dementia, and Alzheimer's Dementia

BACKGROUND

Aging, frontotemporal dementia (FTD), and Alzheimer's dementia (AD) manifest electroencephalography (EEG) alterations, particularly in the beta-to-theta power ratio derived from linear power spectral density (PSD). Given the brain's nonlinear nature, the EEG nonlinear features could provide valuable physiological indicators of aging and cognitive impairment. Multiscale dispersion entropy (MDE) serves as a sensitive nonlinear metric for assessing the information content in EEGs across biologically relevant time scales.

OBJECTIVE

To compare the MDE-derived beta-to-theta entropy ratio with its PSD-based counterpart to detect differences between healthy young and elderly subjects and between different dementia subtypes.

METHODS

Scalp EEG recordings were obtained from two datasets: 1) Aging dataset: 133 healthy young and 65 healthy older adult individuals; and 2) Dementia dataset: 29 age-matched healthy controls (HC), 23 FTD, and 36 AD participants. The beta-to-theta ratios based on MDE vs. PSD were analyzed for both datasets. Finally, the relationships between cognitive performance and the beta-to-theta ratios were explored in HC, FTD, and AD.

RESULTS

In the Aging dataset, older adults had significantly higher beta-to-theta entropy ratios than young adults. In the Dementia dataset, this ratio outperformed the beta-to-theta PSD approach in distinguishing between HC, FTD, and AD. The AD participants had a significantly lower beta-to-theta entropy ratio than FTD, especially in the temporal region, unlike its corresponding PSD-based ratio. The beta-to-theta entropy ratio correlated significantly with cognitive performance.

CONCLUSION

Our study introduces the beta-to-theta entropy ratio using nonlinear MDE for EEG analysis, highlighting its potential as a sensitive biomarker for aging and cognitive impairment.

Neural connectivity patterns explain why adolescents perceive the world as moving slow

That younger individuals perceive the world as moving slower than adults is a familiar phenomenon. Yet, it remains an open question why that is. Using event segmentation theory, electroencephalogram (EEG) beamforming and nonlinear causal relationship estimation using artificial neural network methods, we studied neural activity while adolescent and adult participants segmented a movie. We show when participants were instructed to segment a movie into meaningful units, adolescents partitioned incoming information into fewer encapsulated segments or episodes of longer duration than adults. Importantly, directed communication between medial frontal and lower-level perceptual areas and between occipito-temporal regions in specific neural oscillation spectrums explained behavioral differences between groups. Overall, the study reveals that a different organization of directed communication between brain regions and inefficient transmission of information between brain regions are key to understand why younger people perceive the world as moving slow.

Characterization of antiseizure medications effects on the EEG neurodynamic by fractal dimension

Objectives

An important challenge in epilepsy is to define biomarkers of response to treatment. Many electroencephalography (EEG) methods and indices have been developed mainly using linear methods, e.g., spectral power and individual alpha frequency peak (IAF). However, brain activity is complex and non-linear, hence there is a need to explore EEG neurodynamics using nonlinear approaches. Here, we use the Fractal Dimension (FD), a measure of whole brain signal complexity, to measure the response to anti-seizure therapy in patients with Focal Epilepsy (FE) and compare it with linear methods.

Materials

Twenty-five drug-responder (DR) patients with focal epilepsy were studied before (t1, named DR-t1) and after (t2, named DR-t2) the introduction of the anti-seizure medications (ASMs). DR-t1 and DR-t2 EEG results were compared against 40 age-matched healthy controls (HC).

Methods

EEG data were investigated from two different angles: frequency domain-spectral properties in δ, θ, α, β, and γ bands and the IAF peak, and time-domain-FD as a signature of the nonlinear complexity of the EEG signals. Those features were compared among the three groups.

Results

The δ power differed between DR patients pre and post-ASM and HC (DR-t1 vs. HC, p < 0.01 and DR-t2 vs. HC, p < 0.01). The θ power differed between DR-t1 and DR-t2 (p = 0.015) and between DR-t1 and HC (p = 0.01). The α power, similar to the δ, differed between DR patients pre and post-ASM and HC (DR-t1 vs. HC, p < 0.01 and DR-t2 vs. HC, p < 0.01). The IAF value was lower for DR-t1 than DR-t2 (p = 0.048) and HC (p = 0.042). The FD value was lower in DR-t1 than in DR-t2 (p = 0.015) and HC (p = 0.011). Finally, Bayes Factor analysis showed that FD was 195 times more likely to separate DR-t1 from DR-t2 than IAF and 231 times than θ.

Discussion

FD measured in baseline EEG signals is a non-linear brain measure of complexity more sensitive than EEG power or IAF in detecting a response to ASMs. This likely reflects the non-oscillatory nature of neural activity, which FD better describes.

Conclusion

Our work suggests that FD is a promising measure to monitor the response to ASMs in FE.

Decoding imagined speech with delay differential analysis

Speech decoding from non-invasive EEG signals can achieve relatively high accuracy (70-80%) for strictly delimited classification tasks, but for more complex tasks non-invasive speech decoding typically yields a 20-50% classification accuracy. However, decoder generalization, or how well algorithms perform objectively across datasets, is complicated by the small size and heterogeneity of existing EEG datasets. Furthermore, the limited availability of open access code hampers a comparison between methods. This study explores the application of a novel non-linear method for signal processing, delay differential analysis (DDA), to speech decoding. We provide a systematic evaluation of its performance on two public imagined speech decoding datasets relative to all publicly available deep learning methods. The results support DDA as a compelling alternative or complementary approach to deep learning methods for speech decoding. DDA is a fast and efficient time-domain open-source method that fits data using only few strong features and does not require extensive preprocessing.

The value of linear and non-linear quantitative EEG analysis in paediatric epilepsy surgery: a machine learning approach

Epilepsy surgery is effective for patients with medication-resistant seizures, however 20-40% of them are not seizure free after surgery. Aim of this study is to evaluate the role of linear and non-linear EEG features to predict post-surgical outcome. We included 123 paediatric patients who underwent epilepsy surgery at Bambino Gesù Children Hospital (January 2009-April 2020). All patients had long term video-EEG monitoring. We analysed 1-min scalp interictal EEG (wakefulness and sleep) and extracted 13 linear and non-linear EEG features (power spectral density (PSD), Hjorth, approximate entropy, permutation entropy, Lyapunov and Hurst value). We used a logistic regression (LR) as feature selection process. To quantify the correlation between EEG features and surgical outcome we used an artificial neural network (ANN) model with 18 architectures. LR revealed a significant correlation between PSD of alpha band (sleep), Mobility index (sleep) and the Hurst value (sleep and awake) with outcome. The fifty-four ANN models gave a range of accuracy (46-65%) in predicting outcome. Within the fifty-four ANN models, we found a higher accuracy (64.8% ± 7.6%) in seizure outcome prediction, using features selected by LR. The combination of PSD of alpha band, mobility and the Hurst value positively correlate with good surgical outcome.

Closed-Loop Transcranial Electrical Neurostimulation for Sustained Attention Enhancement: A Pilot Study towards Personalized Intervention Strategies

Sustained attention is pivotal for tasks like studying and working for which focus and low distractions are necessary for peak productivity. This study explores the effectiveness of adaptive transcranial direct current stimulation (tDCS) in either the frontal or parietal region to enhance sustained attention. The research involved ten healthy university students performing the Continuous Performance Task-AX (AX-CPT) while receiving either frontal or parietal tDCS. The study comprised three phases. First, we acquired the electroencephalography (EEG) signal to identify the most suitable metrics related to attention states. Among different spectral and complexity metrics computed on 3 s epochs of EEG, the Fuzzy Entropy and Multiscale Sample Entropy Index of frontal channels were selected. Secondly, we assessed how tDCS at a fixed 1.0 mA current affects attentional performance. Finally, a real-time experiment involving continuous metric monitoring allowed personalized dynamic optimization of the current amplitude and stimulation site (frontal or parietal). The findings reveal statistically significant improvements in mean accuracy (94.04 vs. 90.82%) and reaction times (262.93 vs. 302.03 ms) with the adaptive tDCS compared to a non-stimulation condition. Average reaction times were statistically shorter during adaptive stimulation compared to a fixed current amplitude condition (262.93 vs. 283.56 ms), while mean accuracy stayed similar (94.04 vs. 93.36%, improvement not statistically significant). Despite the limited number of subjects, this work points out the promising potential of adaptive tDCS as a tailored treatment for enhancing sustained attention.

Respiratory modulation of the heart rate: A potential biomarker of cardiorespiratory function in human

In recent decade, wearable digital devices have shown potentials for the discovery of novel biomarkers of humans' physiology and behavior. Heart rate (HR) and respiration rate (RR) are most crucial bio-signals in humans' digital phenotyping research. HR is a continuous and non-invasive proxy to autonomic nervous system and ample evidence pinpoints the critical role of respiratory modulation of cardiac function. In the present study, we recorded longitudinal (7 days, 4.63 ± 1.52) HR and RR of 89 freely behaving human subjects (Female: 39, age 57.28 ± 5.67, Male: 50, age 58.48 ± 6.32) and analyzed their dynamics using linear models and information theoretic measures. While HR's linear and nonlinear characteristics were expressed within the plane of the HR-RR directed flow of information (HR→RR - RR→HR), their dynamics were determined by its RR→HR axis. More importantly, RR→HR quantified the effect of alcohol consumption on individuals' cardiorespiratory function independent of their consumed amount of alcohol, thereby signifying the presence of this habit in their daily life activities. The present findings provided evidence for the critical role of the respiratory modulation of HR, which was previously only studied in non-human animals. These results can contribute to humans' phenotyping research by presenting RR→HR as a digital diagnosis/prognosis marker of humans' cardiorespiratory pathology.

Classification of Sleep Quality and Aging as a Function of Brain Complexity: A Multiband Non-Linear EEG Analysis

Understanding and classifying brain states as a function of sleep quality and age has important implications for developing lifestyle-based interventions involving sleep hygiene. Current studies use an algorithm that captures non-linear features of brain complexity to differentiate awake electroencephalography (EEG) states, as a function of age and sleep quality. Fifty-eight participants were assessed using the Pittsburgh Sleep Quality Inventory (PSQI) and awake resting state EEG. Groups were formed based on age and sleep quality (younger adults n = 24, mean age = 24.7 years, SD = 3.43, good sleepers n = 11; older adults n = 34, mean age = 72.87; SD = 4.18, good sleepers n = 9). Ten non-linear features were extracted from multiband EEG analysis to feed several classifiers followed by a leave-one-out cross-validation. Brain state complexity accurately predicted (i) age in good sleepers, with 75% mean accuracy (across all channels) for lower frequencies (alpha, theta, and delta) and 95% accuracy at specific channels (temporal, parietal); and (ii) sleep quality in older groups with moderate accuracy (70 and 72%) across sub-bands with some regions showing greater differences. It also differentiated younger good sleepers from older poor sleepers with 85% mean accuracy across all sub-bands, and 92% at specific channels. Lower accuracy levels (<50%) were achieved in predicting sleep quality in younger adults. The algorithm discriminated older vs. younger groups excellently and could be used to explore intragroup differences in older adults to predict sleep intervention efficiency depending on their brain complexity.

Age-Related Characteristics of Resting-State Electroencephalographic Signals and the Corresponding Analytic Approaches: A Review

The study of the effects of aging on neural activity in the human brain has attracted considerable attention in neurophysiological, neuropsychiatric, and neurocognitive research, as it is directly linked to an understanding of the neural mechanisms underlying the disruption of the brain structures and functions that lead to age-related pathological disorders. Electroencephalographic (EEG) signals recorded during resting-state conditions have been widely used because of the significant advantage of non-invasive signal acquisition with higher temporal resolution. These advantages include the capability of a variety of linear and nonlinear signal analyses and state-of-the-art machine-learning and deep-learning techniques. Advances in artificial intelligence (AI) can not only reveal the neural mechanisms underlying aging but also enable the assessment of brain age reliably by means of the age-related characteristics of EEG signals. This paper reviews the literature on the age-related features, available analytic methods, large-scale resting-state EEG databases, interpretations of the resulting findings, and recent advances in age-related AI models.

Local contribution to the somatosensory evoked potentials in rat's thalamus

Local Field Potential (LFP), despite its name, often reflects remote activity. Depending on the orientation and synchrony of their sources, both oscillations and more complex waves may passively spread in brain tissue over long distances and be falsely interpreted as local activity at such distant recording sites. Here we show that the whisker-evoked potentials in the thalamic nuclei are of local origin up to around 6 ms post stimulus, but the later (7-15 ms) wave is overshadowed by a negative component reaching from cortex. This component can be analytically removed and local thalamic LFP can be recovered reliably using Current Source Density analysis. We used model-based kernel CSD (kCSD) method which allowed us to study the contribution of local and distant currents to LFP from rat thalamic nuclei and barrel cortex recorded with multiple, non-linear and non-regular multichannel probes. Importantly, we verified that concurrent recordings from the cortex are not essential for reliable thalamic CSD estimation. The proposed framework can be used to analyze LFP from other brain areas and has consequences for general LFP interpretation and analysis.

Local signal variability and functional connectivity: Sensitive measures of the excitation-inhibition ratio?

A novel network version of permutation entropy, the inverted joint permutation entropy (JPEinv), holds potential as non-invasive biomarker of abnormal excitation-inhibition (E-I) ratio in Alzheimer's disease (AD). In this computational modelling study, we test the hypotheses that this metric, and related measures of signal variability and functional connectivity, are sensitive to altered E-I ratios. The E-I ratio in each neural mass of a whole-brain computational network model was systematically varied. We evaluated whether JPEinv, local signal variability (by permutation entropy) and functional connectivity (by weighted symbolic mutual information (wsMI)) were related to E-I ratio, on whole-brain and regional level. The hub disruption index can identify regions primarily affected in terms of functional connectivity strength (or: degree) by the altered E-I ratios. Analyses were performed for a range of coupling strengths, filter and time-delay settings. On whole-brain level, higher E-I ratios were associated with higher functional connectivity (by JPEinv and wsMI) and lower local signal variability. These relationships were nonlinear and depended on the coupling strength, filter and time-delay settings. On regional level, hub-like regions showed a selective decrease in functional degree (by JPEinv and wsMI) upon a lower E-I ratio, and non-hub-like regions showed a selective increase in degree upon a higher E-I ratio. These results suggest that abnormal functional connectivity and signal variability, as previously reported in patients across the AD continuum, can inform us about altered E-I ratios.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-10003-x.

Investigating the interpretability of schizophrenia EEG mechanism through a 3DCNN-based hidden layer features aggregation framework

BACKGROUND AND OBJECTIVE

Electroencephalogram (EEG) signals record brain activity, with growing interest in quantifying neural activity through complexity analysis as a potential biological marker for schizophrenia. Presently, EEG complexity analysis primarily relies on manual feature extraction, which is subjective and yields varied findings in studies involving schizophrenia and healthy controls.

METHODS

This study aims to leverage deep learning methods for enhanced EEG complexity exploration, aiding early schizophrenia screening and diagnosis. Our proposed approach utilizes a three-dimensional Convolutional Neural Network (3DCNN) to extract enhanced data features for early schizophrenia identification and subsequent complexity analysis. Leveraging the spatiotemporal capabilities of 3DCNN, we extract advanced latent features and employ knowledge distillation to reintegrate these features into the original channels, creating feature-enhanced data.

RESULTS

We employ a 10-fold cross-validation strategy, achieving the average accuracies of 99.46% and 98.06% in subject-dependent experiments on Dataset 1(14SZ and 14HC) and Dataset 2 (45SZ and 39HC). The average accuracy for subject-independent is 96.04% and 92.67% on both datasets. Feature extraction and classification are conducted on both the re-aggregated data and the original data. Our results demonstrate that re-aggregated data exhibit superior classification performance and a more stable training process after feature extraction. In the complexity analysis of re-aggregated data, we observe lower entropy features in schizophrenic patients compared to healthy controls, with more pronounced differences in the temporal and frontal lobes. Analyzing Katz's Fractal Dimension (KFD) across three sub-bands of lobe channels reveals the lowest α band KFD value in schizophrenia patients.

CONCLUSIONS

This emphasizes the ability of our method to enhance the discrimination and interpretability in schizophrenia detection and analysis. Our approach enhances the potential for EEG-based schizophrenia diagnosis by leveraging deep learning, offering superior discrimination capabilities and richer interpretive insights.

Corticomuscular cross-recurrence analysis reveals between-limb differences in motor control among individuals with ACL reconstruction

Surgical reconstruction of the anterior cruciate ligament (ACL) and subsequent physical therapy can help athletes return to competition; however, re-injury rates remain disproportionately high due, in part, to lingering biomechanical and neurological factors that are not fully addressed during rehabilitation. Prior reports indicate that individuals exhibit altered electrical activity in both brain and muscle after ACL reconstruction (ACLR). In this investigation, we aimed to extend existing approaches by introducing a novel non-linear analysis of corticomuscular dynamics, which does not assume oscillatory coupling between brain and muscle: Corticomuscular cross-recurrence analysis (CM-cRQA). Our findings indicate that corticomuscular dynamics vary significantly between involved (injured) and uninvolved legs of participants with ACLR during voluntary isometric contractions between the brain and both the vastus medialis and lateralis. This finding points to a potential lingering neural deficit underlying re-injury for athletes after surgical reconstruction, namely the dynamical structure of neuromuscular (brain to quad muscle) coordination, which is significantly asymmetric, between limbs, in those who have ACLR.

Heterogeneous correlate and potential diagnostic biomarker of tinnitus based on nonlinear dynamics of resting-state EEG recordings

Tinnitus is a heterogeneous condition of hearing a rattling sound when there is no auditory stimulus. This rattling sound is associated with abnormal synchronous oscillations in auditory and non-auditory cortical areas. Since tinnitus is a highly heterogeneous condition with no objective detection criteria, it is necessary to search for indicators that can be compared between and within participants for diagnostic purposes. This study introduces heterogeneous though comparable indicators of tinnitus through investigation of spontaneous fluctuations in resting-state brain dynamics. The proposed approach uses nonlinear measures of chaos theory, to detect tinnitus and cross correlation patterns to reflect many of the previously reported neural correlates of tinnitus. These indicators may serve as effective measures of tinnitus risk even at early ages before any symptom is reported. The approach quantifies differences in oscillatory brain dynamics of tinnitus and normal subjects. It demonstrates that the left temporal areas of subjects with tinnitus exhibit larger lyapunov exponent indicating irregularity of brain dynamics in these regions. More complex dynamics is further recognized in tinnitus cases through entropy. We use this evidence to distinguish tinnitus patients from normal participants. Besides, we illustrate that certain anticorrelation patterns appear in these nonlinear measures across temporal and frontal areas in the brain perhaps corresponding to increased/decreased connectivity in certain brain networks and a shift in the balance of excitation and inhibition in tinnitus. Additionally, the main correlations are lost in tinnitus participants compared to control group suggesting involvement of distinct neural mechanisms in generation and persistence of tinnitus.

An Efficient Group Federated Learning Framework for Large-Scale EEG-Based Driver Drowsiness Detection

To avoid traffic accidents, monitoring the driver's electroencephalogram (EEG) signals to assess drowsiness is an effective solution. However, aggregating the personal data of these drivers may lead to insufficient data usage and pose a risk of privacy breaches. To address these issues, a framework called Group Federated Learning (Group-FL) for large-scale driver drowsiness detection is proposed, which can efficiently utilize diverse client data while protecting privacy. First, by arranging the clients into different levels of groups and gradually aggregating their model parameters from low-level groups to high-level groups, communication and time costs are reduced. In addition, to solve the problem of notable variations in EEG signals among different clients, a global-personalized deep neural network is designed. The global model extracts shared features from various clients, while the personalized model extracts fine-grained features from each client and outputs classification results. Finally, to address special issues such as scale/category imbalance and data pollution, three checking modules are designed for adjusting grouping, evaluating client data, and effectively applying personalized models. Through extensive experimentation, the effectiveness of each component within the framework was validated, and a mean accuracy, F1-score, and Area Under Curve (AUC) of 81.0%, 82.0%, and 87.9% was achieved, respectively, on a publicly available dataset comprising 11 subjects.

The Evolution of Neuroimaging Technologies to Evaluate Neural Activity Related to Knee Pain and Injury Risk

PURPOSE OF REVIEW

In this review, we present recent findings and advancements in the use of neuroimaging to evaluate neural activity relative to ACL injury risk and patellofemoral pain. In particular, we describe prior work using fMRI and EEG that demonstrate the value of these techniques as well as the necessity of continued development in this area. Our goal is to support future work by providing guidance for the successful application of neuroimaging techniques that most effectively expose pain and injury mechanisms.

RECENT FINDINGS

Recent studies that utilized both fMRI and EEG indicate that athletes who are at risk for future ACL injury exhibit divergent brain activity both during active lower extremity movement and at rest. Such activity patterns are likely due to alterations to cognitive, visual, and attentional processes that manifest as coordination deficits during naturalistic movement that may result in higher risk of injury. Similarly, in individuals with PFP altered brain activity in a number of key regions is related to subjective pain judgements as well as measures of fear of movement. Although these findings may begin to allow objective pain assessment and identification, continued refinement is needed. One key limitation across both ACL and PFP related work is the restriction of movement during fMRI and EEG data collection, which drastically limits ecological validity. Given the lack of sufficient research using EEG and fMRI within a naturalistic setting, our recommendation is that researchers target the use of mobile, source localized EEG as a primary methodology for exposing neural mechanisms of ACL injury risk and PFP. Our contention is that this method provides an optimal balance of spatial and temporal resolution with ecological validity via naturalistic movement.

Macroscopic resting-state brain dynamics are best described by linear models

It is typically assumed that large networks of neurons exhibit a large repertoire of nonlinear behaviours. Here we challenge this assumption by leveraging mathematical models derived from measurements of local field potentials via intracranial electroencephalography and of whole-brain blood-oxygen-level-dependent brain activity via functional magnetic resonance imaging. We used state-of-the-art linear and nonlinear families of models to describe spontaneous resting-state activity of 700 participants in the Human Connectome Project and 122 participants in the Restoring Active Memory project. We found that linear autoregressive models provide the best fit across both data types and three performance metrics: predictive power, computational complexity and the extent of the residual dynamics unexplained by the model. To explain this observation, we show that microscopic nonlinear dynamics can be counteracted or masked by four factors associated with macroscopic dynamics: averaging over space and over time, which are inherent to aggregated macroscopic brain activity, and observation noise and limited data samples, which stem from technological limitations. We therefore argue that easier-to-interpret linear models can faithfully describe macroscopic brain dynamics during resting-state conditions.

NLDyn - An open source MATLAB toolbox for the univariate and multivariate nonlinear dynamical analysis of physiological data

BACKGROUND AND OBJECTIVE

We present NLDyn, an open-source MATLAB toolbox tailored for in-depth analysis of nonlinear dynamics in biomedical signals. Our objective is to offer a user-friendly yet comprehensive platform for researchers to explore the intricacies of time series data.

METHODS

NLDyn integrates approximately 80 distinct methods, encompassing both univariate and multivariate nonlinear dynamics, setting it apart from existing solutions. This toolbox combines state-of-the-art nonlinear dynamical techniques with advanced multivariate entropy methods, providing users with powerful analytical capabilities. NLDyn enables analyses with or without a sliding window, and users can easily access and customize default parameters.

RESULTS

NLDyn generates results that are both exportable and visually informative, facilitating seamless integration into research and presentations. Its ongoing development ensures it remains at the forefront of nonlinear dynamics analysis.

CONCLUSIONS

NLDyn is a valuable resource for researchers in the biomedical field, offering an intuitive interface and a wide array of nonlinear analysis tools. Its integration of advanced techniques empowers users to gain deeper insights from their data. As we continually refine and expand NLDyn's capabilities, we envision it becoming an indispensable tool for the exploration of complex dynamics in biomedical signals.

Clinical Sensitivity of Fractal Neurodynamics

Among the significant advances in the understanding of the organization of the neuronal networks that coordinate the body and brain, their complex nature is increasingly important, resulting from the interaction between the very large number of constituents strongly organized hierarchically and at the same time with "self-emerging." This awareness drives us to identify the measures that best quantify the "complexity" that accompanies the continuous evolutionary dynamics of the brain. In this chapter, after an introductory section (Sect. 15.1), we examine how the Higuchi fractal dimension is able to perceive physiological processes (15.2), neurological (15.3) and psychiatric (15.4) disorders, and neuromodulation effects (15.5), giving a mention of other methods of measuring neuronal electrical activity in addition to electroencephalography, such as magnetoencephalography and functional magnetic resonance. Conscious that further progress will support a deeper understanding of the temporal course of neuronal activity because of continuous interaction with the environment, we conclude confident that the fractal dimension has begun to uncover important features of the physiology of brain activity and its alterations.

The Three-Lead EEG Sensor: Introducing an EEG-Assisted Depression Diagnosis System Based on Ant Lion Optimization

For depression diagnosis, traditional methods such as interviews and clinical scales have been widely leveraged in the past few decades, but they are subjective, time-consuming, and labor-consuming. With the development of affective computing and Artificial Intelligence (AI) technologies, Electroencephalogram (EEG)-based depression detection methods have emerged. However, previous research has virtually neglected practical application scenarios, as most studies have focused on analyzing and modeling EEG data. Furthermore, EEG data is typically obtained from specialized devices that are large, complex to operate, and poorly ubiquitous. To address these challenges, a wearable three-lead EEG sensor with flexible electrodes was developed to obtain prefrontal-lobe EEG data. Experimental measurements show that the EEG sensor achieves promising performance (background noise of no more than 0.91 μVpp, Signal-to-Noise Ratio (SNR) of 26--48 dB, and electrode-skin contact impedance of less than 1 K Ω). In addition, EEG data from 70 depressed patients and 108 healthy controls were collected using the EEG sensor, and the linear and nonlinear features were extracted. The features were then weighted and selected using the Ant Lion Optimization (ALO) algorithm to improve classification performance. The experimental results show that the k-NN classifier achieves a classification accuracy of 90.70%, specificity of 96.53%, and sensitivity of 81.79%, indicating the promising potential of the three-lead EEG sensor combined with the ALO algorithm and the k-NN classifier for EEG-assisted depression diagnosis.

ComsystanJ: A collection of Fiji/ImageJ2 plugins for nonlinear and complexity analysis in 1D, 2D and 3D

Complex systems such as the global climate, biological organisms, civilisation, technical or social networks exhibit diverse behaviours at various temporal and spatial scales, often characterized by nonlinearity, feedback loops, and emergence. These systems can be characterized by physical quantities such as entropy, information, chaoticity or fractality rather than classical quantities such as time, velocity, energy or temperature. The drawback of these complexity quantities is that their definitions are not always mathematically exact and computational algorithms provide estimates rather than exact values. Typically, evaluations can be cumbersome, necessitating specialized tools. We are therefore introducing ComsystanJ, a novel and user-friendly software suite, providing a comprehensive set of plugins for complex systems analysis, without the need for prior programming knowledge. It is platform independent, end-user friendly and extensible. ComsystanJ combines already known algorithms and newer methods for generalizable analysis of 1D signals, 2D images and 3D volume data including the generation of data sets such as signals and images for testing purposes. It is based on the framework of the open-source image processing software Fiji and ImageJ2. ComsystanJ plugins are macro recordable and are maintained as open-source software. ComsystanJ includes effective surrogate analysis in all dimensions to validate the features calculated by the different algorithms. Future enhancements of the project will include the implementation of parallel computing for image stacks and volumes and the integration of artificial intelligence methods to improve feature recognition and parameter calculation.

How far neuroscience is from understanding brains

The cellular biology of brains is relatively well-understood, but neuroscientists have not yet generated a theory explaining how brains work. Explanations of how neurons collectively operate to produce what brains can do are tentative and incomplete. Without prior assumptions about the brain mechanisms, I attempt here to identify major obstacles to progress in neuroscientific understanding of brains and central nervous systems. Most of the obstacles to our understanding are conceptual. Neuroscience lacks concepts and models rooted in experimental results explaining how neurons interact at all scales. The cerebral cortex is thought to control awake activities, which contrasts with recent experimental results. There is ambiguity distinguishing task-related brain activities from spontaneous activities and organized intrinsic activities. Brains are regarded as driven by external and internal stimuli in contrast to their considerable autonomy. Experimental results are explained by sensory inputs, behavior, and psychological concepts. Time and space are regarded as mutually independent variables for spiking, post-synaptic events, and other measured variables, in contrast to experimental results. Dynamical systems theory and models describing evolution of variables with time as the independent variable are insufficient to account for central nervous system activities. Spatial dynamics may be a practical solution. The general hypothesis that measurements of changes in fundamental brain variables, action potentials, transmitter releases, post-synaptic transmembrane currents, etc., propagating in central nervous systems reveal how they work, carries no additional assumptions. Combinations of current techniques could reveal many aspects of spatial dynamics of spiking, post-synaptic processing, and plasticity in insects and rodents to start with. But problems defining baseline and reference conditions hinder interpretations of the results. Furthermore, the facts that pooling and averaging of data destroy their underlying dynamics imply that single-trial designs and statistics are necessary.

Multivariate multiscale entropy (mMSE) as a tool for understanding the resting-state EEG signal dynamics: the spatial distribution and sex/gender-related differences

BACKGROUND

The study aimed to determine how the resting-state EEG (rsEEG) complexity changes both over time and space (channels). The complexity of rsEEG and its sex/gender differences were examined using the multivariate Multiscale Entropy (mMSE) in 95 healthy adults. Following the probability maps (Giacometti et al. in J Neurosci Methods 229:84-96, 2014), channel sets have been identified that correspond to the functional networks. For each channel set the area under curve (AUC), which represents the total complexity, MaxSlope-the maximum complexity change of the EEG signal at thefine scales (1:4 timescales), and AvgEnt-to the average entropy level at coarse-grained scales (9:12 timescales), respectively, were extracted. To check dynamic changes between the entropy level at the fine and coarse-grained scales, the difference in mMSE between the #9 and #4 timescale (DiffEnt) was also calculated.

RESULTS

We found the highest AUC for the channel sets corresponding to the somatomotor (SMN), dorsolateral network (DAN) and default mode (DMN) whereas the visual network (VN), limbic (LN), and frontoparietal (FPN) network showed the lowest AUC. The largest MaxSlope were in the SMN, DMN, ventral attention network (VAN), LN and FPN, and the smallest in the VN. The SMN and DAN were characterized by the highest and the LN, FPN, and VN by the lowest AvgEnt. The most stable entropy were for the DAN and VN while the LN showed the greatest drop of entropy at the coarse scales. Women, compared to men, showed higher MaxSlope and DiffEnt but lower AvgEnt in all channel sets.

CONCLUSIONS

Novel results of the present study are: (1) an identification of the mMSE features that capture entropy at the fine and coarse timescales in the channel sets corresponding to the main resting-state networks; (2) the sex/gender differences in these features.

Method for quantifying arousal and consciousness in healthy states and severe brain injury via EEG-based measures of corticothalamic physiology

BACKGROUND

Characterization of normal arousal states has been achieved by fitting predictions of corticothalamic neural field theory (NFT) to electroencephalographic (EEG) spectra to yield relevant physiological parameters.

NEW METHOD

A prior fitting method is extended to distinguish conscious and unconscious states in healthy and brain injured subjects by identifying additional parameters and clusters in parameter space.

RESULTS

Fits of NFT predictions to EEG spectra are used to estimate neurophysiological parameters in healthy and brain injured subjects. Spectra are used from healthy subjects in wake and sleep and from patients with unresponsive wakefulness syndrome, in a minimally conscious state (MCS), and emerged from MCS. Subjects cluster into three groups in parameter space: conscious healthy (wake and REM), sleep, and brain injured. These are distinguished by the difference X-Y between corticocortical (X) and corticothalamic (Y) feedbacks, and by mean neural response rates α and β to incoming spikes. X-Y tracks consciousness in healthy individuals, with smaller values in wake/REM than sleep, but cannot distinguish between brain injuries. Parameters α and β differentiate deep sleep from wake/REM and brain injury.

COMPARISON WITH EXISTING METHODS

Other methods typically rely on laborious clinical assessment, manual EEG scoring, or evaluation of measures like Φ from integrated information theory, for which no efficient method exists. In contrast, the present method can be automated on a personal computer.

CONCLUSION

The method provides a means to quantify consciousness and arousal in healthy and brain injured subjects, but does not distinguish subtypes of brain injury.

Early detection of Parkinson's disease: Systematic analysis of the influence of the eyes on quantitative biomarkers in resting state electroencephalography

While resting state electroencephalography (EEG) provides relevant information on pathological changes in Parkinson's disease, most studies focus on the eyes-closed EEG biomarkers. Recent evidence has shown that both eyes-open EEG and reactivity to eyes-opening can also differentiate Parkinson's disease from healthy aging, but no consensus has been reached on a discriminatory capability benchmark. The aim of this study was to determine the resting-state EEG biomarkers suitable for real-time application that can differentiate Parkinson's patients from healthy subjects under both eyes closed and open. For this, we analysed and compared the quantitative EEG analyses of 13 early-stage cognitively normal Parkinson's patients with an age and sex-matched healthy group. We found that Parkinson's disease exhibited abnormal excessive theta activity in eyes-closed, which was reflected by a significantly higher relative theta power, a higher time percentage with a frequency peak in the theta band and a reduced alpha/theta ratio, while Parkinson's patients showed a significantly steeper non-oscillatory spectral slope activity than that of healthy subjects. We also found considerably less alpha and beta reactivity to eyes-opening in Parkinson's disease plus a significant moderate correlation between these EEG-biomarkers and the MDS-UPDRS score, used to assesses the clinical symptoms of Parkinson's Disease. Both EEG recordings with the eyes open and reactivity to eyes-opening provided additional information to the eyes-closed condition. We thus strongly recommend that both eyes open and closed be used in clinical practice recording protocols to promote EEG as a complementary non-invasive screening method for the early detection of Parkinson's disease, which would allow clinicians to design patient-oriented treatment and improve the patient's quality of life.

Entropy and fractal analysis of brain-related neurophysiological signals in Alzheimer's and Parkinson's disease

Brain-related neuronal recordings, such as local field potential, electroencephalogram and magnetoencephalogram, offer the opportunity to study the complexity of the human brain at different spatial and temporal scales. The complex properties of neuronal signals are intrinsically related to the concept of 'scale-free' behavior and irregular dynamic, which cannot be fully described through standard linear methods, but can be measured by nonlinear indexes. A remarkable application of these analysis methods on electrophysiological recordings is the deep comprehension of the pathophysiology of neurodegenerative diseases, that has been shown to be associated to changes in brain activity complexity. In particular, a decrease of global complexity has been associated to Alzheimer's disease, while a local increase of brain signals complexity characterizes Parkinson's disease. Despite the recent proliferation of studies using fractal and entropy-based analysis, the application of these techniques is still far from clinical practice, due to the lack of an agreement about their correct estimation and a conclusive and shared interpretation. Along with the aim of helping towards the realization of a multidisciplinary audience to approach nonlinear methods based on the concepts of fractality and irregularity, this survey describes the implementation and proper employment of the mostly known and applied indexes in the context of Alzheimer's and Parkinson's diseases.

Spectral quantitative and semi-quantitative EEG provide complementary information on the life-long effects of early childhood malnutrition on cognitive decline

Objective

This study compares the complementary information from semi-quantitative EEG (sqEEG) and spectral quantitative EEG (spectral-qEEG) to detect the life-long effects of early childhood malnutrition on the brain.

Methods

Resting-state EEGs (N = 202) from the Barbados Nutrition Study (BNS) were used to examine the effects of protein-energy malnutrition (PEM) on childhood and middle adulthood outcomes. sqEEG analysis was performed on Grand Total EEG (GTE) protocol, and a single latent variable, the semi-quantitative Neurophysiological State (sqNPS) was extracted. A univariate linear mixed-effects (LME) model tested the dependence of sqNPS and nutritional group. sqEEG was compared with scores on the Montreal Cognitive Assessment (MoCA). Stable sparse classifiers (SSC) also measured the predictive power of sqEEG, spectral-qEEG, and a combination of both. Multivariate LME was applied to assess each EEG modality separately and combined under longitudinal settings.

Results

The univariate LME showed highly significant differences between previously malnourished and control groups (p < 0.001); age (p = 0.01) was also significant, with no interaction between group and age detected. Childhood sqNPS (p = 0.02) and adulthood sqNPS (p = 0.003) predicted MoCA scores in adulthood. The SSC demonstrated that spectral-qEEG combined with sqEEG had the highest predictive power (mean AUC 0.92 ± 0.005). Finally, multivariate LME showed that the combined spectral-qEEG+sqEEG models had the highest log-likelihood (-479.7).

Conclusion

This research has extended our prior work with spectral-qEEG and the long-term impact of early childhood malnutrition on the brain. Our findings showed that sqNPS was significantly linked to accelerated cognitive aging at 45-51 years of age. While sqNPS and spectral-qEEG produced comparable results, our study indicated that combining sqNPS and spectral-qEEG yielded better performance than either method alone, suggesting that a multimodal approach could be advantageous for future investigations.

Significance

Based on our findings, a semi-quantitative approach utilizing GTE could be a valuable diagnostic tool for detecting the lasting impacts of childhood malnutrition. Notably, sqEEG has not been previously explored or reported as a biomarker for assessing the longitudinal effects of malnutrition. Furthermore, our observations suggest that sqEEG offers unique features and information not captured by spectral quantitative EEG analysis and could lead to its improvement.

ECoG spiking activity and signal dimension are early predictive measures of epileptogenesis in a translational mouse model of traumatic brain injury

The latency between traumatic brain injury (TBI) and the onset of epilepsy (PTE) represents an opportunity for counteracting epileptogenesis. Antiepileptogenesis trials are hampered by the lack of sensitive biomarkers that allow to enrich patient's population at-risk for PTE. We aimed to assess whether specific ECoG signals predict PTE in a clinically relevant mouse model with ∼60% epilepsy incidence. TBI was provoked in adult CD1 male mice by controlled cortical impact on the left parieto-temporal cortex, then mice were implanted with two perilesional cortical screw electrodes and two similar electrodes in the hemisphere contralateral to the lesion site. Acute seizures and spikes/sharp waves were ECoG-recorded during 1 week post-TBI. These early ECoG events were analyzed according to PTE incidence as assessed by measuring spontaneous recurrent seizures (SRS) at 5 months post-TBI. We found that incidence, number and duration of acute seizures during 3 days post-TBI were similar in PTE mice and mice not developing epilepsy (No SRS mice). Control mice with cortical electrodes (naïve, n = 5) or with electrodes and craniotomy (sham, n = 5) exhibited acute seizures but did not develop epilepsy. The daily number of spikes/sharp waves at the perilesional electrodes was increased similarly in PTE (n = 15) and No SRS (n = 8) mice vs controls (p < 0.05, n = 10) from day 2 post-injury. Differently, the daily number of spikes/sharp waves at both contralateral electrodes showed a progressive increase in PTE mice vs No SRS and control mice. In particular, spikes number was higher in PTE vs No SRS mice (p < 0.05) at 6 and 7 days post-TBI, and this measure predicted epilepsy development with high accuracy (AUC = 0.77, p = 0.03; CI 0.5830-0.9670). The cut-off value was validated in an independent cohort of TBI mice (n = 12). The daily spike number at the contralateral electrodes showed a circadian distribution in PTE mice which was not observed in No SRS mice. Analysis of non-linear dynamics at each electrode site showed changes in dimensionality during 4 days post-TBI. This measure yielded the best discrimination between PTE and No SRS mice (p < 0.01) at the cortical electrodes contralateral to injury. Data show that epileptiform activity contralateral to the lesion site has the the highest predictive value for PTE in this model reinforcing the hypothesis that the hemisphere contralateral to the lesion core may drive epileptogenic networks after TBI.

EEG response of dexmedetomidine during drug induced sleep endoscopy

Introduction

Dexmedetomidine is one of the anesthetics of choice for drug induced sleep endoscopy (DISE), with advantages including limited respiratory depression, analgesia, and decreased incidence of emergence delirium. However, challenges with determining sedation levels and prolonged recovery have limited its usage. An improved understanding of the effect of dexmedetomidine on the level of sedation and the corresponding electroencephalographic (EEG) changes could help overcome these barriers.

Methods

Fifty-one patients received dexmedetomidine sedation with Richmond Agitation-Sedation Scale (RASS) score assessment and continuous EEG monitoring via SedLine for DISE. We constructed a pharmacokinetic model to determine continuous dexmedetomidine blood concentration. From the SedLine, we extracted the patient state index (PSI), and from the EEG we calculated the spectral edge frequency 95% (SEF95) and the correlation dimension (CD), a type of fractal dimension used to assess the complexity of a system. These metrics were subsequently compared against one another and with the dexmedetomidine concentration.

Results

Our pharmacokinetic model yielded a two-compartment model with volumes of 51.8 L and 106.2 L, with clearances of 69.5 and 168.9 L/h, respectively, and a time to effect of 9 min, similar to prior studies. Based on this model, decreasing RASS score, SEF95, CD, and PSI were all significantly associated with increasing dexmedetomidine concentration (p < 0.001, p = 0.006, p < 0.001 respectively). The CD, SEF95, and PSI better captured the effects of increasing dexmedetomidine concentration as compared to the RASS score. Simulating dexmedetomidine concentration based on titration to target levels derived from CD and PSI confirmed commonly used dexmedetomidine infusion dosages.

Conclusion

Dexmedetomidine use for DISE confirmed previous pharmacokinetic models seen with dexmedetomidine. Complex EEG metrics such as PSI and CD, as compared to RASS score and SEF95, better captured changes in brain state from dexmedetomidine and have potential to improve the monitoring of dexmedetomidine sedation.

Effects of Sleep Reactivity on Sleep Macro-Structure, Orderliness, and Cortisol After Stress: A Preliminary Study in Healthy Young Adults

Purpose

To investigate changes and links of stress and high sleep reactivity (H-SR) on the macro-structure and orderliness of sleep and cortisol levels in good sleepers (GS).

Patients and Methods

Sixty-two GS (18-40 years old) were recruited, with 32 in the stress group and 30 in the control group. Each group was further divided into H-SR and low SR subgroups based on the Ford Insomnia Response to Stress Test. All participants completed two nights of polysomnography in a sleep laboratory. Before conducting polysomnography on the second night, the stress group completed the Trier Social Stress Test and saliva was collected.

Results

The duration of NREM sleep stages 1, 2 (N1, N2) and rapid eye movement sleep (REM) decreased, and the values of approximate entropy, sample entropy, fuzzy entropy, and multiscale entropy increased under stress and SR effects. Stress increased rapid eye movement density, and H-SR increased cortisol reactivity.

Conclusion

Stress can damage the sleep and increase cortisol release in GS, especially those with H-SR. N1, N2 and REM sleep are more easily affected, while NREM sleep stage 3 sleep is relatively stable.

Two-channel EEG based diagnosis of panic disorder and major depressive disorder using machine learning and non-linear dynamical methods

The current study aimed to investigate the possibility of rapid and accurate diagnoses of Panic disorder (PD) and Major depressive disorder (MDD) using machine learning. The support vector machine method was applied to 2-channel EEG signals from the frontal lobes (Fp1 and Fp2) of 149 participants to classify PD and MDD patients from healthy individuals using non-linear measures as features. We found significantly lower correlation dimension and Lempel-Ziv complexity in PD patients and MDD patients in the left hemisphere compared to healthy subjects at rest. Most importantly, we obtained a 90% accuracy in classifying MDD patients vs. healthy individuals, a 68% accuracy in classifying PD patients vs. controls, and a 59% classification accuracy between PD and MDD patients. In addition to demonstrating classification performance in a simplified setting, the observed differences in EEG complexity between subject groups suggest altered cortical processing present in the frontal lobes of PD patients that can be captured through non-linear measures. Overall, this study suggests that machine learning and non-linear measures using only 2-channel frontal EEGs are useful for aiding the rapid diagnosis of panic disorder and major depressive disorder.

Inter-individual variability during neurodevelopment: an investigation of linear and nonlinear resting-state EEG features in an age-homogenous group of infants

Electroencephalography measures are of interest in developmental neuroscience as potentially reliable clinical markers of brain function. Features extracted from electroencephalography are most often averaged across individuals in a population with a particular condition and compared statistically to the mean of a typically developing group, or a group with a different condition, to define whether a feature is representative of the populations as a whole. However, there can be large variability within a population, and electroencephalography features often change dramatically with age, making comparisons difficult. Combined with often low numbers of trials and low signal-to-noise ratios in pediatric populations, establishing biomarkers can be difficult in practice. One approach is to identify electroencephalography features that are less variable between individuals and are relatively stable in a healthy population during development. To identify such features in resting-state electroencephalography, which can be readily measured in many populations, we introduce an innovative application of statistical measures of variance for the analysis of resting-state electroencephalography data. Using these statistical measures, we quantified electroencephalography features commonly used to measure brain development-including power, connectivity, phase-amplitude coupling, entropy, and fractal dimension-according to their intersubject variability. Results from 51 6-month-old infants revealed that the complexity measures, including fractal dimension and entropy, followed by connectivity were the least variable features across participants. This stability was found to be greatest in the right parietotemporal region for both complexity feature, but no significant region of interest was found for connectivity feature. This study deepens our understanding of physiological patterns of electroencephalography data in developing brains, provides an example of how statistical measures can be used to analyze variability in resting-state electroencephalography in a homogeneous group of healthy infants, contributes to the establishment of robust electroencephalography biomarkers of neurodevelopment through the application of variance analyses, and reveals that nonlinear measures may be most relevant biomarkers of neurodevelopment.

Kernel-based Nonlinear Manifold Learning for EEG-based Functional Connectivity Analysis and Channel Selection with Application to Alzheimer's Disease

Dynamical, causal, and cross-frequency coupling analysis using the electroencephalogram (EEG) has gained significant attention for diagnosing and characterizing neurological disorders. Selecting important EEG channels is crucial for reducing computational complexity in implementing these methods and improving classification accuracy. In neuroscience, measures of (dis)similarity between EEG channels are often used as functional connectivity (FC) features, and important channels are selected via feature selection. Developing a generic measure of (dis)similarity is important for FC analysis and channel selection. In this study, learning of (dis)similarity information within the EEG is achieved using kernel-based nonlinear manifold learning. The focus is on FC changes and, thereby, EEG channel selection. Isomap and Gaussian Process Latent Variable Model (Isomap-GPLVM) are employed for this purpose. The resulting kernel (dis)similarity matrix is used as a novel measure of linear and nonlinear FC between EEG channels. The analysis of EEG from healthy controls (HC) and patients with mild to moderate Alzheimer's disease (AD) are presented as a case study. Classification results are compared with other commonly used FC measures. Our analysis shows significant differences in FC between bipolar channels of the occipital region and other regions (i.e. parietal, centro-parietal, and fronto-central) between AD and HC groups. Furthermore, our results indicate that FC changes between channels along the fronto-parietal region and the rest of the EEG are important in diagnosing AD. Our results and its relation to functional networks are consistent with those obtained from previous studies using fMRI, resting-state fMRI and EEG.

Sensory-Processing Sensitivity Is Associated with Increased Neural Entropy

BACKGROUND

This study aimed at answering the following research questions: (1) Does the self-reported level of sensory-processing sensitivity (SPS) correlate with complexity, or criticality features of the electroencephalogram (EEG)? (2) Are there significant EEG differences comparing individuals with high and low levels of SPS?

METHODS

One hundred fifteen participants were measured with 64-channel EEG during a task-free resting state. The data were analyzed using criticality theory tools (detrended fluctuation analysis, neuronal avalanche analysis) and complexity measures (sample entropy, Higuchi's fractal dimension). Correlations with the 'Highly Sensitive Person Scale' (HSPS-G) scores were determined. Then, the cohort's lowest and the highest 30% were contrasted as opposites. EEG features were compared between the two groups by applying a Wilcoxon signed-rank test.

RESULTS

During resting with eyes open, HSPS-G scores correlated significantly positively with the sample entropy and Higuchi's fractal dimension (Spearman's ρ = 0.22, p < 0.05). The highly sensitive group revealed higher sample entropy values (1.83 ± 0.10 vs. 1.77 ± 0.13, p = 0.031). The increased sample entropy in the highly sensitive group was most pronounced in the central, temporal, and parietal regions.

CONCLUSION

For the first time, neurophysiological complexity features associated with SPS during a task-free resting state were demonstrated. Evidence is provided that neural processes differ between low- and highly-sensitive persons, whereby the latter displayed increased neural entropy. The findings support the central theoretical assumption of enhanced information processing and could be important for developing biomarkers for clinical diagnostics.

Assessing stroke rehabilitation degree based on quantitative EEG index and nonlinear parameters

The assessment of motor function is critical to the rehabilitation of stroke patients. However, commonly used evaluation methods are based on behavior scoring, which lacks neurological indicators that directly reflect the motor function of the brain. The objective of this study was to investigate whether resting-state EEG indicators could improve stroke rehabilitation evaluation. We recruited 68 participants and recorded their resting-state EEG data. According to Brunnstrom stage, the participants were divided into three groups: severe, moderate, and mild. Ten quantitative electroencephalographic (QEEG) and five non-linear parameters of resting-state EEG were calculated for further analysis. Statistical tests were performed, and the genetic algorithm-support vector machine was used to select the best feature combination for classification. We found the QEEG parameters show significant differences in Delta, Alpha1, Alpha2, DAR, and DTABR (P < 0.05) among the three groups. Regarding nonlinear parameters, ApEn, SampEn, Lz, and C0 showed significant differences (P < 0.05). The optimal feature classification combination accuracy rate reached 85.3%. Our research shows that resting-state EEG indicators could be used for stroke rehabilitation evaluation.

Complex Brain-Heart Mapping in Mental and Physical Stress

OBJECTIVE

The central and autonomic nervous systems are deemed complex dynamic systems, wherein each system as a whole shows features that the individual system sub-components do not. They also continuously interact to maintain body homeostasis and appropriate react to endogenous and exogenous stimuli. Such interactions are comprehensively referred to functional brain-heart interplay (BHI). Nevertheless, it remains uncertain whether this interaction also exhibits complex characteristics, that is, whether the dynamics of the entire nervous system inherently demonstrate complex behavior, or if such complexity is solely a trait of the central and autonomic systems. Here, we performed complexity mapping of the BHI dynamics under mental and physical stress conditions.

METHODS AND PROCEDURES

Electroencephalographic and heart rate variability series were obtained from 56 healthy individuals performing mental arithmetic or cold-pressure tasks, and physiological series were properly combined to derive directional BHI series, whose complexity was quantified through fuzzy entropy.

RESULTS

The experimental results showed that BHI complexity is mainly modulated in the efferent functional direction from the brain to the heart, and mainly targets vagal oscillations during mental stress and sympathovagal oscillations during physical stress.

CONCLUSION

We conclude that the complexity of BHI mapping may provide insightful information on the dynamics of both central and autonomic activity, as well as on their continuous interaction.

CLINICAL IMPACT

This research enhances our comprehension of the reciprocal interactions between central and autonomic systems, potentially paving the way for more accurate diagnoses and targeted treatments of cardiovascular, neurological, and psychiatric disorders.

Evidence of Chaos in Electroencephalogram Signatures of Human Performance: A Systematic Review

(1) Background: Chaos, a feature of nonlinear dynamical systems, is well suited for exploring biological time series, such as heart rates, respiratory records, and particularly electroencephalograms. The primary purpose of this article is to review recent studies using chaos theory and nonlinear dynamical methods to analyze human performance in different brain processes. (2) Methods: Several studies have examined chaos theory and related analytical tools for describing brain dynamics. The present study provides an in-depth analysis of the computational methods that have been proposed to uncover brain dynamics. (3) Results: The evidence from 55 articles suggests that cognitive function is more frequently assessed than other brain functions in studies using chaos theory. The most frequently used techniques for analyzing chaos include the correlation dimension and fractal analysis. Approximate, Kolmogorov and sample entropy account for the largest proportion of entropy algorithms in the reviewed studies. (4) Conclusions: This review provides insights into the notion of the brain as a chaotic system and the successful use of nonlinear methods in neuroscience studies. Additional studies of brain dynamics would aid in improving our understanding of human cognitive performance.

Nonlinear dynamics of EEG responses to unmanned vehicle visual detection with different levels of task difficulty

The main objective of this study was to examine the presence of chaos in the EEG recordings of brain activity under simulated unmanned ground vehicle visual detection scenarios with different levels of task difficulty. One hundred and fifty people participated in the experiment and completed four visual detection task scenarios: (1) change detection, (2) a threat detection task, (3) a dual-task with different change detection task rates, and (4) a dual-task with different threat detection task rates. We used the largest Lyapunov exponent and correlation dimension of the EEG data and performed 0-1 tests on the EEG data. The results revealed a change in the level of nonlinearity in the EEG data corresponding to different levels of cognitive task difficulty. The differences in EEG nonlinearity measures among the studied levels of task difficulty, as well as between a single task scenario and a dual-task scenario, have also been assessed. The results increase our understanding of the nature of unmanned systems' operational requirements.

The time scales of irreversibility in spontaneous brain activity are altered in obsessive compulsive disorder

We study how obsessive-compulsive disorder (OCD) affects the complexity and time-reversal symmetry-breaking (irreversibility) of the brain resting-state activity as measured by magnetoencephalography (MEG). Comparing MEG recordings from OCD patients and age/sex matched control subjects, we find that irreversibility is more concentrated at faster time scales and more uniformly distributed across different channels of the same hemisphere in OCD patients than in control subjects. Furthermore, the interhemispheric asymmetry between homologous areas of OCD patients and controls is also markedly different. Some of these differences were reduced by 1-year of Kundalini Yoga meditation treatment. Taken together, these results suggest that OCD alters the dynamic attractor of the brain's resting state and hint at a possible novel neurophysiological characterization of this psychiatric disorder and how this therapy can possibly modulate brain function.

Stability, change, and reliable individual differences in electroencephalography measures: a lifespan perspective on progress and opportunities

Electroencephalographic (EEG) methods have great potential to serve both basic and clinical science approaches to understand individual differences in human neural function. Importantly, the psychometric properties of EEG data, such as internal consistency and test-retest reliability, constrain their ability to differentiate individuals successfully. Rapid and recent technological and computational advancements in EEG research make it timely to revisit the topic of psychometric reliability in the context of individual difference analyses. Moreover, pediatric and clinical samples provide some of the most salient and urgent opportunities to apply individual difference approaches, but the changes these populations experience over time also provide unique challenges from a psychometric perspective. Here we take a developmental neuroscience perspective to consider progress and new opportunities for parsing the reliability and stability of individual differences in EEG measurements across the lifespan. We first conceptually map the different profiles of measurement reliability expected for different types of individual difference analyses over the lifespan. Next, we summarize and evaluate the state of the field's empirical knowledge and need for testing measurement reliability, both internal consistency and test-retest reliability, across EEG measures of power, event-related potentials, nonlinearity, and functional connectivity across ages. Finally, we highlight how standardized pre-processing software for EEG denoising and empirical metrics of individual data quality may be used to further improve EEG-based individual differences research moving forward. We also include recommendations and resources throughout that individual researchers can implement to improve the utility and reproducibility of individual differences analyses with EEG across the lifespan.