Poincaré analysis of the electroencephalogram during sevoflurane anesthesia

Chaotic nature of the electroencephalogram during shallow and deep anesthesia: From analysis of the Lyapunov exponent

In conscious states, the electrodynamics of the cortex are reported to work near a critical point or phase transition of chaotic dynamics, known as the edge-of-chaos, representing a boundary between stability and chaos. Transitions away from this boundary disrupt cortical information processing and induce a loss of consciousness. The entropy of the electroencephalogram (EEG) is known to decrease as the level of anesthesia deepens. However, whether the chaotic dynamics of electroencephalographic activity shift from this boundary to the side of stability or the side of chaotic enhancement during anesthesia-induced loss of consciousness remains poorly understood. We investigated the chaotic properties of EEGs at two different depths of clinical anesthesia using the maximum Lyapunov exponent, which is mathematically regarded as a formal measure of chaotic nature, using the Rosenstein algorithm. In 14 adult patients, 12 s of electroencephalographic signals were selected during two depths of clinical anesthesia (sevoflurane concentration 2% as relatively deep anesthesia, sevoflurane concentration 0.6% as relatively shallow anesthesia). Lyapunov exponents, correlation dimensions and approximate entropy were calculated from these electroencephalographic signals. As a result, maximum Lyapunov exponent was generally positive during sevoflurane anesthesia, and both maximum Lyapunov exponents and correlation dimensions were significantly greater during deep anesthesia than during shallow anesthesia despite reductions in approximate entropy. The chaotic nature of the EEG might be increased at clinically deeper inhalational anesthesia, despite the decrease in randomness as reflected in the decreased entropy, suggesting a shift to the side of chaotic enhancement under anesthesia.

Assessing consciousness in patients with disorders of consciousness using soft-clustering

Consciousness is something we experience in our everyday life, more especially between the time we wake up in the morning and go to sleep at night, but also during the rapid eye movement (REM) sleep stage. Disorders of consciousness (DoC) are states in which a person's consciousness is damaged, possibly after a traumatic brain injury. Completely locked-in syndrome (CLIS) patients, on the other hand, display covert states of consciousness. Although they appear unconscious, their cognitive functions are mostly intact. Only, they cannot externally display it due to their quadriplegia and inability to speak. Determining these patients' states constitutes a challenging task. The ultimate goal of the approach presented in this paper is to assess these CLIS patients consciousness states. EEG data from DoC patients are used here first, under the assumption that if the proposed approach is able to accurately assess their consciousness states, it will assuredly do so on CLIS patients too. This method combines different sets of features consisting of spectral, complexity and connectivity measures in order to increase the probability of correctly estimating their consciousness levels. The obtained results showed that the proposed approach was able to correctly estimate several DoC patients' consciousness levels. This estimation is intended as a step prior attempting to communicate with them, in order to maximise the efficiency of brain-computer interfaces (BCI)-based communication systems.

Application of Soft-Clustering to Assess Consciousness in a CLIS Patient

Completely locked-in (CLIS) patients are characterized by sufficiently intact cognitive functions, but a complete paralysis that prevents them to interact with their surroundings. On one hand, studies have shown that the ability to communicate plays an important part in these patients' quality of life and prognosis. On the other hand, brain-computer interfaces (BCIs) provide a means for them to communicate using their brain signals. However, one major problem for such patients is the difficulty to determine if they are conscious or not at a specific time. This work aims to combine different sets of features consisting of spectral, complexity and connectivity measures, to increase the probability of correctly estimating CLIS patients' consciousness levels. The proposed approach was tested on data from one CLIS patient, which is particular in the sense that the experimenter was able to point out one time frame Δt during which he was undoubtedly conscious. Results showed that the method presented in this paper was able to detect increases and decreases of the patient's consciousness levels. More specifically, increases were observed during this Δt, corroborating the assertion of the experimenter reporting that the patient was definitely conscious then. Assessing the patients' consciousness is intended as a step prior attempting to communicate with them, in order to maximize the efficiency of BCI-based communication systems.

EEG artifact removal using sub-space decomposition, nonlinear dynamics, stationary wavelet transform and machine learning algorithms

