Detection of quadratic phase coupling by cross-bicoherence and spectral Granger causality in bifrequencies interactions

Quadratic Phase Coupling (QPC) serves as an essential statistical instrument for evaluating nonlinear synchronization within multivariate time series data, especially in signal processing and neuroscience fields. This study explores the precision of QPC detection using numerical estimates derived from cross-bicoherence and bivariate Granger causality within a straightforward, yet noisy, instantaneous multiplier model. It further assesses the impact of accidental statistically significant bifrequency interactions, introducing new metrics such as the ratio of bispectral quadratic phase coupling and the ratio of bivariate Granger causality quadratic phase coupling. Ratios nearing 1 signify a high degree of accuracy in detecting QPC. The coupling strength between interacting channels is identified as a key element that introduces nonlinearities, influencing the signal-to-noise ratio in the output channel. The model is tested across 59 experimental conditions of simulated recordings, with each condition evaluated against six coupling strength values, covering a wide range of carrier frequencies to examine a broad spectrum of scenarios. The findings demonstrate that the bispectral method outperforms bivariate Granger causality, particularly in identifying specific QPC under conditions of very weak couplings and in the presence of noise. The detection of specific QPC is crucial for neuroscience applications aimed at better understanding the temporal and spatial coordination between different brain regions.

Heart Rate Variability as a Tool for Seizure Prediction: A Scoping Review

The most critical burden for People with Epilepsy (PwE) is represented by seizures, the unpredictability of which severely impacts quality of life. The design of real-time warning systems that can detect or even predict ictal events would enhance seizure management, leading to high benefits for PwE and their caregivers. In the past, various research works highlighted that seizure onset is anticipated by significant changes in autonomic cardiac control, which can be assessed through heart rate variability (HRV). This manuscript conducted a scoping review of the literature analyzing HRV-based methods for detecting or predicting ictal events. An initial search on the PubMed database returned 402 papers, 72 of which met the inclusion criteria and were included in the review. These results suggest that seizure detection is more accurate in neonatal and pediatric patients due to more significant autonomic modifications during the ictal transitions. In addition, conventional metrics are often incapable of capturing cardiac autonomic variations and should be replaced with more advanced methodologies, considering non-linear HRV features and machine learning tools for processing them. Finally, studies investigating wearable systems for heart monitoring denoted how HRV constitutes an efficient biomarker for seizure detection in patients presenting significant alterations in autonomic cardiac control during ictal events.

Pairwise and higher-order measures of brain-heart interactions in children with temporal lobe epilepsy

OBJECTIVE

While it is well-known that epilepsy has a clear impact on the activity of both the central nervous system (CNS) and the autonomic nervous system (ANS), its role on the complex interplay between CNS and ANS has not been fully elucidated yet. In this work, pairwise and higher-order predictability measures based on the concepts of Granger causality (GC) and Partial Information Decomposition (PID) were applied on time series of electroencephalographic (EEG) brain wave amplitude and heart rate variability (HRV) in order to investigate directed brain-heart interactions associated with the occurrence of focal epilepsy.

APPROACH

HRV and the envelopes of δ and α EEG activity recorded from ipsilateral (ipsi-EEG) and contralateral (contra-EEG) scalp regions were analyzed in 18 children suffering from temporal lobe epilepsy monitored during pre-ictal, ictal and post-ictal periods. After linear parametric model identification, we compared pairwise GC measures computed between HRV and a single EEG component with PID measures quantifying the unique, redundant and synergistic information transferred from ipsi-EEG and contra-EEG to HRV.

MAIN RESULTS

The analysis of GC revealed a dominance of the information transfer from EEG to HRV and negligible transfer from HRV to EEG, suggesting that CNS activities drive the ANS modulation of the heart rhythm, but did not evidence clear differences between δ and α rhythms, ipsi-EEG and contra-EEG, or pre- and post-ictal periods. On the contrary, PID revealed that epileptic seizures induce a reorganization of the interactions from brain to heart, as the unique predictability of HRV originated from the ipsi-EEG for the δ waves and from the contra-EEG for the α waves in the pre-ictal phase, while these patterns were reversed after the seizure.

