A real-time automated point-process method for the detection and correction of erroneous and ectopic heartbeats

Monitoring surgical nociception using multisensor physiological models

Monitoring nociception, the flow of information associated with harmful stimuli through the nervous system even during unconsciousness, is critical for proper anesthesia care during surgery. Currently, this is done by tracking heart rate and blood pressure by eye. Monitoring objectively a patient's nociceptive state remains a challenge, causing drugs to often be over- or underdosed intraoperatively. Inefficient management of surgical nociception may lead to more complex postoperative pain management and side effects such as postoperative cognitive dysfunction, particularly in elderly patients. We collected a comprehensive and multisensor prospective observational dataset focused on surgical nociception (101 surgeries, 18,582 min, and 49,878 nociceptive stimuli), including annotations of all nociceptive stimuli occurring during surgery and medications administered. Using this dataset, we developed indices of autonomic nervous system activity based on physiologically and statistically rigorous point process representations of cardiac action potentials and sweat gland activity. Next, we constructed highly interpretable supervised and unsupervised models with appropriate inductive biases that quantify surgical nociception throughout surgery. Our models track nociceptive stimuli more accurately than existing nociception monitors. We also demonstrate that the characterizing signature of nociception learned by our models resembles the known physiology of the response to pain. Our work represents an important step toward objective multisensor physiology-based markers of surgical nociception. These markers are derived from an in-depth characterization of nociception as measured during surgery itself rather than using other experimental models as surrogates for surgical nociception.

Publication guidelines for human heart rate and heart rate variability studies in psychophysiology-Part 1: Physiological underpinnings and foundations of measurement

This Committee Report provides methodological, interpretive, and reporting guidance for researchers who use measures of heart rate (HR) and heart rate variability (HRV) in psychophysiological research. We provide brief summaries of best practices in measuring HR and HRV via electrocardiographic and photoplethysmographic signals in laboratory, field (ambulatory), and brain-imaging contexts to address research questions incorporating measures of HR and HRV. The Report emphasizes evidence for the strengths and weaknesses of different recording and derivation methods for measures of HR and HRV. Along with this guidance, the Report reviews what is known about the origin of the heartbeat and its neural control, including factors that produce and influence HRV metrics. The Report concludes with checklists to guide authors in study design and analysis considerations, as well as guidance on the reporting of key methodological details and characteristics of the samples under study. It is expected that rigorous and transparent recording and reporting of HR and HRV measures will strengthen inferences across the many applications of these metrics in psychophysiology. The prior Committee Reports on HR and HRV are several decades old. Since their appearance, technologies for human cardiac and vascular monitoring in laboratory and daily life (i.e., ambulatory) contexts have greatly expanded. This Committee Report was prepared for the Society for Psychophysiological Research to provide updated methodological and interpretive guidance, as well as to summarize best practices for reporting HR and HRV studies in humans.

Transcranial direct current stimulation improves heart rate variability: A systematic review and meta-analysis

BACKGROUND

Heart rate variability (HRV) is a useful tool for evaluating cardiovascular autonomic nervous system (ANS) functions. This systematic review and meta-analysis examined the potential effects of transcranial direct current stimulation (tDCS) protocols on HRV parameters.

METHODS

This study acquired 97 comparisons from 24 qualified studies for data synthesis. Using standardized mean difference (SMD), individual and overall effect sizes were estimated to show differences in HRV variables between active tDCS and sham stimulation conditions. More positive effect size values indicated that active tDCS caused greater increases in HRV than sham stimulation. Furthermore, moderator variable analyses were performed to determine whether changes in HRV variables differed depending on (a) task types (physical stress versus psychological stress versus resting condition), (b) targeted brain regions, (c) stimulation polarity, (d) characteristics of participants, and (e) specific HRV variables. Finally, we used meta-regression analyses to determine whether different tDCS parameters (i.e., the number of tDCS sessions, stimulation duration, and density) were associated with changes in HRV patterns.

RESULTS

The random-effects model meta-analysis showed that tDCS protocols significantly improved HRV variables (SMD = 0.400; P < 0.001). Moreover, for increasing HRV during the physical stress task (SMD = 1.352; P = 0.001), anodal stimulation on the M1 was effective, while combined polarity stimulation on the PFC improved HRV during the psychological stress task (SMD = 0.550; P < 0.001) and resting condition (SMD = 0.192; P = 0.012). Additional moderator variables and meta-regression analyses failed to show that tDCS protocols had positive effects in certain conditions, such as different stimulus polarity, characteristics of participants, specific HRV variables, and tDCS parameters.

CONCLUSION

These findings tentatively suggest that using tDCS protocols to stimulate optimal targeted brain areas may be effective in improving HRV patterns potentially related to cardiovascular ANS functions.

Attention to cardiac sensations enhances the heartbeat-evoked potential during exhalation

Respiration and cardiac activity intricately interact through complex physiological mechanisms. The heartbeat-evoked potential (HEP) is an EEG fluctuation reflecting the cortical processing of cardiac signals. We recently found higher HEP amplitude during exhalation than inhalation during a task involving attention to cardiac sensations. This may have been due to reduced cardiac perception during inhalation and heightened perception during exhalation through attentional mechanisms. To investigate relationships between HEP, attention, and respiration, we introduced an experimental setup that included tasks related to cardiac and respiratory interoceptive and exteroceptive attention. Results revealed HEP amplitude increases during the interoceptive tasks over fronto-central electrodes. When respiratory phases were taken into account, HEP increases were primarily driven by heartbeats recorded during exhalation, specifically during the cardiac interoceptive task, while inhalation had minimal impact. These findings emphasize the role of respiration in cardiac interoceptive attention and could have implications for respiratory interventions to fine-tune cardiac interoception.

Intracortical brain-heart interplay: An EEG model source study of sympathovagal changes

The interplay between cerebral and cardiovascular activity, known as the functional brain-heart interplay (BHI), and its temporal dynamics, have been linked to a plethora of physiological and pathological processes. Various computational models of the brain-heart axis have been proposed to estimate BHI non-invasively by taking advantage of the time resolution offered by electroencephalograph (EEG) signals. However, investigations into the specific intracortical sources responsible for this interplay have been limited, which significantly hampers existing BHI studies. This study proposes an analytical modeling framework for estimating the BHI at the source-brain level. This analysis relies on the low-resolution electromagnetic tomography sources localization from scalp electrophysiological recordings. BHI is then quantified as the functional correlation between the intracortical sources and cardiovascular dynamics. Using this approach, we aimed to evaluate the reliability of BHI estimates derived from source-localized EEG signals as compared with prior findings from neuroimaging methods. The proposed approach is validated using an experimental dataset gathered from 32 healthy individuals who underwent standard sympathovagal elicitation using a cold pressor test. Additional resting state data from 34 healthy individuals has been analysed to assess robustness and reproducibility of the methodology. Experimental results not only confirmed previous findings on activation of brain structures affecting cardiac dynamics (e.g., insula, amygdala, hippocampus, and anterior and mid-cingulate cortices) but also provided insights into the anatomical bases of brain-heart axis. In particular, we show that the bidirectional activity of electrophysiological pathways of functional brain-heart communication increases during cold pressure with respect to resting state, mainly targeting neural oscillations in theδ $$ \delta $$ ,β $$ \beta $$ , andγ $$ \gamma $$ bands. The proposed approach offers new perspectives for the investigation of functional BHI that could also shed light on various pathophysiological conditions.

Integrative neuro-cardiovascular dynamics in response to test anxiety: A brain-heart axis study