Blind source separation (BSS) methods have received a great deal of attention in electroencephalogram (EEG) artifact elimination as they are routine and standard signal processing tools to remove artifacts and reserve desired neural information. On the other hand, a classifier should follow BSS methods to automatically identify artifactual sources and remove them in the following steps. In addition, removing all detected artifactual components leads to loss of information since some desired information related to neural activity leaks to these sources. So, an approach should be employed to detect and suppress the artifacts and reserve neural activity. This study introduces a novel method based on EEG and Poincare planes in the phase space to detect artifactual components estimated by second-order blind identification (SOBI). Artifacts are detected using a mixture of well-known conventional classifiers and were removed employing stationary wavelet transform (SWT) to reserve neural information. The proposed method is a combination of signal processing techniques and machine learning algorithms, including multi-layer perceptron (MLP), K-nearest neighbor (KNN), naïve Bayes, and support vector machine (SVM) which have significant results while applying our proposed method to different scenarios. Simulated, semi-simulated, and real EEG signals are employed to evaluate the proposed method, and several evaluation criteria are calculated. We achieved acceptable results, for example, 98% average accuracy and 97% average sensitivity in artifactual EEG component detection or about 2% as mean square error in EEG reconstruction after artifact removal. Results showed that the proposed method is effective and can be used in future studies as we have considered different real-world scenarios to evaluate it.

Network Analysis of Time Series: Novel Approaches to Network Neuroscience

In the last two decades, there has been an explosion of interest in modeling the brain as a network, where nodes correspond variously to brain regions or neurons, and edges correspond to structural or statistical dependencies between them. This kind of network construction, which preserves spatial, or structural, information while collapsing across time, has become broadly known as "network neuroscience." In this work, we provide an alternative application of network science to neural data: network-based analysis of non-linear time series and review applications of these methods to neural data. Instead of preserving spatial information and collapsing across time, network analysis of time series does the reverse: it collapses spatial information, instead preserving temporally extended dynamics, typically corresponding to evolution through some kind of phase/state-space. This allows researchers to infer a, possibly low-dimensional, "intrinsic manifold" from empirical brain data. We will discuss three methods of constructing networks from nonlinear time series, and how to interpret them in the context of neural data: recurrence networks, visibility networks, and ordinal partition networks. By capturing typically continuous, non-linear dynamics in the form of discrete networks, we show how techniques from network science, non-linear dynamics, and information theory can extract meaningful information distinct from what is normally accessible in standard network neuroscience approaches.

Assessing rheoencephalography dynamics through analysis of the interactions among brain and cardiac networks during general anesthesia

Cerebral blood flow (CBF) reflects the rate of delivery of arterial blood to the brain. Since no nutrients, oxygen or water can be stored in the cranial cavity due to space and pressure restrictions, a continuous perfusion of the brain is critical for survival. Anesthetic procedures are known to affect cerebral hemodynamics, but CBF is only monitored in critical patients due, among others, to the lack of a continuous and affordable bedside monitor for this purpose. A potential solution through bioelectrical impedance technology, also known as rheoencephalography (REG), is proposed, that could fill the existing gap for a low-cost and effective CBF monitoring tool. The underlying hypothesis is that REG signals carry information on CBF that might be recovered by means of the application of advanced signal processing techniques, allowing to track CBF alterations during anesthetic procedures. The analysis of REG signals was based on geometric features extracted from the time domain in the first place, since this is the standard processing strategy for this type of physiological data. Geometric features were tested to distinguish between different anesthetic depths, and they proved to be capable of tracking cerebral hemodynamic changes during anesthesia. Furthermore, an approach based on Poincaré plot features was proposed, where the reconstructed attractors form REG signals showed significant differences between different anesthetic states. This was a key finding, providing an alternative to standard processing of REG signals and supporting the hypothesis that REG signals do carry CBF information. Furthermore, the analysis of cerebral hemodynamics during anesthetic procedures was performed by means of studying causal relationships between global hemodynamics, cerebral hemodynamics and electroencephalogram (EEG) based-parameters. Interactions were detected during anesthetic drug infusion and patient positioning (Trendelenburg positioning and passive leg raise), providing evidence of the causal coupling between hemodynamics and brain activity. The provided alternative of REG signal processing confirmed the hypothesis that REG signals carry information on CBF. The simplicity of the technology, together with its low cost and easily interpretable outcomes, should provide a new opportunity for REG to reach standard clinical practice. Moreover, causal relationships among the hemodynamic physiological signals and brain activity were assessed, suggesting that the inclusion of REG information in depth of anesthesia monitors could be of valuable use to prevent unwanted CBF alterations during anesthetic procedures.

