PyBioS: A freeware computer software for analysis of cardiovascular signals

Recovery of the baroreflex and autonomic modulation after anesthesia with MS-222 in bullfrogs

The time course for recovery after anesthesia is poorly described for tricaine methanesulfonate (MS-222). We suggest that the baroreflex and the heart rate variability (HRV) could be used to index the recovery of the autonomic modulation after anesthesia. We analyzed the recovery profile of behavioral and physiological parameters over time to analyze the progression of recovery after anesthesia of American bullfrogs with MS-222. Mean heart rate stabilized after 17 h, whereas the baroreflex efficiency index took 23 h and the baroreflex operating gain, 29 h. Mean arterial pressure recovered after 26 h. Power spectral density peaked at 23 h and again after 40 h. Baroreflex was a relevant component of the first phase of HRV, while autonomic modulation for resting may take longer than 40 h. We suggest that physiological recovery is a complex phenomenon with multiple progressive phases, and the baroreflex may be a useful tool to observe the first substantial recovery of post-instrumentation capacity for autonomic modulation.

Analysis of fetal heart rate fluctuations in women diagnosed with preeclampsia during the latent phase of labor

Introduction

Fetal heart rate variability (fHRV) is a tool used to investigate the functioning of the fetal autonomic nervous system. Despite the significance of preeclampsia, fHRV during the latent phase of labor has not been extensively studied. This study aimed to evaluate fetal cardiac autonomic activity by using linear and nonlinear indices of fHRV analysis in women diagnosed with preeclampsia without hypertensive treatment during gestation, compared to normotensive women during the latent phase of labor.

Methods

A cross-sectional and exploratory study was conducted among pregnant women in the latent phase of labor, forming three study groups: normotensive or control (C, 38.8 ± 1.3 weeks of pregnancy, n = 22), preeclampsia with moderate features (P, 37.6 ± 1.4 weeks of pregnancy n = 10), and preeclampsia with severe features (SP, 36.9 ± 1.2 weeks of pregnancy, n = 12). None of the participants received anti-hypertensive treatment during their pregnancy. Linear and nonlinear features of beat-to-beat fHRV, including temporal, frequency, symbolic dynamics, and entropy measures, were analyzed to compare normotensive and preeclamptic groups.

Results

Significantly lower values of multiscale entropy (MSE) and short-term complexity index (Ci) were observed in the preeclamptic groups compared to the C group (p < 0.05). Additionally, higher values of SDNN (standard deviation of R-R intervals) and higher values of low-frequency power (LF) were found in the P group compared to the C group.

Conclusion

Our findings indicate that changes in the complexity of fetal heart rate fluctuations may indicate possible disruptions in the autonomic nervous system of fetuses in groups affected by undiagnosed preeclampsia during pregnancy. Reduced complexity and shifts in fetal autonomic cardiac activity could be associated with preeclampsia's pathophysiological mechanisms during the latent phase of labor.

Heart rate dynamics and asymmetry during sympathetic activity stimulation and post-stimulation recovery in ski mountaineers-a pilot exploratory study

There is a lack of studies on non-linear heart rate (HR) variability in athletes. We aimed to assess the usefulness of short-term HR dynamics and asymmetry parameters to evaluate the neural modulation of cardiac activity based on non-stationary RR interval series by studying their changes during sympathetic nervous system activity stimulation (isometric handgrip test) and post-stimulation recovery in professional ski mountaineers. The correlation between the changes in the parameters and the respiratory rate (RespRate) and also the duration of the career was analyzed. Short-term (5 min) and ultra-short-term (1 min) rates of patterns with no variations (0V), number of acceleration runs of length 1 (AR1), and short-term Porta's Index were greater, whereas Guzik's Index (GI) was smaller during sympathetic stimulation compared to rest. GI increased and the number of AR1 decreased during recovery. Greater increases in GI and RMSSD were associated with greater decreases in RespRate during recovery. Greater increases in RespRate from rest to short-term sympathetic stimulation were associated with greater increases in 0V (Max-min method) and AR1 but also with greater decreases in decelerations of short-term variance and accelerations and decelerations of long-term variance. Greater increases in 0V (Max-min method) and number of AR1 during sympathetic stimulation were associated with a shorter career duration. Greater decreases in these parameters during recovery were associated with a longer career duration. Changes in measures of HR dynamics and asymmetry, calculated based on short-term non-stationary RRi time series induced by sympathetic stimulation and post-stimulation recovery, reflected sympathovagal shift and were associated with condition-related alterations in RespRate and career duration in athletes who practice ski mountaineering.

