Discovery of Causal Paths in Cardiorespiratory Parameters: A Time-Independent Approach in Elite Athletes

Dual-task improvement of older adults after treadmill walking combined with blood flow restriction of low occlusion pressure: the effect on the heart-brain axis

OBJECTIVE

This study explored the impact of one session of low-pressure leg blood flow restriction (BFR) during treadmill walking on dual-task performance in older adults using the neurovisceral integration model framework.

METHODS

Twenty-seven older adults participated in 20-min treadmill sessions, either with BFR (100 mmHg cuff pressure on both thighs) or without it (NBFR). Dual-task performance, measured through light-pod tapping while standing on foam, and heart rate variability during treadmill walking were compared.

RESULTS

Following BFR treadmill walking, the reaction time (p = 0.002) and sway area (p = 0.012) of the posture dual-task were significantly reduced. Participants exhibited a lower mean heart rate (p < 0.001) and higher heart rate variability (p = 0.038) during BFR treadmill walking. Notably, BFR also led to band-specific reductions in regional brain activities (theta, alpha, and beta bands, p < 0.05). The topology of the EEG network in the theta and alpha bands became more star-like in the post-test after BFR treadmill walking (p < 0.005).

CONCLUSION

BFR treadmill walking improves dual-task performance in older adults via vagally-mediated network integration with superior neural economy. This approach has the potential to prevent age-related falls by promoting cognitive reserves.

Causal Inference for Hypertension Prediction

Hypertension is a leading cause of cardiovascular disease and premature death worldwide and it puts a heavy burden on the healthcare system. It is, therefore, very important to detect and evaluate hypertension and related cardiovascular events as to for efficient diagnosis, treatment and management. Hypertension can be evaluated with noninvasive cardiac signals, such as electrocardiogram (ECG) and photoplethysmogram (PPG) signals. Most of the previous studies predicted hypertension from ECG and PPG signals with extracted features that are correlated with hypertension. However, correlation is sometimes unreliable and may be affected by confounding factors. In this study, we propose a causal inference based approach to identify feature variables from ECG and PPG signals that are potentially causally related with hypertension. The method of greedy equivalence search was employed to construct the causal graph of features and hypertension. With causal features identified from the causal graph, we used machine learning models to diagnose hypertension. The machine learning classification models achieve great classification performance, among which random forest model has the best classification performance, with accuracy being 0.987, precision being 0.990, recall being 0.981, and F1-score being 0.985. The results show that the causal inference based approach can effectively predict hyper-tension.Clinical relevance- This paper proposes a new hypertension risk prediction method, which uses causality instead of correlation as the feature screening criteria to establish a causal graph of hypertension, which can predict the hypertension more reliably.

Factors Associated with Reduced Heart Rate Variability in the General Japanese Population: The Iwaki Cross-Sectional Research Study

Although many studies have reported factors associated with reduced heart rate variability (HRV) in Western populations, evidence is limited among Asian populations. Therefore, we investigated the factors associated with reduced HRV values in a general Japanese population by measuring HRV among the participants of the Iwaki Health Promotion Project who underwent medical examination in 2019. We performed 90-s HRV measurements in 1065 participants. Of these, we evaluated the coefficient of variation in R–R intervals (CVRR) and standard deviation in R–R intervals (SDNN). Blood was collected under a fasting condition, and investigations of glucose metabolism, lipid metabolism, renal function, liver function, advanced glycation end products, and blood pressure were performed. A multivariate regression analysis of the association between CVRR or SDNN and blood test parameters and blood pressure in 987 participants with adequately completed HRV assessments showed that reduced CVRR or SDNN was associated with higher levels of glycated hemoglobin (HbA1c), glycoalbumin, blood glucose, triglycerides, creatinine, plasma pentosidine, and diastolic blood pressure. In the general Japanese population, higher levels of HbA1c, glycoalbumin, blood glucose, triglycerides, creatinine, plasma pentosidine, and diastolic blood pressure are associated with reduced CVRR or SDNN, which are typical HRV parameters.

Granger causality test with nonlinear neural-network-based methods: Python package and simulation study

BACKGROUND AND OBJECTIVE

Causality defined by Granger in 1969 is a widely used concept, particularly in neuroscience and economics. As there is an increasing interest in nonlinear causality research, a Python package with a neural-network-based causality analysis approach was created. It allows performing causality tests using neural networks based on Long Short-Term Memory (LSTM), Gated Recurrent Unit (GRU), or Multilayer Perceptron (MLP). The aim of this paper is to present the nonlinear method for causality analysis and the created Python package.

METHODS

