Recurrence quantification analysis of heart rate variability during continuous incremental exercise test in obese subjects

Agreement Between Heart Rate Variability - Derived vs. Ventilatory and Lactate Thresholds: A Systematic Review with Meta-Analyses

BACKGROUND

Determining thresholds by measuring blood lactate levels (lactate thresholds) or gas exchange (ventilatory thresholds) that delineate the different exercise intensity domains is crucial for training prescription. This systematic review with meta-analyses aims to assess the overall validity of the first and second heart rate variability - derived threshold (HRVT1 and HRVT2, respectively) by computing global effect sizes for agreement and correlation between HRVTs and reference - lactate and ventilatory (LT-VTs) - thresholds. Furthermore, this review aims to assess the impact of subjects' characteristics, HRV methods, and study protocols on the agreement and correlation between LT-VTs and HRVTs.

METHODS

Systematic computerised searches for studies determining HRVTs during incremental exercise in humans were conducted. The agreements and correlations meta-analyses were conducted using a random-effect model. Causes of heterogeneity were explored by subgroup analysis and meta-regression with subjects' characteristics, incremental exercise protocols, and HRV methods variables. The methodological quality was assessed using QUADAS-2 and STARDHRV tools. The risk of bias was assessed by funnel plots, fail-safe N test, Egger's test of the intercept, and the Begg and Mazumdar rank correlation test.

RESULTS

Fifty included studies (1160 subjects) assessed 314 agreements (95 for HRVT1, 219 for HRVT2) and 246 correlations (82 for HRVT1, 164 for HRVT2) between LT-VTs and HRVTs. The standardized mean differences were trivial between HRVT1 and LT1-VT1 (SMD = 0.08, 95% CI -0.04-0.19, n = 22) and between HRVT2 and LT2-VT2 (SMD = -0.06, 95% CI -0.15-0.03, n = 42). The correlations were very strong between HRVT1 and LT1-VT1 (r = 0.85, 95% CI 0.75-0.91, n = 22), and between HRVT2 and LT2-VT2 (r = 0.85, 95% CI 0.80-0.89, n = 41). Moreover, subjects' characteristics, type of ergometer, or initial and incremental workload had no impact on HRVTs determination.

CONCLUSION

HRVTs showed trivial differences and very strong correlations with LT-VTs and might thus serve as surrogates. These results emphasize the usefulness of HRVTs as promising, accessible, and cost-effective means for exercise and clinical prescription purposes.

Heart Rate Variability Thresholds: Agreement with Established Approaches and Reproducibility in Trained Females and Males

PURPOSE

To determine in trained females and males i) the agreement between the gas exchange threshold (GET), lactate threshold 1 (LT1), and heart rate variability threshold 1 (HRVT1), as well as between the respiratory compensation point (RCP), lactate threshold 2 (LT2), and heart rate variability threshold 2 (HRVT2), and ii) the reproducibility of HRVT1 and HRVT2 during 2-min incremental step protocols.

METHODS

Fifty-seven trained participants (24 females) completed a 2-min step incremental test to task failure. Nineteen participants (eight females) completed a second test to evaluate reproducibility. Gas exchange and ventilatory responses, blood lactate concentration, and RR time series were recorded to assess the oxygen consumption (V̇O 2 ) and heart rate (HR) associated with the GET, RCP, LT1, LT2, HRVT1, and HRVT2.

RESULTS

V̇O 2 -GET versus V̇O 2 -HRVT1 and HR-GET versus HR-HRVT1 were statistically different for females (29.5 ± 4.0 vs 34.6 ± 6.1 mL·kg -1 ·min -1 ; 154 ± 11 vs 166 ± 12 bpm) and for males (33.9 ± 4.2 vs 42.7 ± 4.6 mL·kg -1 ·min -1 ; 145 ± 11 vs 165 ± 9 bpm; P < 0.001). V̇O 2 and HR at HRVT1 were greater than at LT1 ( P < 0.05). V̇O 2 -RCP versus V̇O 2 -HRVT2 and HR-RCP versus HR-HRVT2 were not statistically different for females (40.1 ± 4.7 vs 39.5 ± 6.7 mL·kg -1 ·min -1 ; 177 ± 9 vs 176 ± 9 bpm) and males (48.4 ± 5.4 vs 47.8 ± 4.8 mL·kg -1 ·min -1 ; 176 ± 8 vs 175 ± 9 bpm; P > 0.05). V̇O 2 and HR responses at LT2 were similar to HRVT2 ( P > 0.05). Intraclass correlation coefficient for V̇O 2 -HRVT1, HR-HRVT1, V̇O 2 -HRVT2, and HR-HRVT2 indicated good reproducibility when comparing the two different time points to standard methods.