Electroencephalogram variability analysis for monitoring depth of anesthesia

Objective. In this paper, a new approach of extracting and measuring the variability in electroencephalogram (EEG) was proposed to assess the depth of anesthesia (DOA) under general anesthesia.Approach. The EEG variability (EEGV) was extracted as a fluctuation in time interval that occurs between two local maxima of EEG. Eight parameters related to EEGV were measured in time and frequency domains, and compared with state-of-the-art DOA estimation parameters, including sample entropy, permutation entropy, median frequency and spectral edge frequency of EEG. The area under the receiver-operator characteristics curve (AUC) and Pearson correlation coefficient were used to validate its performance on 56 patients.Main results. Our proposed EEGV-derived parameters yield significant difference for discriminating between awake and anesthesia stages at a significance level of 0.05, as well as improvement in AUC and correlation coefficient on average, which surpasses the conventional features of EEG in detection accuracy of unconscious state and tracking the level of consciousness.Significance. To sum up, EEGV analysis provides a new perspective in quantifying EEG and corresponding parameters are powerful and promising for monitoring DOA under clinical situations.

Poincaré Plot Area of Gamma-Band EEG as a Measure of Emergence From Inhalational General Anesthesia

The Poincaré plot obtained from electroencephalography (EEG) has been used to evaluate the depth of anesthesia. A standalone EEG Analyzer application was developed; raw EEG signals obtained from a bispectral index (BIS) monitor were analyzed using an on-line monitoring system. Correlations between Poincaré plot parameters and other measurements associated with anesthesia depth were evaluated during emergence from inhalational general anesthesia. Of the participants, 20 were adults anesthetized with sevoflurane (adult_{_SEV}), 20 were adults anesthetized with desflurane (adult_{_DES}), and 20 were pediatric patients anesthetized with sevoflurane (ped_{_SEV}). EEG signals were preprocessed through six bandpass digital filters (f0: 0.5–47 Hz, f1: 0.5–8 Hz, f2: 8–13 Hz, f3: 13–20 Hz, f4: 20–30 Hz, and f5: 30–47 Hz). The Poincaré plot-area ratio (PP_AR = PP_{A_fx}/PP_{A_f0}, fx = f1∼f5) was analyzed at five frequency ranges. Regardless of the inhalational anesthetic used, there were strong linear correlations between the logarithm of PP_AR at f5 and BIS (R² = 0.67, 0.79, and 0.71, in the adult_{_SEV}, adult_{_DES}, and ped_{_SEV} groups, respectively). As an additional observation, a part of EMG activity at the gamma range of 30–47 Hz probably influenced the calculations of BIS and PP_{AR_f5} with a non-negligible level. The logarithm of PP_AR in the gamma band was most sensitive to state changes during the emergence process and could provide a new non-proprietary parameter that correlates with changes in BIS during measurement of anesthesia depth.

Consciousness Detection in a Complete Locked-in Syndrome Patient through Multiscale Approach Analysis

Completely locked-in state (CLIS) patients are unable to speak and have lost all muscle movement. From the external view, the internal brain activity of such patients cannot be easily perceived, but CLIS patients are considered to still be conscious and cognitively active. Detecting the current state of consciousness of CLIS patients is non-trivial, and it is difficult to ascertain whether CLIS patients are conscious or not. Thus, it is important to find alternative ways to re-establish communication with these patients during periods of awareness, and one such alternative is through a brain-computer interface (BCI). In this study, multiscale-based methods (multiscale sample entropy, multiscale permutation entropy and multiscale Poincaré plots) were applied to analyze electrocorticogram signals from a CLIS patient to detect the underlying consciousness level. Results from these different methods converge to a specific period of awareness of the CLIS patient in question, coinciding with the period during which the CLIS patient is recorded to have communicated with an experimenter. The aim of the investigation is to propose a methodology that could be used to create reliable communication with CLIS patients.