SIGNIFICANCE

These results highlight the importance of considering higher-order interactions elicited by PID for the study of the neuro-autonomic effects of focal epilepsy, and may have neurophysiological and clinical implications.

The Assessment of Autonomic Nervous System Activity Based on Photoplethysmography in Healthy Young Men

Noninvasive assessment of autonomic nervous system (ANS) activity is of great importance, but the accuracy of the method used, which is primarily based on electrocardiogram-derived heart rate variability (HRV), has long been suspected. We investigated the feasibility of photoplethysmography (PPG) in ANS evaluation. Data of 32 healthy young men under four different ANS activation patterns were recorded: baseline, slow deep breathing (parasympathetic activation), cold pressor test (peripheral sympathetic activation), and mental arithmetic test (cardiac sympathetic activation). We extracted 110 PPG-based features to construct classification models for the four ANS activation patterns. Using interpretable models based on random forest, the main PPG features related to ANS activation were obtained. Results showed that pulse rate variability (PRV) exhibited similar changes to HRV across the different experiments. The four ANS patterns could be better classified using more PPG-based features compared with using HRV or PRV features, for which the classification accuracies were 0.80, 0.56, and 0.57, respectively. Sensitive features of parasympathetic activation included features of nonlinear (sample entropy), frequency, and time domains of PRV. Sensitive features of sympathetic activation were features of the amplitude and frequency domain of PRV of the PPG derivatives. Subsequently, these sensitive PPG-based features were used to fit the improved HRV parameters. The fitting results were acceptable (p < 0.01), which might provide a better method of evaluating ANS activity using PPG.

Characterizing the orientation selectivity in V1 and V4 of macaques by quadratic phase coupling

OBJECTIVE

Orientation selectivity is one of the significant characteristics of neurons in the primary visual cortex (V1). Some neurons in extrastriate visual cortical areas also exhibit certain orientation selectivity. But it is still not well understood that how the orientation selectivity generates. Most previous studies about the orientation selectivity are based on the spike firing rate. However, the spikes are prone to be biased by the detection and sorting algorithms. Then, in this paper, the local field potential (LFP) is adopted to investigate the mechanism of orientation selectivity.

APPROACH

We used the quadratic phase coupling (QPC), which was calculated by wavelet bicoherence, to describe the characteristics of orientation selectivity available in V1 and V4. The raw wideband neural signals were recorded by two chronically implanted multi-electrode arrays, which were placed in V1 and V4 respectively in two macaques performing a selective visual attention task.

MAIN RESULTS

There is a strong correlation between the total bicoherence (TotalBic), which is a quantization for the overall QPC of frequency pairs in gamma band, and the grating orientation. Furthermore, the QPC distribution at the non-preferred orientation is mainly concentrated in the low frequencies (30-40 Hz) of gamma; while the QPC distribution at the preferred orientation concentrates in both the low frequencies and high frequencies (60-80 Hz) of gamma. In addition, the TotalBic of the gamma-band LFP between V1 and V4 varies with the grating orientations, indicating that the QPC is available in the feedforward link and the gamma-band LFP in V1 modulates the QPC in V4.

SIGNIFICANCE

The QPC reflects the orientations of the sinusoidal grating and describes the interaction of gamma-band LFP between different brain regions. Our results suggest that the QPC is an alternative avenue to explore the mechanism for generating orientation selectivity of visual neurons effectively.

An improved method to evaluate heart rate variability based on time-variant cardiorespiratory relation

