Poincaré plot indexes of heart rate variability: Relationships with other nonlinear variables

Heart rate variability helps classify phenotype in systemic sclerosis

We aimed to develop a systemic sclerosis (SSc) subtypes classifier tool to be used at the patient's bedside. We compared the heart rate variability (HRV) at rest (5-min) and in response to orthostatism (5-min) of patients (n = 58) having diffuse (n = 16, dcSSc) and limited (n = 38, lcSSc) cutaneous forms. The HRV was evaluated from the beat-to-beat RR intervals in time-, frequency-, and nonlinear-domains. The dcSSc group differed from the lcSSc group mainly by a higher heart rate (HR) and a lower HRV, in decubitus and orthostatism conditions. Stand-up maneuver lowered HR standard deviation (sd_HR), the major axis length of the fitted ellipse of Poincaré plot of RR intervals (SD2), and the correlation dimension (CorDim) in the dcSSc group while increased these HRV indexes in the lcSSc group (p = 0.004, p = 0.002, and p = 0.004, respectively). We identified the 5 most informative and discriminant HRV variables. We then compared 341 classifying models (1 to 5 variables combinations × 11 classifier algorithms) according to mean squared error, logloss, sensitivity, specificity, precision, accuracy, area under curve of the ROC-curves and F1-score. F1-score ranged from 0.823 for the best 1-variable model to a maximum of 0.947 for the 4-variables best model. Most specific and precise models included sd_HR, SD2, and CorDim. In conclusion, we provided high performance classifying models able to distinguish diffuse from limited cutaneous SSc subtypes easy to perform at the bedside from ECG recording. Models were based on 1 to 5 HRV indexes used as nonlinear markers of autonomic integrated influences on cardiac activity.

Poincaré analysis detects pathological limb loading rate variability in post-anterior cruciate ligament reconstruction individuals

BACKGROUND

Post-ACLR individuals can experience repeated exposure to variable limb loading, which contributes to development of knee osteoarthritis. Variable limb loading can present as loading rate variability (LRV) and is magnified during tasks like fast walking when the system is stressed. Nonlinear measures that evaluate temporal variability have successfully detected changes in gait variability associated with altered motor control, however, appropriately describing and uncovering the nature of gait variability has been challenging. Here, Poincaré analysis, a nonlinear method unique in its ability to capture different aspects of variability, served to uncover and quantify changes in limb LRV. It was hypothesized that post-ACLR individuals' overloaded limbs would quantitatively and graphically demonstrate greater short-term stride-to-stride and long-term limb LRV during fast walking compared to the underloaded and healthy control limbs.

METHODS

Fourteen post-ACLR individuals and fourteen healthy controls completed a walking protocol on an instrumented treadmill where they walked at 1.0 m/s and 1.5 m/s for 5-minutes each. A Welch's test was performed to compare differences in short-term and long-term LRV metrics for the post-ACLR individuals' overloaded and underloaded limbs and the healthy controls' right limbs.

RESULTS

Analyses revealed that the post-ACLR individuals' overloaded limb exhibited significantly greater short-term and long-term values compared to the underloaded and healthy control limbs at 1.5 m/s (p<0.05). Additionally, the loading rate data was widely scattered across the plots for the overloaded limb, indicating greater LRV.

SIGNIFICANCE

Poincaré analysis successfully identified that post-ACLR overloaded limbs exhibited impaired motor control during fast walking based on quantitative and graphical changes in variability. This highlights the clinical applications of Poincaré analysis, with the plots potentially serving as an easy-to-interpret diagnostic tool for pathological limb LRV.

Automated extraction of heart rate variability from magnetoencephalography signals

Magnetoencephalography (MEG) measures magnetic fluctuations in the brain generated by neural processes, some of which, such as cardiac signals, are generally removed as artifacts and discarded. However, heart rate variability (HRV) has long been regarded as a biomarker related to autonomic function, suggesting the cardiac signal in MEG contains valuable information that can provide supplemental health information about a patient. To enable access to these ancillary HRV data, we created an automated extraction tool capable of capturing HRV directly from raw MEG data with artificial intelligence. Five scans were conducted with simultaneous MEG and electrocardiogram (ECG) acquisition, which provides a ground truth metric for assessing our algorithms and data processing pipeline. In addition to directly comparing R-peaks between the MEG and ECG signals, this work explores the variation of the corresponding HRV output in time, frequency, and non-linear domains. After removing outlier intervals and aligning the ECG and derived cardiac MEG signals, the RMSE between the RR-intervals of each was RMSE1 = 2 ms, RMSE2 = 2 ms, RMSE3 = 8 ms, RMSE4 = 4 ms, RMSE5 = 13 ms. The findings indicate that cardiac artifacts from MEG data carry sufficient signal to approximate an individual's HRV metrics.

Post-exercise heart rate recovery and its speed are associated with resting-reactivity cardiovagal modulation in healthy women

The present study sought to expand upon prior investigations of the relationship between post-exercise heart rate recovery (HRR) and cardiovagal resting-reactivity modulation. HRR from 1st to 5th min after maximal exercise test was correlated with a cardiovagal index of heart rate variability (SD1) at resting (supine and orthostatic positions) and its reactivity after the orthostatic stress test in 34 healthy women. Statistical analysis employed non-parametric tests with a p-value set at 5%. HRR, ∆%HRR, and coefficient of HRR (CHRR) at the 3rd and 5th min correlated with SD1 and SD1n (normalized units) in the supine position (rs = 0.36 to 0.47; p =  < 0.01). From the 1st to 5th min, HRR, ∆%HRR, and CHRR correlated with SD1 and SD1n in the orthostatic position (rs = 0.29 to 0.47; p =  ≤ 0.01 to 0.05), except for HRR at 5th min with SD1n (p = 0.06). Following the orthostatic stress test, HRR at 3rd and HRR, %∆HRR at 5th min correlated with ∆absSD1 (rs = 0.28 to 0.35; p = 0.02 to 0.05). All HRR measurements at 1st min correlated with ∆absSD1n (rs = 0.32 to 0.38; p = 0.01 to 0.03), and the CHRR at 1st min correlated with ∆%SD1(rs = 0.37; p = 0.01). After the sample was divided into high and low cardiovagal modulation subgroups, the subgroup with high modulation at rest (supine and orthostatic) and higher cardiovagal reactivity (reduction) showed faster HRR (p =  ≤ 0.01 to 0.05; ES:0.37 to 0.50). HRR throughout the 1st to 5th min positively correlates with cardiovagal modulation in the orthostatic position, and the 3rd and 5th min positively correlate with cardiovagal modulation in both postures at rest. Faster HRR following the maximal exercise test is associated with high resting-reactivity cardiovagal modulation in healthy women.

Feasibility of wearable monitors to detect heart rate variability in children with hand, foot and mouth disease

Hand foot and mouth disease (HFMD) is caused by a variety of enteroviruses, and occurs in large outbreaks in which a small proportion of children deteriorate rapidly with cardiopulmonary failure. Determining which children are likely to deteriorate is difficult and health systems may become overloaded during outbreaks as many children require hospitalization for monitoring. Heart rate variability (HRV) may help distinguish those with more severe diseases but requires simple scalable methods to collect ECG data.We carried out a prospective observational study to examine the feasibility of using wearable devices to measure HRV in 142 children admitted with HFMD at a children's hospital in Vietnam. ECG data were collected in all children. HRV indices calculated were lower in those with enterovirus A71 associated HFMD compared to those with other viral pathogens.HRV analysis collected from wearable devices is feasible in a low and middle income country (LMIC) and may help classify disease severity in HFMD.

Short-Term Heart Rate Variability in Metabolic Syndrome: A Systematic Review and Meta-Analysis

BACKGROUND

Our aim was to determine the differences in short-term heart rate variability (HRV) between patients with metabolic syndrome (MS) and healthy controls.

METHODS

We searched electronic databases for primary works with short-term HRV recordings (≤30 min) that made comparisons between individuals with MS versus healthy controls. This systematic review and meta-analysis (MA) was performed according to PRISMA guidelines and registered at PROSPERO (CRD42022358975).

RESULTS

Twenty-eight articles were included in the qualitative synthesis and nineteen met the criteria for the MA. Patients with MS showed decreased SDNN (-0.36 [-0.44, -0.28], p < 0.001), rMSSD (-7.59 [-9.98, -5.19], p < 0.001), HF (-0.36 [-0.51, -0.20], p < 0.00001) and LF (-0.24 [-0.38, -0.1], p = 0.001). In subsequent subanalyses, we found a decrease in SDNN (-0.99 (-1.45, -0.52], p < 0.001), rMSSD (-10.18 [-16.85, -3.52], p < 0.01) and HF (-1.04 [-1.97, -0.1] p < 0.05) in women. In men, only LF showed a significant lower value (-0.26 [-0.5, -0.02], p < 0.05). We could not perform MA for non-linear variables.