Potential Motor Benefits of Visual Feedback of Error Reduction for Older Adults

This study investigated how visual feedback of virtual error reduction (ER) modified the visuomotor performance of older adults with limited attentional capacity. Error structures of young and older adults during birhythmic force tracking were contrasted when the visualized error size was exact or half of the actual size. As compared with full-size error feedback, ER feedback improved the force tracking symmetry of older adults, but undermined that of young adults. Extended Poincaré analysis revealed that young adults presented greater short-term error variability (mean value of κ-lagged SD1 of the error signal) with ER feedback, which led to a smaller mean value of κ-lagged SD1 of the error signal for older adults. The ER-related task improvement of the older adults was negatively correlated with the size of the tracking errors with real error feedback and positively correlated with ER-related increases in force spectral symmetry and decreases in the mean value of κ-lagged SD1 of the error signal. ER feedback could advance visuomotor tasks for older adults who perform worse with full-size visual feedback by the enhancement of self-efficacy and stabilization of negative internal feedback.

Hand grip strength variability during serial testing as an entropic biomarker of aging: a Poincaré plot analysis

Background

The Poincaré plot method can be used for both qualitative and quantitative assessment of self-similarity in usually periodic functions, hence the idea of applying it to the study of homeostasis of living organisms. From the analysis of numerous scientific data, it can be concluded that hand functionality can be correlated with the state of the human body as a biological system exposed to various forms of ontogenetic stress.

Methods

