Despite the widespread diffusion of nonlinear methods for heart rate variability (HRV) analysis, the presence and the extent to which nonlinear dynamics contribute to short-term HRV are still ...controversial. This work aims at testing the hypothesis that different types of nonlinearity can be observed in HRV depending on the method adopted and on the physiopathological state. Two entropy-based measures of time series complexity (normalized complexity index, NCI) and regularity (information storage, IS), and a measure quantifying deviations from linear correlations in a time series (Gaussian linear contrast, GLC), are applied to short HRV recordings obtained in young (Y) and old (O) healthy subjects and in myocardial infarction (MI) patients monitored in the resting supine position and in the upright position reached through head-up tilt. The method of surrogate data is employed to detect the presence and quantify the contribution of nonlinear dynamics to HRV. We find that the three measures differ both in their variations across groups and conditions and in the percentage and strength of nonlinear HRV dynamics. NCI and IS displayed opposite variations, suggesting more complex dynamics in O and MI compared to Y and less complex dynamics during tilt. The strength of nonlinear dynamics is reduced by tilt using all measures in Y, while only GLC detects a significant strengthening of such dynamics in MI. A large percentage of detected nonlinear dynamics is revealed only by the IS measure in the Y group at rest, with a decrease in O and MI and during T, while NCI and GLC detect lower percentages in all groups and conditions. While these results suggest that distinct dynamic structures may lie beneath short-term HRV in different physiological states and pathological conditions, the strong dependence on the measure adopted and on their implementation suggests that physiological interpretations should be provided with caution.
This work presents a comparison between different approaches for the model-free estimation of information-theoretic measures of the dynamic coupling between short realizations of random processes. ...The measures considered are the mutual information rate (MIR) between two random processes X and Y and the terms of its decomposition evidencing either the individual entropy rates of X and Y and their joint entropy rate, or the transfer entropies from X to Y and from Y to X and the instantaneous information shared by X and Y. All measures are estimated through discretization of the random variables forming the processes, performed either via uniform quantization (binning approach) or rank ordering (permutation approach). The binning and permutation approaches are compared on simulations of two coupled non-identical Hènon systems and on three datasets, including short realizations of cardiorespiratory (CR, heart period and respiration flow), cardiovascular (CV, heart period and systolic arterial pressure), and cerebrovascular (CB, mean arterial pressure and cerebral blood flow velocity) measured in different physiological conditions, i.e., spontaneous vs paced breathing or supine vs upright positions. Our results show that, with careful selection of the estimation parameters (i.e., the embedding dimension and the number of quantization levels for the binning approach), meaningful patterns of the MIR and of its components can be achieved in the analyzed systems. On physiological time series, we found that paced breathing at slow breathing rates induces less complex and more coupled CR dynamics, while postural stress leads to unbalancing of CV interactions with prevalent baroreflex coupling and to less complex pressure dynamics with preserved CB interactions. These results are better highlighted by the permutation approach, thanks to its more parsimonious representation of the discretized dynamic patterns, which allows one to explore interactions with longer memory while limiting the curse of dimensionality.
In this work, we extend to the multivariate case the classical correlation analysis used in the field of network physiology to probe dynamic interactions between organ systems in the human body. To ...this end, we define different correlation-based measures of the multivariate interaction (MI) within and between the brain and body subnetworks of the human physiological network, represented, respectively, by the time series of δ, θ, α, and β electroencephalographic (EEG) wave amplitudes, and of heart rate, respiration amplitude, and pulse arrival time (PAT) variability (η, ρ, π). MI is computed: (i) considering all variables in the two subnetworks to evaluate overall brain-body interactions; (ii) focusing on a single target variable and dissecting its global interaction with all other variables into contributions arising from the same subnetwork and from the other subnetwork; and (iii) considering two variables conditioned to all the others to infer the network topology. The framework is applied to the time series measured from the EEG, electrocardiographic (ECG), respiration, and blood volume pulse (BVP) signals recorded synchronously via wearable sensors in a group of healthy subjects monitored at rest and during mental arithmetic and sustained attention tasks. We find that the human physiological network is highly connected, with predominance of the links internal of each subnetwork (mainly η-ρ and δ-θ, θ-α, α-β), but also statistically significant interactions between the two subnetworks (mainly η-β and η-δ). MI values are often spatially heterogeneous across the scalp and are modulated by the physiological state, as indicated by the decrease of cardiorespiratory interactions during sustained attention and by the increase of brain-heart interactions and of brain-brain interactions at the frontal scalp regions during mental arithmetic. These findings illustrate the complex and multi-faceted structure of interactions manifested within and between different physiological systems and subsystems across different levels of mental stress.
Heart Rate Variability (HRV) and Blood Pressure Variability (BPV) are widely employed tools for characterizing the complex behavior of cardiovascular dynamics. Usually, HRV and BPV analyses are ...carried out through short-term (ST) measurements, which exploit ~five-minute-long recordings. Recent research efforts are focused on reducing the time series length, assessing whether and to what extent Ultra-Short-Term (UST) analysis is capable of extracting information about cardiovascular variability from very short recordings. In this work, we compare ST and UST measures computed on electrocardiographic R-R intervals and systolic arterial pressure time series obtained at rest and during both postural and mental stress. Standard time-domain indices are computed, together with entropy-based measures able to assess the regularity and complexity of cardiovascular dynamics, on time series lasting down to 60 samples, employing either a faster linear parametric estimator or a more reliable but time-consuming model-free method based on nearest neighbor estimates. Our results are evidence that shorter time series down to 120 samples still exhibit an acceptable agreement with the ST reference and can also be exploited to discriminate between stress and rest. Moreover, despite neglecting nonlinearities inherent to short-term cardiovascular dynamics, the faster linear estimator is still capable of detecting differences among the conditions, thus resulting in its suitability to be implemented on wearable devices.
Apnea and other breathing-related disorders have been linked to the development of hypertension or impairments of the cardiovascular, cognitive or metabolic systems. The combined assessment of ...multiple physiological signals acquired during sleep is of fundamental importance for providing additional insights about breathing disorder events and the associated impairments. In this work, we apply information-theoretic measures to describe the joint dynamics of cardiorespiratory physiological processes in a large group of patients reporting repeated episodes of hypopneas, apneas (central, obstructive, mixed) and respiratory effort related arousals (RERAs). We analyze the heart period as the target process and the airflow amplitude as the driver, computing the predictive information, the information storage, the information transfer, the internal information and the cross information, using a fuzzy kernel entropy estimator. The analyses were performed comparing the information measures among segments during, immediately before and after the respiratory event and with control segments. Results highlight a general tendency to decrease of predictive information and information storage of heart period, as well as of cross information and information transfer from respiration to heart period, during the breathing disordered events. The information-theoretic measures also vary according to the breathing disorder, and significant changes of information transfer can be detected during RERAs, suggesting that the latter could represent a risk factor for developing cardiovascular diseases. These findings reflect the impact of different sleep breathing disorders on respiratory sinus arrhythmia, suggesting overall higher complexity of the cardiac dynamics and weaker cardiorespiratory interactions which may have physiological and clinical relevance.
Heart rate variability (HRV; variability of the RR interval of the electrocardiogram) results from the activity of several coexisting control mechanisms, which involve the influence of respiration ...(RESP) and systolic blood pressure (SBP) oscillations operating across multiple temporal scales and changing in different physiological states. In this study, multiscale information decomposition is used to dissect the physiological mechanisms related to the genesis of HRV in 78 young volunteers monitored at rest and during postural and mental stress evoked by head-up tilt (HUT) and mental arithmetics (MA). After representing RR, RESP and SBP at different time scales through a recently proposed method based on multivariate state space models, the joint information transfer T RESP , SBP → RR is decomposed into unique, redundant and synergistic components, describing the strength of baroreflex modulation independent of respiration ( U SBP → RR ), nonbaroreflex ( U RESP → RR ) and baroreflex-mediated ( R RESP , SBP → RR ) respiratory influences, and simultaneous presence of baroreflex and nonbaroreflex respiratory influences ( S RESP , SBP → RR ), respectively. We find that fast (short time scale) HRV oscillations—respiratory sinus arrhythmia—originate from the coexistence of baroreflex and nonbaroreflex (central) mechanisms at rest, with a stronger baroreflex involvement during HUT. Focusing on slower HRV oscillations, the baroreflex origin is dominant and MA leads to its higher involvement. Respiration influences independent on baroreflex are present at long time scales, and are enhanced during HUT.
In this study, an analysis of brain, cardiovascular and respiratory dynamics was conducted combining information-theoretic measures with the Network Physiology paradigm during different levels of ...mental stress. Starting from low invasive recordings of electroencephalographic, electrocardiographic, respiratory, and blood volume pulse signals, the dynamical activity of seven physiological systems was probed with one-second time resolution measuring the time series of the δ , θ , α and β brain wave amplitudes, the cardiac period (RR interval), the respiratory amplitude, and the duration of blood pressure wave propagation (pulse arrival time, PAT). Synchronous 5-min windows of these time series, obtained from 18 subjects during resting wakefulness (REST), mental stress induced by mental arithmetic (MA) and sustained attention induced by serious game (SG), were taken to describe the dynamics of the nodes composing the observed physiological network. Network activity and connectivity were then assessed in the framework of information dynamics computing the new information generated by each node, the information dynamically stored in it, and the information transferred to it from the other network nodes. Moreover, the network topology was investigated using directed measures of conditional information transfer and assessing their statistical significance. We found that all network nodes dynamically produce and store significant amounts of information, with the new information being prevalent in the brain systems and the information storage being prevalent in the peripheral systems. The transition from REST to MA was associated with an increase of the new information produced by the respiratory signal time series (RESP), and that from MA to SG with a decrease of the new information produced by PAT. Each network node received a significant amount of information from the other nodes, with the highest amount transferred to RR and the lowest transferred to δ , θ , α and β . The topology of the physiological network underlying such information transfer was node- and state-dependent, with the peripheral subnetwork showing interactions from RR to PAT and between RESP and RR, PAT consistently across states, the brain subnetwork resulting more connected during MA, and the subnetwork of brain–peripheral interactions involving different brain rhythms in the three states and resulting primarily activated during MA. These results have both physiological relevance as regards the interpretation of central and autonomic effects on cardiovascular and respiratory variability, and practical relevance as regards the identification of features useful for the automatic distinction of different mental states.
Network Physiology is a rapidly growing field of study that aims to understand how physiological systems interact to maintain health. Within the information theory framework the information storage ...(IS) allows to measure the regularity and predictability of a dynamic process under stationarity assumption. However, this assumption does not allow to track over time the transient pathways occurring in the dynamical activity of a physiological system. To address this limitation, we propose a time-varying approach based on the recursive least squares algorithm (RLS) for estimating IS at each time instant, in non-stationary conditions. We tested this approach in simulated time-varying dynamics and in the analysis of electroencephalographic (EEG) signals recorded from healthy volunteers and timed with the heartbeat to investigate brain-heart interactions. In simulations, we show that the proposed approach allows to track both abrupt and slow changes in the information stored in a physiological system. These changes are reflected in its evolution and variability over time. The analysis of brain-heart interactions reveals marked differences across the cardiac cycle phases of the variability of the time-varying IS. On the other hand, the average IS values exhibit a weak modulation over parieto-occiptal areas of the scalp. Our study highlights the importance of developing more advanced methods for measuring IS that account for non-stationarity in physiological systems. The proposed time-varying approach based on RLS represents a useful tool for identifying spatio-temporal dynamics within the neurocardiac system and can contribute to the understanding of brain-heart interactions.
This study introduces a framework for the information-theoretic analysis of brain functional connectivity performed at the level of electroencephalogram (EEG) sources. The framework combines the use ...of common spatial patterns to select the EEG components which maximize the variance between two experimental conditions, simultaneous implementation of vector autoregressive modeling (VAR) with independent component analysis to describe the joint source dynamics and their projection to the scalp, and computation of information dynamics measures (information storage, information transfer, statistically significant network links) from the source VAR parameters. The proposed framework was tested on simulated EEGs obtained mixing source signals generated under different coupling conditions, showing its ability to retrieve source information dynamics from the scalp signals. Then, it was applied to investigate scalp and source brain connectivity in a group of children manifesting episodes of focal and generalized epilepsy; the analysis was performed on EEG signals lasting 5 s, collected in two consecutive windows preceding and one window following each ictal episode. Our results show that generalized seizures are associated with a significant decrease from pre-ictal to post-ictal periods of the information stored in the signals and of the information transferred among them, reflecting reduced self-predictability and causal connectivity at the level of both scalp and source brain dynamics. On the contrary, in the case of focal seizures the scalp EEG activity was not discriminated across conditions by any information measure, while source analysis revealed a tendency of the measures of information transfer to increase just before seizures and to decrease just after seizures. These results suggest that focal epileptic seizures are associated with a reorganization of the topology of EEG brain networks which is only visible analyzing connectivity among the brain sources. Our findings emphasize the importance of EEG modeling approaches able to deal with the adverse effects of volume conduction on brain connectivity analysis, and their potential relevance to the development of strategies for prediction and clinical treatment of epilepsy.
Wearable health devices (WHDs) are rapidly gaining ground in the biomedical field due to their ability to monitor the individual physiological state in everyday life scenarios, while providing a ...comfortable wear experience. This study introduces a novel wearable biomedical device capable of synchronously acquiring electrocardiographic (ECG), photoplethysmographic (PPG), galvanic skin response (GSR) and motion signals. The device has been specifically designed to be worn on a finger, enabling the acquisition of all biosignals directly on the fingertips, offering the significant advantage of being very comfortable and easy to be employed by the users. The simultaneous acquisition of different biosignals allows the extraction of important physiological indices, such as heart rate (HR) and its variability (HRV), pulse arrival time (PAT), GSR level, blood oxygenation level (SpO
), and respiratory rate, as well as motion detection, enabling the assessment of physiological states, together with the detection of potential physical and mental stress conditions. Preliminary measurements have been conducted on healthy subjects using a measurement protocol consisting of resting states (i.e., SUPINE and SIT) alternated with physiological stress conditions (i.e., STAND and WALK). Statistical analyses have been carried out among the distributions of the physiological indices extracted in time, frequency, and information domains, evaluated under different physiological conditions. The results of our analyses demonstrate the capability of the device to detect changes between rest and stress conditions, thereby encouraging its use for assessing individuals' physiological state. Furthermore, the possibility of performing synchronous acquisitions of PPG and ECG signals has allowed us to compare HRV and pulse rate variability (PRV) indices, so as to corroborate the reliability of PRV analysis under stationary physical conditions. Finally, the study confirms the already known limitations of wearable devices during physical activities, suggesting the use of algorithms for motion artifact correction.