CONCLUSIONS

Patients with MS showed changes in time-domain analyses, with lower values in SDNN and rMSSD. Regarding frequency-domain analyses, MS patients showed a decrease in HF and LF When sex was used as a grouping variable, the MA was only possible in one of both sexes (men or women) in rMSSD and LF/HF. Lastly, when data for both men and women were available, subanalyses showed a different behavior compared to mixed analyses for SDNN, HF and LF, which might point towards a different impact of MS in men and women.

Prediction of post-operative atrial fibrillation in patients after cardiac surgery using heart rate variability

PURPOSE

Post-operative atrial fibrillation (PoAF) occurs in ~ 30% of patients after cardiac surgery. The etiology of PoAF is complex, but a disbalance in autonomic systems plays an important role. The goal of this study was to assess whether pre-operative heart rate variability analysis can predict the risk of PoAF.

METHODS

Patients without a history of AF with an indication for cardiac surgery were included. Two-hour ECG recordings one day before surgery was used for the HRV analysis. Univariate and multivariate logistic regression, including all HRV parameters, their combination, and clinical variables, were calculated to find the best predictive model for post-operative AF.

RESULTS

One hundred and thirty-seven patients (33 women) were enrolled in the study. PoAF occurred in 48 patients (35%, AF group); the remaining 89 patients were in the NoAF group. AF patients were significantly older (69.1 ± 8.6 vs. 63.4 ± 10.5 yrs., p = 0.002), and had higher CHA₂DS₂-VASc score (3 ± 1.4 vs. 2.5 ± 1.3, p = 0.01). In the multivariate regression model, parameters independently associated with higher risk of AF were pNN50, TINN, absolute power VLF, LF and HF, total power, SD2, and the Porta index. A combination of clinical variables with HRV parameters in the ROC analysis achieved an AUC of 0.86, a sensitivity of 0.95, and a specificity of 0.57 and was more effective in PoAF prediction than a combination of clinical variables alone.

CONCLUSION

A combination of several HRV parameters is helpful in predicting the risk of PoAF. Attenuation of heart rate variability increases the risk for PoAF.

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.

DHA Supplementation During Pregnancy Enhances Maternal Vagally Mediated Cardiac Autonomic Control in Humans

BACKGROUND

DHA is an essential omega-3 (ω-3; n-3) fatty acid that has well-established benefits for the fetus. DHA also has the potential to influence the health of the mother, but this area is understudied.

OBJECTIVES

The objective of this secondary analysis was to determine if DHA was related to maternal heart rate (HR) and heart rate variability (HRV) metrics in a large cohort of pregnant women.

METHODS

In the parent trial (1R01HD086001) eligible participants (≥18 y old, English speaking, carrying a singleton pregnancy, 12-20 wk of gestation) were randomly assigned to consume 200 mg/d or 800 mg/d DHA for the duration of their pregnancy (n = 300). Weight, blood pressure, and magnetocardiograms (MCGs) were collected at 32 wk and 36 wk of gestation (n = 221). Measures of HR and HRV in time-, frequency-, and nonlinear-domains were determined from the isolated maternal MCG. Treatment group and timepoint were examined as predictors in association with HR and HRV metrics using random-intercept mixed-effects ANOVA unadjusted and adjusted models accounting for weight and dietary DHA intake.

RESULTS

Women receiving the higher dose of DHA (800 mg/d) during pregnancy had lower HR, lower sympathetic index, higher vagally mediated HRV indices, and greater HRV complexity when compared with the women who received the lower dose (200 mg/d; all P < 0.05). All the dose relations remained significant even after controlling for the effect of time, maternal weight, and dietary DHA intake.

CONCLUSIONS

DHA supplementation increases vagal tone in pregnant women. Longitudinal studies examining the potential link between DHA, enhanced vagal tone, and reported reduction in early preterm birth are warranted.

Correlation among Poincare plot and traditional heart rate variability indices in adults with different risk levels of metabolic syndrome: a cross-sectional approach from Southern India

OBJECTIVES

Heart rate variability (HRV) is an important marker of cardiac autonomic modulation. Metabolic syndrome (MetS) can alter cardiac autonomic modulation, raising the risk of cardiovascular disease (CVD). Poincaré plot analysis (PPA) is a robust scatter plot-based depiction of HRV and carries similar information to the traditional HRV measures. However, no prior studies have examined the relationship between PPA and traditional HRV measures among different risk levels of MetS. We evaluated the association between the Poincare plot and traditional heart rate variability indices among adults with different risk levels of MetS.

METHODS

We measured anthropometric data and collected fasting blood samples to diagnose MetS. The MetS risk was assessed in 223 participants based on the number of MetS components and was classified as control (n=64), pre-MetS (n=49), MetS (n=56), and severe MetS (n=54). We calculated the Poincaré plot (PP) and traditional HRV measures from a 5 min HRV recording.

RESULTS

Besides the traditional HRV measures, we found that various HRV indices of PPA showed significant differences among the groups. The severe MetS group had significantly lower S (total HRV), SD1 (short-term HRV), SD2 (long-term HRV), and higher SD2/SD1. The values of S, SD1, SD2, and SD2/SD1 were significantly correlated with most traditional HRV measures.

CONCLUSIONS

We found gradual changes in HRV patterns as lower parasympathetic and higher sympathetic activity alongside the rising number of MetS components. The HRV indices of PPA integrating the benefits of traditional HRV indices distinguish successfully between different risk levels of MetS and control subjects.

Cardiac autonomic modulation in response to postural transition during a virtual reality task in individuals with spinal cord injury: A cross-sectional study

PURPOSE

The postural transition from sitting to standing is a moment of dysautonomic occurrence in individuals with Spinal Cord Injury (SCI). Different tools can be used to minimize this event, such as virtual reality. Thus, we aimed to analyze cardiac autonomic modulation in individuals with SCI during postural transition from the sitting to orthostatism position using a cognitive virtual reality (VR) task.

METHODS

Individuals with and without SCI were positioned on the Easy Stand® device, sitting at rest, at 0° considering the angle between the seat and the floor, elevation at 45°, and orthostatism at 90°, for 5 minutes in each position. Heart rate variability (HRV) measures of sympathovagal balance were collected (heart rate receiver: Polar V800). The groups were subdivided into two groups, one that performed VR as an intervention during the postural angle changes and another group that did not perform VR.

RESULTS

We evaluated 76 individuals, 40 with a medical diagnosis of SCI and 36 who composed the able-bodied control group without SCI, matched by age and sex. The HRV results showed that the SCI group who performed the task in VR demonstrated no significant difference in parasympathetic activation and global variability between the sitting versus 90° positions. There was better sympathovagal balance in SCI and able-bodied control groups who performed the VR task between the sitting versus 90° positions.

CONCLUSION

The use of a VR task seems to contribute to better sympathovagal balance, with the potential to reduce dysautonomia during postural changes.

Autonomic nervous system balance in parturient mares: Spontaneous vs induced delivery

Delivery is not easily predictable in horses and the consequences of dystocia can be serious for both the mare and foal. An induction protocol with low doses of oxytocin has been reported as a safe procedure. This study investigates the effect of induced delivery on at-term mares' sympathetic-vagal balance. Fourteen mares were included and divided into two groups, one subjected to spontaneous delivery (SD), and one to induced delivery (ID). In both groups, an ECG was recorded using an elastic belt with integrated smart textile electrodes. The recording started before the delivery (Basal), continued close to delivery (Pre-delivery) and during delivery (Delivery), and ended after parturition (Placental expulsion). From the ECGs, Heart Rate Variability (HRV) parameters relating to time and frequency domains and non-linear analysis were extrapolated. The HRV analysis was performed both within the same group (IntraGA) and between the two groups (InterGA). In the present study, spontaneous and induced delivery did not appear to differ in autonomic nervous system functioning. In IntraGA analysis, both for SD and ID mares, delivery and placental expulsion periods were parasympathetic dominated since vagal-related HRV parameters increased. Moreover, no differences were found in InterGA comparison between SD and ID mares, except for the pre-delivery period of ID mares, during which both branches of the autonomic nervous system were activated. These results are in line with the literature on parasympathetic dominance during parturition and no change in Heart Rate Variability following exogenous oxytocin administration in parturient mares.

Neuro-Cardio-Autonomic Modulations in Children with Duchenne Muscular Dystrophy

BACKGROUND AND OBJECTIVE

Duchenne muscular dystrophy (DMD) is a degenerative X-linked muscle disease. Death frequently results from complications in cardiopulmonary systems. Preclinical/early diagnosis of cardiac autonomic abnormalities may aid initiate cardioprotective therapy and enhance prognosis.

