Multiscale, multiorgan and multivariate complexity analyses of cardiovascular regulation

Advances in Multivariate and Multiscale Physiological Signal Analysis

Physiological systems are characterized by complex dynamics and nonlinear behaviors due to their intricate structural organization and regulatory mechanisms [...].

Stratified Multivariate Multiscale Dispersion Entropy for Physiological Signal Analysis

Multivariate entropy quantification algorithms are becoming a prominent tool for the extraction of information from multi-channel physiological time-series. However, in the analysis of physiological signals from heterogeneous organ systems, certain channels may overshadow the patterns of others, resulting in information loss. Here, we introduce the framework of Stratified Entropy to prioritize each channels' dynamics based on their allocation to respective strata, leading to a richer description of the multi-channel time-series. As an implementation of the framework, three algorithmic variations of the Stratified Multivariate Multiscale Dispersion Entropy are introduced. These variations and the original algorithm are applied to synthetic time-series, waveform physiological time-series, and derivative physiological data. Based on the synthetic time-series experiments, the variations successfully prioritize channels following their strata allocation while maintaining the low computation time of the original algorithm. In experiments on waveform physiological time-series and derivative physiological data, increased discrimination capacity was noted for multiple strata allocations in the variations when benchmarked to the original algorithm. This suggests improved physiological state monitoring by the variations. Furthermore, our variations can be modified to utilize a priori knowledge for the stratification of channels. Thus, our research provides a novel approach for the extraction of previously inaccessible information from multi-channel time series acquired from heterogeneous systems.

Multiscale partial information decomposition of dynamic processes with short and long-range correlations: theory and application to cardiovascular control

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.

Towards a Contactless Stress Classification Using Thermal Imaging

Thermal cameras capture the infrared radiation emitted from a body in a contactless manner and can provide an indirect estimation of the autonomic nervous system (ANS) dynamics through the regulation of the skin temperature. This study investigates the contribution given by thermal imaging for an effective automatic stress detection with the perspective of a contactless stress recognition system. To this aim, we recorded both ANS correlates (cardiac, electrodermal, and respiratory activity) and thermal images from 25 volunteers under acute stress induced by the Stroop test. We conducted a statistical analysis on the features extracted from each signal, and we implemented subject-independent classifications based on the support vector machine model with an embedded recursive feature elimination algorithm. Particularly, we trained three classifiers using different feature sets: the full set of features, only those derived from the peripheral autonomic correlates, and only those derived from the thermal images. Classification accuracy and feature selection results confirmed the relevant contribution provided by the thermal features in the acute stress detection task. Indeed, a combination of ANS correlates and thermal features achieved 97.37% of accuracy. Moreover, using only thermal features we could still successfully detect stress with an accuracy of 86.84% in a contact-free manner.

Dysregulated Neurovascular Control Underlies Declining Microvascular Functionality in People With Non-alcoholic Fatty Liver Disease (NAFLD) at Risk of Liver Fibrosis

Background/Aims

Increasing evidence shows that non-alcoholic fatty liver disease (NAFLD) is associated with dysregulation of microvascular perfusion independently of established cardio-metabolic risk factors. We investigated whether hepatic manifestations of NAFLD such as liver fibrosis and liver fat are associated with microvascular hemodynamics through dysregulation of neurovascular control.

Methods

Microvascular dilator (post-occlusive reactive hyperemia) and sympathetically mediated constrictor (deep inspiratory breath-hold) responses were measured at the forearm and finger, respectively, using laser Doppler fluximetry. Non-linear complexity-based analysis was used to assess the information content and variability of the resting blood flux (BF) signals, attributable to oscillatory flow-motion activity, and over multiple sampling frequencies.

Results