Acute and short-term beetroot juice nitrate-rich ingestion enhances cardiovascular responses following aerobic exercise in postmenopausal women with arterial hypertension: A triple-blinded randomized controlled trial

BACKGROUND

The increase in blood pressure (BP) levels in the postmenopausal period can be partly explained by the decrease in nitric oxide synthases (NOS).

OBJECTIVE

To investigate the acute and one-week effects of beetroot juice nitrate-rich (BRJ-NO3-rich) ingestion on cardiovascular and autonomic performance in response to submaximal aerobic exercise in postmenopausal women with systemic arterial hypertension (SAH) who are physically inactive.

METHODS

Fourteen postmenopausal women with SAH [mean (SD) age: 59(4) y; BMI (kg/m2): 29.2(3.1)] completed submaximal aerobic exercise bouts after an acute and a one-week intervention with BRJ in a placebo-controlled, randomized, triple-blind, crossover design. Participants ingested either BRJ (800 mg of NO3-) or placebo acutely and drank either BRJ (400 mg of NO3-) or placebo every day for the next six days. After two and ½ hours, they performed a session of aerobic submaximal aerobic exercise, and their systolic BP (SBP) and diastolic BP (DBP), flow-mediated dilation (FMD), heart rate (HR) recovery, and HR variability were measured.

RESULTS

In the post-exercise recovery period, SBP dropped significantly in the BRJ-NO3-rich group (-9.28 mmHg [95%CI: -1.68 to -16.88] ES: -0.65, p = 0.019) compared to placebo after acute ingestion. The FMD values increased after acute BRJ-NO3-rich on post-exercise (3.18 % [0.36 to 5.99] ES: 0.87, p = 0.031). After the one-week intervention, FMD values were higher in the BRJ-NO3-rich group before (4.5 % [1.62 to 7.37] ES: 1.21, p = 0.005) and post-exercise measurements (4.2 % [1.52 to 6.87] ES: 1.22, p = 0.004) vs. placebo. HRV indices with remarkable parasympathetic modulation to heart recovered faster on the BRJ-NO3-rich group than placebo group. No between-group differences were identified in values of HR post-exercise recovery in the 30s, 60s, 120s, 180s, and 300s.

CONCLUSIONS

Acute and short-term BRJ-NO3-rich ingestion may enhance cardiovascular and autonomic behavior in response to aerobic exercise in postmenopausal women diagnosed with SAH.

CLINICAL TRIAL REGISTRY NUMBER

https://clinicaltrials.gov/ct2/show/NCT05384340.

CVRanalysis: a free software for analyzing cardiac, vascular and respiratory interactions

Introduction: Simultaneous beat-to-beat R-R intervals, blood pressure and respiration signals are routinely analyzed for the evaluation of autonomic cardiovascular and cardiorespiratory regulations for research or clinical purposes. The more recognized analyses are i) heart rate variability and cardiac coherence, which provides an evaluation of autonomic nervous system activity and more particularly parasympathetic and sympathetic autonomic arms; ii) blood pressure variability which is mainly linked to sympathetic modulation and myogenic vascular function; iii) baroreflex sensitivity; iv) time-frequency analyses to identify fast modifications of autonomic activity; and more recently, v) time and frequency domain Granger causality analyses were introduced for assessing bidirectional causal links between each considered signal, thus allowing the scrutiny of many physiological regulatory mechanisms. Methods: These analyses are commonly applied in various populations and conditions, including mortality and morbidity predictions, cardiac and respiratory rehabilitation, training and overtraining, diabetes, autonomic status of newborns, anesthesia, or neurophysiological studies. Results: We developed CVRanalysis, a free software to analyze cardiac, vascular and respiratory interactions, with a friendly graphical interface designed to meet laboratory requirements. The main strength of CVRanalysis resides in its wide scope of applications: recordings can arise from beat-to-beat preprocessed data (R-R, systolic, diastolic and mean blood pressure, respiration) or raw data (ECG, continuous blood pressure and respiratory waveforms). It has several tools for beat detection and correction, as well as setting of specific areas or events. In addition to the wide possibility of analyses cited above, the interface is also designed for easy study of large cohorts, including batch mode signal processing to avoid running repetitive operations. Results are displayed as figures or saved in text files that are easily employable in statistical softwares. Conclusion: CVRanalysis is freely available at this website: anslabtools.univ-st-etienne.fr. It has been developed using MATLAB® and works on Windows 64-bit operating systems. The software is a standalone application avoiding to have programming skills and to install MATLAB. The aims of this paper area are to describe the physiological, research and clinical contexts of CVRanalysis, to introduce the methodological approach of the different techniques used, and to show an overview of the software with the aid of screenshots.