METHODS

A cross sectional, prospective study of 38 DMD boys compared with 37 age-matched healthy controls was conducted. Lead II electrocardiography and beat-to-beat blood pressure were recorded to assess heart rate variability (HRV), blood pressure variability (BPV), and baroreceptor sensitivity (BRS) in a standardized environment. Data were analysed and correlated with disease severity and genotype.

RESULTS

In the DMD group, the median age at assessment was 8 years [IQR 7-9 years], the median age at disease onset was 3 years [IQR, 2-6 years], and the mean duration of illness was 4 years [IQR, 2.5-5]. DNA sequencing showed deletions in 34/38 (89.5 %) and duplications in 4/38 (10.5%) patients. The median heart rate in DMD children was significantly higher [101.19 (Range, 94.71-108.49)] /min compared to controls [81 (Range, 76.2-92.76)] /min (p < 0.05). All the assessed HRV and BPV parameters were significantly impaired in DMD cases except for the coefficient of variance of systolic blood pressure. Further, BRS parameters were also significantly reduced in DMD, excluding alpha-LF. A positive correlation was found between alpha HF with age at onset and duration of illness.

CONCLUSION

This study demonstrates a distinct early impairment of neuro-cardio-autonomic regulation in DMD. Simple yet effective non-invasive techniques such as HRV, BPV, and BRS may help identify cardiac dysfunction in a pre-clinical state, paving the way for early cardio-protective therapies and limiting disease progression in DMD patients.

Central and Peripheral Fatigue in Recreational Trail Runners: A Pilot Study

BACKGROUND

Understanding fatigue mechanisms is crucial for exercise performance. However, scientific evidence on non-invasive methods for assessing fatigue in trail running competitions is scarce, especially when vertical kilometer trail running races (VK) are considered. The main purpose of this study was to assess the autonomic nervous system (ANS) activity (i.e., central fatigue) and the state of muscle activation (i.e., peripheral fatigue) before and after a VK competition.

METHODS

A cross-sectional pilot study was performed. After applying inclusion/exclusion criteria, 8 recreational male trail runners (31.63 ± 7.21 yrs, 1.75 m ± 0.05 m, 70.38 ± 5.41 kg, BMI: 22.88 ± 0.48, running experience: 8.0 ± 3.63 yrs, weekly training volume: 58.75 ± 10.35 km) volunteered to participate and were assessed for both central (i.e., via heart rate variability, HRV) and peripheral (via tensiomyography, TMG) fatigue before and after a VK race.

RESULTS

After the VK, resting heart rate, RMSSD (p = 0.01 for both) and SDNN significantly decreased (p = 0.02), while the stress score and the sympathetic-parasympathetic ratio increased (p = 0.01 and p = 0.02, respectively). The TMG analyses suggest that runners already suffered peripheral fatigue before the VK and that 20-30 min are enough for muscular recovery after the race. In summary, our data suggest that participants experienced a pre-competition fatigue status. Further longitudinal studies are necessary to investigate the mechanisms underlying fatigue during trail running races, while training periodization and tapering strategies could play a key role for minimizing pre-competition fatigue status.

Physiological effects of Therapeutic Body Wraps in healthy volunteers: An observational study

INTRODUCTION

The clinical management of severe anxiety is challenging. Along with specific medication, several nonpharmacological treatments exist, among which the Therapeutic Body Wraps (TBW). While TBW is clinically known to be efficient, the absence of objective physiological measurements raised some debates about its harmlessness.

AIMS

This observational study investigated TBW in healthy volunteers.

METHODS

Measurements were done in 26 participants, who experienced TBW, which consisted of tightly wrapping the body in one layer of wet cold sheets and several warm blankets. These were compared to 13 participants, who rested in supine position.

RESULTS

Our results showed increase in the trunk skin temperature by the end of TBW similar to rest. Somatosensory perception as assessed with quantitative sensory testing was stable after both TBW and rest. The heart and breath rates decreased both during TBW and rest. It was accompanied by increase in heart rate variability parameters and decrease in salivary cortisol levels.

DISCUSSION

Our data indicate anxiolytic effect of TBW in healthy participants.

IMPLICATIONS FOR PRACTICE

The TBW is harmless and might be used in patients, who are unable to relax by themselves and/or without anxiolytic medication. Further studies are necessary to investigate physiological response to TBW in clinical population.

Bedside tracking of functional autonomic age in preterm infants

BACKGROUND

Preterm birth predisposes infants to adverse outcomes that, without early intervention, impacts their long-term health. To assist bedside monitoring, we developed a tool to track the autonomic maturation of the preterm by assessing heart rate variability (HRV) changes during intensive care.

METHODS

Electrocardiogram (ECG) recordings were longitudinally recorded in 67 infants (26-38 weeks postmenstrual age (PMA)). Supervised machine learning was used to generate a functional autonomic age (FAA), by combining 50 computed HRV features from successive 5-minute ECG epochs (median of 23 epochs per infant). Performance of the FAA was assessed by correlation to PMA, clinical outcomes and the infant's functional brain age (FBA), an index of maturation derived from the electroencephalogram.

RESULTS

The FAA was strongly correlated to PMA (r = 0.86, 95% CI: 0.83-0.93) with a mean absolute error (MAE) of 1.66 weeks and also accurately estimated FBA (MAE = 1.58 weeks, n = 54 infants). The relationship between PMA and FAA was not confounded by neurodevelopmental outcome (p = 0.18, n = 45), sex (p = 0.88, n = 56), patent ductus arteriosus (p = 0.08, n = 56), IVH (p = 0.63, n = 56) or body weight at birth (p = 0.95, n = 56).

CONCLUSIONS

The FAA, an index derived from the ubiquitous ECG signal, offers direct avenues towards estimating autonomic maturation at the bedside during intensive care monitoring.

IMPACT

The development of a tool to track functional autonomic age in preterm infants based on heart rate variability features in the electrocardiogram provides a rapid and specialized view of autonomic maturation at the bedside. Functional autonomic age is linked closely to postmenstrual age and central nervous system function response, as determined by its relationship to functional brain age from the electroencephalogram. Tracking functional autonomic age during neonatal intensive care unit monitoring offers a unique insight into cardiovascular health in infants born extremely preterm and their maturational trajectories to term age.

The Impact on the Stress-Associated Autonomic Response of Physiotherapy Students Receiving Interferential Current in an Electrotherapy Training Session

Electrical currents are didactic contents widely applied in the training of physiotherapy students, but the treatment is considered a stressful situation for both the patient who receives it and the student who applies it. The aim of this study was to evaluate the stress-associated autonomic response of physiotherapy students receiving interferential current by measuring and analysing heart rate variability. An observational case-control study was conducted. Ninety healthy male volunteers, all physiotherapy degree students, were enrolled while attending laboratory practice during the 2020-2021 academic year. Participants were randomly allocated to a sham electrotherapy group (44 subjects), in which heart rate variability was recorded for 10 min, both at rest and during the application of sham technique on the lower back (10 min), and an electrotherapy group (46 subjects), applying the same procedure with the electrical current flowing. Outcome measures included baseline (seated position) and postintervention (prone position) time domain parameter, diameters of the Poincaré plot 1 and 2, stress score, and sympathetic/parasympathetic ratio. The sham electrotherapy group exhibited significant increases in time domain parameter (p = 0.027) and diameters of the Poincaré plot 1 (p = 0.032), with a small effect size (d ≤ 0.5). The electrotherapy group exhibited significant increases in time domain parameter and diameters of the Poincaré plot 1 and 2 (p < 0.001) and decreases in the stress score and sympathetic/parasympathetic ratio (p < 0.001), with a large effect size (d > 0.8) other than for the time domain parameter (d = 0.42), indicating increased parasympathetic and decreased sympathetic activity. After interventions, there were significant differences between groups in diameters of the Poincaré plot 2 (p < 0.001), stress score (p = 0.01) and sympathetic/parasympathetic ratio (p = 0.003), with moderate effect size (d > 0.5). The application of the interferential current technique produces stress-associated autonomic response characterized by greater parasympathetic activity and decreased sympathetic activity. Further studies are needed to determine possible adverse effects.

Pre-competitive anxiety and autonomic responses in professional U-20 futsal players: Effect of the competition phase and game location

The study examined whether pre-competition anxiety and autonomic responses in elite futsal players change across the playoffs at away vs home venues. Heart rate variability, somatic and cognitive anxiety and self-confidence (by CSAI-2R questionnaire) were evaluated in nine male futsal players from a professional U-20 team, before competitive matches at the quarter-final, semi-final and final stages played at home and away venues. Two-way ANOVA for repeated measures was used, considering factor one the game location and factor two the playoff stage, with significance set at p<0.05. Significant effects of the playoff stage were demonstrated in the Mean of HR (F = 4.643; p = 0.014) and SD2 index (F = 14.83; p=<0.001)(quarter versus final). No difference was found for somatic and cognitive anxiety and self-confidence between the two factors. The results suggest that play at the final stage of the playoff, regardless of the game location, may cause higher physiological stress, demonstrated by elevated HR and decrease of SD2 index.