Measurements were made in 189 adults (113 men) with NAFLD, with (n = 65) and without (n = 124) type 2 diabetes mellitus (T2DM), age = 50.9 ± 11.7 years (mean ± SD). Microvascular dilator and constrictor capacity were both negatively associated with age (r = −0.178, p = 0.014, and r = −0.201, p = 0.007, respectively) and enhanced liver fibrosis (ELF) score (r = −0.155, p = 0.038 and r = −0.418, p < 0.0001, respectively). There was no association with measures of liver fat, obesity or T2DM. Lempel-Ziv complexity (LZC) and sample entropy (SE) of the BF signal measured at the two skin sites were associated negatively with age (p < 0.01 and p < 0.001) and positively with ELF score (p < 0.05 and p < 0.0001). In individuals with an ELF score ≥7.8 the influence of both neurogenic and respiratory flow-motion activity on LZC was up-rated (p < 0.0001).

Conclusion

Altered microvascular network functionality occurs in adults with NAFLD suggesting a mechanistic role for dysregulated neurovascular control in individuals at risk of severe liver fibrosis.

Multiscale, multidomain analysis of microvascular flow dynamics

NEW FINDINGS

What is the topic of this review? We describe a range of techniques in the time, frequency and information domains and their application alone and together for the analysis of blood flux signals acquired using laser Doppler fluximetry. What advances does it highlight? This review highlights the idea of using quantitative measures in different domains and scales to gain a better mechanistic understanding of the complex behaviours in the microcirculation.

ABSTRACT

To date, time- and frequency-domain metrics of signals acquired through laser Doppler fluximetry have been unable to provide consistent and robust measures of the changes that occur in the microcirculation in healthy individuals at rest or in response to a provocation, or in patient cohorts. Recent studies have shown that in many disease states, such as metabolic and cardiovascular disease, there appears to be a reduction in the adaptive capabilities of the microvascular network and a consequent reduction in physiological information content. Here, we introduce non-linear measures for assessing the information content of fluximetry signals and demonstrate how they can yield deeper understanding of network behaviour. In addition, we show how these methods may be adapted to accommodate the multiple time scales modulating blood flow and how they can be used in combination with time- and frequency-domain metrics to discriminate more effectively between the different mechanistic influences on network properties.

Enhanced flow-motion complexity of skin microvascular perfusion in Sherpas and lowlanders during ascent to high altitude

An increased and more effective microvascular perfusion is postulated to play a key role in the physiological adaptation of Sherpa highlanders to the hypobaric hypoxia encountered at high altitude. To investigate this, we used Lempel-Ziv complexity (LZC) analysis to explore the spatiotemporal dynamics of the variability of the skin microvascular blood flux (BF) signals measured at the forearm and finger, in 32 lowlanders (LL) and 46 Sherpa highlanders (SH) during the Xtreme Everest 2 expedition. Measurements were made at baseline (BL) (LL: London 35 m; SH: Kathmandu 1300 m) and at Everest base camp (LL and SH: EBC 5,300 m). We found that BF signal content increased with ascent to EBC in both SH and LL. At both altitudes, LZC of the BF signals was significantly higher in SH, and was related to local slow-wave flow-motion activity over multiple spatial and temporal scales. In SH, BF LZC was also positively associated with LZC of the simultaneously measured tissue oxygenation signals. These data provide robust mechanistic information of microvascular network functionality and flexibility during hypoxic exposure on ascent to high altitude. They demonstrate the importance of a sustained heterogeneity of network perfusion, associated with local vaso-control mechanisms, to effective tissue oxygenation during hypobaric hypoxia.

Multivariate Multiscale Dispersion Entropy of Biomedical Times Series