NLDyn - An open source MATLAB toolbox for the univariate and multivariate nonlinear dynamical analysis of physiological data

BACKGROUND AND OBJECTIVE

We present NLDyn, an open-source MATLAB toolbox tailored for in-depth analysis of nonlinear dynamics in biomedical signals. Our objective is to offer a user-friendly yet comprehensive platform for researchers to explore the intricacies of time series data.

METHODS

NLDyn integrates approximately 80 distinct methods, encompassing both univariate and multivariate nonlinear dynamics, setting it apart from existing solutions. This toolbox combines state-of-the-art nonlinear dynamical techniques with advanced multivariate entropy methods, providing users with powerful analytical capabilities. NLDyn enables analyses with or without a sliding window, and users can easily access and customize default parameters.

RESULTS

NLDyn generates results that are both exportable and visually informative, facilitating seamless integration into research and presentations. Its ongoing development ensures it remains at the forefront of nonlinear dynamics analysis.

CONCLUSIONS

NLDyn is a valuable resource for researchers in the biomedical field, offering an intuitive interface and a wide array of nonlinear analysis tools. Its integration of advanced techniques empowers users to gain deeper insights from their data. As we continually refine and expand NLDyn's capabilities, we envision it becoming an indispensable tool for the exploration of complex dynamics in biomedical signals.

Investigating autonomic nervous system dysfunction among patients with post-COVID condition and prolonged cardiovascular symptoms

Heart Rate Variability (HRV) and arterial pressure (AP) variability and their responses to head-up tilt test (HUTT) were investigated in Post-COVID-19 syndrome (PCS) patients reporting tachycardia and/or postural hypotension. Besides tachycardia, PCS patients also showed attenuation of the following HRV parameters: RMSSD [square root of the mean of the sum of the squares of differences between adjacent normal-to-normal (NN) intervals] from statistical measures; the power of RR (beat-to-beat interval) spectra at HF (high frequency) from the linear method spectral analysis; occurrence of 2UV (two unlike variation) pattern of RR from the nonlinear method symbolic analysis; and the new family of statistics named sample entropy, when compared to control subjects. Basal AP and LF (low frequency) power of systolic AP were similar between PCS patients and control subjects, while 0 V (zero variation) patterns of AP from the nonlinear method symbolic analysis were exacerbated in PCS patients. Despite tachycardia and a decrease in RMSSD, no parameter of HRV changed during HUTT in PCS patients compared to control subjects. PCS patients reassessed after 6 months showed higher HF power of RR spectra and a higher percentage of 2UV pattern of RR. Moreover, the reassessed PCS patients showed a lower occurrence of 0 V patterns of AP, while the HUTT elicited HR (heart rate) and AP responses identical to control subjects. The HRV and AP variability suggest an autonomic dysfunction with sympathetic predominance in PCS patients. In contrast, the lack of responses of HRV and AP variability indices during HUTT indicates a marked impairment of autonomic control. Of note, the reassessment of PCS patients showed that the noxious effect of COVID-19 on autonomic control tended to fade over time.

Changes in the autonomic cardiorespiratory activity in parturient women with severe and moderate features of preeclampsia

Background

Cardiorespiratory coupling (CRC) is a physiological phenomenon that reflects the mutual interaction between the cardiac and respiratory control systems. It is mainly associated with efferent vagal activity from the central autonomic network. Few studies have explored the autonomic changes of CRC in preeclampsia, a critical obstetric complication related to possible autonomic dysfunctions and inflammatory disturbances. This study examined the autonomic mechanisms of CRC in women with severe and moderate preeclampsia and healthy controls by applying nonlinear methods based on information theory, such as mutual information (MI) and Renyi's mutual information (RMI) and the linear and nonlinear analysis of the Pulse-Respiration Quotient (PRQ).

Methods

We studied three groups of parturient women in the third trimester of pregnancy with a clinical diagnosis of preeclampsia without severe symptoms (P, 38.5 ± 1.4 weeks of pregnancy, n=19), preeclampsia with severe symptoms (SP, 37.5 ± 0.9 weeks of pregnancy, n=22), and normotensive control women (C, 39.1 ± 1.3 weeks of pregnancy, n=20). 10-minutes of abdominal electrocardiograms (ECG) and respiratory signals (RESP) were recorded in all the participants. Subsequently, we obtained the maternal beat-to-beat (RR) and breath-to-breath (BB) time series from ECG and RESP, respectively. The CRC between RR and BB was quantified by nonlinear methods based on information theory, such as MI and RMI, along with the analysis of the novel index of PRQ. Subsequently, we computed the mean PRQ (mPRQ) and the normalized permutation entropy (nPermEn_PRQ) from the PRQ time series generated from BB and RR. In addition, we examined the vagal activity in the three groups by the logarithm of the median of the distribution of the absolute values of successive RR differences (logRSA).