Frequency domain analysis of heart rate variability (HRV) is a noninvasive method to evaluate the autonomic nervous system (ANS), but the traditional parameters of HRV, i.e., the power spectra of the high-frequency (HF) and low-frequency bands (LF), cannot estimate the activity of the parasympathetic (PNS) and sympathetic nervous systems (SNS) well. The aim of our study was to provide a corrected method to better distinguish the contributions of the PNS and SNS in the HRV spectrum. Respiration has a gating effect on cardiac vagal efferent activity, which induces respiration-locked heart rate (HR) changes because of the fast effect of the PNS. So the respiration-related heart rate (HR_r) is closely related to PNS activity. In this study, HR was decomposed into HR_r and the respiration-unrelated component (HR_ru) based on empirical mode decomposition (EMD) and the relationship between HR and respiration. Time-frequency analysis of HR_r and HR_ru was defined as HF_r and LF_ru, respectively, with specific adaptive bands for every signal. Two experimental data sets, representing SNS and PNS activation, respectively, were used for efficiency analysis of our method. Our results show that the corrected HRV predicted ANS activity well. HF_r could be an index of PNS activity, LF_ru mainly reflected SNS activity, and LF_ru/HF_r could be more accurate in representing the sympathovagal balance.NEW & NOTEWORTHY This study includes the time-varying relationship between respiration and heart rate in the analysis of heart rate variability. Correction for low-frequency and high-frequency components based on respiration significantly improved evaluation of the sympathetic and parasympathetic nervous systems.

Autonomic impairment of patients in coma with different Glasgow coma score assessed with heart rate variability

PRIMARY OBJECTIVE

The objective of this study is to assess the functional state of the autonomic nervous system in healthy individuals and in individuals in coma using measures of heart rate variability (HRV) and to evaluate its efficiency in predicting mortality.

DESIGN AND METHODS

Retrospective group comparison study of patients in coma classified into two subgroups, according to their Glasgow coma score, with a healthy control group. HRV indices were calculated from 7 min of artefact-free electrocardiograms using the Hilbert-Huang method in the spectral range 0.02-0.6 Hz. A special procedure was applied to avoid confounding factors. Stepwise multiple regression logistic analysis (SMLRA) and ROC analysis evaluated predictions.

RESULTS

Progressive reduction of HRV was confirmed and was associated with deepening of coma and a mortality score model that included three spectral HRV indices of absolute power values of very low, low and very high frequency bands (0.4-0.6 Hz). The SMLRA model showed sensitivity of 95.65%, specificity of 95.83%, positive predictive value of 95.65%, and overall efficiency of 95.74%.

CONCLUSIONS

HRV is a reliable method to assess the integrity of the neural control of the caudal brainstem centres on the hearts of patients in coma and to predict patient mortality.

Assignment of Empirical Mode Decomposition Components and Its Application to Biomedical Signals

Objectives: Empirical mode decomposition (EMD) is a frequently used signal processing approach which adaptively decomposes a signal into a set of narrow-band components known as intrinsic mode functions (IMFs). For multi-trial, multivariate (multiple simultaneous recordings), and multi-subject analyses the number and signal properties of the IMFs can deviate from each other between trials, channels and subjects. A further processing of IMFs, e.g. a simple ensemble averaging, should determine which IMFs of one signal correspond to IMFs from another signal. When the signal properties have similar characteristics, the IMFs are assigned to each other. This problem is known as correspondence problem.

Methods: From the mathematical point of view, in some cases the correspondence problem can be transformed into an assignment problem which can be solved e.g. by the Kuhn-Munkres algorithm (KMA) by which a minimal cost matching can be found. We use the KMA for solving classic assignment problems, i.e. the pairwise correspondence between two sets of IMFs of equal cardinalities, and for pairwise correspondences between two sets of IMFs with different cardinalities representing an unbalanced assignment problem which is a special case of the k-cardinality assignment problem.

Results: A KMA-based approach to solve the correspondence problem was tested by using simulated, heart rate variability (HRV), and EEG data. The KMA-based results of HRV decomposition are compared with those obtained from a hierarchical cluster analysis (state-of-the-art). The major difference between the two approaches is that there is a more consistent assignment pattern using KMA. Integrating KMA into complex analysis concepts enables a comprehensive exploitation of the key advantages of the EMD. This can be demonstrated by non-linear analysis of HRV-related IMFs and by an EMD-based cross-frequency coupling analysis of the EEG data.

