Search

In the study of complex biomedical systems represented by multivariate stochastic processes, such as the cardiovascular and respiratory systems, an issue of great relevance is the description of the system dynamics spanning multiple temporal scales. Recently, the quantification of multiscale complexity based on linear parametric models, incorporating autoregressive coefficients and fractional integration, encompassing short term dynamics and long-range correlations, was extended to multivariate time series. Within this Vector AutoRegressive Fractionally Integrated (VARFI) framework formalized for Gaussian processes, in this work we propose to estimate the Transfer Entropy, or equivalently Granger Causality, in the cardiovascular and respiratory systems. This allows to quantify the information flow and assess directed interactions accounting for the simultaneous presence of short-term dynamics and long-range correlations. The proposed approach is first tested on simulations of benchmark VARFI processes where the transfer entropy could be computed from the known model parameters. Then, it is applied to experimental data consisting of heart period, systolic arterial pressure and respiration time series measured in healthy subjects monitored at rest and during mental and postural stress. Both simulations and real data analysis revealed that the proposed method highlights the dependence of the information transfer on the balance between short-term and longrange correlations in coupled dynamical systems.

The increasing interest in innovative solutions for health and physiological monitoring has recently fostered the development of smaller biomedical devices. These devices are capable of recording an increasingly large number of biosignals simultaneously, while maximizing the user’s comfort. In this study, we have designed and realized a novel wearable multisensor ring-shaped probe that enables synchronous, real-time acquisition of photoplethysmographic (PPG) and galvanic skin response (GSR) signals. The device integrates both the PPG and GSR sensors onto a single probe that can be easily placed on the finger, thereby minimizing the device footprint and overall size. The system enables the extraction of various physiological indices, including heart rate (HR) and its variability, oxygen saturation (SpO2), and GSR levels, as well as their dynamic changes over time, to facilitate the detection of different physiological states, e.g., rest and stress. After a preliminary SpO2 calibration procedure, measurements have been carried out in laboratory on healthy subjects to demonstrate the feasibility of using our system to detect rapid changes in HR, skin conductance, and SpO2 across various physiological conditions (i.e., rest, sudden stress-like situation and breath holding). The early findings encourage the use of the device in daily-life conditions for real-time monitoring of different physiological states.

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 [Formula: see text] and [Formula: see text] and the terms of its decomposition evidencing either the individual entropy rates of [Formula: see text] and [Formula: see text] and their joint entropy rate, or the transfer entropies from [Formula: see text] to [Formula: see text] and from [Formula: see text] to [Formula: see text] and the instantaneous information shared by [Formula: see text] and [Formula: see text]. 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.

Arterial compliance is an important parameter influencing ventricular-arterial coupling, depending on structural and functional mechanics of arteries. In this study, the spontaneous beat-to-beat variability of arterial compliance was investigated in time and frequency domains in thirty-nine young and healthy subjects monitored in the supine resting state and during head-up tilt. Spectral decomposition was applied to retrieve the spectral content of the time series associated to low (LF) and high frequency (HF) oscillatory components. Our results highlight: (i) a decrease of arterial compliance with tilt, in agreement with previous studies; (ii) an increase of the LF power content concurrent with a decrease of the HF power, potentially reflecting changes in vasomotor tone, blood pressure and heart rate variability associated with higher sympathetic activity and vagal withdrawal occurring with tilt.

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.

Many complex systems in physics, biology and engineering are modeled as dynamical networks and described using multivariate time series analysis. Recent developments have shown that the emergent dynamics of a network system are significantly affected by interactions involving multiple network nodes which cannot be described using pairwise links. While these higher-order interactions can be probed using information-theoretic measures, a rigorous framework to describe them in the frequency domain is still lacking. This work presents an approach for the spectral decomposition of multivariate information measures, capable of identifying higher-order synergistic and redundant interactions between oscillatory processes. We show theoretically that synergy and redundancy can coexist at different frequencies among the output signals of a network system and can be detected only using the proposed spectral method. To demonstrate the broad applicability of the framework, we provide parametric and non-parametric data-efficient estimators for the spectral information measures, and employ them to describe multivariate interactions in three complex systems producing rich oscillatory dynamics, namely the human brain, a ring of electronic oscillators, and the global climate system. In these systems, we show that the use of our framework for the spectral decomposition of information measures reveals multivariate and higher-order interactions not detectable in the time domain. Our results are exemplary of how the frequency-specific analysis of multivariate dynamics can aid the implementation of assessment and control strategies in real-world network systems.