Results

The MI and RMI values were significantly lower (p<0.05) in the preeclamptic groups compared to the control group. However, no significant differences were found between the preeclamptic groups. The logRSA and nPermEn_PRQ indices were significantly lower (p<0.05) in SP compared to C and P.

Conclusion

Our data suggest that parturient women with severe and mild preeclampsia may manifest an altered cardiorespiratory coupling compared with normotensive control women. Disrupted CRC in severe preeclampsia could be associated with vagal withdrawal and less complex cardiorespiratory dynamics. The difference in vagal activity between the preeclamptic groups may suggest a further reduction in vagal activity associated with the severity of the disease.

Comparison of Heart Autonomic Control between Hemodynamically Stable and Unstable Patients during Hemodialysis Sessions: A Bayesian Approach

Intradialytic hypotension is a common complication during hemodialysis sessions. The analysis of successive RR interval variability using nonlinear methods represents a promising tool for evaluating the cardiovascular response to acute volemic changes. Thus, the present study aims to compare the variability of successive RR intervals between hemodynamically stable (HS) and unstable (HU) patients during a hemodialysis session, through linear and nonlinear methods. Forty-six chronic kidney disease patients volunteered in this study. Successive RR intervals and blood pressures were recorded throughout the hemodialysis session. Hemodynamic stability was defined based on the delta of systolic blood pressure (higher SBP-lower SBP). The cutoff for hemodynamic stability was defined as 30 mm Hg, and patients were stratified as: HS ([n = 21]: ≤29.9 mm Hg) or HU ([n = 25]: ≥30 mm Hg). Linear methods (low-frequency [LFnu] and high-frequency [HFnu] spectra) and nonlinear methods (multiscale entropy [MSE] for Scales 1-20, and fuzzy entropy) were applied. The area under the MSE curve at Scales 1-5 (MSE_1-5), 6-20 (MSE_6-20), and 1-20 (MSE_1-20) were also used as nonlinear parameters. Frequentist and Bayesian inferences were applied to compare HS and HU patients. The HS patients exhibited a significantly higher LFnu and lower HFnu. For MSE parameters, Scales 3-20 were significantly higher, as well as MSE_1-5, MSE_6-20, and MSE_1-20 in HS, when compared to HU patients (p < 0.05). Regarding Bayesian inference, the spectral parameters demonstrated an anecdotal (65.9%) posterior probability favoring the alternative hypothesis, while MSE exhibited moderate to very strong probability (79.4 to 96.3%) at Scales 3-20, and MSE_1-5, MSE_6-20, and MSE_1-20. HS patients exhibited a higher heart-rate complexity than HU patients. In addition, the MSE demonstrated a greater potential than spectral methods to differentiate variability patterns in successive RR intervals.

Induced Relaxation Enhances the Cardiorespiratory Dynamics in COVID-19 Survivors

Most COVID-19 survivors report experiencing at least one persistent symptom after recovery, including sympathovagal imbalance. Relaxation techniques based on slow-paced breathing have proven to be beneficial for cardiovascular and respiratory dynamics in healthy subjects and patients with various diseases. Therefore, the present study aimed to explore the cardiorespiratory dynamics by linear and nonlinear analysis of photoplethysmographic and respiratory time series on COVID-19 survivors under a psychophysiological assessment that includes slow-paced breathing. We analyzed photoplethysmographic and respiratory signals of 49 COVID-19 survivors to assess breathing rate variability (BRV), pulse rate variability (PRV), and pulse-respiration quotient (PRQ) during a psychophysiological assessment. Additionally, a comorbidity-based analysis was conducted to evaluate group changes. Our results indicate that all BRV indices significantly differed when performing slow-paced breathing. Nonlinear parameters of PRV were more appropriate for identifying changes in breathing patterns than linear indices. Furthermore, the mean and standard deviation of PRQ exhibited a significant increase while sample and fuzzy entropies decreased during diaphragmatic breathing. Thus, our findings suggest that slow-paced breathing may improve the cardiorespiratory dynamics of COVID-19 survivors in the short term by enhancing cardiorespiratory coupling via increased vagal activity.