Using DeepLabCut as a Real-Time and Markerless Tool for Cardiac Physiology Assessment in Zebrafish

Simple Summary

With the advancement of existing technology, artificial intelligence is widely applied in various fields of research, including cardiovascular studies. In this study, we explored the feasibility of conducting a markerless cardiac physiology assessment in zebrafish embryos by using DeepLabCut (DLC), a deep learning tool for motion analysis. Several cardiac parameters, such as heart rate, diastolic–systolic volumes (EDV/ESV), stroke volume, cardiac output, shortening fraction, and ejection fraction were obtained by the DLC-trained model and then compared to the previous published methods, Time Series Analysis and Kymograph. This new method has several advantages, having full automation, precise detection, and real-time labelling. This network was also trained to analyze zebrafish with cardiovascular defects (pericardial edema) induced by chemical treatments with ethanol and ponatinib. It was revealed that the heart rate, EDV/ESV, stroke volume, and cardiac output from both the ethanol and ponatinib groups displayed significant reductions compared with the control. Hopefully, this trained DLC network can contribute to a better understanding and investigation of the existing cardiovascular system and abnormalities.

Abstract

DeepLabCut (DLC) is a deep learning-based tool initially invented for markerless pose estimation in mammals. In this study, we explored the possibility of adopting this tool for conducting markerless cardiac physiology assessment in an important aquatic toxicology model of zebrafish (Danio rerio). Initially, high-definition videography was applied to capture heartbeat information at a frame rate of 30 frames per second (fps). Next, 20 videos from different individuals were used to perform convolutional neural network training by labeling the heart chamber (ventricle) with eight landmarks. Using Residual Network (ResNet) 152, a neural network with 152 convolutional neural network layers with 500,000 iterations, we successfully obtained a trained model that can track the heart chamber in a real-time manner. Later, we validated DLC performance with the previously published ImageJ Time Series Analysis (TSA) and Kymograph (KYM) methods. We also evaluated DLC performance by challenging experimental animals with ethanol and ponatinib to induce cardiac abnormality and heartbeat irregularity. The results showed that DLC is more accurate than the TSA method in several parameters tested. The DLC-trained model also detected the ventricle of zebrafish embryos even in the occurrence of heart abnormalities, such as pericardial edema. We believe that this tool is beneficial for research studies, especially for cardiac physiology assessment in zebrafish embryos.

Hydration Influence on the Autonomic Recovery of the Coronary Diseases Patient: Geometric Indices Analysis