Objective. In this work, an analytical framework for the multiscale analysis of multivariate Gaussian processes is presented, whereby the computation of Partial Information Decomposition measures is achieved accounting for the simultaneous presence of short-term dynamics and long-range correlations. Approach. We consider physiological time series mapping the activity of the cardiac, vascular and respiratory systems in the field of Network Physiology. In this context, the multiscale representation of transfer entropy within the network of interactions among Systolic arterial pressure (S), respiration (R) and heart period (H), as well as the decomposition into unique, redundant and synergistic contributions, is obtained using a Vector AutoRegressive Fractionally Integrated (VARFI) framework for Gaussian processes. This novel approach allows to quantify the directed information flow accounting for the simultaneous presence of short-term dynamics and long-range correlations among the analyzed processes. Additionally, it provides analytical expressions for the computation of the information measures, by exploiting the theory of state space models. The approach is first illustrated in simulated VARFI processes and then applied to H, S and R time series measured in healthy subjects monitored at rest and during mental and postural stress. Main Results. We demonstrate the ability of the VARFI modeling approach to account for the coexistence of short-term and long-range correlations in the study of multivariate processes. Physiologically, we show that postural stress induces larger redundant and synergistic effects from S and R to H at short time scales, while mental stress induces larger information transfer from S to H at longer time scales, thus evidencing the different nature of the two stressors. Significance. The proposed methodology allows to extract useful information about the dependence of the information transfer on the balance between short-term and long-range correlations in coupled dynamical systems, which cannot be observed using standard methods that do not consider long-range correlations.

The amount of information exchanged per unit of time between two dynamic processes is an important concept for the analysis of complex systems. Theoretical formulations and data-efficient estimators have been recently introduced for this quantity, known as the mutual information rate (MIR), allowing its continuous-time computation for event-based data sets measured as realizations of coupled point processes. This work presents the implementation of MIR for point process applications in Network Physiology and cardiovascular variability, which typically feature short and noisy experimental time series. We assess the bias of MIR estimated for uncoupled point processes in the frame of surrogate data, and we compensate it by introducing a corrected MIR (cMIR) measure designed to return zero values when the two processes do not exchange information. The method is first tested extensively in synthetic point processes including a physiologically-based model of the heartbeat dynamics and the blood pressure propagation times, where we show the ability of cMIR to compensate the negative bias of MIR and return statistically significant values even for weakly coupled processes. The method is then assessed in real point-process data measured from healthy subjects during different physiological conditions, showing that cMIR between heartbeat and pressure propagation times increases significantly during postural stress, though not during mental stress. These results document that cMIR reflects physiological mechanisms of cardiovascular variability related to the joint neural autonomic modulation of heart rate and arterial compliance.

While the standard network description of complex systems is based on quantifying the link between pairs of system units, higher-order interactions (HOIs) involving three or more units often play a major role in governing the collective network behavior. This work introduces a new approach to quantify pairwise and HOIs for multivariate rhythmic processes interacting across multiple time scales. We define the so-called O-information rate (OIR) as a new metric to assess HOIs for multivariate time series, and present a framework to decompose the OIR into measures quantifying Granger-causal and instantaneous influences, as well as to expand all measures in the frequency domain. The framework exploits the spectral representation of vector autoregressive and state space models to assess the synergistic and redundant interaction among groups of processes, both in specific bands of interest and in the time domain after whole-band integration. Validation of the framework on simulated networks illustrates how the spectral OIR can highlight redundant and synergistic HOIs emerging at specific frequencies, which cannot be detected using time-domain measures. The applications to physiological networks described by heart period, arterial pressure and respiration variability measured in healthy subjects during a protocol of paced breathing, and to brain networks described by electrocorticographic signals acquired in an animal experiment during anesthesia, document the capability of our approach to identify informational circuits relevant to well-defined cardiovascular oscillations and brain rhythms and related to specific physiological mechanisms involving autonomic control and altered consciousness. The proposed framework allows a hierarchically-organized evaluation of timeand frequency-domain interactions in dynamic networks mapped by multivariate time series, and its high flexibility and scalability make it suitable for the investigation of networks beyond pairwise interactions in neuroscience, physiology and many other fields.