The created functions with the autoregressive (AR) and Generalized Radial Basis Functions (GRBF) neural network models were tested on simulated signals in two cases: with nonlinear dependency and with absence of causality from Y to X signal. The train-test split (70/30) was used. Errors obtained on the test set were compared using the Wilcoxon signed-rank test to determine the presence of the causality. For the chosen model, the proposed method of study the change of causality over time was presented.

RESULTS

In the case when X was a polynomial of Y, nonlinear methods were able to detect the causality, while the AR model did not manage to indicate it. The best results (in terms of the prediction accuracy) were obtained for the MLP for the lag of 150 (MSE equal to 0.011, compared to 0.041 and 0.036 for AR and GRBF, respectively). When there was no causality between the signals, none of the proposed and AR models did indicate false causality, while it was detected by GRBF models in one case. Only the proposed models gave the expected results in each of the tested scenarios.

CONCLUSIONS

The proposed method appeared to be superior to the compared methods. They were able to detect non-linear causality, make accurate forecasting and not indicate false causality. The created package enables easy usage of neural networks to study the causal relationship between signals. The neural-networks-based approach is a suitable method that allows the detection of a nonlinear causal relationship, which cannot be detected by the classical Granger method. Unlike other similar tools, the package allows for the study of changes in causality over time.

Discovering Ecological Relationships in Flowing Freshwater Ecosystems

Knowledge of ecological responses to changes in the environment is vital to design appropriate measures for conserving biodiversity. Experimental studies are the standard to identify ecological cause-effect relationships, but their results do not necessarily translate to field situations. Deriving ecological cause-effect relationships from observational field data is, however, challenging due to potential confounding influences of unmeasured variables. Here, we present a causal discovery algorithm designed to reveal ecological relationships in rivers and streams from observational data. Our algorithm (a) takes into account the spatial structure of the river network, (b) reveals the complete network of ecological relationships, and (c) shows the directions of these relationships. We apply our algorithm to data collected in the US state of Ohio to better understand causes of reductions in fish and invertebrate community integrity. We found that nitrogen is a key variable underlying fish and invertebrate community integrity in Ohio, likely negatively impacting both. We also found that fish and community integrity are each linked to one physical habitat quality variable. Our algorithm further revealed a split between physical habitat quality and water quality variables, indicating that causal relations between these groups of variables are likely absent. Our approach is able to reveal networks of ecological relationships in rivers and streams based on observational data, without the need to formulate a priori hypotheses. This is an asset particularly for diagnostic assessments of the ecological state and potential causes of biodiversity impairment in rivers and streams.

Respiratory Activity during Exercise: A Feasibility Study on Transition Point Estimation Using Impedance Pneumography

The current diagnostic procedures for assessing physiological response to exercise comprise blood lactates measurements, ergospirometry, and electrocardiography. The first is not continuous, the second requires specialized equipment distorting natural breathing, and the last is indirect. Therefore, we decided to perform the feasibility study with impedance pneumography as an alternative technique. We attempted to determine points in respiratory-related signals, acquired during stress test conditions, that suggest a transition similar to the gas exchange threshold. In addition, we analyzed whether or not respiratory activity reaches steady states during graded exercise. Forty-four students (35 females), practicing sports on different levels, performed a graded exercise test until exhaustion on cycloergometer. Eventually, the results from 34 of them were used. The data were acquired with Pneumonitor 2. The signals demonstrated that the steady state phenomenon is not as evident as for heart rate. The results indicated respiratory rate approaches show the transition point at the earliest (more than 6 min before the end of the exercise test on average), and the tidal volume ones at the latest (less than 5 min). A combination gave intermediate findings. The results showed the impedance pneumography appears reasonable for the transition point estimation, but this should be further studied with the reference.

Changes in Short-Term and Ultra-Short Term Heart Rate, Respiratory Rate, and Time-Domain Heart Rate Variability Parameters during Sympathetic Nervous System Activity Stimulation in Elite Modern Pentathlonists-A Pilot Study

Monitoring of markers reflecting cardiac autonomic activity before and during stressful situations may be useful for identifying the physiological state of an athlete and may have medical or performance implications. The study aimed to determine group and individual changes in short-term (5 min) and ultra-short-term (1 min) heart rate (HR), respiratory rate (RespRate), and time-domain heart rate variability (HRV) parameters during sympathetic nervous system activity (SNSa) stimulation among professional endurance athletes. Electrocardiographic recordings were performed in stable measurement conditions (Baseline) and during SNSa stimulation via isometric handgrip in 12 elite modern pentathlonists. Significant increases in short-term HR and decreases in time-domain HRV parameters with no changes in RespRate were observed during SNSa stimulation. Significant differences were observed between Baseline (all minutes) and the last (i.e., 5th) minute of SNSa stimulation for ultra-short-term parameters. Analysis of intra-individual changes revealed some heterogeneity in responses. The study provides baseline responses of HR, RespRate, and time-domain HRV parameters to SNSa stimulation among elite pentathlonists, which may be useful for identifying abnormal responses among fatigued or injured (e.g., concussed) athletes. More attention to individual analysis seems to be necessary when assessing physiological responses to sympathetic stimuli in professional endurance athletes.