Test anxiety (TA), a recognized form of social anxiety, is the most prominent cause of anxiety among students and, if left unmanaged, can escalate to psychiatric disorders. TA profoundly impacts both central and autonomic nervous systems, presenting as a dual manifestation of cognitive and autonomic components. While limited studies have explored the physiological underpinnings of TA, none have directly investigated the intricate interplay between the CNS and ANS in this context. In this study, we introduce a non-invasive, integrated neuro-cardiovascular approach to comprehensively characterize the physiological responses of 27 healthy subjects subjected to test anxiety induced via a simulated exam scenario. Our experimental findings highlight that an isolated analysis of electroencephalographic and heart rate variability data fails to capture the intricate information provided by a brain-heart axis assessment, which incorporates an analysis of the dynamic interaction between the brain and heart. With respect to resting state, the simulated examination induced a decrease in the neural control onto heartbeat dynamics at all frequencies, while the studying condition induced a decrease in the ascending heart-to-brain interplay at EEG oscillations up to 12Hz. This underscores the significance of adopting a multisystem perspective in understanding the complex and especially functional directional mechanisms underlying test anxiety.

A Unified Framework for Investigating Aperiodic and Periodic Components in the Hearbeat Dynamics Spectrum: a Feasibility Study

Heart Rate Variability (HRV) series is a widely used, non-invasive, and easy-to-acquire time-resolved signal for evaluating autonomic control on cardiovascular activity. Despite the recognition that heartbeat dynamics contains both periodic and aperiodic components, the majority of HRV modeling studies concentrate on only one component. On the one hand, there are models based on self-similarity and 1/f behavior that focus on the aperiodic component; on the other hand, there is the conventional division of the spectral domain into narrow-band oscillations, which considers HRV as a combination of periodic components. Taking inspiration from a recent parametrization of EEG power spectra, here we evaluate the applicability of a unified modeling framework to quantitatively assess heartbeat dynamics spectra as a mixture of aperiodic and periodic components. The proposed model is applied on publicly-available, real HRV series collected during postural changes from 10 healthy subjects. Results show that the proposed modeling effectively characterizes different experimental conditions and may complement HRV standard analysis defined in the frequency domain.

Dynamic fluctuations in ascending heart-to-brain communication under mental stress

Dynamical information exchange between central and autonomic nervous systems, as referred to functional brain-heart interplay, occurs during emotional and physical arousal. It is well documented that physical and mental stress lead to sympathetic activation. Nevertheless, the role of autonomic inputs in nervous system-wise communication under mental stress is yet unknown. In this study, we estimated the causal and bidirectional neural modulations between electroencephalogram (EEG) oscillations and peripheral sympathetic and parasympathetic activities using a recently proposed computational framework for a functional brain-heart interplay assessment, namely the sympathovagal synthetic data generation model. Mental stress was elicited in 37 healthy volunteers by increasing their cognitive demands throughout three tasks associated with increased stress levels. Stress elicitation induced an increased variability in sympathovagal markers, as well as increased variability in the directional brain-heart interplay. The observed heart-to-brain interplay was primarily from sympathetic activity targeting a wide range of EEG oscillations, whereas variability in the efferent direction seemed mainly related to EEG oscillations in the γ band. These findings extend current knowledge on stress physiology, which mainly referred to top-down neural dynamics. Our results suggest that mental stress may not cause an increase in sympathetic activity exclusively as it initiates a dynamic fluctuation within brain-body networks including bidirectional interactions at a brain-heart level. We conclude that directional brain-heart interplay measurements may provide suitable biomarkers for a quantitative stress assessment and bodily feedback may modulate the perceived stress caused by increased cognitive demand.

Characterization of autonomic states by complex sympathetic and parasympathetic dynamics

Assessment of heartbeat dynamics provides a promising framework for non-invasive monitoring of cardiovascular and autonomic states. Nevertheless, the non-specificity of such measurements among clinical populations and healthy conditions associated with different autonomic states severely limits their applicability and exploitation in naturalistic conditions. This limitation arises especially when pathological or postural change-related sympathetic hyperactivity is compared to autonomic changes across age and experimental conditions. In this frame, we investigate the intrinsic irregularity and complexity of cardiac sympathetic and vagal activity series in different populations, which are associated with different cardiac autonomic dynamics. Sample entropy, fuzzy entropy, and distribution entropy are calculated on the recently proposed sympathetic and parasympathetic activity indices (SAI and PAI) series, which are derived from publicly available heartbeat series of congestive heart failure patients, elderly and young subjects watching a movie in the supine position, and healthy subjects undergoing slow postural changes. Results show statistically significant differences between pathological/old subjects and young subjects in the resting state and during slow tilt, with interesting trends in SAI- and PAI-related entropy values. Moreover, while CHF patients and healthy subjects in upright position show the higher cardiac sympathetic activity, elderly and young subjects in resting state showed higher vagal activity. We conclude that quantification of intrinsic cardiac complexity from sympathetic and vagal dynamics may provide new physiology insights and improve on the non-specificity of heartbeat-derived biomarkers.

Modeling brain-heart interactions from Poincaré plot-derived measures of sympathetic-vagal activity

Recent studies suggest that the interaction between the brain and heart plays a key role in cognitive processes, and measuring these interactions is crucial for understanding the interaction between the central and autonomic nervous systems. However, studying this bidirectional interplay presents methodological challenges, and there is still much room for exploration. This paper presents a new computational method called the Poincaré Sympathetic-Vagal Synthetic Data Generation Model (PSV-SDG) for estimating brain-heart interactions. The PSV-SDG combines EEG and cardiac sympathetic-vagal dynamics to provide time-varying and bidirectional estimators of mutual interplay. The method is grounded in the Poincaré plot, a heart rate variability method to estimate sympathetic-vagal activity that can account for potential non-linearities. This algorithm offers a new approach and computational tool for functional assessment of the interplay between EEG and cardiac sympathetic-vagal activity. The method is implemented in MATLAB under an open-source license. • A new brain-heart interaction modeling approach is proposed. • The modeling is based on coupled synthetic data generators of EEG and heart rate series. • Sympathetic and vagal activities are gathered from Poincaré plot geometry.

A two-step pre-processing tool to remove Gaussian and ectopic noise for heart rate variability analysis

Artifacts in the Electrocardiogram (ECG) degrade the quality of the recorded signal and are not conducive to heart rate variability (HRV) analysis. The two types of noise most often found in ECG recordings are technical and physiological artifacts. Current preprocessing methods primarily attend to ectopic beats but do not consider technical issues that affect the ECG. A secondary aim of this study was to investigate the effect of increasing increments of artifacts on 24 of the most used HRV measures. A two-step preprocessing approach for denoising HRV is introduced which targets each type of noise separately. First, the technical artifacts in the ECG are eliminated by applying complete ensemble empirical mode decomposition with adaptive noise. The second step removes physiological artifacts from the HRV signal using a combination filter of single dependent rank order mean and an adaptive filtering algorithm. The performance of the two-step pre-processing tool showed a high correlation coefficient of 0.846 and RMSE value of 7.69 × 10^-5 for 6% of added ectopic beats and 6 dB Gaussian noise. All HRV measures studied except HF peak and LF peak are significantly affected by both types of noise. Frequency measures of Total power, HF power, and LF power and fragmentation measures; PAS, PIP, and PSS are the most sensitive to both types of noise.

Detection and categorization of severe cardiac disorders based solely on heart period measurements