Conclusions: The successful application to HRV and EEG analysis demonstrates that our solutions can be used for automated EMD-based processing concepts for biomedical signals.

The neural bases of ictal tachycardia in temporal lobe seizures

OBJECTIVE

Due to limited information from scalp electroencephalographic (EEG) recordings, brain areas driving changes in cardiac rhythm during Temporal lobe (TL) seizures are not clearly identified. Using stereotactic EEG (SEEG) recordings, we aimed at identifying which of the brain regions involved in autonomic control trigger ictal tachycardia.

METHODS

The neural activity of several mesial temporal lobe structures including amygdala, hippocampus, insula, and lateral temporal lobe recorded with SEEG were collected during 37 TL seizures in 9 patients, using indices based on High Frequency Activity (HFA). R-R intervals (RR) monitoring and time-frequency spectral analysis were performed to assess parasympathetic (High frequency power (HF)) and sympathetic (Low frequency/High frequency (LF/HF) ratio) reactivities.

RESULTS

Tachycardia was associated with a significant increase in LF/HF ratio and decrease in HF. Autonomic cardiac changes were accompanied by simultaneous SEEG signal changes with an increase in seizure-related HFA in anterior hippocampal formation and amygdala, but not in insula.

CONCLUSION

In our sample, TL seizures are thus accompanied by an early decrease in parasympathetic control of cardiac rhythm and by an increase of sympathetic tone, concomitant to seizure activity in anterior hippocampus and amygdala.

SIGNIFICANCE

These results support a pivotal role of hippocampus and amygdala in tachycardia occurring during TL seizures.

Cerebral and neural regulation of cardiovascular activity during mental stress

Background

Mental arithmetic has been verified inducing cerebral and cardiovascular responses. However, the mechanism and sequential responses are still ambiguous. This study aims to reveal the mechanism of cardiovascular and autonomic responses and the related scalp positions that regulate the autonomic nerves system (ANS) during MA task.

Methods

34 healthy male subjects aged between 19 and 27 years old (mean age 23.6 ± 2.3 years) were recruited in. Electrocardiogram, impedance cardiography, beat-to-beat blood pressure and electroencephalography were measured simultaneously and continuously during the experiments. And the analysis of time–frequency, approximate entropy and Pearson correlation coefficient were adopted. For statistical comparison, paired t test is utilized in the study.

Results

The results showed that mental arithmetic task increased heart rate (from 72.35 ± 1.88 to 80.38 ± 2.34), blood pressure (systolic blood pressure: from 112.09 ± 3.23 to 126.79 ± 3.44; diastolic blood pressure: from 74.15 ± 1.93 to 81.20 ± 1.97), and cardiac output (from 8.71 ± 0.30 to 9.68 ± 0.35), and the mental arithmetic induced physiological responses could be divided into two stages, the first stage (10–110 s) and late stage (150–250 s). The high frequency power component (HF) of HRV decreased during MA, but the normalized low frequency power component (nLF) and LF/HF ratio of HRV increased only at the late stage. Moreover, during first stage, the correlations between approximate entropy of electroencephalography at Fp2, Fz, F4, F7 and the corresponding time–frequency results of HF were significant. During the late stage, the correlations between approximate entropy of electroencephalography at Fp2, Fz, C3, C4 and the corresponding nLF was significant.

Conclusions

Our results demonstrated that (1) mental stress induces time-dependent ANS activity and cardiovascular response. (2) Parasympathetic activity is lower during mental arithmetic task, but sympathetic nerve is activated only during late stage of mental arithmetic task. (3) Brain influences the cardiac activity through prefrontal and temporal cortex with the activation of ANS during mental arithmetic.

Comparative study of methods for solving the correspondence problem in EMD applications