A Pilot Study of the Reliability and Agreement of Heart Rate, Respiratory Rate and Short-Term Heart Rate Variability in Elite Modern Pentathlon Athletes

Research on reliability of heart rate variability (HRV) parameters in athletes has received increasing attention. The aims of this study were to examine the inter-day reliability of short-term (5 min) and ultra-short-term (1 min) heart rate (HR), respiratory rate (RespRate) and HRV parameters, agreement between short-term and ultra-short-term parameters, and association between differences in HR, RespRate and HRV parameters in elite modern pentathletes. Electrocardiographic recordings were performed in stable measurement conditions with a week interval between tests. Relative reliability was evaluated by intra-class correlation coefficients, absolute reliability was evaluated by within-subject coefficient of variation, and agreement was evaluated using Bland-Altman (BA) plot with limits of agreement and defined a priori maximum acceptable difference. Short-term HR, RespRate, log transformed (ln) root mean square of successive normal-to-normal interval differences (lnRMSSD), ln high frequency (lnHF) and SD2/SD1 HRV indices and ultra-short-term HR, RespRate and lnRMSSD presented acceptable, satisfactory inter-day reliability. Although there were no significant differences between short-term and ultra-short-term HR, RespRate and lnRMSSD, no parameter showed acceptable differences with BA plots. Differences in time-domain and non-linear HRV parameters were more correlated with differences in HR than with differences in RespRate. Inverse results were observed for frequency-domain parameters. Short-term HR, RespRate, lnRMSSD, lnHF, and SD2/SD1 and ultra-short-term HR, RespRate and lnRMSSD could be used as reliable parameters in endurance athletes. However, practitioners should interpret changes in HRV parameters with regard to concomitant differences in HR and RespRate and caution should be taken before considering 5 min and 1 min parameters as interchangeable.

Effect of Acute Hypoxia on Cardiorespiratory Coherence in Male Runners

Understanding the mechanisms of oxygen supply regulation, which involves the respiratory and cardiovascular systems, during human adaptation to intense physical activity, accompanied by hypoxemia, is important for the management of a training process. The objectives of this study were to investigate the cardiorespiratory coherence (CRC) changes in the low-frequency band in response to hypoxic exposure and to verify a dependence of these changes upon sports qualification level in athletes. Twenty male runners aged 17-25 years were exposed to acute normobaric hypoxia (10% O2) for 10 min. Respiration, gas exchange, and heart rate were measured at baseline, during hypoxia, and after the exposure. To evaluate cardiorespiratory coupling, squared coherence was calculated based on 5-s averaged time series of heart and respiratory rhythms. Based on sports qualification level achieved over 4 years after the experimental testing, athletes were retrospectively divided into two groups, one high level (HLG, n = 10) and the other middle level (MLG, n = 10). No differences in anthropometric traits were observed between the groups. In the pooled group, acute hypoxia significantly increased CRC at frequencies 0.030-0.045 Hz and 0.075 Hz. In response to hypoxia, oxygen consumption decreased in HLG, and carbon dioxide production and ventilation increased in MLG. At 0.070-0.080 Hz frequencies in hypoxia, the CRC in HLG was higher than in MLG. Thus, highly qualified athletes enhance intersystem integration in response to hypoxia. This finding can be a physiological sign for the prognosis of qualification level in runners.

From Real-World Patient Data to Individualized Treatment Effects Using Machine Learning: Current and Future Methods to Address Underlying Challenges

Clinical decision making needs to be supported by evidence that treatments are beneficial to individual patients. Although randomized control trials (RCTs) are the gold standard for testing and introducing new drugs, due to the focus on specific questions with respect to establishing efficacy and safety vs. standard treatment, they do not provide a full characterization of the heterogeneity in the final intended treatment population. Conversely, real-world observational data, such as electronic health records (EHRs), contain large amounts of clinical information about heterogeneous patients and their response to treatments. In this paper, we introduce the main opportunities and challenges in using observational data for training machine learning methods to estimate individualized treatment effects and make treatment recommendations. We describe the modeling choices of the state-of-the-art machine learning methods for causal inference, developed for estimating treatment effects both in the cross-section and longitudinal settings. Additionally, we highlight future research directions that could lead to achieving the full potential of leveraging EHRs and machine learning for making individualized treatment recommendations. We also discuss how experimental data from RCTs and Pharmacometric and Quantitative Systems Pharmacology approaches can be used to not only improve machine learning methods, but also provide ways for validating them. These future research directions will require us to collaborate across the scientific disciplines to incorporate models based on RCTs and known disease processes, physiology, and pharmacology into these machine learning models based on EHRs to fully optimize the opportunity these data present.