Cardiac disorders are common conditions associated with a high mortality rate. Due to their potential for causing serious symptoms, it is desirable to constantly monitor cardiac status using an accessible device such as a smartwatch. While electrocardiograms (ECGs) can make the detailed diagnosis of cardiac disorders, the examination is typically performed only once a year for each individual during health checkups, and it requires expert medical practitioners to make comprehensive judgments. Here we describe a newly developed automated system for alerting individuals about cardiac disorders solely by measuring a series of heart periods. For this purpose, we examined two metrics of heart rate variability (HRV) and analyzed 1-day ECG recordings of more than 1,000 subjects in total. We found that a metric of local variation was more efficient than conventional HRV metrics for alerting cardiac disorders, and furthermore, that a newly introduced metric of local-global variation resulted in superior capacity for discriminating between premature contraction and atrial fibrillation. Even with a 1-minute recording of heart periods, our new detection system had a diagnostic performance even better than that of the conventional analysis method applied to a 1-day recording.

Comprehensive HRV estimation pipeline in Python using Neurokit2: Application to sleep physiology

NeuroKit2 is a Python Toolbox for Neurophysiological Signal Processing. The presented method is an adaptation of NeuroKit2 to simplify and automate computation of the various mathematical estimates of heart rate variability (HRV) or similar time series. By default, the present approach accepts as input electrocardiogram's R-R intervals (RRIs) or peak times, i.e., timestamp of each consecutive R peak, but the RRIs or peak times can also stem from other biosensors such as photoplethysmography (PPGs) or represent more general kinds of biological or non-biological time series oscillations. The data may be derived from a single or several sources such as conventional univariate heart rate time series or intermittently weakly coupled fetal and maternal heart rate data. The method describes preprocessing and computation of an output of 124 HRV measures including measures with a dynamic, time-series-specific optimal time delay-based complexity estimation with a user-definable time window length. I also provide an additional layer of HRV estimation looking at the temporal fluctuations of the HRV estimates themselves, an approach not yet widely used in the field, yet showing promise (doi: 10.3389/fphys.2017.01112). To demonstrate the application of the methodology, I present an approach to studying the dynamic relationships between sleep state architecture and multi-dimensional HRV metrics in 31 subjects. NeuroKit2's documentation is extensive. Here, I attempted to simplify things summarizing all you need to produce the most extensive HRV estimation output available to date as open source and all in one place. The presented Jupyter notebooks allow the user to run HRV analyses quickly and at scale on univariate or multivariate time-series data. I gratefully acknowledge the excellent support from the NeuroKit team.•Univariate or multivariate time series input; ingestion, preprocessing, and computation of 124 HRV metrics.•Estimation of intra- and inter-individual higher order temporal fluctuations of HRV metrics.•Application to a sleep dataset recorded using Apple Watch and expert sleep labeling.

A new framework for modeling the bidirectional interplay between brain oscillations and cardiac sympathovagal activity

The study of functional brain-heart interplay (BHI) aims to describe the dynamical interactions between central and peripheral autonomic nervous systems. Here, we introduce the Sympathovagal Synthetic Data Generation Model, which constitutes a new computational framework for the assessment of functional BHI. The model estimates the bidirectional interplay with novel quantifiers of cardiac sympathovagal activity gathered from Laguerre expansions of RR series (from the ECG), as an alternative to the classical spectral analysis. The main features of the model are time-varying coupling coefficients linking Electroencephalography (EEG) oscillations and cardiac sympathetic or parasympathetic activity, for either ascending or descending direction of the information flow. In this proof-of-concept study, functional BHI is quantified in the direction from-heart-to-brain, on data from 16 human volunteers undergoing a cold-pressor test. Results show that thermal stress induces heart-to-brain functional interplay originating from sympathetic and parasympathetic activities and sustaining EEG oscillations mainly in the δ and γ bands. The proposed computational framework could provide a viable tool for the functional assessment of the causal interplay between cortical and cardiac sympathovagal dynamics.

Cardiac sympathetic-vagal activity initiates a functional brain-body response to emotional arousal

We investigate the temporal dynamics of brain and cardiac activities in healthy subjects who underwent an emotional elicitation through videos. We demonstrate that, within the first few seconds, emotional stimuli modulate heartbeat activity, which in turn stimulates an emotion intensity (arousal)–specific cortical response. The emotional processing is then sustained by a bidirectional brain–heart interplay, where the perceived arousal level modulates the amplitude of ascending heart-to-brain neural information flow. These findings may constitute fundamental knowledge linking neurophysiology and psychiatric disorders, including the link between depressive symptoms and cardiovascular disorders.

A century-long debate on bodily states and emotions persists. While the involvement of bodily activity in emotion physiology is widely recognized, the specificity and causal role of such activity related to brain dynamics has not yet been demonstrated. We hypothesize that the peripheral neural control on cardiovascular activity prompts and sustains brain dynamics during an emotional experience, so these afferent inputs are processed by the brain by triggering a concurrent efferent information transfer to the body. To this end, we investigated the functional brain–heart interplay under emotion elicitation in publicly available data from 62 healthy subjects using a computational model based on synthetic data generation of electroencephalography and electrocardiography signals. Our findings show that sympathovagal activity plays a leading and causal role in initiating the emotional response, in which ascending modulations from vagal activity precede neural dynamics and correlate to the reported level of arousal. The subsequent dynamic interplay observed between the central and autonomic nervous systems sustains the processing of emotional arousal. These findings should be particularly revealing for the psychophysiology and neuroscience of emotions.

Validation of Continuous Monitoring System for Epileptic Users in Outpatient Settings

Epilepsy is a chronic disease with a significant social impact, given that the patients and their families often live conditioned by the possibility of an epileptic seizure and its possible consequences, such as accidents, injuries, or even sudden unexplained death. In this context, ambulatory monitoring allows the collection of biomedical data about the patients’ health, thus gaining more knowledge about the physiological state and daily activities of each patient in a more personalized manner. For this reason, this article proposes a novel monitoring system composed of different sensors capable of synchronously recording electrocardiogram (ECG), photoplethysmogram (PPG), and ear electroencephalogram (EEG) signals and storing them for further processing and analysis in a microSD card. This system can be used in a static and/or ambulatory way, providing information about the health state through features extracted from the ear EEG signal and the calculation of the heart rate variability (HRV) and pulse travel time (PTT). The different applied processing techniques to improve the quality of these signals are described in this work. A novel algorithm used to compute HRV and PTT robustly and accurately in ambulatory settings is also described. The developed device has also been validated and compared with other commercial systems obtaining similar results. In this way, based on the quality of the obtained signals and the low variability of the computed parameters, even in ambulatory conditions, the developed device can potentially serve as a support tool for clinical decision-taking stages.

Neural Correlates of Non-ordinary States of Consciousness in Pranayama Practitioners: The Role of Slow Nasal Breathing