Purpose: To evaluate the recovery period of autonomic modulation, through geometric indices of heart rate variability (HRV), on coronary artery disease (CAD) patients submitted to a cardiovascular rehabilitation session (CR), associated with hydration. Methods: Thirty male participants of a CR program, diagnosed with CAD were submitted to the control (CP) and hydration protocol (HP) characterized by a CR session. Only during HP were the participants given 8 equal portions of water. The water amount was determined through the hydric loss measured at the CP. During the protocols, the heart rate was measured beat-by-beat at rest (5-10 minutes[M1]) and at recovery (0-5 minutes [M2], 5-10 minutes [M3], 15-20 minutes [M4], 25-30 minutes [M5], 40-45 minutes [M6], 55-60 minutes [M7]) for the HRV analysis, performed by the geometric indices: TINN, RRTRI, SD1, SD2 and SD1/SD2 ratio. Results: Statistically significant differences were observed between the protocols (SD1, p_value = 0.022), moments (TINN, p_value = 0.001; SD1, p_value = 0.019; SD2, p_value = 0.001; SD1/SD2, p_value = 0.001) and moments vs. protocol interaction (SD1, p_value = 0.019). The SD1 index pointed to acceleration of parasympathetic recovery in the first minutes after exercising (HP recovery after M3 [86.07 ± 32.31%] vs. CP recovery after M5[86.43 ± 24.56]) and increase in global variability (TINN-HP remained increased in longer, until M5 (M1 83.10 ± 55.76 ms to M5 116.82 ± 67.54 ms) vs. CP that remained increased for a short time, until M2 (M1 77.93 ± 68.56 ms to M2 134.82 ± 56.08 ms). Conclusions: In CAD patients, hydration promoted a more efficient recovery on parasympathetic autonomic modulation and increased the global HRV in the recovery period.

Reduction of Cardiac Autonomic Modulation and Increased Sympathetic Activity by Heart Rate Variability in Patients With Long COVID

Although several clinical manifestations of persistent long coronavirus disease (COVID-19) have been documented, their effects on the cardiovascular and autonomic nervous system over the long term remain unclear. Thus, we examined the presence of alterations in cardiac autonomic functioning in individuals with long-term manifestations. The study was conducted from October 2020 to May 2021, and an autonomic assessment was performed to collect heart rate data for the heart rate variability (HRV) analysis. The study participants were divided into the long COVID clinical group, the intragroup, which included patients who were hospitalized, and those who were not hospitalized and were symptomatic for different periods (≤3, >3, ≤6, and >6 months), with and without dyspnoea. The control group, the intergroup, comprised of COVID-free individuals. Our results demonstrated that the long COVID clinical group showed reduced HRV compared with the COVID-19-uninfected control group. Patients aged 23–59 years developed COVID symptoms within 30 days after infection, whose diagnosis was confirmed by serologic or reverse-transcription polymerase chain reaction (swab) tests, were included in the study. A total of 155 patients with long COVID [95 women (61.29%), mean age 43.88 ± 10.88 years and 60 men (38.71%), mean age 43.93 ± 10.11 years] and 94 controls [61 women (64.89%), mean age 40.83 ± 6.31 and 33 men (35.11%), mean age 40.69 ± 6.35 years] were included. The intragroup and intergroup comparisons revealed a reduction in global HRV, increased sympathetic modulation influence, and a decrease in parasympathetic modulation in long COVID. The intragroup showed normal sympathovagal balance, while the intergroup showed reduced sympathovagal balance. Our findings indicate that long COVID leads to sympathetic excitation influence and parasympathetic reduction. The excitation can increase the heart rate and blood pressure and predispose to cardiovascular complications. Short-term HRV analysis showed good reproducibility to verify the cardiac autonomic involvement.

A comparative study of heart rate variability and physical fitness in women with moderate and severe fibromyalgia

The aim of this study was to compare the physical fitness and cardiac autonomic activity among women with moderate and severe fibromyalgia (FM) and healthy women. This study included 35 women with FM (age: 46.2±8.9 years) and 17 healthy women (age: 44.3±9.9 years). Participants with FM were divided into moderate FM (n=15) and severe FM (n=20) according to the total score obtained in FM impact questionnaire. The heart rate variability was monitored using a portable cardiac monitor with participants resting in supine position during 10 min. Thereafter, the participants performed the chair sit and reach test, the chair stand test, and the 6-min walk test to measure the lower-body flexibility, lower-body muscle strength, and cardiorespiratory fitness, respectively. The lower-body muscle strength and cardiorespiratory fitness were both reduced in moderate and severe FM compared to healthy women (P<0.01), with greater reduction in severe FM when compared to moderate FM (P<0.05). In addition, the parasympathetic indexes of heart rate variability were all similarly decreased in both moderate and severe FM, when compared to healthy women (P<0.05). The cardiac parasympathetic activity is similarly decreased in women with both moderate and severe FM in comparison to healthy women, despite a greater physical deconditioning in severe FM.

Reliability of Symbolic Analysis of Heart Rate Variability and Its Changes During Sympathetic Stimulation in Elite Modern Pentathlon Athletes: A Pilot Study

Background and Purpose

Most studies on heart rate variability (HRV) in professional athletes concerned linear, time-, and frequency-domain indices, and there is lack of studies on non-linear parameters in this group. The study aimed to determine the inter-day reliability, and group-related and individual changes of short-term symbolic dynamics (SymDyn) measures during sympathetic nervous system activity (SNSa) stimulation among elite modern pentathletes.

Methods

Short-term electrocardiographic recordings were performed in stable measurement conditions with a 7-day interval between tests. SNSa stimulation via isometric handgrip strength test was conducted on the second day of study. The occurrence rate of patterns without variations (0V), with one variation (1V), two like (2LV), and two unlike variations (2UV) obtained using three approaches (the Max–min, the σ, and the Equal-probability methods) were analyzed. Relative and absolute reliability were evaluated.

Results

All SymDyn indices obtained using the Max–min method, 0V, and 2UV obtained using the σ method, 2UV obtained using the Equal-probability method presented acceptable inter-day reliability (the intraclass correlation coefficient between .91 and .99, Cohen’s d between −.08 and .10, the within-subject coefficient of variation between 4% and 22%). 2LV, 2UV, and 0V obtained using the Max–min and σ methods significantly decreased and increased, respectively, during SNSa stimulation—such changes were noted for all athletes. There was no significant association between differences in SymDyn parameters and respiratory rate in stable conditions and while comparing stable conditions and SNSa stimulation.

Conclusion

SymDyn indices may be used as reliable non-respiratory-associated parameters in laboratory settings to detect autonomic nervous system (ANS) activity modulations in elite endurance athletes. These findings provide a potential solution for addressing the confounding influence of respiration frequency on HRV-derived inferences of cardiac autonomic function. For this reason, SymDyn may prove to be preferable for field-based monitoring where measurements are unsupervised.

Cardiac Autonomic Modulation in Subjects with Amyotrophic Lateral Sclerosis (ALS) during an Upper Limb Virtual Reality Task: A Prospective Control Trial

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease. As a result of the rapid progression and severity of the disease, people with ALS experience loss of functionality and independence. Furthermore, it has already been described presence of autonomic dysfunction. Despite the increasing use of virtual reality (VR) in the treatment of different diseases, the use of virtual reality environment as an intervention program for ALS patients is innovative. The benefits and limitations have not yet been proven. Our objective was to evaluate the autonomic function of individuals with amyotrophic lateral sclerosis throughout the virtual reality task. The analysis of autonomic function was completed before, during, and after the virtual reality task using the upper limbs; also, all steps lasted ten minutes in a sitting position. Heart rate variability (HRV) was taken via the Polar® RS800CX cardiofrequencymeter. The following questionnaire was enforced: Amyotrophic Lateral Sclerosis Functional Rating Scale-Revised (ALSFRS) and Fatigue Severity Scale (FSS). Different types of HRV were revealed for the groups, indicating that the ALS group has reduced HRV, with most of the representative indices of the sympathetic nervous system. Besides, the physiological process of reducing parasympathetic activity from rest to VR activity (vagal withdrawal), with reduction in HF (ms²) and an increase in HR from rest to activity, and a further increase throughout recovery, with withdrawal of sympathetic nervous system, occurs just for the control group (CG), with no alterations between rest, activity, and recovery in individuals with ALS. We could conclude that patients with ALS have the reduction of HRV with the sympathetic predominance when equated to the healthy CG. Besides that, the ALS individuals have no capability to adapt the autonomic nervous system when likened to the CG during therapy based on VR and their recovery.

Driver Monitoring of Automated Vehicles by Classification of Driver Drowsiness Using a Deep Convolutional Neural Network Trained by Scalograms of ECG Signals

Driver drowsiness is one of the leading causes of traffic accidents. This paper proposes a new method for classifying driver drowsiness using deep convolution neural networks trained by wavelet scalogram images of electrocardiogram (ECG) signals. Three different classes were defined for drowsiness based on video observation of driving tests performed in a simulator for manual and automated modes. The Bayesian optimization method is employed to optimize the hyperparameters of the designed neural networks, such as the learning rate and the number of neurons in every layer. To assess the results of the deep network method, heart rate variability (HRV) data is derived from the ECG signals, some features are extracted from this data, and finally, random forest and k-nearest neighbors (KNN) classifiers are used as two traditional methods to classify the drowsiness levels. Results show that the trained deep network achieves balanced accuracies of about 77% and 79% in the manual and automated modes, respectively. However, the best obtained balanced accuracies using traditional methods are about 62% and 64%. We conclude that designed deep networks working with wavelet scalogram images of ECG signals significantly outperform KNN and random forest classifiers which are trained on HRV-based features.

Does the ventricle limit cardiac contraction rate in the anoxic turtle (Trachemys scripta)? II. In vivo and in vitro assessment of the prevalence of cardiac arrhythmia and atrioventricular block

Previous studies have reported evidence of atrio-ventricular (AV) block in the oxygen-limited Trachemys scripta heart. However, if cardiac arrhythmia occurs in live turtles during prolonged anoxia exposure remains unknown. Here, we compare the effects of prolonged anoxic submergence and subsequent reoxygenation on cardiac electrical activity through in vivo electrocardiogram (ECG) recordings of 21 °C- and 5 °C-acclimated turtles to assess the prevalence of cardiac arrhythmia. Additionally, to elucidate the influence of extracellular conditions on the prominence of cardiac arrhythmia, we exposed spontaneously contracting T. scripta right atrium and electrically coupled ventricle strip preparations to extracellular conditions that sequentially and additively approximated the shift from the normoxic to anoxic extracellular condition of warm- and cold-acclimated turtles. Cardiac arrhythmia was prominent in 21 °C anoxic turtles. Arrhythmia was qualitatively evidenced by groupings of contractions in pairs and trios and quantified by an increased coefficient of variation of the RR interval. Similarly, exposure to combined anoxia, acidosis, and hyperkalemia induced arrhythmia in vitro that was not counteracted by hypercalcemia or combined hypercalcemia and heightened adrenergic stimulation. By comparison, cold acclimation primed the turtle heart to be resilient to cardiac arrhythmia. Although cardiac irregularities were present intermittently, no change in the variation of the RR interval occurred in vivo with prolonged anoxia exposure at 5 °C. Moreover, the in vitro studies at 5 °C highlighted the importance of adrenergic stimulation in counteracting AV block. Finally, at both acclimation temperatures, cardiac arrhythmia and irregularities ceased upon reoxygenation, indicating that the T. scripta heart recovers from anoxia-induced disruptions to cardiac excitation.

•

Cardiac arrhythmia was prominent in 21 °C anoxic turtles.

•

Cold acclimation primes the turtle heart to be resilient to the cardiac arrhythmia induced by prolonged anoxic submergence.

•

Adrenergic stimulation counteracts atrioventricular block at 5 °C.

•

The turtle heart recovers from anoxia-induced disruptions to cardiac electrical activity.

Autism Spectrum Disorder in Children Is Not Associated With Abnormal Autonomic Nervous System Function: Hypothesis and Theory

The quest to understand the pathophysiology of autism spectrum disorder (ASD) has led to extensive literature that purports to provide evidence for autonomic dysfunction based on heart rate and heart rate variability (HRV), in particular respiratory sinus arrhythmia (RSA), a measure of parasympathetic functioning. Many studies conclude that autism is associated with vagal withdrawal and sympathetic hyperactivation based on HRV and electrodermal analyses. We will argue that a critical analysis of the data leads to the hypothesis that autonomic nervous system dysfunction is not a dominant feature of autism. Most children with ASD have normal parasympathetic baseline values and normal autonomic responses to social stimuli. The existing HRV and electrodermal data cannot lead to the conclusion of an over-excitation of the sympathetic nervous system. A small subgroup of ASD children in experimental settings has relatively low RSA values and relatively high heart rates. The data suggest that this is likely associated with a relatively high level of anxiety during study conditions, associated with co-morbidities such as constipation, or due to the use of psychoactive medication. Many studies interpret their data to conform with a preferred hypothesis of autonomic dysfunction as a trait of autism, related to the polyvagal theory, but the HRV evidence is to the contrary. HRV analysis may identify children with ASD having autonomic dysfunction due to co-morbidities.

Assessing rheoencephalography dynamics through analysis of the interactions among brain and cardiac networks during general anesthesia

Cerebral blood flow (CBF) reflects the rate of delivery of arterial blood to the brain. Since no nutrients, oxygen or water can be stored in the cranial cavity due to space and pressure restrictions, a continuous perfusion of the brain is critical for survival. Anesthetic procedures are known to affect cerebral hemodynamics, but CBF is only monitored in critical patients due, among others, to the lack of a continuous and affordable bedside monitor for this purpose. A potential solution through bioelectrical impedance technology, also known as rheoencephalography (REG), is proposed, that could fill the existing gap for a low-cost and effective CBF monitoring tool. The underlying hypothesis is that REG signals carry information on CBF that might be recovered by means of the application of advanced signal processing techniques, allowing to track CBF alterations during anesthetic procedures. The analysis of REG signals was based on geometric features extracted from the time domain in the first place, since this is the standard processing strategy for this type of physiological data. Geometric features were tested to distinguish between different anesthetic depths, and they proved to be capable of tracking cerebral hemodynamic changes during anesthesia. Furthermore, an approach based on Poincaré plot features was proposed, where the reconstructed attractors form REG signals showed significant differences between different anesthetic states. This was a key finding, providing an alternative to standard processing of REG signals and supporting the hypothesis that REG signals do carry CBF information. Furthermore, the analysis of cerebral hemodynamics during anesthetic procedures was performed by means of studying causal relationships between global hemodynamics, cerebral hemodynamics and electroencephalogram (EEG) based-parameters. Interactions were detected during anesthetic drug infusion and patient positioning (Trendelenburg positioning and passive leg raise), providing evidence of the causal coupling between hemodynamics and brain activity. The provided alternative of REG signal processing confirmed the hypothesis that REG signals carry information on CBF. The simplicity of the technology, together with its low cost and easily interpretable outcomes, should provide a new opportunity for REG to reach standard clinical practice. Moreover, causal relationships among the hemodynamic physiological signals and brain activity were assessed, suggesting that the inclusion of REG information in depth of anesthesia monitors could be of valuable use to prevent unwanted CBF alterations during anesthetic procedures.

What Does Sleeping Brain Tell About Stress? A Pilot Functional Near-Infrared Spectroscopy Study Into Stress-Related Cortical Hemodynamic Features During Sleep

People with mental stress often experience disturbed sleep, suggesting stress-related abnormalities in brain activity during sleep. However, no study has looked at the physiological oscillations in brain hemodynamics during sleep in relation to stress. In this pilot study, we aimed to explore the relationships between bedtime stress and the hemodynamics in the prefrontal cortex during the first sleep cycle. We tracked the stress biomarkers, salivary cortisol, and secretory immunoglobulin A (sIgA) on a daily basis and utilized the days of lower levels of measured stress as natural controls to the days of higher levels of measured stress. Cortical hemodynamics was measured using a cutting-edge wearable functional near-infrared spectroscopy (fNIRS) system. Time-domain, frequency-domain features as well as nonlinear features were derived from the cleaned hemodynamic signals. We proposed an original ensemble algorithm to generate an average importance score for each feature based on the assessment of six statistical and machine learning techniques. With all channels counted in, the top five most referred feature types are Hurst exponent, mean, the ratio of the major/minor axis standard deviation of the Poincaré plot of the signal, statistical complexity, and crest factor. The left rostral prefrontal cortex (RLPFC) was the most relevant sub-region. Significantly strong correlations were found between the hemodynamic features derived at this sub-region and all three stress indicators. The dorsolateral prefrontal cortex (DLPFC) is also a relevant cortical area. The areas of mid-DLPFC and caudal-DLPFC both demonstrated significant and moderate association to all three stress indicators. No relevance was found in the ventrolateral prefrontal cortex. The preliminary results shed light on the possible role of the RLPCF, especially the left RLPCF, in processing stress during sleep. In addition, our findings echoed the previous stress studies conducted during wake time and provides supplementary evidence on the relevance of the dorsolateral prefrontal cortex in stress responses during sleep. This pilot study serves as a proof-of-concept for a new research paradigm to stress research and identified exciting opportunities for future studies.

Detection of Changes on Parameters Related to Heart Rate Variability after Applying Current Interferential Therapy in Subjects with Non-Specific Low Back Pain

Interferential current therapy (ICT) is an electrotherapeutic intervention that combines the advantages of high permeability from middle frequency currents and efficient tissue stimulation from low frequency currents, delivering the maximum current with high tissue permeability. The aim was to evaluate the effects of ICT on heart rate variability (HRV) and on pain perception in patients with non-specific chronic low back pain (NSCLBP). In the study, 49 patients with NSCLBP were randomly divided into an experimental (EG) and a sham group (SG). All participants received a single intervention, ICT, or simulated intervention. Outcome measures including baseline (sit-down position) and postintervention (prone position) pain, heart rate (HR), time domain parameter (rMSSD), diameters of the Poincaré plot (SD1, SD2), stress score (SS), and sympathetic/parasympathetic (S/PS) ratio were investigated. In both groups, significant statistical differences were found in perceived pain and in all HRV parameters except in HRmax. Between-group comparisons showed statistically significant differences in all variables except for HRmin and HRmean in favor of the experimental group. These changes reported an increase in parasympathetic activity (rMSSD) (p < 0.05) and a decrease in sympathetic activity (increase in SD2 and decrease in SS) (p < 0.001) and perceived pain (p < 0.001), with a greater size effect (η2 = 0.44) in favor of the experimental group. In conclusion, a single session of ICT can shift the autonomic balance towards increase parasympathetic dominance and decrease the sympathetic dominance and intensity of pain perceived by patients with NSCLBP.

Psychophysiological dynamics of emotional reactivity: Interindividual reactivity characterization and prediction by a machine learning approach

The fast reaction of the autonomic nervous system (ANS) to an emotional challenge (EC) is the result of a functional coupling between parasympathetic (PNS) and sympathetic (SNS) branches. This coupling can be characterized by measures of cross-correlations between electrodermal activity (EDA) (under the influence of the SNS) and the RR interval (the interval between R peaks) (under the influence of the PNS and the SNS). Significant interindividual variability has previously been reported in SNS-PNS coupling in emotional situations, and the present study aimed to identify interindividual cross-correlation variability in ANS reactivity. We therefore studied EDA and the RR interval in 62 healthy subjects, recorded during a 24-minute EC. A Gaussian Mixture Model was used to cluster tonic EDA-RR cross-correlations during the EC. This identified two clusters that were characterized by significant or non-significant cross-correlations (SCC and NCC clusters, respectively). The SCC cluster reported higher negative emotion after the EC, while the NCC cluster reported higher scores on the Center for Epidemiologic Studies-Depression scale. The latter finding suggests that NCC is a pathological mood pattern with altered negative perception. Furthermore, a machine learning model that included three parameters indexing the functionality of both branches of the ANS, measured at baseline, predicted cluster membership. Our results are a first step in detecting dysfunctional ANS reactivity in general population.

Heart Rate Variability in Psychology: A Review of HRV Indices and an Analysis Tutorial

The use of heart rate variability (HRV) in research has been greatly popularized over the past decades due to the ease and affordability of HRV collection, coupled with its clinical relevance and significant relationships with psychophysiological constructs and psychopathological disorders. Despite the wide use of electrocardiograms (ECG) in research and advancements in sensor technology, the analytical approach and steps applied to obtain HRV measures can be seen as complex. Thus, this poses a challenge to users who may not have the adequate background knowledge to obtain the HRV indices reliably. To maximize the impact of HRV-related research and its reproducibility, parallel advances in users' understanding of the indices and the standardization of analysis pipelines in its utility will be crucial. This paper addresses this gap and aims to provide an overview of the most up-to-date and commonly used HRV indices, as well as common research areas in which these indices have proven to be very useful, particularly in psychology. In addition, we also provide a step-by-step guide on how to perform HRV analysis using an integrative neurophysiological toolkit, NeuroKit2.

Recent Research for Unobtrusive Atrial Fibrillation Detection Methods Based on Cardiac Dynamics Signals: A Survey

Atrial fibrillation (AF) is the most common cardiac arrhythmia. It tends to cause multiple cardiac conditions, such as cerebral artery blockage, stroke, and heart failure. The morbidity and mortality of AF have been progressively increasing over the past few decades, which has raised widespread concern about unobtrusive AF detection in routine life. The up-to-date non-invasive AF detection methods include electrocardiogram (ECG) signals and cardiac dynamics signals, such as the ballistocardiogram (BCG) signal, the seismocardiogram (SCG) signal and the photoplethysmogram (PPG) signal. Cardiac dynamics signals can be collected by cushions, mattresses, fabrics, or even cameras, which is more suitable for long-term monitoring. Therefore, methods for AF detection by cardiac dynamics signals bring about extensive attention for recent research. This paper reviews the current unobtrusive AF detection methods based on the three cardiac dynamics signals, summarized as data acquisition and preprocessing, feature extraction and selection, classification and diagnosis. In addition, the drawbacks and limitations of the existing methods are analyzed, and the challenges in future work are discussed.

Life events are associated with elevated heart rate and reduced heart complexity to acute psychological stress

The current study examined whether the exposure to life events and reported impact of life events are associated with altered cardiac reactivity to an acute psychological stressor. Participants (N = 69) completed the Life Experience Survey (LES) and Positive and Negative Affect Schedule (PANAS) and undertook a standardized social-evaluative stress task. Cardiac activity was measured via heart rate and non-linear heart rate variability (HRV) indices Sample Entropy, SD1, SD2 and SD1/SD2 ratio. Heart rate and non-linear HRV were measured before, during and after stress exposure. Findings suggest higher heart rate reactivity in individuals reporting higher number and impact of negative and total life events. Decreases in Sample Entropy were evident for number as well as impact of life events. No associations were found for SD1, SD2 and SD1/SD2 ratio. Findings suggest that life-events are associated with elevated heart rate and diminished heart rate complexity in response to acute stress.

Physiological Reactions in the Therapist and Turn-Taking during Online Psychotherapy with Children and Adolescents with Autism Spectrum Disorder

This study aims to analyze the relationship between the sociocognitive skills of a group of children and adolescents with autism spectrum disorder (ASD) at verbal level 1, the variability of the therapist's heart rate (HRV), and the conversational turn-taking during online psychotherapy sessions. Initially, we assessed the intelligence, narrative, and behavioral characteristics of the participants. We videotaped the online sessions and recorded the therapist's HRV via a smart wireless sensor. Finally, we analyzed the video sessions using an observation system and the therapist's HRV using the Poincaré technique. The results show that the patients' communicative intention was related to their narrative, intellectual and social competencies. Furthermore, the turn-taking between the therapist and the participant was associated with the patient's emotional and behavioral difficulties. On the other side, the therapist's heart rate variability (HRV) was related to the synchrony between the therapist and the participant with more significant stress on the therapist, when he shared and expanded the conversation with the patient, and when the patient broadened and shared the conversation with the therapist.

Heart Rate Variability Monitoring during Interferential Current Application in the Lower Back Area: A Cross-Sectional Study

Vasovagal reactions may occur occasionally during electrical stimulation using interferential current (IFC). The purpose of this study was to examine variations in autonomic activity during the application of IFC in asymptomatic participants by analysis of their heart rate variability (HRV). Seventy-three male volunteers were randomly assigned to a placebo group (n = 36; HRV was documented for 10 min, both at rest and during a placebo intervention) and an intervention group (n = 37; HRV was documented for 10 min in two conditions labelled as (1) rest and (2) application of IFC technique on the lumbar segment). The diameters of the Poincaré plot (SD1, SD2), stress score (SS), and the ratio between sympathetic and parasympathetic activity (S/PS) were measured. After interventions, differences amongst the placebo group and the IFC group were found in SD2 (p < 0.001), SS (p = 0.01) and S/PS ratio (p = 0.003). The IFC technique was associated with increased parasympathetic modulation, which could induce a vasovagal reaction. Monitorization of adverse reactions should be implemented during the application of IFC technique. HRV indicators might have a part in prevention of vasovagal reactions. Further studies in patients with lumbar pain are needed to explore possible differences in HRV responses due to the presence of chronic pain.

Reduced heart-rate variability and increased risk of hypertension-a prospective study of the ELSA-Brasil

Although autonomic disturbances are known to precede hypertension onset, the risks underlying different ranges of blood pressure and impaired cardiac autonomic modulation are still unknown. This study aimed to identify the risk of hypertension incidence related to low heart-rate variability profile in normotensive blood pressure subcategories: normal (<120/80 mmHg) and prehypertension (120/80-139/89 mmHg) in a 4-year follow-up. 7665 participants free of hypertension at baseline were examined. They were allocated into one of two groups (<P25 and ≥P25) based on age-specific value distributions of heart-rate variability ranges presented in the literature. The relationship between heart-rate variability parameters (SDNN, RMSSD, VLF, LF, and HF) and hypertension incidence was analyzed using Poisson regression. The reference for the independent variable in these regressions was the dataset with all HRV values ≥P25. After full adjustments for age, sex, ethnicity, educational level, body mass index, smoking status, glucose, insulin, dyslipidemia, alcohol consumption, and physical activity, low values for all heart-rate variability indices were significantly associated with an increased relative risk of developing hypertension in the overall sample. Individuals with normal BP levels presented an increased risk for low SDNN, VLF, and LF values, whereas the risk in prehypertensive participants was only verified for low VLF. In conclusion, cardio autonomic disturbances precede the incidence of clinical hypertension and are potentially involved in its pathophysiological basis and progression. Even individuals with normal blood pressure may be at an increased risk for future hypertension in the presence of reduced cardio autonomic control.

Resting Heart Rate Variability Is Associated With Subsequent Orthostatic Hypotension: Comparison Between Healthy Older People and Patients With Rapid Eye Movement Sleep Behavior Disorder

Background: Orthostatic hypotension (OH) caused by autonomic dysfunction is a common symptom in older people and patients with idiopathic rapid eye movement sleep behavior disorder (iRBD). The orthostatic challenge test is a standard autonomic function test that measures a decrease of blood pressure during a postural change from supine to standing positions. Although previous studies have reported that changes in heart rate variability (HRV) are associated with autonomic dysfunction, no study has investigated the relationship between HRV before standing and the occurrence of OH in an orthostatic challenge test. This study aims to examine the connection between HRV in the supine position and the occurrence of OH in an orthostatic challenge test. Methods: We measured the electrocardiograms of patients with iRBD and healthy older people during an orthostatic challenge test, in which the supine and standing positions were held for 15 min, respectively. The subjects were divided into three groups: healthy controls (HC), OH-negative iRBD [OH (-) iRBD], and OH-positive iRBD [OH (+) iRBD]. HRV measured in the supine position during the test were calculated by time-domain analysis and Poincaré plots to evaluate the autonomic dysfunction. Results: Forty-two HC, 12 OH (-) iRBD, and nine OH (+) iRBD subjects were included. HRV indices in the OH (-) and the OH (+) iRBD groups were significantly smaller than those in the HC group. The multivariate logistic regression analysis for OH identification for the iRBD groups showed the model whose inputs were the HRV indices, i.e., standard deviation 2 (SD2) and the percentage of adjacent intervals that varied by more than 50 ms (pNN50), had a receiver operating characteristic curve with area under the curve of 0.840, the sensitivity to OH (+) of 1.000, and the specificity to OH (-) of 0.583 (p = 0.023). Conclusions: This study showed the possibility that short-term HRV indices in the supine position would predict subsequent OH in iRBD patients. Our results are of clinical importance in terms of showing the possibility that OH can be predicted using only HRV in the supine position without an orthostatic challenge test, which would improve the efficiency and safety of testing.

Assessment of autonomous nerve system through non-linear heart rate variability outcomes in sedentary healthy adults

Background

Heart rate variability (HRV) is a psycho-physiological phenomenon with broad health implications. Different data analysis methods have been used to assess the autonomic nervous system activity, but the validation of new indexes that accurately describe its balance through non-invasive methods (i.e., HRV analysis) is of clinical interest. This study aimed: (i) to evaluate the association of the Stress Score (SS) and the Sympathetic/Parasympathetic Ratio (S/PS) with time domain and frequency domain analysis of HRV, and (ii) to set reference values of SS and S/PS in sedentary healthy adults.

Methods

A total of 156 sedentary healthy adults (38.4 ± 15.57 years old, 81 women), aged were involved in this study. HRV was measured for 15 min in a supine position at rest. SS and S/PS were calculated from the non-linear HRV analyses based on Poincare Plot.

Results

Stress Score showed a non-linear negative power-law relationship with SDNN (β = -0.969; R ² = 0.963; P < 0.001), RMSSD (β = -0.867; R ² = 0.722; P < 0.001), high frequency (β = -0.834; R ² = 0.752; P =< 0.001), low frequency (β = -0.627; R ² = 0.330; P < 0.001), SD1 (β = -0.867; R ² = 0.722; P < 0.001) and SD2 (β = -1.000; R ² > 0.999; P < 0.001). There was observed a negative cubic relationship between SS with PNN50 (β = -1.972; R ² = 0.644; P < 0.001). A linear regression model was conducted between SS with Ratio Low/High Frequency (β = 0.026; R ² < 0.001; P = 0.750). Non-linear power-law regression models were built between S/PS and SDNN (β = -0.990; R ² = 0.981; P < 0.001), RMSSD (β = -0.973; R ² = 0.939; P < 0.001), high frequency (β = -0.928; R ² = 0.970; P < 0.001), low frequency (β = -2.344; R ² = 0.557; P < 0.001), SD1 (β = -0.973; R ² = 0.939; P < 0.001) and SD2 (β = -0.611; R ² = 0.908; P < 0.001). A non-linear negative regression model was built between S/PS and PNN50 (β = -3.412; R ² = 0.868; P < 0.001). A linear regression model was conducted between S/PS and SD2/SD1 (β = 0.075; R ² = 0.006; P < 0.001).

Conclusion

Our results support the use of SS as a sympathetic activity marker, and S/PS as an indicator of the sympathetic and parasympathetic activity of the autonomic nervous system in sedentary healthy adults.

Cardiac autonomic recovery following traditional and augmented remote ischemic preconditioning

PURPOSE

While the possible ergogenic benefits of remote ischemic preconditioning (RIPC) make it an attractive training modality, the mechanisms of action remain unclear. Alterations in neural tone have been demonstrated in conjunction with circulatory occlusion, yet investigation of the autonomic nervous system following RIPC treatment has received little attention. We sought to characterize alterations in autonomic balance to both RIPC and augmented RIPC (RIPC_aug) performed while cycling, using acute and sustained autonomic indices.

METHODS

Thirteen participants (8M:5F) recorded baseline waking heart rate variability (HRV) for 5 days prior to treatment. Participants then completed control exercise (CON), RIPC, and RIPC_aug interventions in a randomized cross-over design. Cardiovascular measurements were recorded immediately before and after each intervention at rest, and during an orthostatic challenge. Waking HRV was repeated the morning after each intervention.

RESULTS

RIPC resulted in acutely reduced resting heart rates (HR) (∆ - 4 ± 6 bpm, P = 0.02) and suppressed HR 30 s following the orthostatic challenge compared to CON (64 ± 10 vs 74 ± 9 bpm, P = 0.003). RIPC_aug yielded elevated HRs compared to CON and RIPC prior to (P = 0.003) and during the orthostatic challenge (P = 0.002). RIPC_aug reduced LnSDNN (Baseline 4.39 ± 0.27; CON 4.44 ± 0.39; RIPC 4.41 ± 0.34; RIPC_aug 4.22 ± 0.29, P = 0.02) and LnHfa power (Baseline 7.82 ± 0.54; CON 7.73 ± 1.11; RIPC 7.89 ± 0.78; RIPC_aug 7.23 ± 0.87, P = 0.04) the morning after treatment compared to all other conditions.

CONCLUSIONS

Our data suggest that RIPC may influence HR acutely, possibly through a reduction in cardiac sympathetic activity, and that RIPC_aug reduces HRV through cardiac vagal withdrawal or increased cardiac sympathetic modulation, with alterations persisting until the following morning. These findings imply a dose-response relationship with potential for optimization of performance.

Aircraft Pilots Workload Analysis: Heart Rate Variability Objective Measures and NASA-Task Load Index Subjective Evaluation

Workload and fatigue of aircraft pilots represent an argument of great interest in the framework of human factors and a pivotal point to be considered in aviation safety. 75% of aircraft accidents are related to human errors that, in most cases, are due to high level of mental workload and fatigue. There exist several subjective or objective metrics to quantify the pilots’ workload level, with both linear and nonlinear relationships reported in the literature. The main research objective of the present work is to analyze the relationships between objective and subjective workload measurements by looking for a correlation between metrics belonging to the subjective and biometric rating methods. More particularly, the Heart Rate Variability (HRV) is used for the objective analysis, whereas the NASA-TLX questionnaire is the tool chosen for the subjective evaluation of the workload. Two different flight scenarios were considered for the studies: the take-off phase with the initial climb and the final approach phase with the landing. A Maneuver Error Index (MEI) is also introduced to evaluate the pilot flight performance according to mission requirements. Both qualitative and quantitative correlation analyses were performed among the MEI, subjective and objective measurements. Monotonic relationships were found within the HRV indexes, and a nonlinear relationship is proposed among NASA-TLX and HRV indexes. These findings suggest that the relationship between workload, biometric data, and performance indexes are characterized by intricate patterns of nonlinear relationships.

Datasets for Cognitive Load Inference Using Wearable Sensors and Psychological Traits

This study introduces two datasets for multimodal research on cognitive load inference and personality traits. Different to other datasets in Affective Computing, which disregard participants’ personality traits or focus only on emotions, stress, or cognitive load from one specific task, the participants in our experiments performed seven different tasks in total. In the first dataset, 23 participants played a varying difficulty (easy, medium, and hard) game on a smartphone. In the second dataset, 23 participants performed six psychological tasks on a PC, again with varying difficulty. In both experiments, the participants filled personality trait questionnaires and marked their perceived cognitive load using NASA-TLX after each task. Additionally, the participants’ physiological response was recorded using a wrist device measuring heart rate, beat-to-beat intervals, galvanic skin response, skin temperature, and three-axis acceleration. The datasets allow multimodal study of physiological responses of individuals in relation to their personality and cognitive load. Various analyses of relationships between personality traits, subjective cognitive load (i.e., NASA-TLX), and objective cognitive load (i.e., task difficulty) are presented. Additionally, baseline machine learning models for recognizing task difficulty are presented, including a multitask learning (MTL) neural network that outperforms single-task neural network by simultaneously learning from the two datasets. The datasets are publicly available to advance the field of cognitive load inference using commercially available devices.

Nonlinear Methods Most Applied to Heart-Rate Time Series: A Review

The heart-rate dynamics are one of the most analyzed physiological interactions. Many mathematical methods were proposed to evaluate heart-rate variability. These methods have been successfully applied in research to expand knowledge concerning the cardiovascular dynamics in healthy as well as in pathological conditions. Notwithstanding, they are still far from clinical practice. In this paper, we aim to review the nonlinear methods most used to assess heart-rate dynamics. We focused on methods based on concepts of chaos, fractality, and complexity: Poincaré plot, recurrence plot analysis, fractal dimension (and the correlation dimension), detrended fluctuation analysis, Hurst exponent, Lyapunov exponent entropies (Shannon, conditional, approximate, sample entropy, and multiscale entropy), and symbolic dynamics. We present the description of the methods along with their most notable applications.

Geometric indexes of heart rate variability in healthy individuals exposed to long-term air pollution

The aim of this study was to investigate the autonomic modulation of heart rate in healthy individuals exposed to long-term air pollution through geometric methods. We analyzed data from 109 healthy adults aged 18 to 49, divided into three groups according to the exposure time: period 0 to 15 years of exposure (n = 29), more than 15 years of exposure (n = 31), and control group (n = 49). For the analysis of heart rate variability (HRV), heart rate was recorded beat-to-beat for 20 min in the sitting position. The RR intervals were transformed into geometric indexes, and from them, we calculated the RRTri (triangular index), TINN (triangle interpolation of histogram of intervals NN), and Poincaré plot (SD1, SD2, and SD1/SD2). Significantly lower values were observed in the group of individuals exposed to air pollution for more than 15 years compared with the group of individuals exposed to air pollution for a period of 0-15 years and those not exposed for the RRTri (11.5 vs 13.8 vs 14.0), SD1 (16.4 vs 20.5 vs 20.6), SD2 (60.5 vs 68.1 vs 72.5), and SD1/SD2 (0.27 vs 0.34 vs 0.31), with the effect of this difference being considered large (RRTri), medium (SD1, SD1/SD2), and small (SD2). TINN was not significantly different among groups (198.2 vs 223.1 vs 233.6). Healthy individuals exposed to air pollution for more than 15 years present an autonomic imbalance, characterized by lower parasympathetic modulation and overall HRV.

Preictal Time Assessment using Heart Rate Variability Features in Drug-resistant Epilepsy Patients

Epileptic seizures are by default associated with the occurrence of EEG changes as a consequence of alterations in brain activity seconds to hours before the seizure onset. Additionally, the influence of the autonomic nervous system may also be reflected in the electrocardiogram (ECG) trace. In both biosignals (EEG and ECG), differences have been reported between normal epochs, known as interictal periods and the interval preceding seizures, preictal period. However, the existence of a preictal state indicating the transition between epileptic brain states has not yet been clinically defined. In fact, some studies report differences in preictal location and duration among patients and among seizures occurring in the same patient as well. Based on the above, a study was designed in order to investigate the existence of a preictal interval specific for each seizure, using heart rate variability (HRV) features. Time and frequency domain features (linear and non-linear) were extracted from ECG data acquired in 37 drug-resistant epilepsy patients, comprised in EPILEPSIAE database. A total of 209 seizures occurring in the temporal lobe were analysed. The existence of a transition period before the seizure onset was inspected using a linear discriminant analysis classifier. The overall best performance (88.04%±12.30% of accuracy) was obtained by combining RRMean, NN50 and SD₂ features when discriminating the 50 min of interictal located farthest from seizure onset from the nearest preictal 50 min.

Brain and Body Emotional Responses: Multimodal Approximation for Valence Classification

In order to develop more precise and functional affective applications, it is necessary to achieve a balance between the psychology and the engineering applied to emotions. Signals from the central and peripheral nervous systems have been used for emotion recognition purposes, however, their operation and the relationship between them remains unknown. In this context, in the present work, we have tried to approach the study of the psychobiology of both systems in order to generate a computational model for the recognition of emotions in the dimension of valence. To this end, the electroencephalography (EEG) signal, electrocardiography (ECG) signal and skin temperature of 24 subjects have been studied. Each methodology has been evaluated individually, finding characteristic patterns of positive and negative emotions in each of them. After feature selection of each methodology, the results of the classification showed that, although the classification of emotions is possible at both central and peripheral levels, the multimodal approach did not improve the results obtained through the EEG alone. In addition, differences have been observed between cerebral and peripheral responses in the processing of emotions by separating the sample by sex; though, the differences between men and women were only notable at the peripheral nervous system level.