ActionPytential: An open source tool for analyzing and visualizing cardiac action potential data

The action potential forms the basis of cardiac pacemaking, conduction, and contraction. Action potentials can be recorded from numerous preparation types, including ventricular or atrial trabecules, Purkinje fibers, isolated cardiac myocytes. Numerous techniques are also available as well, such as the conventional microelectrode and the single-cell current clamp techniques, optical mapping, or in silico modeling. With such a vast array of electrophysiological methods comes an array of available hardware and software solutions. In this work, we present a software with an intuitive graphical user interface, ActionPytential, that enables the analysis of any type of cardiac action potential, regardless of acquisition method or tissue type. In most available software tools, the analysis of continuous (gap-free) recordings often requires manual user interaction to segment the individual action potentials. We provide an automated solution for this, both for slow-response and for externally paced action potentials. As of now, ActionPytential calculates 34 parameters from each action potential. The most often utilized ones, including amplitude, maximal rate of depolarization, and action potential duration values, were validated on 1200 action potentials from human, dog, rabbit, guinea pig, and rat cardiac preparations. We also provide new parameters that were previously only measurable manually, including the position and the depth of the notch in potentials showing a spike-and-dome morphology. Further notable features include a Butterworth-type low-pass filter, the averaging of multiple potentials, automated corrections for baseline drifting, aided manual analysis, high-quality plots, and batch processing for any number of potentials. ActionPytential is available for all major platforms (Windows, MacOS, GNU + Linux, BSD).

Heart rate fragmentation is impaired in type 2 diabetes mellitus patients

AIMS

To evaluate the heart rate fragmentation (HRF) of type 2 diabetes mellitus (T2DM) patients and its relationship with heart rate variability (HRV) indices.

METHODS

One hundred sixty-four men, aged 47-57 years were retrospectively analyzed from a database. Participants were T2DM (n = 82) and apparently healthy (n = 82). R-R interval time series recorded by electrocardiogram were collected at the supine position for 10 to 15 min. From HRF, the percentage of inflection points (PIP), percentage of words with zero, one, two, or three inflections points (W₀, W₁, W₂, W₃), and percentage with only type hard, soft, or mixed inflections points type (W^H, W^S, W^M) were quantified.

RESULTS

T2DM presented higher PIP, W^S, W^M and W₃, while W^H and W₁ was lower compared with healthy (p < 0.05). Moreover, a positive moderate correlation was found between W^H and root mean square of the successive R-R differences (RMSSD) and high frequency (HF) indices. In contrast, a negative moderate correlation was found between W^S and W^M with RMSSD and HF indices.

CONCLUSIONS

T2DM have increased fragmentation patterns, and words grouped by inflection type are more closely related to HRV. The HRF approach might be useful to assess heart rate dynamic abnormalities in males with type 2 diabetes.

Multiscale asymmetry reveals changes in the maternal short-term heart rate dynamics of preeclamptic women during labor

BACKGROUND

It is known that acceleration and deceleration patterns in heart rate variability (HRV) are asymmetrically distributed in healthy subjects. Accordingly, novel approaches for assessing the asymmetrical properties of HRV, such as the multiscale asymmetry (MSA), have been applied in the perinatal field.

OBJECTIVE

To study the asymmetry of accelerations and decelerations of maternal short-term cardiac dynamics of thirty-six normotensive and preeclamptic women during labor/nonlabor by MSA analysis.

METHODS

The RR interval time series obtained from these participants were classified into four groups: normotensive (control) without labor C-NL, n= 10; control with labor C-L, n= 10; and two preeclamptic groups with absence or presence of labor P-NL, n= 6; and P-L, n= 10, respectively. Multiscale indices of heart rate asymmetry (HRA) such as Porta (P%), Guzik (G%) and Ehlers (E) were used to explore the changes of HRA in the normotensive and preeclamptic groups in the presence or absence of labor.

RESULTS

The main result of this study shows that preeclamptic women manifest decreased magnitude of decelerations of heart rate dynamics compared to normotensive women indicated by G% and E. We speculate that a lower cardiac parasympathetic response may be manifested in preeclamptic women during labor/nonlabor compared to normotensive women.

CONCLUSIONS

These observations represented a new insight into the autonomic cardiovascular regulation in preeclampsia, which could contribute to the perinatal field in the future.

Heart Rate Variability in Individuals with Down Syndrome: A Scoping Review with Methodological Considerations