The modulatory effect of nasal respiration on integrative brain functions and hence consciousness has recently been unambiguously demonstrated. This effect is sustained by the olfactory epithelium mechanical sensitivity complemented by the existence of massive projections between the olfactory bulb and the prefrontal cortex. However, studies on slow nasal breathing (SNB) in the context of contemplative practices have sustained the fundamental role of respiratory vagal stimulation, with little attention to the contribution of the olfactory epithelium mechanical stimulation. This study aims at disentangling the effects of olfactory epithelium stimulation (proper of nasal breathing) from those related to respiratory vagal stimulation (common to slow nasal and mouth breathing). We investigated the psychophysiological (cardio-respiratory and electroencephalographic parameters) and phenomenological (perceived state of consciousness) aftereffects of SNB (epithelium mechanical – 2.5 breaths/min) in 12 experienced meditators. We compared the nasal breathing aftereffects with those observed after a session of mouth breathing at the same respiratory rate and with those related to a resting state condition. SNB induced (1) slowing of electroencephalography (EEG) activities (delta-theta bands) in prefrontal regions, (2) a widespread increase of theta and high-beta connectivity complemented by an increase of phase-amplitude coupling between the two bands in prefrontal and posterior regions belonging to the Default Mode Network, (3) an increase of high-beta networks small-worldness. (4) a higher perception of being in a non-ordinary state of consciousness. The emerging scenario strongly suggests that the effects of SNB, beyond the relative contribution of vagal stimulation, are mainly ascribable to olfactory epithelium stimulation. In conclusion, slow Pranayama breathing modulates brain activity and hence subjective experience up to the point of inducing a non-ordinary state of consciousness.

Functional assessment of bidirectional cortical and peripheral neural control on heartbeat dynamics: a brain-heart study on thermal stress

The study of functional brain-heart interplay (BHI) from non-invasive recordings has gained much interest in recent years. Previous endeavors aimed at understanding how the two dynamical systems exchange information, providing novel holistic biomarkers and important insights on essential cognitive aspects and neural system functioning. However, the interplay between cardiac sympathovagal and cortical oscillations still has much room for further investigation. In this study, we introduce a new computational framework for a functional BHI assessment, namely the Sympatho-Vagal Synthetic Data Generation Model, combining cortical (electroencephalography, EEG) and peripheral (cardiac sympathovagal) neural dynamics. The causal, bidirectional neural control on heartbeat dynamics was quantified on data gathered from 26 human volunteers undergoing a cold-pressor test. Results show that thermal stress induces heart-to-brain functional interplay sustained by EEG oscillations in the delta and gamma bands, primarily originating from sympathetic activity, whereas brain-to-heart interplay originates over central brain regions through sympathovagal control. The proposed methodology provides a viable computational tool for the functional assessment of the causal interplay between cortical and cardiac neural control.

Functional brain-heart interplay extends to the multifractal domain

The study of functional brain-heart interplay has provided meaningful insights in cardiology and neuroscience. Regarding biosignal processing, this interplay involves predominantly neural and heartbeat linear dynamics expressed via time and frequency domain-related features. However, the dynamics of central and autonomous nervous systems show nonlinear and multifractal behaviours, and the extent to which this behaviour influences brain-heart interactions is currently unknown. Here, we report a novel signal processing framework aimed at quantifying nonlinear functional brain-heart interplay in the non-Gaussian and multifractal domains that combines electroencephalography (EEG) and heart rate variability series. This framework relies on a maximal information coefficient analysis between nonlinear multiscale features derived from EEG spectra and from an inhomogeneous point-process model for heartbeat dynamics. Experimental results were gathered from 24 healthy volunteers during a resting state and a cold pressor test, revealing that synchronous changes between brain and heartbeat multifractal spectra occur at higher EEG frequency bands and through nonlinear/complex cardiovascular control. We conclude that significant bodily, sympathovagal changes such as those elicited by cold-pressure stimuli affect the functional brain-heart interplay beyond second-order statistics, thus extending it to multifractal dynamics. These results provide a platform to define novel nervous-system-targeted biomarkers. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Intensification of functional neural control on heartbeat dynamics in subclinical depression

Subclinical depression (dysphoria) is a common condition that may increase the risk of major depression and leads to impaired quality of life and severe comorbid somatic diseases. Despite its prevalence, specific biological markers are unknown; consequently, the identification of dysphoria currently relies exclusively on subjective clinical scores and structured interviews. Based on recent neurocardiology studies that link brain and cardiovascular disorders, it was hypothesized that multi-system biomarkers of brain-body interplay may effectively characterize dysphoria. Thus, an ad hoc computational technique was developed to quantify the functional bidirectional brain-heart interplay. Accordingly, 32-channel electroencephalographic and heart rate variability series were obtained from 24 young dysphoric adults and 36 healthy controls. All participants were females of a similar age, and results were obtained during a 5-min resting state. The experimental results suggest that a specific feature of dysphoria is linked to an augmented functional central-autonomic control to the heart, which originates from central, frontopolar, and occipital oscillations and acts through cardiovascular sympathovagal activity. These results enable further development of a large set of novel biomarkers for mood disorders based on comprehensive brain-body measurements.

Toward brain-heart computer interfaces: a study on the classification of upper limb movements using multisystem directional estimates

OBJECTIVE

Brain-computer interfaces (BCI) exploit computational features from brain signals to perform a given task. Despite recent neurophysiology and clinical findings indicating the crucial role of functional interplay between brain and cardiovascular dynamics in locomotion, heartbeat information remains to be included in common BCI systems. In this study, we exploit the multidimensional features of directional and functional interplay between electroencephalographic and heartbeat spectra to classify upper limb movements into three classes.

APPROACH

We gathered data from 26 healthy volunteers that performed 90 movements; the data were processed using a recently proposed framework for brain-heart interplay (BHI) assessment based on synthetic physiological data generation. Extracted BHI features were employed to classify, through sequential forward selection scheme and k-nearest neighbors algorithm, among resting state and three classes of movements according to the kind of interaction with objects.

MAIN RESULTS

The results demonstrated that the proposed brain-heart computer interface (BHCI) system could distinguish between rest and movement classes automatically with an average 90% of accuracy.

SIGNIFICANCE

Further, this study provides neurophysiology insights indicating the crucial role of functional interplay originating at the cortical level onto the heart in the upper limb neural control. The inclusion of functional BHI insights might substantially improve the neuroscientific knowledge about motor control, and this may lead to advanced BHCI systems performances.

Measuring Heart Rate Variability in Free-Living Conditions Using Consumer-Grade Photoplethysmography: Validation Study

Background

Heart rate variability (HRV) is used to assess cardiac health and autonomic nervous system capabilities. With the growing popularity of commercially available wearable technologies, the opportunity to unobtrusively measure HRV via photoplethysmography (PPG) is an attractive alternative to electrocardiogram (ECG), which serves as the gold standard. PPG measures blood flow within the vasculature using color intensity. However, PPG does not directly measure HRV; it measures pulse rate variability (PRV). Previous studies comparing consumer-grade PRV with HRV have demonstrated mixed results in short durations of activity under controlled conditions. Further research is required to determine the efficacy of PRV to estimate HRV under free-living conditions.

Objective

This study aims to compare PRV estimates obtained from a consumer-grade PPG sensor with HRV measurements from a portable ECG during unsupervised free-living conditions, including sleep, and examine factors influencing estimation, including measurement conditions and simple editing methods to limit motion artifacts.

Methods

A total of 10 healthy adults were recruited. Data from a Microsoft Band 2 and a Shimmer3 ECG unit were recorded simultaneously using a smartphone. Participants wore the devices for >90 min during typical day-to-day activities and while sleeping. After filtering, ECG data were processed using a combination of discrete wavelet transforms and peak-finding methods to identify R-R intervals. P-P intervals were edited for deletion using methods based on outlier detection and by removing sections affected by motion artifacts. Common HRV metrics were compared, including mean N-N, SD of N-N intervals, percentage of subsequent differences >50 ms (pNN50), root mean square of successive differences, low-frequency power (LF), and high-frequency power. Validity was assessed using root mean square error (RMSE) and Pearson correlation coefficient (R2).