Due to the non-linearity of numerous physiological recordings, non-linear analysis of multi-channel signals has been extensively used in biomedical engineering and neuroscience. Multivariate multiscale sample entropy (MSE–mvMSE) is a popular non-linear metric to quantify the irregularity of multi-channel time series. However, mvMSE has two main drawbacks: (1) the entropy values obtained by the original algorithm of mvMSE are either undefined or unreliable for short signals (300 sample points); and (2) the computation of mvMSE for signals with a large number of channels requires the storage of a huge number of elements. To deal with these problems and improve the stability of mvMSE, we introduce multivariate multiscale dispersion entropy (MDE–mvMDE), as an extension of our recently developed MDE, to quantify the complexity of multivariate time series. We assess mvMDE, in comparison with the state-of-the-art and most widespread multivariate approaches, namely, mvMSE and multivariate multiscale fuzzy entropy (mvMFE), on multi-channel noise signals, bivariate autoregressive processes, and three biomedical datasets. The results show that mvMDE takes into account dependencies in patterns across both the time and spatial domains. The mvMDE, mvMSE, and mvMFE methods are consistent in that they lead to similar conclusions about the underlying physiological conditions. However, the proposed mvMDE discriminates various physiological states of the biomedical recordings better than mvMSE and mvMFE. In addition, for both the short and long time series, the mvMDE-based results are noticeably more stable than the mvMSE- and mvMFE-based ones. For short multivariate time series, mvMDE, unlike mvMSE, does not result in undefined values. Furthermore, mvMDE is faster than mvMFE and mvMSE and also needs to store a considerably smaller number of elements. Due to its ability to detect different kinds of dynamics of multivariate signals, mvMDE has great potential to analyse various signals.

Vagal contributions to fetal heart rate variability: an omics approach

OBJECTIVE

Fetal heart rate variability (fHRV) is an important indicator of health and disease, yet its physiological origins, neural contributions, in particular, are not well understood. We aimed to develop novel experimental and data analytical approaches to identify fHRV measures reflecting the vagus nerve contributions to fHRV.

APPROACH

In near-term ovine fetuses, a comprehensive set of 46 fHRV measures was computed from fetal pre-cordial electrocardiogram recorded during surgery and 72 h later without (n  =  24) and with intra-surgical bilateral cervical vagotomy (n  =  15).

MAIN RESULTS

The fetal heart rate did not change due to vagotomy. We identify fHRV measures specific to the vagal modulation of fHRV: multiscale time irreversibility asymmetry index (AsymI), detrended fluctuation analysis (DFA) α ₁, Kullback-Leibler permutation entropy (KLPE) and scale-dependent Lyapunov exponent slope (SDLE α).

SIGNIFICANCE

We provide a systematic delineation of vagal contributions to fHRV across signal-analytical domains which should be relevant for the emerging field of bioelectronic medicine and the deciphering of the 'vagus code'. Our findings also have clinical significance for in utero monitoring of fetal health during surgery. Key points •Fetal surgery causes a complex pattern of changes in heart rate variability measures with an overall reduction of complexity or variability. •At 72 h after surgery, many of the HRV measures recover and this recovery is delayed by an intrasurgical cervical bilateral vagotomy. •We identify HRV pattern representing complete vagal withdrawal that can be understood as part of 'HRV code', rather than any single HRV measure. •We identify HRV biomarkers of recovery from fetal surgery and discuss the effect of anticholinergic medication on this recovery.

Towards understanding the complexity of cardiovascular oscillations: Insights from information theory

Cardiovascular complexity is a feature of healthy physiological regulation, which stems from the simultaneous activity of several cardiovascular reflexes and other non-reflex physiological mechanisms. It is manifested in the rich dynamics characterizing the spontaneous heart rate and blood pressure variability (HRV and BPV). The present study faces the challenge of disclosing the origin of short-term HRV and BPV from the statistical perspective offered by information theory. To dissect the physiological mechanisms giving rise to cardiovascular complexity in different conditions, measures of predictive information, information storage, information transfer and information modification were applied to the beat-to-beat variability of heart period (HP), systolic arterial pressure (SAP) and respiratory volume signal recorded non-invasively in 61 healthy young subjects at supine rest and during head-up tilt (HUT) and mental arithmetics (MA). Information decomposition enabled to assess simultaneously several expected and newly inferred physiological phenomena, including: (i) the decreased complexity of HP during HUT and the increased complexity of SAP during MA; (ii) the suppressed cardiorespiratory information transfer, related to weakened respiratory sinus arrhythmia, under both challenges; (iii) the altered balance of the information transferred along the two arms of the cardiovascular loop during HUT, with larger baroreflex involvement and smaller feedforward mechanical effects; and (iv) an increased importance of direct respiratory effects on SAP during HUT, and on both HP and SAP during MA. We demonstrate that a decomposition of the information contained in cardiovascular oscillations can reveal subtle changes in system dynamics and improve our understanding of the complexity changes during physiological challenges.