Individuals with Down syndrome (DS) present similar heart rate variability (HRV) parameters at rest but different responses to selected movement maneuvers in comparison to individuals without DS, which indicates reduced vagal regulation. The present study undertakes a scoping review of research on HRV in individuals with DS, with special attention paid to the compliance of the studies with standards and methodological paper guidelines for HRV assessment and interpretation. A review was performed using PubMed, Web of Science and CINAHL databases to search for English language publications from 1996 to 2020 with the MESH terms "heart rate variability" and "down syndrome", with the additional inclusion criteria of including only human participants and empirical investigations. From 74 studies, 15 were included in the review. None of the reviewed studies met the recommendations laid out by the standards and guidelines for providing the acquisition of RR intervals and necessary details on HRV analysis. Since authors publishing papers on this research topic do not adhere to the prescribed standards and guidelines when constructing the methodology, results of the research papers on the topic are not directly comparable. Authors need to design the study methodology more robustly by following the aforementioned standards, guidelines and recommendations.

Heart rate variability as a biomarker in patients with Chronic Chagas Cardiomyopathy with or without concomitant digestive involvement and its relationship with the Rassi score

Background

Dysautonomia plays an ancillary role in the pathogenesis of Chronic Chagas Cardiomyopathy (CCC), but is the key factor causing digestive organic involvement. We investigated the ability of heart rate variability (HRV) for death risk stratification in CCC and compared alterations of HRV in patients with isolated CCC and in those with the mixed form (CCC + digestive involvement). Thirty-one patients with CCC were classified into three risk groups (low, intermediate and high) according to their Rassi score. A single-lead ECG was recorded for a period of 10–20 min, RR series were generated and 31 HRV indices were calculated. The HRV was compared among the three risk groups and regarding the associated digestive involvement. Four machine learning models were created to predict the risk class of patients.

Results

Phase entropy is decreased and the percentage of inflection points is increased in patients from the high-, compared to the low-risk group. Fourteen patients had the mixed form, showing decreased triangular interpolation of the RR histogram and absolute power at the low-frequency band. The best predictive risk model was obtained by the support vector machine algorithm (overall F1-score of 0.61).

Conclusions

The mixed form of Chagas' disease showed a decrease in the slow HRV components. The worst prognosis in CCC is associated with increased heart rate fragmentation. The combination of HRV indices enhanced the accuracy of risk stratification. In patients with the mixed form of Chagas disease, a higher degree of sympathetic autonomic denervation may be associated with parasympathetic impairment.

Analysis of Short-Term Heart Rate Asymmetry in High-Performance Athletes and Non-Athletes

Heart rate asymmetry (HRA) refers to how asymmetrically the acceleration and deceleration patterns in heartbeat fluctuations are distributed. There is limited evidence regarding HRA changes in athletes and their association with autonomic regulation. This study aimed to compare the short-term HRA of high-performance athletes and non-athletes during an autonomic function test by calculating relevant HRA measures. This exploratory study obtained beat-to-beat RR interval time series from 15 high-performance athletes and 12 non-athletes during a standardized autonomic function test. This test includes rest, postural change, controlled respiration, prolonged orthostatism, exercise, and recovery phases. The following HRA parameters were computed from the RR time series for both groups: asymmetric spread index (ASI), slope index (SI), Porta’s index (PI), Guzik’s index (GI), and Ehlers’ index (EI). We found significant differences (p < 0.01) in the mean value of several HRA parameters between athletes and non-athletes and across the autonomic function test phases, mainly in postural change and recovery phases. Our results indicate that high-performance athletes manifest a higher number and magnitude of cardiac decelerations than non-athletes after an orthostatic challenge, as indicated by GI and EI. In addition, lower HRA was found in athletes in the recovery phase than in non-athletes, as indicated by ASI.

Bitter Orange (Citrus aurantium L.) Intake Before Submaximal Aerobic Exercise Is Safe for Cardiovascular and Autonomic Systems in Healthy Males: A Randomized Trial

Background

There are still no studies of the cardiovascular safety of the isolated use of Citrus aurantium in aerobic submaximal exercise.

Objective

To evaluate the effect of C. aurantium supplementation on the recovery of cardiorespiratory and autonomic parameters following a session of submaximal aerobic exercise.

Methods

Twelve healthy male adults achieved a crossover, randomized, double-blind, and placebo-controlled trial. C. aurantium (600 mg, p-synephrine at 30% amount [180 mg]) or placebo (600 mg of starch) were ingested 90 min before evaluation in randomized days. We evaluated systolic blood pressure (SBP), diastolic blood pressure (DBP), pulse pressure (PP), mean arterial pressure (MAP), heart rate (HR) and, HR variability indexes at Rest and during 60 min of recovery from exercise.