Results

Data sets for 10 days and 9 corresponding nights were acquired. The mean RMSE was 182 ms (SD 48) during the day and 158 ms (SD 67) at night. R2 ranged from 0.00 to 0.66, with 2 of 19 (2 nights) trials considered moderate, 7 of 19 (2 days, 5 nights) fair, and 10 of 19 (8 days, 2 nights) poor. Deleting sections thought to be affected by motion artifacts had a minimal impact on the accuracy of PRV measures. Significant HRV and PRV differences were found for LF during the day and R-R, SDNN, pNN50, and LF at night. For 8 of the 9 matched day and night data sets, R2 values were higher at night (P=.08). P-P intervals were less sensitive to rapid R-R interval changes.

Conclusions

Owing to overall poor concurrent validity and inconsistency among participant data, PRV was found to be a poor surrogate for HRV under free-living conditions. These findings suggest that free-living HRV measurements would benefit from examining alternate sensing methods, such as multiwavelength PPG and wearable ECG.

Intrinsic Complexity of Sympathetic and Parasympathetic Dynamics from HRV series: a Preliminary Study on Postural Changes

The analysis of complex heartbeat dynamics has been widely used to characterize heartbeat autonomic control in healthy and pathological conditions. However, underlying physiological correlates of complexity measurements from heart rate variability (HRV) series have not been identified yet. To this extent, we investigated intrinsic irregularity and complexity of cardiac sympathetic and vagal activity time series during postural changes. We exploited our recently proposed HRV-based, time-varying Sympathetic and Parasympathetic Activity Indices (SAI and PAI) and performed Sample Entropy, Fuzzy Entropy, and Distribution Entropy calculations on publicly-available heartbeat series gathered from 10 healthy subjects undergoing resting state and passive slow tilt sessions. Results show significantly higher entropy values during the upright position than resting state in both SAI and PAI series. We conclude that an increase in HRV complexity resulting from postural changes may derive from sympathetic and vagal activities with higher complex dynamics.

An Inhomogeneous Point-process Model for the Assessment of the Brain-to-Heart Functional Interplay: a Pilot Study

We propose a novel computational framework for the estimation of functional directional brain-to-heart interplay in an instantaneous fashion. The framework is based on inhomogeneous point-process models for human heartbeat dynamics and employs inverse-Gaussian probability density functions characterizing the timing of R-peak events. The instantaneous estimation of the functional directional coupling is based on the definition of point-process transfer entropy, which is here retrieved from heart rate variability (HRV) and Electroencephalography (EEG) power spectral series gathered from 12 healthy subjects undergoing significant sympathovagal changes induced by a cold-pressor test. Results suggest that EEG oscillations dynamically influence heartbeat dynamics with specific time delays in the 30-60s and 90-120s ranges, and through a functional activity over specific cortical regions.

Placebo Analgesia, Nocebo Hyperalgesia, and the Cardiovascular System: A Qualitative Systematic Review

Background: Placebo/nocebo effects involve the autonomic nervous system, including cardiac activity, but studies have reported inconsistent findings on how cardiac activity is modulated following a placebo/nocebo effect. However, no systematic review has been conducted to provide a clear picture of cardiac placebo responses. Objective: The main goal of the present study is to review the effects of placebo analgesia and nocebo hyperalgesia on cardiac activity as measured by blood pressure, heart rate, and heart rate variability. Methods: Using several Boolean keyword combinations, the PubMed, EMBASE, PsycINFO, Cochrane Review Library, and ISI Web of Knowledge databases were searched until January 5, 2020, to find studies that analyzed blood pressure, heart rate, or heart rate variability indexes following a placebo analgesic/nocebo hyperalgesic effect. Results: Nineteen studies were found, with some reporting more than one index of cardiac activity; eight studies were on blood pressure, 14 studies on heart rate, and six on heart rate variability. No reliable association between placebo/nocebo effects and blood pressure or heart rate was found. However, placebo effects reduced, and nocebo effects increased low-frequency heart rate variability, and heart rate variability significantly predicted placebo effects in two studies. Conclusion: Placebo/nocebo effects can have reliable effects on heart rate variability, but not on heart rate and blood pressure.

Impact of sex and depressed mood on the central regulation of cardiac autonomic function

Cardiac autonomic dysregulation has been implicated in the comorbidity of major psychiatric disorders and cardiovascular disease, potentially through dysregulation of physiological responses to negative stressful stimuli (here, shortened to stress response). Further, sex differences in these comorbidities are substantial. Here, we tested the hypothesis that mood- and sex-dependent alterations in brain circuitry implicated in the regulation of the stress response are associated with reduced peripheral parasympathetic activity during negative emotional arousal. Fifty subjects (28 females) including healthy controls and individuals with major depression, bipolar psychosis and schizophrenia were evaluated. Functional magnetic resonance imaging and physiology (cardiac pulse) data were acquired during a mild visual stress reactivity challenge. Associations between changes in activity and functional connectivity of the stress response circuitry and variations in cardiovagal activity [normalized high frequency power of heart rate variability (HFn)] were evaluated using GLM analyses, including interactions with depressed mood and sex across disorders. Our results revealed that in women with high depressed mood, lower cardiovagal activity in response to negative affective stimuli was associated with greater activation of hypothalamus and right amygdala and reduced connectivity between hypothalamus and right orbitofrontal cortex, amygdala, and hippocampus. No significant associations were observed in women with low levels of depressed mood or men. Our results revealed mood- and sex-dependent interactions in the central regulation of cardiac autonomic activity in response to negative affective stimuli. These findings provide a potential pathophysiological mechanism for previously observed sex differences in the comorbidity of major depression and cardiovascular disease.

Quantifying Functional Links between Brain and Heartbeat Dynamics in the Multifractal Domain: a Preliminary Analysis

Quantification of brain-heart interplay (BHI) has mainly been performed in the time and frequency domains. However, such functional interactions are likely to involve nonlinear dynamics associated with the two systems. To this extent, in this preliminary study we investigate the functional coupling between multifractal properties of Electroencephalography (EEG) and Heart Rate Variability (HRV) series using a channel- and time scale-wise maximal information coefficient analysis. Experimental results were gathered from 24 healthy volunteers undergoing a resting state and a cold-pressure test, and suggest that significant changes between the two experimental conditions might be associated with nonlinear quantifiers of the multifractal spectrum. Particularly, major brain-heart functional coupling was associated with the secondorder cumulant of the multifractal spectrum. We conclude that a functional nonlinear relationship between brain- and heartbeat-related multifractal sprectra exist, with higher values associated with the resting state.

Assessing Autonomic Function from Electrodermal Activity and Heart Rate Variability During Cold-Pressor Test and Emotional Challenge

Standard functional assessment of autonomic nervous system (ANS) activity on cardiovascular control relies on spectral analysis of heart rate variability (HRV) series. However, difficulties in obtaining a reliable measure of sympathetic activity from HRV spectra limits the exploitation of sympatho-vagal metrics. On the other hand, measures of electrodermal activity (EDA) have been demonstrated to provide a reliable quantifier of sympathetic dynamics. In this study we propose novel indices of phasic autonomic regulation mechanisms by combining HRV and EDA correlates and thoroughly investigating their time-varying dynamics. HRV and EDA series were gathered from 26 healthy subjects during a cold-pressor test and emotional stimuli. Instantaneous linear and nonlinear (bispectral) estimates of vagal dynamics were obtained from HRV through inhomogeneous point-process models, and combined with a sensitive maker of sympathetic tone from EDA spectral power. A wavelet decomposition analysis was applied to estimate phasic components of the proposed sympatho-vagal indices. Results show significant statistical differences for the proposed indices between the cold-pressor elicitation and previous resting state. Furthermore, an accuracy of 73.08% was achieved for the automatic emotional valence recognition. The proposed nonlinear processing of phasic ANS markers brings novel insights on autonomic functioning that can be exploited in the field of affective computing and psychophysiology.

