- Department of Anesthesiology, Kyoto Prefectural University of Medicine, Kyoto, Japan
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 (PPAR = PPA_fx/PPA_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 PPAR at f5 and BIS (R2 = 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 PPAR_f5 with a non-negligible level. The logarithm of PPAR 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.
Introduction
The need to maintain a proper depth of general anesthesia (GA) during surgery is an important aspect of anesthesiology. An underdose of anesthetics increases the risk of intraoperative awareness, which may cause long-term psychological problems for the patient (Leslie et al., 2010; Bischoff and Rundshagen, 2011). An overdose of anesthetics may lead to postoperative neurocognitive dysfunction (Cole and Kharasch, 2018). Among various technologies used to optimize the depth of anesthesia, indices for monitoring processed electroencephalograms (EEG), such as bispectral index (BIS), and patient state index (PSI), are popular and have been used in various clinical studies (Fahy and Chau, 2018). Although these anesthesia depth monitoring methods are not standardized, their clinical usefulness is unquestionable as they are widely used by anesthesiologists worldwide. However, most of these devices use proprietary algorithms to measure the depth of anesthesia; moreover, clinical studies using these monitors have demonstrated controversial outcomes. Some trials reported that a BIS protocol reduced the incidence of intraoperative awareness (Ekman et al., 2004; Myles et al., 2004), while others failed to show the superiority of BIS for prevention of intraoperative awareness (Avidan et al., 2008, 2011; Lewis et al., 2019). Similarly, several clinical studies have concluded that the advantages of EEG-guided management include a lower incidence of delirium (Chan et al., 2013; Radtke et al., 2013; Luo and Zou, 2018; MacKenzie et al., 2018), whereas more recent studies failed to show the superiority of EEG-guided management (Wildes et al., 2019; Tang et al., 2020). These controversial outcomes have been explained by several factors (Abbott and Pearse, 2019; Whitlock and Avidan, 2020). Carrying out research using parameters for which the calculation algorithms are unknown raises concerns that simple comparisons cannot be performed among studies. Thus, there is a need to develop a multifaceted and non-proprietary algorithm for the evaluation of anesthesia depth.
There are several notable research approaches regarding the quantification of anesthesia depth, such as cardiorespiratory interactions in distinguishing awake from anesthetized states (Musizza et al., 2007; Kenwright et al., 2015), anesthesia-induced alterations of functional connectivity across the cortex (potentially important for both consciousness and anesthesia) (Chang et al., 2016; Li et al., 2019), and recurrence quantification analysis (Becker et al., 2010). Also, various methods based on entropy analysis have been reported (Bruhn et al., 2000b; Jiang et al., 2015; Su et al., 2016). Some of them seem useful for identifying burst suppression observed in the deep anesthesia phase (Liang et al., 2015). On the other hand, few reports of EEG parameters show the effectiveness in the shallow anesthesia phase leading to the arousal stage. As the BIS monitor combines multiple algorithms to calculate the BIS value for anesthesia depth, it is unlikely that there will be a universal index with a single EEG parameter covering all anesthesia phases at the different anesthesia depth levels.
One approach involves the Poincaré plot, which has been used to analyze various physiological signals (Tulppo et al., 1996; Carrasco et al., 2001; Golińska, 2013). The Poincaré plot is a type of recurrence plot used to quantify self-similarity in processes, usually periodic functions; it can be used to distinguish chaos from randomness by embedding a dataset into a higher-dimensional state space (Golińska, 2013). The Poincaré plot generated from EEG signals has shown correlations with the spectral edge frequency below which 95% of a given signal’s power (SEF95) is located during inhalational anesthesia (Hayashi et al., 2015a,b). However, this correlation between SEF95 and Poincaré plot indices, similar to the correlation between BIS and SEF95, was observed only when the depth of anesthesia is somewhat deep and dissociated during the awakening process from GA. Therefore, in this study, we focused on the emergence from GA, and we thought it is crucial to analyze the higher frequency band of EEG to capture the arousal state. We applied Poincaré plot analysis to processed EEG signals with finite impulse response (FIR) bandpass filters to investigate systematic changes to the Poincaré plot during GA. The Poincaré plot analysis applied to the bandpass-filtered EEG is positioned between time- and frequency domain analysis. We think it is a new method that may incorporate time- and frequency-domain analysis characteristics in univariate time-series analysis. Here, we report that the Poincaré plot-area ratio (PPAR) of gamma-band EEG constitutes a new independent parameter with sensitivity for state changes from anesthesia to arousal during the emergence process.
Materials and Methods
Anesthesia Management and Data Acquisition
All experiment protocols involving humans were conducted in accordance with the principles of the Declaration of Helsinki. The current study was approved (No. ERB-C-1074-2) by the Institutional Review Board for Human Experiments at the Kyoto Prefectural University of Medicine (IRB of KPUM), and for this non-interventional and noninvasive retrospective observational study, informed patient consent was waived by the IRB of KPUM; patients were provided an opt-out option, of which they were notified in the preoperative anesthesia clinic. In our facility, the use of a BIS monitor is routine for adult and pediatric patients who undergo surgery involving GA. The anesthesiologists in charge of management did not receive notice of the study and planned the anesthesia methods for scheduled surgeries in accordance with our facility’s standard care protocol, without any feedback regarding the on-line analysis of processed EEG signals. Patients were not premedicated before anesthesia induction, in accordance with our facility’s standard protocol. An anesthetic gas monitor (IntelliVue G5 Anesthesia Gas Module, Philips, Amsterdam, Netherlands) was used for the measurement of end-tidal sevoflurane (etSEV) and desflurane (etDES) concentrations. The end-tidal anesthetic gas concentration (etAG) was automatically recorded at 1-min intervals on a data server, then retrieved after the completion of anesthesia management. Based on the etAG data for each minute, etAG values at 20-s intervals were obtained by spline interpolation. In pediatric patients, GA was induced with an inhalational mixed gas, whereas in adult patients, rapid induction with propofol was employed (because of the irritating effects of desflurane on the respiratory tract during slow induction, all anesthesiologists in our facility refrain from the use of desflurane in children). After anesthesia induction, rocuronium (0.8–1.0 mg × kg–1) was intravenously administered, and tracheal intubation was conducted. The anesthesia was maintained with approximately 2.5% sevoflurane or 6% desflurane, small doses of fentanyl (1 μg × kg–1 per dose), and additional maintenance doses rocuronium (0.2 mg × kg–1 at intervals of 20–30 min). The timing of sugammadex administration was significantly earlier in the adult_SEV group than in the adult_DES group (p < 0.05) (Table 1). The timing of the end of surgery was significantly earlier in the adult_SEV group than in either adult_DES (p < 0.05) or ped_SEV groups (p < 0.01) (Table 1). In this EEG analysis, we focused on the 10-min period beginning around the end of surgery and ending with emergence from anesthesia; the study of EEG during the emergence process is important for preventing accidental intraoperative awakening. We analyzed the changes in correlations over time of various Poincaré plot parameters with etAG, EEG parameters such as observed BIS, SEF95, and total power (TP), and electromyography (EMG) parameter EMGlow.
Poincaré Plot Indices
A BIS Quatro sensor was mounted on the frontal region, in accordance with the manufacturer’s recommendations. The digitized EEG packets with a sampling frequency of 128 Hz were obtained through the serial output of the BIS monitor that sent a packet of sixteen sets of EEG μV data (32 bits) and eight packets per second (128 Hz); 8 s of EEG yielded 1,024 data points. Besides, a packet containing processed EEG valuables such as EMGlow (absolute power in the 70–110 Hz range, and, and values in decibel [dB] with respect to 0.0001 μV2) was obtained through the serial output of the BIS monitor once every second. For each 8-s data set, raw EEG signals were reconstructed from each digitized EEG packet and processed through six settings of FIR bandpass 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 (Figure 1A). For each FIR bandpass-filtered EEG dataset, the SEF95, TP, and power spectrum were calculated. The Poincaré plot was constructed from the FIR-filtered EEG using 8-s epochs of the EEG signal (1,024 data points). The calculation of Poincaré plot parameters using FIR-filtered EEG A pair of EEG voltages with 1/128-s time lag (the shortest time lag under 128 Hz sampling rate of BIS monitor’s EEG) was plotted in the XY plane, first at a specific time (X: x(k)) and then after a time delay (Y: x(k+1)). To analyze the distribution of EEG patterns in the Poincaré plot, the standard deviation (SD) of the EEG voltage dispersion was measured along and perpendicular to the diagonal line of identity (SD1 and SD2, respectively, Figure 2). SD1 and SD2 represent the minor and the major semi-axes of this fitted ellipse. SD1 is the standard deviation of the distances of points from axis 1 and determines the width of the ellipse (short-term variability), while SD2 equals the standard deviations from axis 2 and length of the ellipse (long-term variability). The equations for the Poincaré plot parameters are as follows (Brennan et al., 2001; Golińska, 2013; Khandoker et al., 2013):
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g001.jpg)
Figure 1. EEG Analyzer software as a real-time monitor of the Poincaré plot. The EEG Analyzer (downloadable for free use from the blog site Science to Medicine, EEG Analyzer ver. 54_GP, http://www.anesth-kpum.org/blog_ts/?p=3169) was developed using the open-source Processing 3.5 software package with the Apache Commons Mathematics Library (version 3.6.1, Apache Software Foundation, Forrest Hill, MD, United States) and a Java Virtual Machine (Oracle, Redwood Shores, CA, United States) with Java class libraries (javax.swing, java.AWT, and java.io packages). (A) Flow chart of on-line Poincaré plot analysis of processed EEG data. EEG waves obtained from the BIS VISTA A-3000 were passed through six bandpass filters for fast Fourier transform frequency analysis. Then, processed EEG data at each filtered range was subjected to Poincaré plot analysis during general anesthesia management. (B) The EEG Analyzer (connected to a VISTA A-3000 BIS monitor to collect EEG packets through an RS-232 interface) displays BIS, SEF95, total power, frequency spectrum, and the Poincaré plot of bandpass-filtered EEG data. Digitized EEG packets with a sampling frequency of 128 Hz were obtained through the serial output of the BIS monitor. Eight seconds of EEG analysis yields 1,024 data points, resulting in a Poincaré plot with 1,024 data points. The Fc-PIS shown on the screen was not yet optimized. BIS, bispectral index; EEG, electroencephalogram; Fc-PIS, frequency analysis-cooperated Poincaré plot-integrated score; SEF95, spectral edge frequency below 95% of the power of a given signal is located.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g002.jpg)
Figure 2. Poincaré plot: SD1/SD2 and Poincaré plot area (SD1 × SD2 × π). A pair of EEG voltages with 1/128-s time lag was plotted in the XY plane, first at a specific time (X: x(k)) and then after a time delay (Y: x(k+1)). SD1 represents the standard deviation of the Poincaré plot perpendicular to the line-of-identity, while SD2 represents the standard deviation of the Poincaré plot along the line-of-identity. EEG, electroencephalogram.
where SD(xk) is a standard deviation of the time series xk.
The SD1/SD2 ratio, which characterizes the sharpness of the scattered pattern, has been reportedly used for estimation of anesthesia depth (Hayashi et al., 2015a,b). However, in the present study, the Poincaré plot area (PPA) was calculated with the equation SD1 × SD2 × π (Data Sheet 1, Supplementary Method, the example codes of Poincaré plot parameters–calculations with sample data using Python in Jupyter Notebook, and Processing). Our novel approach comprised dividing the PPA of each frequency range by the PPA of the f0 range (PPA_f0) and defining the result as the PPAR (PPAR = PPA_fx/PPA_f0, fx = f1∼f5). The value of BIS was simultaneously collected from digital packets received from the BIS monitor, then recorded at 3-s intervals into the output data file.
Prior to this study, we created a software application named EEG Analyzer (downloadable for free use from our blog site, Science to Medicine, EEG Analyzer ver. 54_GP,1) (Figure 1B), through which raw EEG signals are obtained from a VISTA A-3000 BIS monitor (Application revision 3.22, Medtronic, Minneapolis, MN, United States) in the legacy mode using an RS-232 interface to a personal computer (Surface Pro 4, Microsoft Co., Redmond, WA, United States). The Processing 3.5 software package (ver. 3.5.3,2, MIT Media Lab, Massachusetts Institute of Technology, Cambridge, MA, United States) was used with the Apache Commons Mathematics Library (version 3.6.1,3, Apache Software Foundation, Forrest Hill, MD, United States) and a Java Virtual Machine (Oracle, Redwood Shores, CA, United States) with Java class libraries (javax.swing, java.AWT, and java.io packages). Through Processing’s integrated development environment, the program code for the EEG analysis was compiled using a Java Virtual Machine to build an execution file as a standalone application that functions in both Microsoft Windows 10 (Microsoft Co.) and Mac OS X (Apple, Inc. Cupertino, CA, United States).
Data Processing and Statistics
For parametric regression with curve fitting function and creation of graphs, Microsoft Excel for Mac (ver. 16.16.5, Microsoft Corp., Redmond, WA, United States) and RINEARN Graph 3D (ver. 5.6,4, RINEARN, Kyoto, Japan) were used, respectively. Changes in various EEG parameters between two-time points, 10 min before emergence (EM–10) and at the time of emergence (EM0), were statistically compared using paired t-tests with InStat 3 (GraphPad Software, San Diego, CA, United States). Root mean square error (RMSE) and coefficients of determination (R2) were calculated as estimators for the regression analyses. For group comparisons of patients’ background characteristics, one-way analysis of variance (ANOVA) was first used to evaluate variation among group means. Only when the variation among group means was greater than expected by chance (p < 0.05), Student-Newman-Keuls multiple comparisons test was performed to evaluate the mean variation between the two groups. p < 0.05 was considered significant. For time-course comparisons, Student t-tests (paired, two tail p-value) were used. The p-values < 0.05 were considered significant.
Results
Time Course of etAG, BIS, and SEF95 From Hypnosis to Emergence
Data were analyzed from 20 adults anesthetized with sevoflurane (adult_SEV group), 20 adults anesthetized with desflurane (adult_DES group), and 20 pediatric patients (aged 1–10 years) anesthetized with sevoflurane (ped_SEV group). All 60 patients from the three groups were also included in a combined group (all_combined group) (Table 1). All of the patients had the American Society of Anesthesiologists’ physical status values of 1–2 without any neurological diseases (Supplementary Table 1).
In the time-course from the hypnotic state to arousal during the 10-min process of emergence from GA (Figure 3), etSEV in the adult_SEV group, etDES in the adult_DES group, and etSEV in the ped_SEV group decreased from 0.97 ± 0.32, 3.71 ± 0.84, and 1.34 ± 0.75 to 0.17 ± 0.10, 0.68 ± 0.42, and 0.27 ± 0.16, respectively (p < 0.05 for all three groups). BIS values of the three groups increased from 52.9 ± 8.9, 41.9 ± 8.8, and 65.4 ± 14.9 to 80.5 ± 13.8, 84.7 ± 16.4, and 88.4 ± 9.1, respectively (p < 0.05 for all three groups). SEF95 increased from 16.8 ± 2.5 Hz, 13.6 ± 1.9 Hz, and 20.8 ± 4.5 Hz to 21.7 ± 4.4 Hz, 20.3 ± 4.6 Hz, and 24.1 ± 3.3 Hz, respectively (p < 0.05 for all three groups). RelativeBetaRatio (RBR), a subparameter of BIS, tended to increase (p < 0.05). EMGlow, which corresponds to the EMG power of 70–110 Hz, also increased for all three groups (p < 0.05). In the all_combined group, significant differences (p < 0.05) between EM–10 and EM0 were detected for etAG, BIS, SEF95, RBR, and EMGlow (Figure 3). Supplementary Tables 2, 3 show the detailed measurements in each of the four groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g003.jpg)
Figure 3. Time-course changes of parameters in the adult_SEV, adult_DES, ped_SEV, and all_combined groups. etAG, BIS, SEF95 RelativeBetaRatio (BetaRatio), and EMGlow of EEG. Plots of the time course from the hypnotic condition to the awake state of (1) etAG, (2) BIS, (3) SEF95, and the BIS subparameters (4) RelativeBetaRatio (BetaRatio), and (5) EMGlow. adult_SEV, adults anesthetized with sevoflurane (n = 20); adult_DES, adults anesthetized with desflurane (n = 20); ped_SEV, pediatric patients anesthetized with sevoflurane (n = 20); all_combined, all patients from the three groups (n = 60). Data are shown as mean ± SD. ∗p < 0.05, ∗∗p < 0.01, and ∗∗∗p < 0.001 by Student’s paired t-test between the values at EM–10 and EM0. BIS, bispectral index; EM–10, 10 min before emergence; EM0, at the time of emergence; etAG, end-tidal anesthetic gas concentration; SEF95, spectral edge frequency below which 95% of the power of a given signal is located; EMGlow, a BIS monitor-derived electromyography parameter corresponding to absolute power in the 70–110 Hz range, and values in decibel (dB) with respect to 0.0001 μV2.
Time Course of TP and Poincaré Plot Parameters From Hypnosis to Arousal
The time-course changes of TP were analyzed at each frequency band. In each frequency band, TP decreased over time, although small increases were observed in the last few minutes of TPf5 in the adult_SEV group, as well as in the last few minutes of TPf4 and TPf5 in the adult_DES group (p < 0.05) (Supplementary Figure 1A). Next, Poincaré plot parameters (SD1 and SD2), which were obtained from the Poincaré plot-analysis of a pair of FIR-filtered EEG μV with a 1/128-s time lag, were analyzed at each frequency band. The time-course changes of SD1/SD2 depended strongly on the frequency band: gradual increases of SD1/SD2 were observed in f0, f3, f4, and f5, while gradual decreases were observed in f1 and f2 (Supplementary Figure 1B). However, the overall mean value changes of SD1/SD2 across all frequency bands were small compared to their standard deviations.
PPA, which is calculated from SD2 × SD2 × π, changed dynamically during the emergence process (Supplementary Figure 1C). The absolute values of PPA in the ped_SEV group were approximately five times larger than those values in the adult groups. The changes in PPA were not consistent over time: the values of PPA_f4 and PPA_f5 suddenly increased over the last 2–3 min in all three groups (p < 0.05, PPA_f5). A similar increase in PPA_f0 was observed during the final phase of emergence. Next, the calculated absolute values of PPA were significantly different among the three groups. Therefore, the adjustment was performed: the ratios of PPA at each frequency range to PPA_f0 were calculated as PPAR (PPAR of each frequency-range/PPAR of the f0: 0.5–47 Hz range) throughout the time course from the hypnotic condition to the awake state. The results showed that the PPAR values for all three groups fit into the same plotting range. Drastic changes in PPAR were detected among f1–f5 across all three groups (Figure 4). PPAR_f2 and PPAR_f3 significantly decreased (p < 0.05 for all three groups), while PPAR_f5 significantly increased (p < 0.01 for all three groups). The time-course changes of PPAR_f5 demonstrated notable changes 2–3 min before emergence: although the values of PPAR_f5 remained near 0 for the first 7–8 min, the values began to increase in the final few minutes before emergence. The mean changes in TP, SD1/SD2, PPA, and PPAR values in each frequency range over time are summarized in Figure 5. Changes of parameters in f5, such as PPA_f5 and PPAR_f5, were remarkable among all five hierarchical frequency bands. In the all_combined group, for etAG, significant differences (p < 0.05) were detected in both SD1/SD2 and PPAR at all five hierarchical frequency bands (Figure 4 and Supplementary Table 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g004.jpg)
Figure 4. Time-course changes of parameters in the adult_SEV, adult_DES, ped_SEV, and all_combined groups. (1–5) Plots of the time course of the PPAR (Poincaré plot area of each frequency-range/Poincaré plot area of the f0: 0.5–47 Hz range) throughout the time course from the hypnotic condition to the awake state. (6) log10(PPAR_f5) and PIS. (7) Poincaré plot-integrated score (PIS). PIS = 25 × log10(PPAR_f5) + 112.5. adult_SEV, adults anesthetized with sevoflurane (n = 20); adult_DES, adults anesthetized with desflurane (n = 20); ped_SEV, pediatric patients anesthetized with sevoflurane (n = 20); all_combined group, all patients from the three groups (n = 60). Data are shown as mean ± SD. ∗p < 0.05, ∗∗p < 0.01, and ∗∗∗p < 0.001 by Student’s paired t-test between the values at EM–10 and EM0. EM–10, 10 min before emergence; EM0, at the time of emergence; PPAR, Poincaré plot-area ratio.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g005.jpg)
Figure 5. Time-course changes of mean values of EEG parameters in the adult_SEV, adult_DES, and ped_SEV groups, and etAG in the all_combined group. Mean values of (1) total power, (2) SD1/SD2 of the Poincaré plot, (3) Poincaré plot area (PPA), and (4) Poincaré plot-area ratio (PPAR: PPA of each frequency-range/PPA value of the f0:0.5–47 Hz range) throughout the time course from the hypnotic condition to the awake state. Data are shown as means. adult_SEV, adults anesthetized with sevoflurane (n = 20); adult_DES, adults anesthetized with desflurane (n = 20); ped_SEV, pediatric patients anesthetized with sevoflurane (n = 20); all_combined group, all patients from the three groups (n = 60). Plots of the time course range from the hypnotic condition (EM–10) to the awake state (EM0). EM–10, 10 min before emergence; EM0, at the time of emergence; PPAR, Poincaré plot-area ratio.
Correlations of SD1/SD2 With etAG and BIS, and Correlations of PPAR With etAG and BIS
Next, using 4,000 data points obtained from 20 patients, the relationships of SD1/SD2 with etAG and BIS were evaluated (Supplementary Figure 2A). Differences in the distributions of SD1/SD2, compared with etAG and BIS, are present in each frequency range plot. Analysis of the relationship between SD1/SD2 and etAG showed that when etAG decreased, SD1/SD2 tended to be dispersed across all plots. Overall, both linear and logarithmic correlations between SD1/SD2 of six frequency bands and BIS in all three groups were low (R2 ≤ 0.5). The relationships of PPAR with etAG and BIS were also evaluated (Supplementary Figure 2B). Scatter plots were created to assess PPAR at five frequency bands, compared with etAG. In plots between PPAR_f5 and etAG, hyperbolic shapes with large variation at low etAG values were observed in all three groups. Regarding the linear correlation between PPAR_f1 and BIS, the R2 values in all three groups were less than 0.44. PPAR_f2 and PPAR_f3 showed biphasic changes against BIS in all three groups. Logarithmic correlations between PPAR_f4 and BIS showed R2 values of 0.36, 0.59, and 0.27 for adult_SEV, adult_DES, and ped_SEV groups. The relationships between PPAR_f5 and BIS in the three groups are summarized in Figure 6. Regarding the logarithmic correlations between PPAR_f5 and BIS, R2 values were 0.67, 0.79, and 0.71, respectively [Figure 6A(3)]. Regarding the logarithmic correlation between PPAR_f5 and BIS in the all_combined group, R2 was 0.78 [Figure 6B(3)].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g006.jpg)
Figure 6. Relationship between Poincaré plot-area ratio at f5:30–47 Hz (PPAR_f5), BIS, and etAG. Scatter plots display correlations between (1) etAG (%) and BIS, (2) PPAR_f5 and etAG (%), and (3) PPAR_f5 and BIS (A) in the adult_SEV, adult_DES, and ped_SEV groups and (B) in the all_combined group. (4) 3D-scatter plot between end-tidal concentrations (%MAC), PPAR_f5, and BIS in the all_combined group. (C) Semi-logarithmic scatter plots between Poincaré plot-integrated score (PIS) and BIS in the adult_SEV, adult_DES, ped_SEV, and all_combined groups. PIS = 25 × log10(PPAR_f5)+112.5. Minimum alveolar concentration at 50% (MAC) values of sevoflurane and desflurane were calculated by using 1.71 and 7.25, respectively. adult_SEV, adult patients anesthetized with sevoflurane (n = 20); adult_DES, adult patients anesthetized with desflurane (n = 20); ped_SEV, pediatric patients anesthetized with sevoflurane (n = 20); all_combined, all patients from the three groups (n = 60). RINEARN Graph 3D (free software ver 5.6, https://www.rinearn.com, RINEARN, Kyoto, Japan) was used for the creation of graphs. BIS, bispectral index; etAG, end-tidal anesthetic gas; etDES, end-tidal sevoflurane concentration; etSEV, end-tidal sevoflurane concentration; PIS, Poincaré plot-integrated score; PPAR_f5, Poincaré plot-area ratio at f5:30–47 Hz; R2, coefficient of determination.
Correlations of log10(PPAR_f5) With BIS
During the emergence process, the all_combined group exhibited a gradual reduction of the average etAG from 0.61 ± 0.30% to 0.12 ± 0.08%; furthermore, the average BIS gradually increased from 53.4 ± 14.7 to 84.6 ± 13.6 [Figure 3(1, 2)]. Because the relationship between PPAR_f5 and BIS fit well to a logarithmic curve with high R2 values in all three groups, the relationship between log10(PPAR_f5) and BIS was fitted using linear regression (Figure 6C). Thus, the depth of anesthesia captured by the logarithmic value of PPAR_f5, or by BIS, changed in a nearly linear manner from the state of anesthesia to emergence. Therefore, BIS is an index that can be easily understood by anesthesiologists because it converts the depth of anesthesia into a score ranging from 0 to 100; the logarithmic value of PPAR_f5 can also be easily converted to a 0–100 range. From the data of our 60 cases, PPAR_f5 = 0.3162 (logarithmic value −0.5) was defined as Poincaré plot-integrated score (PIS) = 100, and PPAR_f5 = 0.0031623 (logarithmic value −2.5) was defined as PIS = 50. The linear regression conversion formula for calculating the PIS value from PPAR_f5 then became PIS = 25 × log10 (PPAR_f5)+112.5. Figure 4(6, 7) show the time-course changes of log10(PPAR_f5) and calculated PIS values in the four groups. R2 between PIS and BIS in the three groups, which was between 0.67 and 0.79 (Figure 6C), is mathematically identical to the R2 between the original PPAR_f5 and BIS [Figure 6A(3)]. In the all_combined group, the regression equation BIS = 1.12 × PIS −11.68 was obtained (R2 = 0.78) (Figure 6C).
Correlations of EMGlow With BIS and PPAR_f5
While, during the emergence process, the all_combined group exhibited a gradual increase of the average BIS from 53.4 ± 14.7 to 84.6 ± 13.6 [Figure 3(2)], the average EMGlow gradually increased from 2943.5 ± 465.7 to 5063.9 ± 1082.0 (Figure 3(5), Supplementary Table 2B). R2 of the logarithmic regressions between EMGlow and BIS in the three groups were between 0.24 and 0.54 [Figure 7(1)], and R2 of the exponential regressions between EMGlow and PPAR_f5 were between 0.39 and 0.64 [Figure 7(2)]. EMGlow is an absolute power at the range of 70–110 Hz. These results implied that, although the R2 of the regressions between EMGlow and PPAR_f5 was not so great as that between BIS and PPAR_f5, a part of EMG activity at the gamma range of 30–47 Hz probably influenced the calculations of BIS and PPAR_f5 with a non-negligible level.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g007.jpg)
Figure 7. Relationships of EMGlow with BIS and Poincaré plot parameter PPAR_f5. Scatter plots display correlations between (1) EMGlow and BIS and (2) EMGlow and PPAR_f5 in the adult_SEV, adult_DES, and ped_SEV groups and in the all_combined group. adult_SEV, adult patients anesthetized with sevoflurane (n = 20); adult_DES, adult patients anesthetized with desflurane (n = 20); ped_SEV, pediatric patients anesthetized with sevoflurane (n = 20); all_combined, all patients from the three groups (n = 60). RINEARN Graph 3D (free software ver 5.6, https://www.rinearn.com, RINEARN, Kyoto, Japan) was used for the creation of graphs. EMGlow, a BIS monitor-derived electromyography parameter corresponding to absolute power in the 70–110 Hz range, and values in decibel (dB) with respect to 0.0001 μV2; BIS, bispectral index; PPAR_f5, Poincaré plot-area ratio at f5:30–47 Hz; R2, coefficient of determination.
Case Examples: Video Clips of the Poincaré Plot, EEG Spectrum, and Spectrogram
In the case-by-case plots of BIS and PIS (Supplementary Figure 3), the overall time-course change of calculated PIS fit well to the change of measured BIS in each case, despite slight differences in correlations detected by RMSE and R2 among the three groups (Supplementary Table 3). One representative case from each group of the three patient groups (A: adult_SEV, B: adult_DES, and C: ped_SEV) is shown as a movie clip (Supplementary Video Clips) with spectrograms, time trends of BIS, PIS, and etSEV or etDES at 10 min before emergence from GA. Corresponding spectrograms of the above three cases from each group are shown with the corresponding time trends of BIS, PIS, etSEV, or etDES (Figure 8). Theta-to-alpha oscillations were observed before the dissipation points, followed by the emergence of beta-to-gamma oscillations with abrupt increases in BIS and PIS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627088/fphys-12-627088-HTML/image_m/fphys-12-627088-g008.jpg)
Figure 8. Spectrogram from the hypnotic state to emergence from general anesthesia. Spectrograms of three representative cases [(A) adult_SEV, (B) adult_DES, and (C) ped_SEV] from each of the three patient groups with time trends of BIS, PIS, etSEV, and etDES at 10 min before emergence from general anesthesia. Theta-to-alpha oscillations were observed before the dissipation points, after which beta-to-gamma oscillations were observed with abrupt increases of BIS and PIS. adult_SEV, adult patients anesthetized with sevoflurane; adult_DES, adult patients anesthetized with desflurane; ped_SEV, pediatric patients anesthetized with sevoflurane; BIS, bispectral index; etSEV, end-tidal sevoflurane concentration; etDES, end-tidal sevoflurane concentration; PIS, Poincaré plot-integrated score.
Discussion
Electroencephalography signals represent the total accumulation of action potentials of each nerve cell in the cerebral cortex. A frequency-dependent method is generally used to reduce the complexity of those data. The Poincaré plot is another mathematical approach for the evaluation of those complex data. There is no difference in the neurobiological basis that serves as the starting point. This study developed an EEG analyzer that can display multiple aspects of EEG signals, including unprocessed EEG waves, power spectra, spectrograms, and the Poincaré plot with accompanying parameters, on a single screen. As a result, the PPA of the Poincaré plot reflected the factors of both EEG power and complexity; therefore, PPA changed dynamically during the process of emergence from GA. The PPAR indicates the power and complexity of a specific frequency range of EEG, relative to those attributes of the full EEG frequency range. Notably, regardless of the anesthetic type and patient age, PPAR_f5 uniformly increased in the last few minutes before emergence across all three groups. The logarithmic regression results indicated that PPAR_f5 exhibited positive correlations with BIS, with high R2 values. These results suggested that, during anesthesia depth change, EEG signals contain more high-frequency gamma-wave components and become more complex with time; these tendencies were reflected in the trend of PPAR_f5, although we acknowledge an undefined theoretical foundation in the Poincaré plot analysis for processed EEG.
Because the logarithmic regression analysis between PPAR_f5 and BIS showed a positive correlation with BIS, with high R2 values, the logarithmic value of PPAR_f5 changed linearly, similar to the linear increase of BIS during the investigated 10-min period. Thus, for the period in which the change in depth of anesthesia is difficult to perceive, as the patient does not exhibit active behavior, the logarithm of PPAR_f5 could serve as a simple, independent non-proprietary indicator of GA depth. There were slight differences in regression model error between adult and pediatric patients, as well as between sevoflurane and desflurane. Many recent reports have indicated the need for new monitoring algorithms or calibration to better adjust for age-specific EEG characteristics in both pediatric (Fuentes et al., 2008; Sciusco et al., 2017; Beekoo et al., 2019; Rigouzzo et al., 2019; Wang et al., 2019) and older patients (Purdon P. et al., 2015; Ni et al., 2019; Kreuzer et al., 2020). It may be necessary to incorporate these regulatory mechanisms to obtain better correlations with anesthesia depth.
Gamma-band EEG oscillations above 40 Hz have been recognized as a critical marker of the conscious state (Desmedt and Tomberg, 1994; Franken et al., 1994; Menon et al., 1996; Gross and Gotman, 1999). The BIS proprietary algorithm includes subparameters such as burst suppression ratio, RBR, and sync fast slow (SFS), with multiple regression equations (Sigl and Chamoun, 1994; Rampil, 1998). RBR is defined as the logarithmic ratio of gamma-range spectral power (P30–47 Hz) to the spectral power of the 11–20 Hz frequency band (P11–20 Hz). Several studies have indicated that BIS values are highly dependent on RBR when the value of BIS is greater than 60 (Morimoto et al., 2004; Lee et al., 2019). SFS is the logarithm of the ratio of the bispectral power in the 40–47 Hz waveband to the bispectral power in the 0.5–47 Hz band (Rampil, 1998). At surgical levels of isoflurane anesthesia, BIS and SFS are closely correlated, as are BIS and SEF95 (Morimoto et al., 2004). A good correlation between BIS, SEF95, and SD1/SD2 of the Poincaré plot was also reported at surgical levels of sevoflurane inhalational anesthesia (Hayashi et al., 2015a). However, these strong correlations among BIS, SEF95, and SD1/SD2 of the Poincaré plot have probably been only observed at the surgical level of anesthesia, but not at the shallow anesthetic level during the arousal phase.
Previous reports mentioned that increased BIS was related to an EMG signal associated with muscle relaxation in the awake state (Bruhn et al., 2000a; Sleigh et al., 2001; Messner et al., 2003; Schuller et al., 2015), or reversal of muscle relaxation by sugammadex in the hypnotic state (Aho et al., 2012). Therefore, EMG power may be suitable as an index of the BIS range. Including the non-negligible influence of frontal EMG on the high-frequency band of frontal EEG, the BIS values generally reflect changes in gamma-frequency EEG components, which probably contaminates with EMG at the level that cannot be ignored. In addition to BIS, the spectral entropy monitor uses EEG irregularity and computes state entropy over the frequency range of 0.8–32 Hz, reflecting the EEG-dominant part of the spectrum (Ellerkmann et al., 2004). Additionally, response entropy was computed over the frequency band of 0.8–47 Hz, including both the EEG- and EMG-dominant parts of the recorded spectrum Ellerkmann et al., 2004). Our main result—PPA_f5 or logarithmic PPA_f5 had the best correlation with BIS during the emergence process—seems reasonable, considering that both BIS and spectral entropy use measurements of gamma-frequency EEG as critical components of subparameters that represent anesthetic depth. Although one report demonstrated that the average change in EMG signal was approximately one-tenth of the magnitude of the EEG signal when the submental EMG signal was used to estimate the frontal EMG signal (Sleigh et al., 2001), it is necessary to understand that EEG in the gamma region utilized by BIS and our PPAR_f5 are highly affected by the non-negligible power of EMG.
A limitation of this study was that it focused solely on the emergence process from GA; the study did not include a parametric analysis involving patients in a more profound state of GA. Because the BIS uses at least three subparameters of anesthesia state, the combination of SEF95 with the PPAR could be used to create a new depth-of-anesthesia monitor with a known algorithm. We acknowledge that the finding associated with PPAR_f5 may only be applicable for anesthetic phases ranging between light anesthesia level and awareness. In this study, poor correlations were observed between etAG and anesthetic depth parameters, such as BIS and PPAR_f5. Because we used a semi-closed circuit for GA, measured end-tidal anesthetic concentration would not reflect the “true” concentration. Pharmacokinetics simulation of putative brain concentration of anesthetic during emergence from anesthesia will sure lead to more appropriate assessment of anesthetic action on EEG changes, although we could not handle it for technical reasons in this study. The influence of EMG from the forehead also requires further examination through case-control studies because the effects of EMG on BIS have been evaluated in awake volunteers (Messner et al., 2003; Schuller et al., 2015). The use of standard parameters derived from a non-proprietary algorithm to indicate anesthetic depth could clarify outcomes of clinical trials that have targeted prevention of awareness or delirium by using different parameters derived from specific proprietary algorithms (Abbott and Pearse, 2019; Whitlock and Avidan, 2020). Another limitation of this study involved its observational design; thus, we could not eliminate the influences of noxious stimuli from surgical procedures near the end of surgery. Variation in opioid (fentanyl) dosage may also have affected patient awareness. Further analyses (e.g., controlled trials) could enhance the clinical utility of evaluating the PPAR as an indicator of anesthetic depth. Ultimately, we aim to propose an open protocol for an anesthetic depth algorithm that can be used by anesthesiologists.
The clinical use of EEG-based indices for anesthetic depth monitoring remains controversial due to the lack of a neurophysiological interaction between drug-specific neurophysiological signatures and the meaning of the indices, as noted in a recent educational review (Purdon P.L. et al., 2015). We currently agree with the importance of evaluation of unprocessed EEG data; spectrograms have received greater emphasis in terms of the interpretation of anesthetic-induced brain states defined by drug-specific neurophysiological signatures. Furthermore, the relationship between age and postoperative cognitive impairment is becoming an important consideration for measuring anesthesia depth (Purdon P. et al., 2015; Ni et al., 2019; Kreuzer et al., 2020). Among the various components of EEG-guided anesthesia, our results showed an important aspect, such that Poincaré plot parameters obtained from EEG signals processed through gamma-band frequency filters revealed a notable change during emergence from GA, although it cannot be denied the influence of EMG power in this EEG range. In conclusion, our results suggest that under inhalational GA, the logarithm of PPAR_f5 (PIS) could serve as a non-proprietary anesthetic depth parameter independent of the proprietary BIS algorithm.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by The Institutional Review Board for Human Experiments at the Kyoto Prefectural University of Medicine. Written informed consent from the participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.
Author Contributions
KH was the first author and performed the data collection and data analysis. AK, KA, MK, and MS performed manuscript preparation. TS was the principal investigator responsible for study conception and design, data analysis, and manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
Support was provided solely from departmental sources.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Kazuko Hayashi, M.D., Ph.D. of the Department of Anesthesiology, Chubu Medical Center, Kyoto, Japan for technical advice about Poincaré plot analysis. We thank Richard Lipkin, Ph.D., and Ryan Chastain-Gross, Ph.D., from Edanz Group (https://en-author-services.edanzgroup.com/ac) for editing a draft of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.627088/full#supplementary-material
Supplementary Figure 1 | Total EEG power and Poincaré plot parameters throughout the time course from the hypnotic condition to the awake state in anesthetized patients. (A) Total preprocessed EEG power after passing through one of five bandpass filters, (B) SD1/SD2 of the Poincaré plot of the preprocessed EEG signal after passing through one of six bandpass filters, (C) Area of the Poincaré plot (PPA) of the processed EEG signal after passing through one of six bandpass filters.
Supplementary Figure 2 | Relationships of Poincaré plot parameters with BIS and etAG. (A) The relationships of Poincaré plot SD1/SD2 with BIS and etAG. (B) The relationships of Poincaré plot Poincaré plot area ratio (PPAR) with BIS and etAG.
Supplementary Figure 3 | Case-by-case plots of time-course change of PIS and BIS, and linear regressions between PIS and BIS. (A) Adult sevoflurane (adult_SEV) group. (B) Adult sevoflurane (adult_DES) group. (C) Adult sevoflurane (ped_SEV) group.
Supplementary Table 1 | Detailed information about patients and anesthesia management.
Supplementary Table 2 | Comparison of changes in EEG parameters. (A) The comparison between three groups. (B) The comparison in the all_combined group.
Supplementary Table 3 | Evaluation of linear regression of the three groups.
Supplementary Data Sheet 1 | The example codes of Poincaré plot parameters – calculations with sample data using Python in Jupyter Notebook, and Processing.
Supplementary Video Clips (Separate files, three clips) | one case from each group, 10 min movie clips:
Three representative cases from each group are shown as movie clips (35–36 s), which were displayed by decoding from the recorded binary files under “REVIEW mode” of the EEG analyzer and illustrate a 10 min period from the hypnotic state to arousal. The movie clips include the screen of the EEG analyzer (ver. 49_SR), including the Poincaré plot of processed EEG through specific bandpass filters and the EEG power spectrum. (Please ignore the Fc-PIS values on the screen because the algorithm for the calculation of Fc-PIS value was not optimized when we recorded the analyses. Fc-PIS: frequency analysis-cooperated Poincaré plot-integrated score).
Supplementary Video 1 | Supplementary_Video_1_adult_SEv_case4_ EEGanalyzer.mp4, EEG analyzer screen (10 min) of adult_SEV case #4.
Supplementary Video 2 | Supplementary_Video_2_adult_DES_ case2_EEGanalyzer.mp4, EEG analyzer screen (10 min) of adult_DES case #2.
Supplementary Video 3 | Supplementary_Video_3_ped_SEV_case5_ EEGanalyzer.mp4, EEG analyzer screen (10 min) of ped_SEV case #5.
Footnotes
- ^ http://www.anesth-kpum.org/blog_ts/?p=3169
- ^ https://processing.org
- ^ https://commons.apache.org
- ^ https://www.rinearn.com
References
Abbott, T. E. F., and Pearse, R. M. (2019). Depth of anesthesia and postoperative delirium. JAMA 321, 459–460. doi: 10.1001/jama.2019.0164
Aho, A. J., Kamata, K., Yli-Hankala, A., Lyytikäinen, L. P., Kulkas, A., and Jäntti, V. (2012). Elevated BIS and entropy values after sugammadex or neostigmine: an electroencephalographic or electromyographic phenomenon? Acta Anaesthesiol. Scand. 56, 465–473. doi: 10.1111/j.1399-6576.2011.02647.x
Avidan, M. S., Zhang, L., Burnside, B. A., Finkel, K. J., Searleman, A. C., Selvidge, J. A., et al. (2008). Anesthesia awareness and the bispectral index. N. Engl. J. Med. 358, 1097–1108. doi: 10.1056/NEJMoa0707361
Avidan, M. S., Jacobsohn, E., Glick, D., Burnside, B. A., Zhang, L., Villafranca, A., et al. (2011). Prevention of intraoperative awareness in a high-risk surgical population. N. Engl. J Med. 365, 591–600. doi: 10.1056/NEJMoa1100403
Becker, K., Schneider, G., Eder, M., Ranft, A., Kochs, E. F., Zieglgänsberger, W., et al. (2010). Anaesthesia monitoring by recurrence quantification analysis of EEG data. PLoS One 5:e8876. doi: 10.1371/journal.pone.0008876
Beekoo, D., Yuan, K., Dai, S., Chen, L., Di, M., Wang, S., et al. (2019). Analyzing electroencephalography (EEG) waves provides a reliable tool to assess the depth of sevoflurane anesthesia in pediatric patients. Med. Sci. Monit. 25, 4035–4040. doi: 10.12659/MSM.915640
Bischoff, P., and Rundshagen, I. (2011). Awareness under general anesthesia. Dtsch Arztebl. Int. 108, 1–7. doi: 10.3238/arztebl.2011.0001
Brennan, M., Palaniswami, M., and Kamen, P. (2001). Do existing measures of Poincaré plot geometry reflect nonlinear features of heart rate variability. IEEE Trans. Biomed. Eg. 48, 1342–1347. doi: 10.1109/10.9593304
Bruhn, J., Bouillon, T. W., and Shafer, S. L. (2000a). Electromyographic activity falsely elevates the bispectral index. Anesthesiology 92, 1485–1487. doi: 10.1097/00000542-200005000-200005042
Bruhn, J., Röpcke, H., and Hoeft, A. (2000b). Approximate entropy as an electroencephalographic measure of anesthetic drug effect during desflurane anesthesia. Anesthesiology 92, 715–726. doi: 10.1097/00000542-200003000-00016
Carrasco, S., Gaitan, M. J., Gonzalez, R., and Yanez, O. (2001). Correlation among Poincaré plot indexes and time and frequency domain measures of heart rate variability. J. Med. Eng. Technol. 25, 240–248. doi: 10.1080/03091900110086651
Chan, M. T., Cheng, B. C., Lee, T. M., and Gin, T., Coda Trial, and Group. (2013). BIS-guided anesthesia decreases postoperative delirium and cognitive decline. J. Neurosurg. Anesthesiol. 25, 33–42. doi: 10.1097/ANA.0b013e3182712fba
Chang, C., Raven, E. P., and Duyn, J. H. (2016). Brain-heart interactions: challenges and opportunities with functional magnetic resonance imaging at ultra-high field. Philos. Trans. A Math. Phys. Eng. Sci. 374:20150188. doi: 10.1098/rsta.2015.0188
Cole, D. J., and Kharasch, E. D. (2018). Postoperative brain function: toward a better understanding and the American Society of Anesthesiologists perioperative brain health initiative. Anesthesiology 129, 861–863. doi: 10.1097/ALN.0000000000002085
Desmedt, J. E., and Tomberg, C. (1994). Transient phase-locking of 40 Hz electrical oscillations in prefrontal and parietal human cortex reflects the process of conscious somatic perception. Neurosci. Lett. 168, 126–129. doi: 10.1016/0304-3940(94)90432-4
Ekman, A., Lindholm, M. L., Lennmarken, C., and Sandin, R. (2004). Reduction in the incidence of awareness using BIS monitoring. Acta Anaesthesiol. Scand. 48, 20–26. doi: 10.1111/j.1399-6576.2004.00260.x
Ellerkmann, R. K., Liermann, V. M., Alves, T. M., Wenningmann, I., Kreuer, S., Wilhelm, W., et al. (2004). Spectral entropy and bispectral index as measures of the electroencephalographic effects of sevoflurane. Anesthesiology 101, 1275–1282. doi: 10.1097/00000542-200412000-00006
Fahy, B. G., and Chau, D. F. (2018). The technology of processed electroencephalogram monitoring devices for assessment of depth of anesthesia. Anesth. Analg. 126, 111–117. doi: 10.1213/ANE.0000000000002331
Franken, P., Dijk, D. J., Toblerm, I., and Borbely, A. A. (1994). High-frequency components of the rat electrocorticogram are modulated by the vigilance states. Neurosci. Lett. 167, 89–92. doi: 10.1016/0304-3940(94)91034-0
Fuentes, R., Cortinez, L. I., Struys, M. M., Delfino, A., and Munoz, H. (2008). The dynamic relationship between end-tidal sevoflurane concentrations, bispectral index, and cerebral state index in children. Anesth. Analg. 107, 1573–1578. doi: 10.1213/ane.0b013e318181ef88
Golińska, A. K. (2013). Poincaré plots in analysis of selected biomedical signals. Stud. Logic Grammar Rhetoric. 35, 117–127. doi: 10.2478/slgr-2013-0031
Gross, D. W., and Gotman, J. (1999). Correlation of high-frequency oscillations with the sleep-wake cycle and cognitive activity in humans. Neuroscience 94, 1005–1018. doi: 10.1016/s0306-4522(99)00343-7
Hayashi, K., Mukai, N., and Sawa, T. (2015a). Poincaré analysis of the electroencephalogram during sevoflurane anesthesia. Clin. Neurophysiol. 126, 404–411. doi: 10.1016/j.clinph.2014.04.019
Hayashi, K., Yamada, T., and Sawa, T. (2015b). Comparative study of Poincaré plot analysis using short electroencephalogram signals during anaesthesia with spectral edge frequency 95 and bispectral index. Anaesthesia 70, 310–317. doi: 10.1111/anae.12885
Jiang, G. J. A., Fan, S., Abbod, M. F., Huang, H., Lan, J., Tsai, F., et al. (2015). Sample entropy analysis of EEG signals via artificial neural networks to model patients’ consciousness level based on anesthesiologists experience. Biomed. Res. Int. 2015:343478. doi: 10.1155/2015/343478
Kenwright, D. A., Bernjak, A., Draegni, T., Dzeroski, S., Entwistle, M., Horvat, M., et al. (2015). The discriminatory value of cardiorespiratory interactions in distinguishing awake from anaesthetised states: a randomised observational study. Anaesthesia 70, 1356–1368. doi: 10.1111/anae.13208
Khandoker, A. H., Karmakar, C., Brennan, M., Palaniswami, M., and Voss, A. (2013). “Quantitative Poincaré Plot,” in Poincaré Plot Methods for Heart Rate Variability Analysis, eds A. H. Khandoker, C. Karmakar, M. Brennan, M. Palaniswami, and A. Voss (New York, NY: Springer Science+Business Media), 13–24. doi: 10.1007/978-1-4614-7375-6
Kreuzer, M., Stern, M. A., Hight, D., Berger, S., Schneider, G., Sleigh, J. W., et al. (2020). Spectral and entropic features are altered by age in the electroencephalogram in patients under sevoflurane anesthesia. Anesthesiology 132, 1003–1016. doi: 10.1097/ALN.0000000000003182
Lee, H. C., Ryu, H., Yoon, S. B., Yang, S. M., Oh, H., et al. (2019). Data driven investigation of bispectral index algorithm. Sci. Rep. 9:13769. doi: 10.1038/s41598-019-50391-x
Leslie, K., Chan, M. T., Myles, P. S., Forbes, A., and McCulloch, T. J. (2010). Posttraumatic stress disorder in aware patients from the B-aware trial. Anesth. Analg. 110, 823–828. doi: 10.1213/ANE.0b013e3181b8b6ca
Lewis, S. R., Pritchard, M. W., Fawcett, L. J., and Punjasawadwong, Y. (2019). Bispectral index for improving intraoperative awareness and early postoperative recovery in adults. Cochrane Database Syst. Rev. 9:CD003843. doi: 10.1002/14651858.CD003843
Li, D., Vlisides, P. E., Kelz, M. B., Avidan, M. S., and Mashour, G. A., ReCCognition Study, et al. (2019). Dynamic cortical connectivity during general anesthesia in healthy volunteers. Anesthesiology 130, 870–884. doi: 10.1097/ALN.0000000000002656
Liang, Z., Wang, Y., Sun, X., Li, D., Voss, L. J., Sleigh, J. W., et al. (2015). EEG entropy measures in anesthesia. Front. Comput. Neurosci. 9:16. doi: 10.3389/fncom.2015.00016
Luo, C., and Zou, W. (2018). Cerebral monitoring of anaesthesia on reducing cognitive dysfunction and postoperative delirium: a systematic review. J. Int. Med. Res. 46, 4100–4110. doi: 10.1177/0300060518786406
MacKenzie, K. K., Britt-Spells, A. M., Sands, L. P., and Leung, J. M. (2018). Processed electroencephalogram monitoring and postoperative delirium: a systematic review and meta-analysis. Anesthesiology 129, 417–427. doi: 10.1097/ALN.0000000000002323
Menon, V., Freeman, W. J., Cutillo, B. A., Desmond, J. E., Ward, M. F., Bressler, S. L., et al. (1996). Spatio-temporal correlations in human gamma band electrocorticograms. Electroencephalogr. Clin. Neurophysiol. 98, 89–102. doi: 10.1016/0013-4694(95)00206-5
Messner, M., Beese, U., Romstock, J., Dinkel, M., and Tschaikowsky, K. (2003). The bispectral index declines during neuromuscular block in fully awake persons. Anesth. Analg. 97, 488–491. doi: 10.1213/01.ane.0000072741.78244.c0
Morimoto, Y., Hagihira, S., Koizumi, Y., Ishida, K., Matsumoto, M., and Sakabe, T. (2004). The relationship between bispectral index and electroencephalographic parameters during isoflurane anesthesia. Anesth. Analg. 98, 1336–1340. doi: 10.1213/01.ane.0000105867.17108.b6
Musizza, B., Stefanovska, A., McClintock, P. V. E., Palus, M., Petrovcic, J., Ribaric, S., et al. (2007). Interactions between cardiac, respiratory and EEG-delta oscillations in rats during anaesthesia. J. Physiol. 580, 315–326. doi: 10.1113/jphysiol.2006.126748
Myles, P. S., Leslie, K., McNeil, J., Forbes, A., and Chan, M. T. (2004). Bispectral index monitoring to prevent awareness during anaesthesia: the B-Aware randomised controlled trial. Lancet 363, 1757–1763. doi: 10.1016/S0140-6736(04)16300-9
Ni, K., Cooter, M., Gupta, D. K., Thomas, J., Hopkins, T. J., Miller, T. E., et al. (2019). Paradox of age: older patients receive higher age-adjusted minimum alveolar concentration fractions of volatile anaesthetics yet display higher bispectral index values. Brit. J. Anaesth. 123, 288–297. doi: 10.1016/j.bja.2019.05.040
Purdon, P., Pavone, K. J., Akeju, O., Smith, A. C., Sampson, A. L., Lee, J., et al. (2015). The ageing brain: age-dependent changes in the electroencephalogram during propofol and sevoflurane general anaesthesia. Brit. J. Anaesth. 115, i46–i57. doi: 10.1093/bja/aev213
Purdon, P. L., Sampson, A., Pavone, K. J., and Brown, E. N. (2015). Clinical electroencephalography for anesthesiologists. Part I: background and basic signatures. Anesthesiology 123, 937–960. doi: 10.1097/ALN.0000000000000841
Radtke, F. M., Franck, M., Lendner, J., Krüger, S., Wernecke, K. D., and Spies, C. D. (2013). Monitoring depth of anaesthesia in a randomized trial decreases the rate of postoperative delirium but not postoperative cognitive dysfunction. Br. J. Anaesth. 110, i98–i105. doi: 10.1093/bja/aet055
Rampil, I. J. (1998). A primer for EEG signal processing in anesthesia. Anesthesiology 89, 980–1002. doi: 10.1097/00000542-199810000-00023
Rigouzzo, A., Khoy-Ear, L., Laude, D., Louvet, N., Moutard, M. L., et al. (2019). EEG profiles during general anesthesia in children: A comparative study between sevoflurane and propofol. Paediatr. Anaesth. 29, 250–257. doi: 10.1111/pan.13579
Sigl, J. C., and Chamoun, N. G. (1994). An introduction to bispectral analysis for the electroencephalogram. J. Clin. Monit. 10, 392–404. doi: 10.1007/BF01618421
Schuller, P., Newell, S., Strickland, P., and Barry, J. (2015). Response of bispectral index to neuromuscular block in awake volunteers. Brit. J. Anaesth. 115, i95–i103. doi: 10.1093/bja/aev072
Sciusco, A., Standing, J. F., Sheng, Y., Raimondo, P., Cinnella, G., and Dambrosio, M. (2017). Effect of age on the performance of bispectral and entropy indices during sevoflurane pediatric anesthesia: a pharmacometric study. Paediatr. Anaesth. 27, 399–408. doi: 10.1111/pan.13086
Sleigh, J. W., Steyn-Ross, D. A., Steyn-Ross, M. L., Williams, M. L., and Smith, P. (2001). Comparison of changes in electroencephalographic measures during induction of general anaesthesia: influence of the gamma frequency band and electromyogram signal. Br. J. Anaesth. 86, 50–58. doi: 10.1093/bja/86.1.50
Su, C., Liang, Z., Li, X., Li, D., Li, Y., and Ursino, M. (2016). A comparison of multiscale permutation entropy measures in on-line depth of anesthesia monitoring. PLoS One 11:e0164104. doi: 10.1371/journal.pone.0164104
Tang, C. J., Jin, Z., Sands, L. P., Pleasants, D., Tabatabai, S., Hong, Y., et al. (2020). ADAPT-2: a randomized clinical trial to reduce intraoperative EEG suppression in older surgical patients undergoing major noncardiac surgery. Anesth. Analg. 131, 1228–1236. doi: 10.1213/ANE.0000000000004713
Tulppo, M. P., Makikallio, T. H., Takala, T. E., Seppanen, T., and Huikuri, H. V. (1996). Quantitative beat-to-beat analysis of heart rate dynamics during exercise. Am. J. Physiol. 271, H244–H252. doi: 10.1152/ajpheart.1996.271.1.H244
Wang, F., Zhang, J., Yu, J., Tian, M., Cui, X., and Wu, A. (2019). Variation of bispectral index in children aged 1-12 years under propofol anesthesia: an observational study. BMC Anesthesiol. 19:145. doi: 10.1186/s12871-019-0815-6
Whitlock, E. L., and Avidan, M. S. (2020). Three blind mice: a tail of discordant trials. Br. J. Anaesth. 124, 121–125. doi: 10.1016/j.bja.2019.09.035
Wildes, T. S., Mickle, A. M., Abdallah, A. B., Maybrier, H. R., Oberhaus, J., Budelier, T. P., et al. (2019). Effect of electroencephalography-guided anesthetic administration on postoperative delirium among older adults undergoing major surgery: The ENGAGES randomized clinical trial. JAMA 321, 473–483. doi: 10.1001/jama.2018.22005
Keywords: anesthesia, electroencephalography, emergence, Poincaré plot, anesthesia depth
Citation: Hayase K, Kainuma A, Akiyama K, Kinoshita M, Shibasaki M and Sawa T (2021) Poincaré Plot Area of Gamma-Band EEG as a Measure of Emergence From Inhalational General Anesthesia. Front. Physiol. 12:627088. doi: 10.3389/fphys.2021.627088
Received: 10 November 2020; Accepted: 05 January 2021;
Published: 09 February 2021.
Edited by:
Ahsan H. Khandoker, Khalifa University, United Arab EmiratesReviewed by:
Satoshi Hagihira, Osaka University, JapanWajid Aziz, University of Azad Jammu and Kashmir, Pakistan
Copyright © 2021 Hayase, Kainuma, Akiyama, Kinoshita, Shibasaki and Sawa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Teiji Sawa, anesth@koto.kpu-m.ac.jp