Slow 0.1 Hz Breathing and Body Posture Induced Perturbations of RRI and Respiratory Signal Complexity and Cardiorespiratory Coupling

Objective: We explored the physiological background of the non-linear operating mode of cardiorespiratory oscillators as the fundamental question of cardiorespiratory homeodynamics and as a prerequisite for the understanding of neurocardiovascular diseases. We investigated 20 healthy human subjects for changes using electrocardiac RR interval (RRI) and respiratory signal (Resp) Detrended Fluctuation Analysis (DFA, α1RRI, α2RRI, α1Resp, α2Resp), Multiple Scaling Entropy (MSERRI1-4, MSERRI5-10, MSEResp1-4, MSEResp5-10), spectral coherence (CohRRI-Resp), cross DFA (ρ1 and ρ2) and cross MSE (XMSE1-4 and XMSE5-10) indices in four physiological conditions: supine with spontaneous breathing, standing with spontaneous breathing, supine with 0.1 Hz breathing and standing with 0.1 Hz breathing. Main results: Standing is primarily characterized by the change of RRI parameters, insensitivity to change with respiratory parameters, decrease of CohRRI-Resp and insensitivity to change of in ρ1, ρ2, XMSE1-4, and XMSE5-10. Slow breathing in supine position was characterized by the change of the linear and non-linear parameters of both signals, reflecting the dominant vagal RRI modulation and the impact of slow 0.1 Hz breathing on Resp parameters. CohRRI-Resp did not change with respect to supine position, while ρ1 increased. Slow breathing in standing reflected the qualitatively specific state of autonomic regulation with striking impact on both cardiac and respiratory parameters, with specific patterns of cardiorespiratory coupling. Significance: Our results show that cardiac and respiratory short term and long term complexity parameters have different, state dependent patterns. Sympathovagal non-linear interactions are dependent on the pattern of their activation, having different scaling properties when individually activated with respect to the state of their joint activation. All investigated states induced a change of α1 vs. α2 relationship, which can be accurately expressed by the proposed measure-inter-fractal angle θ. Short scale (α1 vs. MSE1-4) and long scale (α2 vs. MSE5-10) complexity measures had reciprocal interrelation in standing with 0.1 Hz breathing, with specific cardiorespiratory coupling pattern (ρ1 vs. XMSE1-4). These results support the hypothesis of hierarchical organization of cardiorespiratory complexity mechanisms and their recruitment in ascendant manner with respect to the increase of behavioral challenge complexity. Specific and comprehensive cardiorespiratory regulation in standing with 0.1 Hz breathing suggests this state as the potentially most beneficial maneuver for cardiorespiratory conditioning.

Cardiorespiratory Temporal Causal Links and the Differences by Sport or Lack Thereof

Fitness level, fatigue and adaptation are important factors for determining the optimal training schedule and predicting future performance. We think that adding analysis of the mutual relationships between cardiac and respiratory activity enables better athlete profiling and feedback for improving training. Therefore, the main objectives were (1) to apply several methods for temporal causality analysis to cardiorespiratory data; (2) to establish causal links between the signals; and (3) to determine how parameterized connections differed across various subgroups. One hundred elite athletes (31 female) and a control group of 20 healthy students (6 female) took part in the study. All were asked to follow a protocol comprising two 5-min sessions of free breathing - once supine, once standing. The data were collected using Pneumonitor 2. Respiratory-related curves were obtained through impedance pneumography, along with a single-lead ECG. Several signals (e.g., tidal volume, instantaneous respiratory rate, and instantaneous heart rate) were derived and stored as: (1) raw data down-sampled to 25Hz; (2) further down-sampled to 2.5Hz; and (3) beat-by-beat sequences. Granger causality frameworks (pairwise-conditional, spectral or extended), along with Time Series Models with Independent Noise (TiMINo), were studied. The connections enabling the best distinctions were found using recursive feature elimination with a random forest kernel. Temporal causal links are the most evident between tidal volume and instantaneous heart rate signals. Predictions of the “effect” variable were improved by adding preceding “cause” samples, by medians of 20.3% for supine and 14.2% for standing body positions. Parameterized causal link structures and directions distinguish athletes from non-athletes with 83.3% accuracy on average. They may also be used to supplement standard analysis and enable classification into groups exhibiting different static and dynamic components during performance. Physiological markers of training may be extended to include cardiorespiratory data, and causality analysis may improve the resolution of training profiling and the precision of outcome prediction.