Heart Period (H) results from the activity of several coexisting control mechanisms, involving Systolic Arterial Pressure (S) and Respiration (R), which operate across multiple time scales encompassing not only short-term dynamics but also long-range correlations. In this work, multiscale representation of Transfer Entropy (TE) and of its decomposition in the network of these three interacting processes is obtained by extending the multivariate approach based on linear parametric VAR models to the Vector AutoRegressive Fractionally Integrated (VARFI) framework for Gaussian processes. This approach allows to dissect the different contributions to cardiac dynamics accounting for the simultaneous presence of short and long term dynamics. The proposed method is first tested on simulations of a benchmark VARFI model and then applied to experimental data consisting of H, S and R time series measured in healthy subjects monitored at rest and during mental and postural stress. The results reveal that the proposed method can highlight the dependence of the information transfer on the balance between short-term and long-range correlations in coupled dynamical systems.

Background and objective: Recently, various algorithms have been introduced using wrist-worn photo-plethysmography (PPG) to provide high accuracy of instantaneous heart rate (HR) estimation, including during high-intensity exercise. Most studies focus on using acceleration and/or gyroscope signals for the motion artifact (MA) reference, which attenuates or cancels out noise from the MA-corrupted PPG signals. We aim to open and pave the path to find an appropriate MA reference selection for MA cancelation in PPG.Methods: We investigated how the acceleration and gyroscope reference signals correlate with the MAs of the distorted PPG signals and derived both mathematically and experimentally an adaptive MA reference selection approach. We applied our algorithm to five state-of-the-art (SOTA) methods for the performance evaluation. In addition, we compared the four MA reference selection approaches, i.e. with acceleration signal only, with gyroscope signal only, with both signals, and using our proposed adaptive selection.Results: When applied to 47 PPG recordings acquired during intensive physical exercise from two differ-ent datasets, our proposed adaptive MA reference selection method provided higher accuracy than the other MA selection approaches for all five SOTA methods.Conclusion: Our proposed adaptive MA reference selection approach can be used in other MA cancelation methods and reduces the HR estimation error.Significance: We believe that this study helps researchers to address acceleration and gyroscope signals as accurate MA references, which eventually improves the overall performance for estimating HRs through the various algorithms developed by research groups.

While cross-spectral and information-theoretic approaches are widely used for the multivariate analysis of physiological time series, their combined utilization is far less developed in the literature. This study introduces a framework for the spectral decomposition of multivariate information measures, which provides frequency-specific quantifications of the information shared between a target and two source time series and of its expansion into amounts related to how the sources contribute to the target dynamics with unique, redundant and synergistic information. The framework is illustrated in simulations of linearly interacting stochastic processes, showing how it allows us to retrieve amounts of information shared by the processes within specific frequency bands which are otherwise not detectable by time-domain information measures, as well as coupling features which are not detectable by spectral measures. Then, it is applied to the time series of heart period, systolic and diastolic arterial pressure and respiration variability measured in healthy subjects monitored in the resting supine position and during head-up tilt. We show that the spectral measures of unique, redundant and synergistic information shared by these variability series, integrated within specific frequency bands of physiological interest and reflect the mechanisms of short-term regulation of cardiovascular and cardiorespiratory oscillations and their alterations induced by the postural stress. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

We present a framework for the linear parametric analysis of pairwise interactions in bivariate time series in the time and frequency domains, which allows the evaluation of total, causal and instantaneous interactions and connects time- and frequency-domain measures. The framework is applied to physiological time series to investigate the cerebrovascular regulation from the variability of mean cerebral blood flow velocity (CBFV) and mean arterial pressure (MAP), and the cardiovascular regulation from the variability of heart period (HP) and systolic arterial pressure (SAP). We analyze time series acquired at rest and during the early and late phase of head-up tilt in subjects developing orthostatic syncope in response to prolonged postural stress, and in healthy controls. The spectral measures of total, causal and instantaneous coupling between HP and SAP, and between MAP and CBFV, are averaged in the low-frequency band of the spectrum to focus on specific rhythms, and over all frequencies to get time-domain measures. The analysis of cardiovascular interactions indicates that postural stress induces baroreflex involvement, and its prolongation induces baroreflex dysregulation in syncope subjects. The analysis of cerebrovascular interactions indicates that the postural stress enhances the total coupling between MAP and CBFV, and challenges cerebral autoregulation in syncope subjects, while the strong sympathetic activation elicited by prolonged postural stress in healthy controls may determine an increased coupling from CBFV to MAP during late tilt. These results document that the combination of time-domain and spectral measures allows us to obtain an integrated view of cardiovascular and cerebrovascular regulation in healthy and diseased subjects.