Uncovering complex central autonomic networks at rest: a functional magnetic resonance imaging study on complex cardiovascular oscillations

This study aims to uncover brain areas that are functionally linked to complex cardiovascular oscillations in resting-state conditions. Multi-session functional magnetic resonance imaging (fMRI) and cardiovascular data were gathered from 34 healthy volunteers recruited within the human connectome project (the '100-unrelated subjects' release). Group-wise multi-level fMRI analyses in conjunction with complex instantaneous heartbeat correlates (entropy and Lyapunov exponent) revealed the existence of a specialized brain network, i.e. a complex central autonomic network (CCAN), reflecting what we refer to as complex autonomic control of the heart. Our results reveal CCAN areas comprised the paracingulate and cingulate gyri, temporal gyrus, frontal orbital cortex, planum temporale, temporal fusiform, superior and middle frontal gyri, lateral occipital cortex, angular gyrus, precuneous cortex, frontal pole, intracalcarine and supracalcarine cortices, parahippocampal gyrus and left hippocampus. The CCAN visible at rest does not include the insular cortex, thalamus, putamen, amygdala and right caudate, which are classical CAN regions peculiar to sympatho-vagal control. Our results also suggest that the CCAN is mainly involved in complex vagal control mechanisms, with possible links with emotional processing networks.

Characterizing autonomic response to arousing visual-auditory multi-modal task in Autism Spectrum Disorder (ASD)

Sensory abnormalities are widespread in Autism Spectrum Disorder (ASD). However, their definition is still quite subjective and vague. Here we propose a novel approach for characterization of Autonomic Nervous System responses to sensory stimulation based on electrocardiogram (ECG) assessment. In particular, we develop a preliminary study where autonomic responses of both autistic (ASD = 5) and neurotypical (NT = 5) participants have been evaluated in terms of changes in responsiveness to repeated stimuli. Autonomic control has been estimated via high-frequency heart rate variability (HF-HRV) and low-frequency HRV (LF-HRV). Results show significant differences among groups for the HRV measures (p value = 0.0158), supported by expected changes of HF (p value = 0.0079) and LF (p value = 0.0079) trends over stimulations. We thus conclude that an overall decrease in autonomic arousal can give important insights for devising new habituation metrics in NT and ASD individuals.

Resting State Neural Correlates of Cardiac Sympathetic Dynamics in Healthy Subjects

Recent advances in functional Magnetic Resonance Imaging (fMRI) research have uncovered the existence of the central autonomic network (CAN), which comprises brain regions whose activity correlates with autonomic nervous system dynamics. By exploiting the spectral paradigm of heartbeat dynamics, cortical and sub-cortical areas functionally linked to vagal activity have been identified. However, due to methodological limitations, functional neural correlates of cardiac sympathetic dynamics remain uncharacterized. To this extent, we exploit the high spatiotemporal resolution of fMRI data from the Human Connectome Project to study the CAN activity by correlating a recently proposed instantaneous characterization of sympathetic activity (the sympathetic activity index - SAI) from heartbeat dynamics. SAI estimates are embedded into the probabilistic modeling of inhomogeneous point-processes, and are derived from a combination of disentangling coefficients linked to the orthonormal Laguerre functions. By analyzing resting state recordings from 34 young healthy people, we obtain positive correlations between instantaneous SAI estimates and a number of brain regions including frontal pole, insular cortex, frontal and temporal gyri, lateral occipital cortex, paracingulate and cingulate gyri, precuneus and temporal fusiform cortices, as well as thalamus, caudate nucleus, putamen, brain-stem, hippocampus, amygdala, and nucleus accumbens. Our findings significantly extend current knowledge on the CAN, opening new avenues in the characterization of healthy and pathological states in humans.

Noise Detection in Electrocardiography Signal for Robust Heart Rate Variability Analysis: A Deep Learning Approach

Heart rate variability (HRV) analysis is widely used to assess the sympathetic and parasympathetic tones. However, the quality of the derived HRV features is heavily dependent on the accuracy of QRS detection. Noisy electrocardiography (ECG) signals, such as those measured by wearable ECG patches, can lead to inaccuracies in the QRS detection and significantly impair the HRV analysis. Hence, it is critical to employ noise detection algorithms to identify the corrupted segments of the ECG signal and discard them from the analysis. This paper proposes a convolutional neural network to distinguish between usable and unusable ECG segments where usability is defined based on the accuracy of QRS detection. The results indicate that the proposed method has significantly lower error rates compared to both the baseline method (HRV analysis on the noisy signals) and a noise detection method based on four ECG signal quality indices and a support vector machines classifier.

A Bayesian Filtering Approach for Tracking Arousal From Binary and Continuous Skin Conductance Features

OBJECTIVE

Neuroanatomical structures within the cortical and sub-cortical brain regions process emotion and cause subsequent variations in signals such as skin conductance and electrocardiography. The signals often encode information in their continuous-valued amplitudes or waves as well as in their underlying impulsive events. We propose to track psychological arousal from this hybrid source of skin conductance information.

METHODS

We present a point process state-space method in tandem with Bayesian filtering for determining a continuous-valued arousal state from skin conductance measurements. To perform state estimation, we relate arousal to binary- and continuous-valued observations derived from the phasic and tonic parts of a skin conductance signal, and recover model parameters using expectation-maximization. We evaluate our model on both synthetic and two different experimental data sets. Stress was artificially induced in the first experimental data set and the second comprised of a fear conditioning experiment.

RESULTS

Results on the first data set indicate high levels of arousal during exposure to cognitive stress and low arousal during relaxation. Plausible results are also obtained in the fear conditioning data set consistent with previous skin conductance studies in similar experimental contexts.

CONCLUSION

The state-space approach-which does not rely on external classification labels-is able to continuously track an arousal level from skin conductance features.

SIGNIFICANCE

The method is a promising arousal estimation scheme utilizing only skin conductance. The approach could find applications in wearable monitoring and the study of neuropsychiatric conditions such as post-traumatic stress disorder.

Functional Linear and Nonlinear Brain–Heart Interplay during Emotional Video Elicitation: A Maximum Information Coefficient Study

Brain and heart continuously interact through anatomical and biochemical connections. Although several brain regions are known to be involved in the autonomic control, the functional brain–heart interplay (BHI) during emotional processing is not fully characterized yet. To this aim, we investigate BHI during emotional elicitation in healthy subjects. The functional linear and nonlinear couplings are quantified using the maximum information coefficient calculated between time-varying electroencephalography (EEG) power spectra within the canonical bands (δ,θ,α,β and γ), and time-varying low-frequency and high-frequency powers from heartbeat dynamics. Experimental data were gathered from 30 healthy volunteers whose emotions were elicited through pleasant and unpleasant high-arousing videos. Results demonstrate that functional BHI increases during videos with respect to a resting state through EEG oscillations not including the γ band (>30 Hz). Functional linear coupling seems associated with a high-arousing positive elicitation, with preferred EEG oscillations in the θ band ([4,8) Hz) especially over the left-temporal and parietal cortices. Differential functional nonlinear coupling between emotional valence seems to mainly occur through EEG oscillations in the δ,θ,α bands and sympathovagal dynamics, as well as through δ,α,β oscillations and parasympathetic activity mainly over the right hemisphere. Functional BHI through δ and α oscillations over the prefrontal region seems primarily nonlinear. This study provides novel insights on synchronous heartbeat and cortical dynamics during emotional video elicitation, also suggesting that a nonlinear analysis is needed to fully characterize functional BHI.