An adaptive technique for multiscale approximate entropy (MAEbin) threshold (r) selection: application to heart rate variability (HRV) and systolic blood pressure variability (SBPV) under postural stress

Multiscale approximate entropy (MAE) is used to quantify the complexity of a time series as a function of time scale τ. Approximate entropy (ApEn) tolerance threshold selection 'r' is based on either: (1) arbitrary selection in the recommended range (0.1-0.25) times standard deviation of time series (2) or finding maximum ApEn (ApEnmax) i.e., the point where self-matches start to prevail over other matches and choosing the corresponding 'r' (rmax) as threshold (3) or computing rchon by empirically finding the relation between rmax, SD1/SD2 ratio and N using curve fitting, where, SD1 and SD2 are short-term and long-term variability of a time series respectively. None of these methods is gold standard for selection of 'r'. In our previous study [1], an adaptive procedure for selection of 'r' is proposed for approximate entropy (ApEn). In this paper, this is extended to multiple time scales using MAEbin and multiscale cross-MAEbin (XMAEbin). We applied this to simulations i.e. 50 realizations (n = 50) of random number series, fractional Brownian motion (fBm) and MIX (P) [1] series of data length of N = 300 and short term recordings of HRV and SBPV performed under postural stress from supine to standing. MAEbin and XMAEbin analysis was performed on laboratory recorded data of 50 healthy young subjects experiencing postural stress from supine to upright. The study showed that (i) ApEnbin of HRV is more than SBPV in supine position but is lower than SBPV in upright position (ii) ApEnbin of HRV decreases from supine i.e. 1.7324 ± 0.112 (mean ± SD) to upright 1.4916 ± 0.108 due to vagal inhibition (iii) ApEnbin of SBPV increases from supine i.e. 1.5535 ± 0.098 to upright i.e. 1.6241 ± 0.101 due sympathetic activation (iv) individual and cross complexities of RRi and systolic blood pressure (SBP) series depend on time scale under consideration (v) XMAEbin calculated using ApEnmax is correlated with cross-MAE calculated using ApEn (0.1-0.26) in steps of 0.02 at each time scale in supine and upright position and is concluded that ApEn0.26 has highest correlation at most scales (vi) choice of 'r' is critical in interpreting interactions between RRi and SBP and in ascertaining true complexity of the individual RRi and SBP series.

Network connectivity modulates power spectrum scale invariance

Measures of complexity are sensitive in detecting disease, which has made them attractive candidates for diagnostic biomarkers; one complexity measure that has shown promise in fMRI is power spectrum scale invariance (PSSI). Even if scale-free features of neuroimaging turn out to be diagnostically useful, however, their underlying neurobiological basis is poorly understood. Using modeling and simulations of a schematic prefrontal-limbic meso-circuit, with excitatory and inhibitory networks of nodes, we present here a framework for how network density within a control system can affect the complexity of signal outputs. Our model demonstrates that scale-free behavior, similar to that observed in fMRI PSSI data, can be obtained for sufficiently large networks in a context as simple as a linear stochastic system of differential equations, although the scale-free range improves when introducing more realistic, nonlinear behavior in the system. PSSI values (reflective of complexity) vary as a function of both input type (excitatory, inhibitory) and input density (mean number of long-range connections, or strength), independent of their node-specific geometric distribution. Signals show pink noise (1/f) behavior when excitatory and inhibitory influences are balanced. As excitatory inputs are increased and decreased, signals shift towards white and brown noise, respectively. As inhibitory inputs are increased and decreased, signals shift towards brown and white noise, respectively. The results hold qualitatively at the hemodynamic scale, which we modeled by introducing a neurovascular component. Comparing hemodynamic simulation results to fMRI PSSI results from 96 individuals across a wide spectrum of anxiety-levels, we show how our model can generate concrete and testable hypotheses for understanding how connectivity affects regulation of meso-circuits in the brain.

Comparison of heart and respiratory rate variability measures using an intermittent incremental submaximal exercise model