CONCLUSIONS

Whereas HRVT2 is a valid and reproducible estimate of the RCP/LT2, current approaches for HRVT1 estimation did not show good agreement with outcomes at GET and LT1.

Heart Rate Variability-Derived Thresholds for Exercise Intensity Prescription in Endurance Sports: A Systematic Review of Interrelations and Agreement with Different Ventilatory and Blood Lactate Thresholds

BACKGROUND

Exercise intensities are prescribed using specific intensity zones (moderate, heavy, and severe) determined by a 'lower' and a 'higher' threshold. Typically, ventilatory (VT) or blood lactate thresholds (LT), and critical power/speed concepts (CP/CS) are used. Various heart rate variability-derived thresholds (HRVTs) using different HRV indices may constitute applicable alternatives, but a systematic review of the proximity of HRVTs to established threshold concepts is lacking.

OBJECTIVE

This systematic review aims to provide an overview of studies that determined HRVTs during endurance exercise in healthy adults in comparison with a reference VT and/or LT concept.

METHODS

A systematic literature search for studies determining HRVTs in healthy individuals during endurance exercise and comparing them with VTs or LTs was conducted in Scopus, PubMed and Web of Science (until January 2022). Studies claiming to describe similar physiological boundaries to delineate moderate from heavy (HRVTlow vs. VTlow and/or LTlow), and heavy from severe intensity zone (HRVThigh vs. VThigh and/or LThigh) were grouped and their results synthesized.

RESULTS

Twenty-seven included studies (461 participants) showed a mean difference in relative HR between HRVTlow and VTlow of - 0.6%bpm in weighted means and 0.02%bpm between HRVTlow and LTlow. Bias between HR at HRVTlow and VTlow was 1 bpm (limits of agreement (LoA): - 10.9 to 12.8 bpm) and 2.7 bpm (LoA: - 20.4 to 25.8 bpm) between HRVTlow and LTlow. Mean difference in HR between HRVThigh and VThigh was 0.3%bpm in weighted means and 2.9%bpm between HRVThigh and LThigh while bias between HR at HRVThigh and VThigh was - 4 bpm (LoA: - 17.9 to 9.9 bpm) and 2.5 bpm (LoA: - 12.1 to 17.1 bpm) between HRVThigh and LThigh.

CONCLUSION

HRVTlow seems to be a promising approach for the determination of a 'lower' threshold comparable to VTlow and potentially for HRVThigh compared to VThigh, although the latter needs further empirical evaluation. LoA for both intensity zone boundaries indicates bias of HRVTs on an individual level. Taken together, HRVTs can be a promising alternative for prescribing exercise intensity in healthy, male athletes undertaking endurance activities but due to the heterogeneity of study design, threshold concepts, standardization, and lack of female participants, further research is necessary to draw more robust and nuanced conclusions.

Smart Devices for Health and Wellness Applied to Tele-Exercise: An Overview of New Trends and Technologies Such as IoT and AI

This descriptive article explores the use of smart devices for health and wellness in the context of telehealth, highlighting rapidly evolving technologies such as the Internet of Things (IoT) and Artificial Intelligence (AI). Key innovations, benefits, challenges, and opportunities related to the adoption of these technologies are outlined. The article provides a descriptive and accessible approach to understanding the evolution and impact of smart devices in the tele-exercise reality. Nowadays, technological advances provide solutions that were unthinkable just a few years ago. The habits of the general population have also changed over the past few years. Hence, there is a need to investigate this issue and draw the attention of the scientific community to this topic by describing the benefits and challenges associated with each topic. If individuals no longer go to exercise, the exercise must go to their homes instead.

Automatic Detection of Aerobic Threshold through Recurrence Quantification Analysis of Heart Rate Time Series