A robust algorithm for heart rate variability time series artefact correction using novel beat classification

Purpose: Heart rate variability is a commonly used measurement to evaluate functioning of autonomic nervous system, psychophysiological stress, and exercise intensity and recovery. HRV measurements contain artefacts such as extra, missed or misaligned beat detections, which can produce significant distortion on HRV parameters. In this paper, a robust automatic method for artefact detection from HRV time series is proposed. Methods: The proposed detection method is based on time-varying thresholds estimated from distribution of successive RR-interval differences combined with a novel beat classification scheme. The method is validated using simulated extra, missed and misaligned beat detections as well as real artefacts such as atrial and ventricular ectopic beats. Results: The sensitivity of the algorithm to detect simulated missed/extra beats was 100%. The sensitivity to detect real atrial and ventricular ectopic beats was 96.96%, the corresponding specificity being 99.94%. The mean error in HRV parameters after correction was <2% for missed and extra beats as well as for misaligned beats generated with large displacement factors. Misaligned beats with smallest displacement factor were the most difficult to detect and resulted in largest HRV parameter errors after correction, largest errors being <8%. Conclusions: The HRV artefact correction algorithm presented in this study provided comparable specificity and better sensitivity to detect ectopic beats as compared to state-of-the-art algorithms. The proposed algorithm detects abnormal beats with high accuracy, is relatively easy to implement, and secures reliable HRV analysis by reducing the effect of possible artefacts to tolerable level.

Central modulation of parasympathetic outflow is impaired in de novo Parkinson's disease patients

Task- and stimulus-based neuroimaging studies have begun to unveil the central autonomic network which modulates autonomic nervous system activity. In the present study, we aimed to evaluate the central autonomic network without the bias constituted by the use of a task. Additionally, we assessed whether this circuitry presents signs of dysregulation in the early stages of Parkinson’s disease (PD), a condition which may be associated with dysautonomia. We combined heart-rate-variability based methods for time-varying assessments of the autonomic nervous system outflow with resting-state fMRI in 14 healthy controls and 14 de novo PD patients, evaluating the correlations between fMRI time-series and the instantaneous high-frequency component of the heart-rate-variability power spectrum, a marker of parasympathetic outflow. In control subjects, the high-frequency component of the heart-rate-variability power spectrum was significantly anti-correlated with fMRI time-series in several cortical, subcortical and brainstem regions. This complex central network was not detectable in PD patients. In between-group analysis, we found that in healthy controls the brain activation related to the high-frequency component of the heart-rate-variability power spectrum was significantly less than in PD patients in the mid and anterior cingulum, sensorimotor cortex and supplementary motor area, insula and temporal lobe, prefrontal cortex, hippocampus and in a region encompassing posterior cingulum, precuneus and parieto-occipital cortex. Our results indicate that the complex central network which modulates parasympathetic outflow in the resting state is impaired in the early clinical stages of PD.

Resting-state brain correlates of instantaneous autonomic outflow

A prominent pathway of brain-heart interaction is represented by autonomic nervous system (ANS) heartbeat modulation. While within-brain resting state networks have been the object of intense functional Magnetic Resonance Imaging (fMRI) research, technological and methodological limitations have hampered research on the central correlates of cardiovascular control dynamics. Here we combine the high temporal and spatial resolution as well as data volume afforded by the Human Connectome Project with a probabilistic model of heartbeat dynamics to characterize central correlates of sympathetic and parasympathetic ANS activity at rest. We demonstrate an involvement of a number of brain regions such as the Insular cortex, Frontal Gyrus, Lateral Occipital Cortex, Paracingulate and Cingulate Gyrus and Precuneous Cortex, as well as subcortical structures (Thalamus, Putamen, Pallidum, Brain-Stem, Hippocampus, Amygdala, and Right Caudate) in the modulation of ANS-mediated cardiovascular control, possibly indicating a broader definition of the central autonomic network (CAN). Our findings provide a basis for an informed neurobiological interpretation of the numerous studies which employ HRV-related measures as standalone biomarkers in health and disease.

Heart Rate Variability During Exercise: A Comparison of Artefact Correction Methods

Giles, DA and Draper, N. Heart rate variability during exercise: a comparison of artefact correction methods. J Strength Cond Res 32(3): 726-735, 2018-There is a need for standard practice in the collection and processing of RR interval data recorded using heart rate monitors (HRMs) in research. This article assessed the validity of RR intervals and heart rate variability (HRV) data obtained using an HRM during incremental exercise and artefact correction methods. Eighteen participants completed an active orthostatic test and incremental running V[Combining Dot Above]O2max test, while simultaneous recordings using a Polar V800 HRM and an electrocardiogram were made. Artefacts were corrected by deletion; degree zero, linear, cubic, and spline interpolation; and Kubios HRV software. Agreement was assessed using percentage bias, effect size (ES), intraclass correlation coefficients (ICC), and Bland-Altman limits of agreement (LoA). Artefacts increased relative to exercise intensity, to a peak of 4.46% during 80-100% V[Combining Dot Above]O2max. Correction of RR intervals was necessary with unacceptably increased bias, LoA, and ES and reduced ICC in all but resting recordings. All correction methods resulted in data with reduced percentage bias and ES for resting and <60% V[Combining Dot Above]O2max recordings. However, at >60% V[Combining Dot Above]O2max, even with correction, large amounts of variation were present in HRV measures of root mean square of the successive difference of intervals, low-to-high frequency ratio, Poincaré dispersion perpendicular to the axis (SD1), and sample entropy. Linear interpolation produced RR intervals with the lowest bias and ES. However, caution should be given to HRV parameters at high exercise intensities, as large amounts of variation were still present. Recommendations for minimizing artefacts are discussed, along with guidelines for their identification, correction, and reporting.

The Role of Baroreflex Sensitivity in Acute Hypotensive Episodes Prediction in the Intensive Care Unit

A life threatening condition in Intensive Care Unit (ICU) is the Acute Hypotensive Episode (AHE). Patients experiencing an AHE may suffer from irreversible organ damage associated with increased mortality. Predicting the onset of AHE could be of pivotal importance to establish appropriate and timely interventions. We propose a method that, using waveforms widely acquired in ICU, like Arterial Blood Pressure (ABP) and Electrocardiogram (ECG), will extract features relative to the cardiac system to predict whether or not a patient will experience a hypotensive episode. Specifically, we want to assess if there are hidden patterns in the dynamics of baroreflex able to improve the prediction of AHEs. We will investigate the predictive power of features related to the baroreflex by performing classifications with and without them. Results are obtained using 17 classifiers belonging to different model families: classification trees, Support Vector Machines (SVMs), K-Nearest Neighbors (KNNs) replicated with different set of hyper-parameters and logistic regression. On average, the use of baroreflex features in the AHE prediction process increases the Area Under the Curve (AUC) by 10%.

Optimizing the learning rate for adaptive estimation of neural encoding models