To better understand the alterations in cardiorespiratory variability during exercise, the present study characterized the patterns of change in heart rate variability (HRV), respiratory rate variability (RRV), and combined cardiorespiratory variability (HRV-RRV) during an intermittent incremental submaximal exercise model. Six males and six females completed a submaximal exercise protocol consisting of an initial baseline resting period followed by three 10-min bouts of exercise at 20%, 40%, and 60% of maximal aerobic capacity (V̇O2max). The R-R interval and interbreath interval variability were measured at baseline rest and throughout the submaximal exercise. A group of 93 HRV, 83 RRV, and 28 HRV-RRV measures of variability were tracked over time through a windowed analysis using a 5-min window size and 30-s window step. A total of 91 HRV measures were able to detect the presence of exercise, whereas only 46 RRV and 3 HRV-RRV measures were able to detect the same stimulus. Moreover, there was a loss of overall HRV and RRV, loss of complexity of HRV and RRV, and loss of parasympathetic modulation of HRV (up to 40% V̇O2max) with exercise. Conflicting changes in scale-invariant structure of HRV and RRV with increases in exercise intensity were also observed. In summary, in this simultaneous evaluation of HRV and RRV, we found more consistent changes across HRV metrics compared with RRV and HRV-RRV.

In the spotlight: biomedical signal processing--a well established discipline or a paradigm to promising integrated visions?

Biomedical signals carry fundamental information about the nature and the status of the living systems under study. A proper processing of these signals obtains useful physiological and clinical information. A closer integration between signal processing and modeling of the relevant biological systems is capable to directly attribute pathophysiological meaning to the parameters obtained from the processing and vice versa. Another issue of great interest in which BSP plays an important role is the Brain-Computer Interface (BCI) or Brain-Machine Interface (BMI) where fast and reliable signal processing approaches are fundamental for a practical implementation. The physiological mechanisms underlying these heart rate variability findings are supposed to be related to stochastic processes at the cellular level, to influence of respiration on the heart rate, and to the interactions of the multiple feedback loops regulating the cardiovascular system. Another important area in which BSP plays a pivotal role is the "computational genomics and proteomics." It is true that "traditional" biomedical engineering studies biomedical signals which are obtained at the level of the major physiological systems.

The virtual liver: a multidisciplinary, multilevel challenge for systems biology

The liver is the central metabolic organ in human physiology, with functions that are fundamentally important to the detoxification of xenobiotics (drugs), the maintenance of homeostasis of numerous blood metabolites, and the production of mediators of the acute phase response. Liver toxicity, whether actual or implied is the reason for the failure of a significant proportion of many promising novel medicines that consequently never reach the market, and diseases such as atherosclerosis, diabetes, and fatty liver diseases, that are a major burden on current health resources, are directly linked to functional and structural disorders of the liver. This article presents the concepts and approaches underpinning one of the most exciting and ambitious modeling projects in the field of systems biology and systems medicine. This major multidisciplinary research program is aimed at developing a whole-organ model of the human liver, representing its central physiological functions under normal and pathological conditions The model will be composed of a larger battery of interconnected submodels representing liver anatomy and physiology, integrating processes across hierarchical levels in space, time, and structural organization. In this article, we outline the general architecture of the liver model and present first step taken to reach this ambitious goal.

A mouse model of high trait anxiety shows reduced heart rate variability that can be reversed by anxiolytic drug treatment