Results

Citrus aurantium ingestion accelerated the reduction in SBP after exercise, anticipated the return of vagal modulation of the heart after exercise via the HF (ms2), pNN50 (%), and 2 UV% indices. Moreover, rushed the output of sympathetic modulation after exercise via the 0V% index. No unfavorable cardiovascular effects were achieved for HR, DBP, PP, and MAP parameters.

Conclusions

Citrus aurantium was shown to be safe for the cardiovascular and autonomic systems alongside submaximal aerobic exercise in healthy males.

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.

Correlation between heart rate variability and polysomnography-derived scores of obstructive sleep apnea

Obstructive sleep apnea (OSA) is one of the most common sleep disorders and affects nearly a billion people worldwide. Furthermore, it is estimated that many patients with OSA are underdiagnosed, which contributes to the development of comorbidities, such as cardiac autonomic imbalance, leading to high cardiac risk. Heart rate variability (HRV) is a non-invasive, widely used approach to evaluating neural control of the heart. This study evaluates the relationship between HRV indices and the presence and severity of OSA. We hypothesize that HRV, especially the nonlinear methods, can serve as an easy-to-collect marker for OSA early risk stratification. Polysomnography (PSG) exams of 157 patients were classified into four groups: OSA-free (N = 26), OSA-mild (N = 39), OSA-moderate (N = 37), and OSA-severe (N = 55). The electrocardiogram was extracted from the PSG recordings, and a 15-min beat-by-beat series of RR intervals were generated every hour during the first 6 h of sleep. Linear and nonlinear HRV approaches were employed to calculate 32 indices of HRV. Specifically, time- and frequency-domain, symbolic analysis, entropy measures, heart rate fragmentation, acceleration and deceleration capacities, asymmetry measures, and fractal analysis. Results with indices of sympathovagal balance provided support to reinforce previous knowledge that patients with OSA have sympathetic overactivity. Nonlinear indices showed that HRV dynamics of patients with OSA display a loss of physiologic complexity that could contribute to their higher risk of development of cardiovascular disease. Moreover, many HRV indices were found to be linked with clinical scores of PSG. Therefore, a complete set of HRV indices, especially the ones obtained by the nonlinear approaches, can bring valuable information about the presence and severity of OSA, suggesting that HRV can be helpful for in a quick diagnosis of OSA, and supporting early interventions that could potentially reduce the development of comorbidities.

Differences in the Asymmetry of Beat-to-Beat Fetal Heart Rate Accelerations and Decelerations at Preterm and Term Active Labor

The fetal autonomic nervous system responds to uterine contractions during active labor as identified by changes in the accelerations and decelerations of fetal heart rate (FHR). Thus, this exploratory study aimed to characterize the asymmetry differences of beat-to-beat FHR accelerations and decelerations in preterm and term fetuses during active labor. In an observational study, we analyzed 10 min of fetal R-R series collected from women during active preterm labor (32–36 weeks of pregnancy, n = 17) and active term labor (38–40 weeks of pregnancy, n = 27). These data were used to calculate the Deceleration Reserve (DR), which is a novel parameter that quantifies the asymmetry of the average acceleration and deceleration capacity of the heart. In addition, relevant multiscale asymmetric indices of FHR were also computed. Lower values of DR, calculated with the input parameters of T = 50 and s = 10, were associated with labor occurring at the preterm condition (p = 0.0131). Multiscale asymmetry indices also confirmed significant (p < 0.05) differences in the asymmetry of FHR. Fetuses during moderate premature labor may experience more decaying R-R trends and a lower magnitude of decelerations compared to term fetuses. These differences of FHR dynamics might be related to the immaturity of the fetal cardiac autonomic nervous system as identified by this system response to the intense uterine activity at active labor.

EntropyHub: An open-source toolkit for entropic time series analysis