During exercise with increasing intensity, the human body transforms energy with mechanisms dependent upon actual requirements. Three phases of the body's energy utilization are recognized, characterized by different metabolic processes, and separated by two threshold points, called aerobic (AerT) and anaerobic threshold (AnT). These thresholds occur at determined values of exercise intensity(workload) and can change among individuals. They are considered indicators of exercise capacities and are useful in the personalization of physical activity plans. They are usually detected by ventilatory or metabolic variables and require expensive equipment and invasive measurements. Recently, particular attention has focused on AerT, which is a parameter especially useful in the overweight and obese population to determine the best amount of exercise intensity for weight loss and increasing physical fitness. The aim of study is to propose a new procedure to automatically identify AerT using the analysis of recurrences (RQA) relying only on Heart rate time series, acquired from a cohort of young athletes during a sub-maximal incremental exercise test (Cardiopulmonary Exercise Test, CPET) on a cycle ergometer. We found that the minima of determinism, an RQA feature calculated from the Recurrence Quantification by Epochs (RQE) approach, identify the time points where generic metabolic transitions occur. Among these transitions, a criterion based on the maximum convexity of the determinism minima allows to detect the first metabolic threshold. The ordinary least products regression analysis shows that values of the oxygen consumption VO₂, heart rate (HR), and Workload correspondent to the AerT estimated by RQA are strongly correlated with the one estimated by CPET (r > 0.64). Mean percentage differences are <2% for both HR and VO₂ and <11% for Workload. The Technical Error for HR at AerT is <8%; intraclass correlation coefficients values are moderate (≥0.66) for all variables at AerT. This system thus represents a useful method to detect AerT relying only on heart rate time series, and once validated for different activities, in future, can be easily implemented in applications acquiring data from portable heart rate monitors.

Detecting Metabolic Thresholds from Nonlinear Analysis of Heart Rate Time Series: A Review

Heart rate time series are widely used to characterize physiological states and athletic performance. Among the main indicators of metabolic and physiological states, the detection of metabolic thresholds is an important tool in establishing training protocols in both sport and clinical fields. This paper reviews the most common methods, applied to heart rate (HR) time series, aiming to detect metabolic thresholds. These methodologies have been largely used to assess energy metabolism and to identify the appropriate intensity of physical exercise which can reduce body weight and improve physical fitness. Specifically, we focused on the main nonlinear signal evaluation methods using HR to identify metabolic thresholds with the purpose of identifying a method which can represent a useful tool for the real-time settings of wearable devices in sport activities. While the advantages and disadvantages of each method, and the possible applications, are presented, this review confirms that the nonlinear analysis of HR time series represents a solid, robust and noninvasive approach to assess metabolic thresholds.

Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers

BACKGROUND

Human motor systems contain nonlinear features. The purpose of this study was to evaluate the geometric structure of attractors and analyze recurrence in two different pirouettes (jazz and classic) performed by 15 professional dancers.

METHODS

The kinematics of the body's center of mass (CoM) and knee of the supporting leg (LKNE) during the pirouette were measured using the Vicon system. A time series of selected points were resampled, normalized, and randomly reordered. Then, every second time series was flipped to be combined with other time series and make a long time series out of the repetitions of a single task. The attractors were reconstructed, and the convex hull volumes (CHV) were counted for the CoM and LKNE for each pirouette in each direction. Recurrence quantification analysis (RQA) was used to extract additional information.

RESULTS

The CHVs calculated for the LKNE were significantly lower for the jazz pirouette. All RQA measures had the highest values for LKNE along the mediolateral axis for the jazz pirouette. This result underscores the high determinism, high motion recurrence, and complexity of this maneuver.

CONCLUSIONS

The findings offer new insight into the evaluation of the approximation of homogeneity in motion control. A high determinism indicates a highly stable and predictive motion trajectory.

Personalized Metabolic Avatar: A Data Driven Model of Metabolism for Weight Variation Forecasting and Diet Plan Evaluation

Development of predictive computational models of metabolism through mechanistic models is complex and resource demanding, and their personalization remains challenging. Data-driven models of human metabolism would constitute a reliable, fast, and continuously updating model for predictive analytics. Wearable devices, such as smart bands and impedance balances, allow the real time and remote monitoring of physiological parameters, providing for a flux of data carrying information on user metabolism. Here, we developed a data-driven model of end-user metabolism, the Personalized Metabolic Avatar (PMA), to estimate its personalized reactions to diets. PMA consists of a gated recurrent unit (GRU) deep learning model trained to forecast personalized weight variations according to macronutrient composition and daily energy balance. The model can perform simulations and evaluation of diet plans, allowing the definition of tailored goals for achieving ideal weight. This approach can provide the correct clues to empower citizens with scientific knowledge, augmenting their self-awareness with the aim to achieve long-lasting results in pursuing a healthy lifestyle.