Increasing evidence suggests that specific physiological measures may serve as biomarkers for successful treatment to alleviate symptoms of pathological anxiety. Studies of autonomic function investigating parameters such as heart rate (HR), HR variability and blood pressure (BP) indicated that HR variability is consistently reduced in anxious patients, whereas HR and BP data show inconsistent results. Therefore, HR and HR variability were measured under various emotionally challenging conditions in a mouse model of high innate anxiety (high anxiety behaviour; HAB) vs. control normal anxiety-like behaviour (NAB) mice. Baseline HR, HR variability and activity did not differ between mouse lines. However, after cued Pavlovian fear conditioning, both elevated tachycardia and increased fear responses were observed in HAB mice compared to NAB mice upon re-exposure to the conditioning stimulus serving as the emotional stressor. When retention of conditioned fear was tested in the home cage, HAB mice again displayed higher fear responses than NAB mice, while the HR responses were similar. Conversely, in both experimental settings HAB mice consistently exhibited reduced HR variability. Repeated administration of the anxiolytic NK1 receptor antagonist L-822429 lowered the conditioned fear response and shifted HR dynamics in HAB mice to a more regular pattern, similar to that in NAB mice. Additional receiver-operating characteristic (ROC) analysis demonstrated the high specificity and sensitivity of HR variability to distinguish between normal and high anxiety trait. These findings indicate that assessment of autonomic response in addition to freezing might be a useful indicator of the efficacy of novel anxiolytic treatments.

Review and classification of variability analysis techniques with clinical applications

Analysis of patterns of variation of time-series, termed variability analysis, represents a rapidly evolving discipline with increasing applications in different fields of science. In medicine and in particular critical care, efforts have focussed on evaluating the clinical utility of variability. However, the growth and complexity of techniques applicable to this field have made interpretation and understanding of variability more challenging. Our objective is to provide an updated review of variability analysis techniques suitable for clinical applications. We review more than 70 variability techniques, providing for each technique a brief description of the underlying theory and assumptions, together with a summary of clinical applications. We propose a revised classification for the domains of variability techniques, which include statistical, geometric, energetic, informational, and invariant. We discuss the process of calculation, often necessitating a mathematical transform of the time-series. Our aims are to summarize a broad literature, promote a shared vocabulary that would improve the exchange of ideas, and the analyses of the results between different studies. We conclude with challenges for the evolving science of variability analysis.

Quantification of compensatory processes of postnatal hypoxia in newborn piglets applying short-term nonlinear dynamics analysis

Background

Newborn mammals suffering from moderate hypoxia during or after birth are able to compensate a transitory lack of oxygen by adapting their vital functions. Exposure to hypoxia leads to an increase in the sympathetic tone causing cardio-respiratory response, peripheral vasoconstriction and vasodilatation in privileged organs like the heart and brain. However, there is only limited information available about the time and intensity changes of the underlying complex processes controlled by the autonomic nervous system.

Methods

In this study an animal model involving seven piglets was used to examine an induced state of circulatory redistribution caused by moderate oxygen deficit. In addition to the main focus on the complex dynamics occurring during sustained normocapnic hypoxia, the development of autonomic regulation after induced reoxygenation had been analysed. For this purpose, we first introduced a new algorithm to prove stationary conditions in short-term time series. Then we investigated a multitude of indices from heart rate and blood pressure variability and from bivariate interactions, also analysing respiration signals, to quantify the complexity of vegetative oscillations influenced by hypoxia.

Results

The results demonstrated that normocapnic hypoxia causes an initial increase in cardiovascular complexity and variability, which decreases during moderate hypoxia lasting one hour (p < 0.004). After reoxygenation, cardiovascular complexity parameters returned to pre-hypoxic values (p < 0.003), however not respiratory-related complexity parameters.

Conclusions

In conclusion, indices from linear and nonlinear dynamics reflect considerable temporal changes of complexity in autonomous cardio-respiratory regulation due to normocapnic hypoxia shortly after birth. These findings might be suitable for non-invasive clinical monitoring of hypoxia-induced changes of autonomic regulation in newborn humans.

The effect of orthostatic stress on multiscale entropy of heart rate and blood pressure