An increasing number of studies across many research fields from biomedical engineering to finance are employing measures of entropy to quantify the regularity, variability or randomness of time series and image data. Entropy, as it relates to information theory and dynamical systems theory, can be estimated in many ways, with newly developed methods being continuously introduced in the scientific literature. Despite the growing interest in entropic time series and image analysis, there is a shortage of validated, open-source software tools that enable researchers to apply these methods. To date, packages for performing entropy analysis are often run using graphical user interfaces, lack the necessary supporting documentation, or do not include functions for more advanced entropy methods, such as cross-entropy, multiscale cross-entropy or bidimensional entropy. In light of this, this paper introduces EntropyHub, an open-source toolkit for performing entropic time series analysis in MATLAB, Python and Julia. EntropyHub (version 0.1) provides an extensive range of more than forty functions for estimating cross-, multiscale, multiscale cross-, and bidimensional entropy, each including a number of keyword arguments that allows the user to specify multiple parameters in the entropy calculation. Instructions for installation, descriptions of function syntax, and examples of use are fully detailed in the supporting documentation, available on the EntropyHub website- www.EntropyHub.xyz. Compatible with Windows, Mac and Linux operating systems, EntropyHub is hosted on GitHub, as well as the native package repository for MATLAB, Python and Julia, respectively. The goal of EntropyHub is to integrate the many established entropy methods into one complete resource, providing tools that make advanced entropic time series analysis straightforward and reproducible.

Differences in heart rate variability and body composition in breast cancer survivors and women without cancer

The aim of this study was to explore cardiac autonomic changes assessed by linear and nonlinear indexes of heart rate variability (HRV) and body composition modifications in breast cancer survivors and cancer-free control women. Women who were breast cancer survivors (BCS, n = 27) and without cancer with similar characteristics (Control, n = 31) were recruited for this study. We calculated some relevant linear and nonlinear parameters of 5 min of RR interval time series such as mean RR interval (RRave), the corrected Poincaré index (cSD1/SD2), the sample entropy (SampEn), the long-term fractal scaling exponent (α2) and 2UV from symbolic dynamics. Additionally, we indirectly assessed body composition measures such as body weight, fat mass, visceral fat rating (VFR), normalized VRF (nVFR), muscle mass, metabolic age, and total body water. We found that diverse HRV indexes and only one body composition measure showed statistical differences (p < 0.05) between the BCS and Control groups. RRave: 729 (648-802) vs. 795 (713-852) ms; cSD2/SD1: 3.4 (2.7-5.0) vs. 2.9 (2.3-3.5); SampEn: 1.5 (1.3-1.8) vs. 1.7 (1.5-1.8); α2: 0.6 (0.3-0.6) vs. 0.5 (0.4-0.5); 2UV: 7.1 (4.3-11.5) vs. 10.8 (6.4-15.7) and nVFR 0.12 (0.11-0.13) vs. 0.10 (0.08-0.12) points/kg, respectively. The nVFR was strongly significantly correlated with several indexes of HRV only in the BCS group.Our findings suggest that BCS exhibit lower parasympathetic cardiac activity and changes in HRV patterns compared to Controls. A concomitant increase of visceral fat, among other factors, may contribute to cardiac autonomic disturbances and changes in HRV patterns in BCS.

CEPS: An Open Access MATLAB Graphical User Interface (GUI) for the Analysis of Complexity and Entropy in Physiological Signals

BACKGROUND

We developed CEPS as an open access MATLAB^® GUI (graphical user interface) for the analysis of Complexity and Entropy in Physiological Signals (CEPS), and demonstrate its use with an example data set that shows the effects of paced breathing (PB) on variability of heart, pulse and respiration rates. CEPS is also sufficiently adaptable to be used for other time series physiological data such as EEG (electroencephalography), postural sway or temperature measurements.

METHODS

Data were collected from a convenience sample of nine healthy adults in a pilot for a larger study investigating the effects on vagal tone of breathing paced at various different rates, part of a development programme for a home training stress reduction system.

RESULTS

The current version of CEPS focuses on those complexity and entropy measures that appear most frequently in the literature, together with some recently introduced entropy measures which may have advantages over those that are more established. Ten methods of estimating data complexity are currently included, and some 28 entropy measures. The GUI also includes a section for data pre-processing and standard ancillary methods to enable parameter estimation of embedding dimension m and time delay τ ('tau') where required. The software is freely available under version 3 of the GNU Lesser General Public License (LGPLv3) for non-commercial users. CEPS can be downloaded from Bitbucket. In our illustration on PB, most complexity and entropy measures decreased significantly in response to breathing at 7 breaths per minute, differentiating more clearly than conventional linear, time- and frequency-domain measures between breathing states. In contrast, Higuchi fractal dimension increased during paced breathing.

CONCLUSIONS

We have developed CEPS software as a physiological data visualiser able to integrate state of the art techniques. The interface is designed for clinical research and has a structure designed for integrating new tools. The aim is to strengthen collaboration between clinicians and the biomedical community, as demonstrated here by using CEPS to analyse various physiological responses to paced breathing.