Detection of hearing losses (HL) via transient-evoked otoacoustic emissions: towards an automatic classification

Transiently evoked otoacoustic emissions (TEOAEs) are routinely used in the hearing assessment of the auditory periphery. The major contribution of TEOAEs is the early detection of hearing losses in neonates, children, and adults. The evaluation of TEOAE responses by specific signal decomposition techniques offers numerous advantages for current and future research. One methodology, based on recurrence quantification analysis (RQA), can identify adult subjects presenting sensorineural hearing impairments. In two previous papers, the RQA-based approach was succesfully applied in identifying and classifying cases presenting noise and age related hearing losses. The current work investigates further two aspects of the previously proposed RQA-based analysis for hearing loss detection: (i) the reliability of a Training set built from different numbers of ears with normal hearing, and (ii) the threshold set of values of the key hearing loss detecting parameter RAD2D. Results: The Training set built from 158 healthy ears was found to be quite reliable and a similar but slightly minor performance was observed for the training set of 118 normal subjects, used in the past; the proposed ROC-curve method, optimizing the values of RAD2D, shows improved sensibility and specificity in one class discrimination. Conclusions: A complete and simplified procedure, based on the combined use of the traditional TEOAE reproducibility value and on values from the RQA-based RAD2D parameter, is proposed as an improved automatic classifier, in terms of sensitivity and specificity, for different types of hearing losses.

Unsupervised Clustering of Heartbeat Dynamics Allows for Real Time and Personalized Improvement in Cardiovascular Fitness

VO₂max index has a significant impact on overall health. Its estimation through wearables notifies the user of his level of fitness but cannot provide a detailed analysis of the time intervals in which heartbeat dynamics are changed and/or fatigue is emerging. Here, we developed a multiple modality biosignal processing method to investigate running sessions to characterize in real time heartbeat dynamics in response to external energy demand. We isolated dynamic regimes whose fraction increases with the VO₂max and with the emergence of neuromuscular fatigue. This analysis can be extremely valuable by providing personalized feedback about the user’s fitness level improvement that can be realized by developing personalized exercise plans aimed to target a contextual increase in the dynamic regime fraction related to VO₂max increase, at the expense of the dynamic regime fraction related to the emergence of fatigue. These strategies can ultimately result in the reduction in cardiovascular risk.

Physical Education on the Beach: An Alternative Way to Improve Primary School Children's Skill- and Health-Related Outcomes during the COVID-19 Pandemic

The COVID-19 restrictions could preclude children from participating in physical education (PE) interventions. This study aimed to evaluate the efficacy of a PE intervention conducted on the beach on children’s skill- and health-related outcomes, as a possible alternative PE intervention that could be also applied during the COVID-19 pandemic. The study involved 106 primary school children, randomly assigned to the traditional indoor (TI) intervention or to the experimental outdoor (EO) intervention. The intervention period lasted 4 months and consisted of two 1-h sessions per week. Intervention was conducted just before the beginning of the COVID-19 pandemic. Children’s anthropometric parameters (height, weight, BMI, body fat percentage, and abdominal circumference), fitness parameter (VO_2peak), health parameters (resting heart rate, and systolic and diastolic blood pressure), gross motor coordination, and physical activity level were assessed before and after intervention. Both groups significantly improved fitness and motor coordination but worsened some anthropometric parameters (weight, abdominal circumference) after the intervention period. The EO group showed a higher increase of gross motor coordination than the TI group. Results of this study demonstrated that children benefited from a well-structured PE intervention conducted in the natural environment of the beach improving physical fitness and gross motor coordination. Therefore, planning outdoor PE interventions could be an alternative and safe way to encourage and implement physical activity at school during the particular period of COVID-19 pandemic.

Recurrence quantification analysis of heart rate variability to detect both ventilatory thresholds