Cardiovascular control acts over multiple time scales, which introduces a significant amount of complexity to heart rate and blood pressure time series. Multiscale entropy (MSE) analysis has been developed to quantify the complexity of a time series over multiple time scales. In previous studies, MSE analyses identified impaired cardiovascular control and increased cardiovascular risk in various pathological conditions. Despite the increasing acceptance of the MSE technique in clinical research, information underpinning the involvement of the autonomic nervous system in the MSE of heart rate and blood pressure is lacking. The objective of this study is to investigate the effect of orthostatic challenge on the MSE of heart rate and blood pressure variability (HRV, BPV) and the correlation between MSE (complexity measures) and traditional linear (time and frequency domain) measures. MSE analysis of HRV and BPV was performed in 28 healthy young subjects on 1000 consecutive heart beats in the supine and standing positions. Sample entropy values were assessed on scales of 1-10. We found that MSE of heart rate and blood pressure signals is sensitive to changes in autonomic balance caused by postural change from the supine to the standing position. The effect of orthostatic challenge on heart rate and blood pressure complexity depended on the time scale under investigation. Entropy values did not correlate with the mean values of heart rate and blood pressure and showed only weak correlations with linear HRV and BPV measures. In conclusion, the MSE analysis of heart rate and blood pressure provides a sensitive tool to detect changes in autonomic balance as induced by postural change.

Biomedical signal and image processing

Generally, physiological modeling and biomedical signal processing constitute two important paradigms of biomedical engineering (BME): their fundamental concepts are taught starting from undergraduate studies and are more completely dealt with in the last years of graduate curricula, as well as in Ph.D. courses. Traditionally, these two cultural aspects were separated, with the first one more oriented to physiological issues and how to model them and the second one more dedicated to the development of processing tools or algorithms to enhance useful information from clinical data. A practical consequence was that those who did models did not do signal processing and vice versa. However, in recent years,the need for closer integration between signal processing and modeling of the relevant biological systems emerged very clearly [1], [2]. This is not only true for training purposes(i.e., to properly prepare the new professional members of BME) but also for the development of newly conceived research projects in which the integration between biomedical signal and image processing (BSIP) and modeling plays a crucial role. Just to give simple examples, topics such as brain–computer machine or interfaces,neuroengineering, nonlinear dynamical analysis of the cardiovascular (CV) system,integration of sensory-motor characteristics aimed at the building of advanced prostheses and rehabilitation tools, and wearable devices for vital sign monitoring and others do require an intelligent fusion of modeling and signal processing competences that are certainly peculiar of our discipline of BME.

Electrocardiologic and related methods of non-invasive detection and risk stratification in myocardial ischemia: state of the art and perspectives

Background: Electrocardiographic methods still provide the bulk of cardiovascular diagnostics. Cardiac ischemia is associated with typical alterations in cardiac biosignals that have to be measured, analyzed by mathematical algorithms and allegorized for further clinical diagnostics. The fast growing fields of biomedical engineering and applied sciences are intensely focused on generating new approaches to cardiac biosignal analysis for diagnosis and risk stratification in myocardial ischemia.

Objectives: To present and review the state of the art in and new approaches to electrocardiologic methods for non-invasive detection and risk stratification in coronary artery disease (CAD) and myocardial ischemia; secondarily, to explore the future perspectives of these methods.

Methods: In follow-up to the Expert Discussion at the 2008 Workshop on "Biosignal Analysis" of the German Society of Biomedical Engineering in Potsdam, Germany, we comprehensively searched the pertinent literature and databases and compiled the results into this review. Then, we categorized the state-of-the-art methods and selected new approaches based on their applications in detection and risk stratification of myocardial ischemia. Finally, we compared the pros and cons of the methods and explored their future potentials for cardiology.

Results: Resting ECG, particularly suited for detecting ST-elevation myocardial infarctions, and exercise ECG, for the diagnosis of stable CAD, are state-of-the-art methods. New exercise-free methods for detecting stable CAD include cardiogoniometry (CGM); methods for detecting acute coronary syndrome without ST elevation are Body Surface Potential Mapping, functional imaging and CGM. Heart rate variability and blood pressure variability analyses, microvolt T-wave alternans and signal-averaged ECG mainly serve in detecting and stratifying the risk for lethal arrythmias in patients with myocardial ischemia or previous myocardial infarctions. Telemedicine and ambient-assisted living support the electrocardiological monitoring of at-risk patients.

Conclusions: There are many promising methods for the exercise-free, non-invasive detection of CAD and myocardial ischemia in the stable and acute phases. In the coming years, these new methods will help enhance state-of-the-art procedures in routine diagnostics. The future can expect that equally novel methods for risk stratification and telemedicine will transition into clinical routine.