One of the most challenging problems in the study of complex dynamical systems is to find the statistical interdependencies among the system components. Granger causality (GC) represents one of the most employed approaches, based on modeling the system dynamics with a linear vector autoregressive (VAR) model and on evaluating the information flow between two processes in terms of prediction error variances. In its most advanced setting, GC analysis is performed through a state-space (SS) representation of the VAR model that allows to compute both conditional and unconditional forms of GC by solving only one regression problem. While this problem is typically solved through Ordinary Least Square (OLS) estimation, a viable alternative is to use Artificial Neural Networks (ANNs) implemented in a simple structure with one input and one output layer and trained in a way such that the weights matrix corresponds to the matrix of VAR parameters. In this work, we introduce an ANN combined with SS models for the computation of GC. The ANN is trained through the Stochastic Gradient Descent L1 (SGD-L1) algorithm, and a cumulative penalty inspired from penalized regression is applied to the network weights to encourage sparsity. Simulating networks of coupled Gaussian systems, we show how the combination of ANNs and SGD-L1 allows to mitigate the strong reduction in accuracy of OLS identification in settings of low ratio between number of time series points and of VAR parameters. We also report how the performances in GC estimation are influenced by the number of iterations of gradient descent and by the learning rate used for training the ANN. We recommend using some specific combinations for these parameters to optimize the performance of GC estimation. Then, the performances of ANN and OLS are compared in terms of GC magnitude and statistical significance to highlight the potential of the new approach to reconstruct causal coupling strength and network topology even in challenging conditions of data paucity. The results highlight the importance of of a proper selection of regularization parameter which determines the degree of sparsity in the estimated network. Furthermore, we apply the two approaches to real data scenarios, to study the physiological network of brain and peripheral interactions in humans under different conditions of rest and mental stress, and the effects of the newly emerged concept of remote synchronization on the information exchanged in a ring of electronic oscillators. The results highlight how ANNs provide a mesoscopic description of the information exchanged in networks of multiple interacting physiological systems, preserving the most active causal interactions between cardiovascular, respiratory and brain systems. Moreover, ANNs can reconstruct the flow of directed information in a ring of oscillators whose statistical properties can be related to those of physiological networks.

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 delta, theta, alpha, and beta 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 heart rate-respiration and delta-theta, theta-alpha, alpha-beta), but also statistically significant interactions between the two subnetworks (mainly heart rate-beta and heart rate-delta). 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.

The objective of this study is to present a new method to assess in the frequency domain the directed interactions between the spontaneous variability of systolic arterial pressure (SAP) and heart period (HP) from their linear model representation, and to apply it for studying the baroreflex control of arterial pressure in healthy physiological states and after acute myocardial infarction (AMI). The method is based on pole decomposition of the model transfer function and on the following evaluation of causal measures of coupling and gain from the poles associated to low frequency (0.04−0.15 Hz) oscillatory components. It is compared with traditional non-causal approaches for the spectral analysis of the baroreflex gain, and with causal approaches based on the directed coherence, in a group of AMI patients and in Young and Old healthy controls studied at rest and during head-up tilt. Analysis of feedforward interactions from RR to SAP is also performed. Our results support the importance of using local causal approaches to quantify separately baroreflex and feedforward interactions between RR and SAP, allowing both to confirm known physiological behaviors (e.g., weaker baroreflex effectiveness and lower sensitivity after AMI) and to suggest novel findings (e.g., preserved low frequency baroreflex response to postural stress after AMI). We also find that the postural stress alters feedforward interactions selectively across groups, being related to decreased coupling only in Young and to increased gain mostly in AMI. These results have a significance for the clinical assessment of the baroreflex and the physiological evaluation of cardiovascular interactions.