Aims of this study were: to verify if Recurrence Quantification Analysis (RQA) of Heart Rate Variability (HRV) time series could determine both ventilatory thresholds in individuals with different fitness levels, and to assess the validity of RQA method compared to gas-exchange method (GE). The two thresholds were estimated in thirty young individuals during incremental exercise on cycle-ergometer: Heart rate (HR), Oxygen consumption (VO₂) and Workload were measured by the two methods (RQA and GE). Repeated measures ANOVA was used to assess main effects of methods and methods-by-groups interaction effects for HR, VO₂ and Workload at aerobic (AerT) and anaerobic (AnT) thresholds. Validity of RQA at both thresholds was assessed for HR, VO₂ and Workload by Ordinary Least Products (OLP) regression, Typical Percentage Error (TE), Intraclass Correlation Coefficients (ICC) and the Bland Altman plots. No methods-by-groups interaction effects were detected for HR, VO₂ and Workload at AerT and AnT. The OLP analysis showed that at both thresholds RQA and GE methods had very strong correlations (r >0.8) in all variables (HR, VO₂ and Workload). Slope and intercept values always included the 1 and the 0, respectively. At AerT the TE ranged from 4.02% (5.48 bpm) to 10.47% (8.53 Watts) (HR and Workload, respectively) and in all variables ICC values were excellent (≥0.85). At AnT the TE ranged from 2.53% (3.98 bpm) to 6.64% (7.81 Watts) (HR and Workload, respectively) and in all variables ICC values were excellent (≥0.90). Therefore, RQA of HRV time series is a new valid approach to determine both ventilatory thresholds in individuals with different physical fitness levels, it can be used when gas analysis is not possible or not convenient.

Self-organization of whole-gene expression through coordinated chromatin structural transition

The human DNA molecule is a 2-m-long polymer collapsed into the micrometer space of the cell nucleus. This simple consideration rules out any "Maxwell demon"-like explanation of regulation in which a single regulatory molecule (e.g., a transcription factor) finds autonomously its way to the particular target gene whose expression must be repressed or enhanced. A gene-by-gene regulation is still more contrasting with the physical reality when in the presence of cell state transitions involving the contemporary expression change of thousands of genes. This state of affair asks for a statistical mechanics inspired approach where specificity arises from a selective unfolding of chromatin driving the rewiring of gene expression pattern. The arising of "expression waves" marking state transitions related to chromatin structural reorganization through self-organized critical control of whole-genome expression will be described in the present paper. We adopt as a model system the gene expression time course of a cancer cell (MCF-7) population exposed to an efficient stimulus causing a state transition in comparison with an ineffective stimulus. The obtained results will be put into the perspective of biological adaptive systems living on the edge of chaos.

Detection of Age-Related Hearing Losses (ARHL) via Transient-Evoked Otoacoustic Emissions

PURPOSE

The objective of the study was to identify subjects presenting hearing deficits, specifically age-related hearing losses (ARHL), via objective assessment methodologies.

MATERIALS AND METHODS

Initially, 259 subjects (165 men, 94 women) were enrolled in the study. After the application of inclusion criteria, the final number was reduced to 88 subjects (49.8 ± 19.1 ys) subdivided into 64 normal and 83 ARHL cases. The subjects were assessed with traditional audiometry tests and with transiently evoked otoacoustic emissions (TEOAEs). Since each ear has its own acoustic signature, the TEOAE analyses were conducted in terms of ears and not subjects. The TEOAE data were processed by traditional and recurrence quantification analyses, leading to the estimation of the WWR (whole waveform reproducibility) and the new RAD2D (2-dimensional radius) parameters. A plot of WWR vs RAD2D was used to optimize the classification of the cases presenting ARHL.

RESULTS

By using a WWR value of 70% as a classifier, the sensitivity of TEOAEs was estimated as 75.9% and the specificity as 89.1%. By using the RAD2D parameter (with a cut-off value of 1.78), a sensitivity value of 80.7% and a specificity value of 71.9% were obtained. When both parameters were used, a sensitivity value of 85.5% and a specificity value of 92.2% were estimated. In the latter classification paradigm, the number of false negatives decreased from 20 to 12 out of 83 ears (14%).

CONCLUSION

In adult hearing screening assessments, the proposed method optimizes the identification of subjects with a hearing impairment correlated to the presence of age-related hearing loss.