Closed-loop neurotechnologies often need to adaptively learn an encoding model that relates the neural activity to the brain state, and is used for brain state decoding. The speed and accuracy of adaptive learning algorithms are critically affected by the learning rate, which dictates how fast model parameters are updated based on new observations. Despite the importance of the learning rate, currently an analytical approach for its selection is largely lacking and existing signal processing methods vastly tune it empirically or heuristically. Here, we develop a novel analytical calibration algorithm for optimal selection of the learning rate in adaptive Bayesian filters. We formulate the problem through a fundamental trade-off that learning rate introduces between the steady-state error and the convergence time of the estimated model parameters. We derive explicit functions that predict the effect of learning rate on error and convergence time. Using these functions, our calibration algorithm can keep the steady-state parameter error covariance smaller than a desired upper-bound while minimizing the convergence time, or keep the convergence time faster than a desired value while minimizing the error. We derive the algorithm both for discrete-valued spikes modeled as point processes nonlinearly dependent on the brain state, and for continuous-valued neural recordings modeled as Gaussian processes linearly dependent on the brain state. Using extensive closed-loop simulations, we show that the analytical solution of the calibration algorithm accurately predicts the effect of learning rate on parameter error and convergence time. Moreover, the calibration algorithm allows for fast and accurate learning of the encoding model and for fast convergence of decoding to accurate performance. Finally, larger learning rates result in inaccurate encoding models and decoders, and smaller learning rates delay their convergence. The calibration algorithm provides a novel analytical approach to predictably achieve a desired level of error and convergence time in adaptive learning, with application to closed-loop neurotechnologies and other signal processing domains.

Measures of sympathetic and parasympathetic autonomic outflow from heartbeat dynamics

Reliable and effective noninvasive measures of sympathetic and parasympathetic peripheral outflow are of crucial importance in cardiovascular physiology. Although many techniques have been proposed to take up this long-lasting challenge, none has proposed a satisfying discrimination of the dynamics of the two separate branches. Spectral analysis of heart rate variability is the most currently used technique for such assessment. Despite its widespread use, it has been demonstrated that the subdivision in the low-frequency (LF) and high-frequency (HF) bands does not fully reflect separate influences of the sympathetic and parasympathetic branches, respectively, mainly due to their simultaneous action in the LF. Two novel heartbeat-derived autonomic measures, the sympathetic activity index (SAI) and parasympathetic activity index (PAI), are proposed to separately assess the time-varying autonomic nervous system synergic functions. Their efficacy is validated in landmark autonomic maneuvers generally employed in clinical settings. The novel measures move beyond the classical frequency domain paradigm through identification of a set of coefficients associated with a proper combination of Laguerre base functions. The resulting measures were compared with the traditional LF and HF power. A total of 236 ECG recordings were analyzed for validation, including autonomic outflow changes elicited by procedures of different nature and temporal variation, such as postural changes, lower body negative pressure, and handgrip tests. The proposed SAI-PAI measures consistently outperform traditional frequency-domain indexes in tracking expected instantaneous autonomic variations, both vagal and sympathetic, and may aid clinical decision making, showing reduced intersubject variability and physiologically plausible dynamics. NEW & NOTEWORTHY While it is possible to obtain reliable estimates of parasympathetic activity from the ECG, a satisfying method to disentangle the sympathetic component from HRV has not been proposed yet. To overcome this long-lasting limitation, we propose two novel HRV-based indexes, the sympathetic and parasympathetic activity indexes.

Mortality Prediction in Severe Congestive Heart Failure Patients with Multifractal Point-Process Modeling of Heartbeat Dynamics

Multifractal analysis of human heartbeat dynamics has been demonstrated to provide promising markers of Congestive Heart Failure (CHF). Yet, it crucially builds on the interpolation of RR intervals series, which has been generically performed with limited links to CHF pathophysiology. We devise a novel methodology estimating multifractal autonomic dynamics from heartbeat-derived series defined in the continuous time. We hypothesize that markers estimated from our novel framework are also effective for mortality prediction in severe CHF. We merge multifractal analysis within a methodological framework based on inhomogeneous point process models of heartbeat dynamics. Specifically, wavelet coefficients and wavelet leaders are computed over measures extracted from instantaneous statistics of probability density functions characterizing and predicting the time until the next heartbeat event occurs. The proposed approach is tested on data from 94 CHF patients, aiming at predicting survivor and non-survivor individuals as determined after a 4 years follow up. Instantaneous markers of vagal and sympatho-vagal dynamics display power-law scaling for a large range of scales, from s to s. Using standard SVM algorithms, the proposed inhomogeneous point-process representation based multifractal analysis achieved the best CHF mortality prediction accuracy of 79.11 % (sensitivity 90.48%, specificity 67.74%). Our results suggest that heartbeat scaling and multifractal properties in CHF patients are not generated at the sinus-node level, but rather by the intrinsic action of vagal short-term control and of sympatho-vagal fluctuations associated with circadian cardiovascular control, especially within the VLF band. These markers might provide critical information in devising a clinical tool for individualized prediction of survivor and non-survivor CHF patients.

Assessing mood symptoms through heartbeat dynamics: An HRV study on cardiosurgical patients

BACKGROUND

Heart Rate Variability (HRV) is reduced both in depression and in coronary heart disease (CHD) suggesting common pathophysiological mechanisms for the two disorders. Within CHD, cardiac surgery patients (CSP) with postoperative depression are at greater risk of adverse cardiac events. Therefore, CSP would especially benefit from depression early diagnosis. Here we tested whether HRV-multi-feature analysis discriminates CSP with or without depression and provides an effective estimation of symptoms severity.

METHODS

Thirty-one patients admitted to cardiac rehabilitation after first-time cardiac surgery were recruited. Depressive symptoms were assessed with the Center for Epidemiologic Studies Depression Scale (CES-D). HRV features in time, frequency, and nonlinear domains were extracted from 5-min-ECG recordings at rest and used as predictors of "least absolute shrinkage and selection" (LASSO) operator regression model to estimate patients' CES-D score and to predict depressive state.

RESULTS

The model significantly predicted the CES-D score in all subjects (the total explained variance of CES-D score was 89.93%). Also it discriminated depressed and non-depressed CSP with 86.75% accuracy. Seven of the ten most informative metrics belonged to non-linear-domain.

LIMITATIONS

A higher number of patients evaluated also with a structured clinical interview would help to generalize the present findings.

DISCUSSION

To our knowledge this is the first study using a multi-feature approach to evaluate depression in CSP. The high informative power of HRV-nonlinear metrics suggests their possible pathophysiological role both in depression and in CHD. The high-accuracy of the algorithm at single-subject level opens to its translational use as screening tool in clinical practice.

Nonlinear analysis of heart rate variability for the assessment of Dysphoria

Dysphoric patients show symptoms associated with Major Depression, although within a narrowed symptomatology spectrum. In prevailing practice, clinicians assess Dysphoria through psychological scores and questionnaires exclusively, therefore without taking into account objective biomarkers. In this study, we investigated heartbeat linear and nonlinear dynamics aiming to an objective assessment of Dysphoria. To this end, we derived standard and nonlinear measures from heart rate variability (HRV) series gathered from dysphoric (n=14) and nondysphoric (n=17) undergraduate students during 5 minutes of resting state. We performed both statistical and pattern recognition analyses in order to discern the two groups. Results showed significant group-wise differences in HRV nonlinear metrics exclusively, suggesting a crucial role of nonlinear sympatho-vagal dynamics in Dysphoria. Furthermore, we achieved a classification accuracy of 77.52% for the automatic identification of Dysphoria at a single-subject level.