We used the Poincaré plot method to analyze the variability of hand grip strength (HGS), as an entropic biomarker of aging, during 60 repetitive tests of the dominant and nondominant hand, in young and older healthy subjects. An observational cross-sectional study was performed on 80 young adults (18–22 years old, mean age 20.01 years) and 80 older people (65–69 years old, mean age 67.13 years), with a sex ratio of 1:1 for both groups. For statistical analysis, we applied univariate descriptive statistics and inferential statistics (Shapiro–Wilk test, Mann–Whitney U-test for independent large samples, with the determination of the effect size coefficient r, and simple linear regression. We calculated the effect of fatigue and the Poincaré indices SD₁, SD₂, SD₁/SD₂ and the area of the fitting ellipse (AFE) for the test values of each subject.

Results

The analysis of the differences between groups revealed statistically significant results for most HGS-derived indices (p ≤ 0.05), and the magnitude of the differences indicated, in most situations, a large effect size (r > 0.5). Our results demonstrate that the proposed repetitive HGS testing indicates relevant differences between young and older healthy subjects. Through the mathematical modeling of data and the application of the concept of entropy, we provide arguments supporting this new design of HGS testing.

Conclusions

Our results indicate that the variability of HGS during serial testing, which reflects complex repetitive biomechanical functions, represents an efficient indicator for differentiation between young and older hand function patterns from an entropic perspective. In practical terms, the variability of HGS, evaluated by the new serial testing design, can be considered an attractive and relatively simple biomarker to use for gerontological studies.

Poincaré plot analysis of cerebral blood flow signals: Feature extraction and classification methods for apnea detection

OBJECTIVE

Rheoencephalography is a simple and inexpensive technique for cerebral blood flow assessment, however, it is not used in clinical practice since its correlation to clinical conditions has not yet been extensively proved. The present study investigates the ability of Poincaré Plot descriptors from rheoencephalography signals to detect apneas in volunteers.

METHODS

A group of 16 subjects participated in the study. Rheoencephalography data from baseline and apnea periods were recorded and Poincaré Plot descriptors were extracted from the reconstructed attractors with different time lags (τ). Among the set of extracted features, those presenting significant differences between baseline and apnea recordings were used as inputs to four different classifiers to optimize the apnea detection.

RESULTS

Three features showed significant differences between apnea and baseline signals: the Poincaré Plot ratio (SDratio), its correlation (R) and the Complex Correlation Measure (CCM). Those differences were optimized for time lags smaller than those recommended in previous works for other biomedical signals, all of them being lower than the threshold established by the position of the inflection point in the CCM curves. The classifier showing the best performance was the classification tree, with 81% accuracy and an area under the curve of the receiver operating characteristic of 0.927. This performance was obtained using a single input parameter, either SDratio or R.

CONCLUSIONS

Poincaré Plot features extracted from the attractors of rheoencephalographic signals were able to track cerebral blood flow changes provoked by breath holding. Even though further validation with independent datasets is needed, those results suggest that nonlinear analysis of rheoencephalography might be a useful approach to assess the correlation of cerebral impedance with clinical changes.

Quasi-Periodicities Detection Using Phase-Rectified Signal Averaging in EEG Signals as a Depth of Anesthesia Monitor

Phase-rectified signal averaging (PRSA) has been known to be a useful method to detect periodicities in non-stationary biological signals. Determination of quasi-periodicities in electroencephalogram (EEG) is a candidate for quantifying the changes in the depth of anesthesia (DOA). In this paper, DOA monitoring capacity of periodicities detected using PRSA was quantified by assessing EEG signals collected from 56 patients during surgery. The method is compared with sample entropy (SampEn), detrended fluctuation analysis (DFA), and permutation entropy (PE). The performance of quasi-periodicities defined by deceleration capacity and acceleration capacity was tested using the area under the receiver operating characteristic curve (AUC) and Pearson correlation coefficient. During the surgery, a significant difference ( ) in the quasi-periodicities was observed among three different stages under general anesthesia. There is a larger mean AUC and correlation coefficient of quasi-periodicities compared with SampEn, DFA, and PE using expert assessment of conscious level and bispectral index as the gold standard, respectively. Quasi-periodicities detected using PRSA in EEG signals are a powerful monitor of DOA and perform more accurate and robust results compared with SampEn, DFA, and PE. The results do provide a valuable reference to researchers in the field of clinical applications.

Assessment of sedation-analgesia by means of Poincaré analysis of the electroencephalogram

Monitoring the levels of sedation-analgesia may be helpful for managing patient stress on minimally invasive medical procedures. Monitors based on EEG analysis and designed to assess general anesthesia cannot distinguish reliably between a light and deep sedation. In this work, the Poincaré plot is used as a nonlinear technique applied to EEG signals in order to characterize the levels of sedation-analgesia, according to observed categorical responses that were evaluated by means of Ramsay Sedation Scale (RSS). To study the effect of high frequencies due to EMG activity, three different frequency ranges (FR1=0.5-110 Hz, FR2=0.5-30 Hz and FR3=30-110 Hz) were considered. Indexes from power spectral analysis and plasma concentration of propofol and remifentanil were also compared with the bispectral index BIS. An adaptive Neurofuzzy Inference System was applied to model the interaction of the best indexes with respect to RSS score for each analysis, and leave-one-out cross validation method was used. The ability of the indexes to describe the level of sedation-analgesia, according with the RSS score, was evaluated using the prediction probability (Pk). The results showed that the ratio SD1/SD2FR3 contains useful information about the sedation level, and SD1FR2 and SD2FR2 had the best performance classifying response to noxious stimuli. Models including parameters from Poincaré plot emerge as a good estimator of sedation-analgesia levels.

Improved spectrum analysis in EEG for measure of depth of anesthesia based on phase-rectified signal averaging

The definition of the depth of anesthesia (DOA) is still controversial and its measurement is not completely standardized in modern anesthesia. Power spectral analysis is an important method for feature detection in electroencephalogram (EEG) signals. Several spectral parameters derived from EEG have been proposed for measuring DOA in clinical applications. In the present paper, an improved method based on phase-rectified signal averaging (PRSA) is designed to improve the predictive accuracy of relative alpha and beta power, a frequency band power ratio, total power, median frequency (MF), spectral edge frequency 95 (SEF95), and spectral entropy for assessing anesthetic drug effects. Fifty-six patients undergoing general anesthesia in an operating theatre are studied. All EEG signals are continuously recorded from the awake state to the end of the recovery state and then filtered using multivariate empirical mode decomposition (MEMD). All parameters are evaluated using the commercial bispectral index (BIS) and expert assessment of conscious level (EACL), respectively. The ability to predict DOA is estimated using the area under the receiver-operator characteristics curve (AUC). All indicators based on the improved method can clearly discriminate the conscious state from the anesthetized state after filtration (p  <  0.05). A significantly larger mean AUC (p  <  0.05) shows that the improved method performs better than the conventional method to measure the DOA in most circumstances. Especially for raw EEG contaminated by artifacts, when the BIS index is used to indicate the consciousness level, the improvement is 7.37% (p  <  0.05), 9.04% (p  <  0.05), 18.46% (p  <  0.05), 27.73% (p  <  0.05), 14.65% (p  <  0.05), 2.52%, 5.38% and 6.24% (p  <  0.05) for relative alpha and beta power, power ratio, total power, MF, SEF, RE and SE, respectively. However, when the EACL is used to indicate the consciousness level, the improvement is 3.30% (p  <  0.05), 16.69% (p  <  0.05), 15.08% (p  <  0.05), 34.83% (p  <  0.05), 27.78% (p  <  0.05), 5.89% (p  <  0.05), 26.05% (p  <  0.05) and 23.42% (p  <  0.05). Spectral parameters derived from PRSA are more useful to measure the DOA in noisy cases.

Abnormal Motor Activity and Thermoregulation in a Schizophrenia Rat Model for Translational Science

Background

Schizophrenia is accompanied by altered motor activity and abnormal thermoregulation; therefore, the presence of these symptoms can enhance the face validity of a schizophrenia animal model. The goal was to characterize these parameters in freely moving condition of a new substrain of rats showing several schizophrenia-related alterations.

Methods

Male Wistar rats were used: the new substrain housed individually (for four weeks) and treated subchronically with ketamine, and naive animals without any manipulations. Adult animals were implanted with E-Mitter transponders intraabdominally to record body temperature and locomotor activity continuously. The circadian rhythm of these parameters and the acute effects of changes in light conditions were analyzed under undisturbed circumstances, and the effects of different interventions (handling, bed changing or intraperitoneal vehicle injection) were also determined.

Results

Decreased motor activity with fragmented pattern was observed in the new substrain. However, these animals had higher body temperature during the active phase, and they showed wider range of its alterations, too. The changes in light conditions and different interventions produced blunted hyperactivity and altered body temperature responses in the new substrain. Poincaré plot analysis of body temperature revealed enhanced short- and long-term variabilities during the active phase compared to the inactive phase in both groups. Furthermore, the new substrain showed increased short- and long-term variabilities with lower degree of asymmetry suggesting autonomic dysregulation.

Conclusions

In summary, the new substrain with schizophrenia-related phenomena showed disturbed motor activity and thermoregulation suggesting that these objectively determined parameters can be biomarkers in translational research.

EEG Signals Analysis Using Multiscale Entropy for Depth of Anesthesia Monitoring during Surgery through Artificial Neural Networks

In order to build a reliable index to monitor the depth of anesthesia (DOA), many algorithms have been proposed in recent years, one of which is sample entropy (SampEn), a commonly used and important tool to measure the regularity of data series. However, SampEn only estimates the complexity of signals on one time scale. In this study, a new approach is introduced using multiscale entropy (MSE) considering the structure information over different time scales. The entropy values over different time scales calculated through MSE are applied as the input data to train an artificial neural network (ANN) model using bispectral index (BIS) or expert assessment of conscious level (EACL) as the target. To test the performance of the new index's sensitivity to artifacts, we compared the results before and after filtration by multivariate empirical mode decomposition (MEMD). The new approach via ANN is utilized in real EEG signals collected from 26 patients before and after filtering by MEMD, respectively; the results show that is a higher correlation between index from the proposed approach and the gold standard compared with SampEn. Moreover, the proposed approach is more structurally robust to noise and artifacts which indicates that it can be used for monitoring the DOA more accurately.