We address the correspondence problem which arises when applying empirical mode decomposition (EMD) to multi-trial and multi-subject data. EMD decomposes a signal into a set of narrow-band components named intrinsic mode functions (IMFs). The number of IMFs and their signal properties can be different between trials, channels and subjects. In order to assign IMFs with similar characteristics to each other, we compare two assignment methods, unbalanced assignment and k-cardinality assignment and two clustering algorithms, namely hierarchical clustering and density-based spatial clustering of applications with noise based on heart rate variability data of children with temporal lobe epilepsy.

Discussion of "Computational Electrocardiography: Revisiting Holter ECG Monitoring"

This article is part of a For-Discussion-Section of Methods of Information in Medicine about the paper "Computational Electrocardiography: Revisiting Holter ECG Monitoring" written by Thomas M. Deserno and Nikolaus Marx. It is introduced by an editorial. This article contains the combined commentaries invited to independently comment on the paper of Deserno and Marx. In subsequent issues the discussion can continue through letters to the editor.

Epileptic seizure onset detection based on EEG and ECG data fusion

This paper presents a novel method for seizure onset detection using fused information extracted from multichannel electroencephalogram (EEG) and single-channel electrocardiogram (ECG). In existing seizure detectors, the analysis of the nonlinear and nonstationary ECG signal is limited to the time-domain or frequency-domain. In this work, heart rate variability (HRV) extracted from ECG is analyzed using a Matching-Pursuit (MP) and Wigner-Ville Distribution (WVD) algorithm in order to effectively extract meaningful HRV features representative of seizure and nonseizure states. The EEG analysis relies on a common spatial pattern (CSP) based feature enhancement stage that enables better discrimination between seizure and nonseizure features. The EEG-based detector uses logical operators to pool SVM seizure onset detections made independently across different EEG spectral bands. Two fusion systems are adopted. In the first system, EEG-based and ECG-based decisions are directly fused to obtain a final decision. The second fusion system adopts an override option that allows for the EEG-based decision to override the fusion-based decision in the event that the detector observes a string of EEG-based seizure decisions. The proposed detectors exhibit an improved performance, with respect to sensitivity and detection latency, compared with the state-of-the-art detectors. Experimental results demonstrate that the second detector achieves a sensitivity of 100%, detection latency of 2.6s, and a specificity of 99.91% for the MAJ fusion case.

Advantages of signal-adaptive approaches for the nonlinear, time-variant analysis of heart rate variability of children with temporal lobe epilepsy

Major aim of our study is to demonstrate that signal-adaptive approaches improve the nonlinear and time-variant analysis of heart rate variability (HRV) of children with temporal lobe epilepsy (TLE). Nonlinear HRV analyses are frequently applied in epileptic patients. As HRV is characterized by components with oscillatory properties frequency-selective methods are in the focus, whereby application of nonlinear analysis to linear filtered signals seems to be doubtful. Signal-adaptive methods that preserve nonlinear properties and utilize only the signal data for an automatic computation of the result could benefit to nonlinear analysis of HRV. Combinations of (1) the signal-adaptive Matched Gabor Transform with time-variant nonlinear bispectral analysis and of (2) signal-adaptive Empirical Mode Decomposition methods with time-variant nonlinear stability analysis are investigated with regard to their application in the analysis of specific HRV components (respiratory sinus arrhythmia and Mayer wave associated low-frequency HRV components) of 18 children with TLE. Changes of timing and coordination of both HRV components during preictal, ictal and postictal periods occur which can be better quantified by advanced signal-adaptive methods. Both approaches contribute with specific importance to the analysis.

Convergent Cross Mapping: Basic concept, influence of estimation parameters and practical application

In neuroscience, data are typically generated from neural network activity. Complex interactions between measured time series are involved, and nothing or only little is known about the underlying dynamic system. Convergent Cross Mapping (CCM) provides the possibility to investigate nonlinear causal interactions between time series by using nonlinear state space reconstruction. Aim of this study is to investigate the general applicability, and to show potentials and limitation of CCM. Influence of estimation parameters could be demonstrated by means of simulated data, whereas interval-based application of CCM on real data could be adapted for the investigation of interactions between heart rate and specific EEG components of children with temporal lobe epilepsy.