Hostname: page-component-7b9c58cd5d-wdhn8 Total loading time: 0 Render date: 2025-03-15T12:29:35.610Z Has data issue: false hasContentIssue false

Non-linear analysis of heart rate variability for evaluating the acute effects of caffeinated beverages in young adults

Published online by Cambridge University Press:  09 June 2020

Serife G. Caliskan*
Affiliation:
Department of Physics, Science and Art Faculty, Aydın Adnan Menderes University, Aydın, Turkey
Mehmet D. Bilgin
Affiliation:
Department of Biophysics, School of Medicine, Aydın Adnan Menderes University, Aydın, Turkey
*
Author for correspondence: Serife G. Caliskan, PhD, Department of Physics, Aydın Adnan Menderes University Science and Art Faculty, Merkez Kampus, Aytepe, Efeler, Aydın09010, Turkey. Tel: +90 256 212 84 98, Fax: +90 256 213 53 79. E-mail: sgokcecaliskan@gmail.com
Rights & Permissions [Opens in a new window]

Abstract

Caffeinated beverages are the most consumed substances in the world. High rate of uptake of these beverages leads to various cardiovascular disorders ranging from palpitations to coronary failure. The objective of the study is to ascertain how the complexity parameters of heart rate variability are affected by acute consumption of caffeinated beverages in young adults.

Electrocardiogram measurements were performed before consuming drinks. After consuming the drinks, measurements were done again at 30 minutes and 60 minutes. Heart rate variability signals were acquired from electrocardiogram signals. Also, the signals were reconstructed in the phase space and largest Lyapunov exponent, correlation dimension, approximate entropy, and detrended fluctuation analysis values were calculated.

Heart rate increased for energy drink and cola groups but not in coffee group. Non-linear parameter values of energy drink, coffee, and cola group are increased within 60 minutes after drink consumption. This change is statistically significant just for energy drink group.

Energy drink consumption increases the complexity of the cardiovascular system in young adults significantly. Coffee and cola consumption have no significant effect on the non-linear parameters of heart rate variability.

Type
Original Article
Copyright
© The Author(s), 2020. Published by Cambridge University Press

Introduction

Caffeine which is a methylxanthine compound is probably the most consumed stimulant in the world. It increases the autonomic nervous system activity and blocks adenosine receptors (A1 and A2) by releasing the plasma catecholamines.Reference Gonzaga, Vanderlei, Gomes and Valenti1 This effect induces tachycardia and hypertension and also increases respiratory rate, bronchodilatation, lipolysis, and diuresis.Reference Monda, Viggiano and Vicidomini2 It is a dominant ingredient in many beverages, such as coffee, tea, energy drinks, and cola-type drinks. These are the most popular beverages that are consumed uncontrollably by young people for the increasing effect upon attention, arousal, and physical performance. In the last decades, there are increasing reports of especially energy drinks related with cardiovascular events, seizures, and death.Reference Wolk, Ganetsky and Babu3 The effects of these beverages must be investigated in detail because of the high-rate consumption and occasional acute adverse health effects. There are numerous studies which investigated the acute effects of caffeinated beverages on haemodynamic parameters, such as blood pressure and heart rate variability.Reference Gonzaga, Vanderlei, Gomes and Valenti1,Reference Steinke, Lanfear, Dhanapal and Kalus4Reference Ammar, Song, Kluger and White6 For instance, Hajsadeghi et alReference Hajsadeghi, Mohammadpour and Manteghi7 reported a significant heart rate decline and Peveler et alReference Peveler, Sanders, Marczinski and Holmer8 showed a significant increase in systolic blood pressure after the consumption of energy drink. Consequently, a number of studies have shown the adverse effects of caffeinated beverages on cardiovascular system, but the dynamics and the mechanism are still ambiguous.

Linear and non-linear methods are being used for a better understanding of physiological signals and systems. Time-domain and frequency-domain analysis which are linear methods can only present first- and second-order statistics of the signals so they can give limited information about the system. However, non-linear methods reveal and characterise the dynamics and the complexity of the systems in detail. For instance, the physiological interaction between the heart period and autonomic nervous system is being determined by the analysis of heart rate variability using linear and non-linear methods.Reference Voss, Schulz, Schroeder, Baumert and Caminal9 However, analysing heart rate variability through time and frequency domains focus only on sympathovagal balance assessment which is inadequate to describe the complex time behaviour of heart rhythm.Reference Delliaux, Delaforge, Deharo and Chaumet10 Since the cardiovascular system is complex and do not show regular periodicity, it cannot be completely explained by linear methods. To overcome this adequacy, non-linear methods are being used extensively. For instance, Koichubekov et al,Reference Koichubekov, Korshukov, Omarbekova, Riklefs, Sorokina and Mkhitaryan11 Pivatelli et al,Reference Pivatelli, dos Santos and Fernandes12 and Neves et alReference Neves, Takahashi and do Santos-Hiss13 analysed heart rate variability by evaluating the non-linear indices, such as largest Lyapunov exponent, approximate entropy and correlation dimension in patients with hypertension, stable angina, and coronary disease, respectively. Also, Tsai et al reported that detrended fluctuation analysis shown better predictive ability for clinical outcomes when compared to traditional heart rate variability analysis methods.Reference Tsai, Ma and Lin14 In the literature, there are some studies which investigated heart rate variability time series with detrended fluctuation analysis and reported the accuracy of the method.Reference Papaioannou, Vlachopoulos, Ioakeimidis, Alexopoulos and Stefanadis15,Reference Moraes, Silva and Massetti16 To our knowledge, there are no reports in the literature on non-linear analysis of heart rate variability used to investigate the acute effects of caffeinated beverages on cardiovascular system. For this reason, our study will be the first for this subject in the literature.

However, due to the fact that the effects on the human body are not fully known and the interest in these beverages has increased, the need for further study in this issue has arisen. In this study, it is aimed to ascertain how the complexity parameters of heart rate variability are affected by acute consumption of caffeinated beverages. In this respect, expressive information can be obtained about the system’s dynamic structure in more detail.

Materials and methods

Participants

Fourty-eight participants (20 female, 28 male) aged 18–24 years (mean 19.20 ± 0.48) were selected from the students of Aydın Adnan Menderes University. They were identified by applicant acceptance questionnaire which includes information about age, weight, height, body mass index, and caffeine consumption habit. The participants with a body mass index greater than 30 kg/m2, a systemic disease, current alcoholism, smokers, and who has previous adverse reactions to caffeinated beverages were excluded from the study. All subjects were caffeine-naïve and grouped as energy drink (n = 12), coffee (n = 12), cola (n = 12), and control (n = 12). This study conforms to the Declaration of Helsinki and was approved by the Institutional Medical Ethics Committee of Adnan Menderes University (Protocol number: 2017/1211). Informed consent was obtained from all individual participants included in the study.

Experimental protocol

All experiments were performed at the same time of the day for avoiding the different responses due to circadian rhythm. The room was isolated from external factors with the room temperature of about 20°C and relative humidity between 40 and 60%. All participants were instructed to abstain from alcohol and caffeine-containing foods for at least 72 hours prior to the test and to come with an overnight (12 hours) fasting.

Electrocardiogram signals were recorded before and after the consumption of the drinks in 30 minutes and 60 minutes. Coffee group consumed 473 ml of Nescafe classic, energy drink group consumed 473 ml of Red Bull energy drink, cola group consumed 473 ml of Coca Cola, and control group consumed the same amount of water within 5 minutes. These drinks contain 3784 mg, 1514 mg, and 404 mg caffeine, respectively.Reference Papakonstantinou, Kechribari and Sotirakoglou17,Reference Voskoboinik, Koh and Kistler18

Data acquisition

Electrocardiogram recordings were obtained from the participants for 5 minutes in rest and sitting position. Disposable Ag–AgCl electrodes were placed following the Einthoven triangle configuration. Electrocardiogram measurement was performed with sampling frequency of 200 Hz using ECG100C unit and BIOPAC Acqknowledge acquisition software (Biopac System Inc., Santa Barbara, CA) connected to a personal computer. We analysed the 5-minute electrocardiogram signal to extract the heart rate variability signal. All analyses were performed according to the standards set by Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology.19

For acquiring the heart rate variability time series, Pan and Thompkins algorithm was used.Reference Pan and Thompkins20 R peaks were detected automatically and R–R intervals are obtained. The time duration between two consecutive R peaks is termed as RR interval (tRR).Reference Acharya, Faust and Sree21 The sample size of the heart rate variability time series is 1000. The heart rate is determined as:

$Heart\,Rate = {{60} \over {{t_{RR}}}}.$

Non-linear analysis

Phase space reconstruction

For non-linear analysis, heart rate variability time series was reconstructed in the phase space using data’s time-delayed copies. Embedding dimension and time-delay values were calculated using Cao’s methodReference Cao22 and mutual information methodReference Fraser and Swinney23 for reconstruction process. Largest Lyapunov exponent, correlation dimension, and approximate entropy parameters were calculated from the attractor which was reconstructed in the phase space. For this purpose, we used MATLAB (The Mathworks Inc., USA). They represent the complexity and irregularity of the system which was investigated. These were explained in the study of Caliskan et alReference Caliskan, Polatli and Bilgin24 in detail.

Detrended fluctuation analysis

Detrended fluctuation analysis is used to evaluate the fractal behaviour of non-stationary RR dynamics by eliminating extrinsic trends to remove spurious long-term correlations.Reference Tsai, Lin and Ho25 It allows us to have a differentiated view on the correlative structure of variability caused by physiological processes in time series.Reference Gronwald and Hoos26 It is a suitable method for the evaluation of heart rate variability which is a short and non-stationary time series data since it is a modification of the root mean square analysis.Reference Peng, Havlin, Stanley and Goldberger27 In this analysis, first of all linear and polynomial fitted trends which are called the external trends are removed from the integrated time series. Detrended integrated time series was added up in individual scales for calculating the detrended fluctuations. Then, fluctuations are plotted logarithmically against timescales. This graph is consisted of two distinct curves separated by one point. The fractal correlation characters of the time series are obtained from the slopes (α1 and α2 exponent) of the log–log graph. If α is between 0 and 0.5, the large and small values of the time series are more likely to alternate. While α = 0.5 corresponds to white noise, α = 1.5 corresponds to Brownian noise, and if α is between 0.5 and 1.0 then there are positive correlations present in the time series.Reference Hardstone, Poil and Schiavone28 For healthy subjects, α1 tends to 1 (unity) and is bigger than α2 value.Reference Moraes, Silva and Massetti16

In this study, heart rate variability time series was analysed by detrended fluctuation analysis. Both short-term (α1) and long-term (α2) exponents were calculated in order to have detailed information about the fractal characters of the system. Detrended fluctuation analysis parameters were calculated using Kubios HRV Standard (Kubios Heart Rate Variability analysis software version 3.1.0; Kupio, Finland).

Statistical analysis

Statistical analysis of the data was performed using SPSS Statistics for Windows, Version 22.0. (SPSS Inc., Chicago, IL, USA). For comparing the groups, first of all Shapiro–Wilk normality test was applied to determine the normality of the distribution. If the result of Shapiro–Wilk normality test was bigger than 0.05, the data were normally distributed and were analysed using one-way analysis of variance with Bonferroni post-test, but if the result was less than 0.05 the data were not normally distributed and were analysed using Kruskal–Wallis non-parametric test. Repeated one-way analysis of variance test was used for the comparison of different measurement time data of the same group. Descriptive statistics were presented as mean ± standard error of the mean. A p-value of less than 0.05 was considered to indicate statistical significance.

Results

Table 1 displays the characteristics of the groups. There was no significant difference in age, height, weight, and body mass index values between the groups. Furthermore, 60-minute heart rate values of the energy drink and cola groups were increased with regard to pre-drink heart rate values. Only, the increase in energy drink group was statistically significant (p = 0.045). There was no change for coffee group. Heart rate values of control group did not change significantly across time.

Table 1. The characteristics of the participants

Non-linear analysis results of heart rate variability for experimental groups are shown in Table 2. Heart rate variability time series was reconstructed in the phase space in order to investigate the dynamics of the system in detail. Embedding dimension and time-delay values were determined as 3 and 2, respectively. Using the embedding dimension and time-delay values, heart rate variability time series was reconstructed in the phase space. Largest Lyapunov exponent and correlation dimension values of energy drink group at 60 minutes after consumptions were significantly increased when compared to pre-drink largest Lyapunov exponent values (Figs 1 and 2) (p = 0.016, p = 0.003, respectively). These non-linear parameters were increased slightly for coffee and cola groups but these changes were not statistically significant. Also, they did not change for control group across time. Similarly, approximate entropy, detrended fluctuation analysis −α1, and detrended fluctuation analysis −α2 values were increased within 60 minutes of consumption when compared to pre-drink values (Table 2). However, the differences were not statistically significant. There is no statistically significant alteration for the non-linear parameters at 30 minutes of consumption for all experimental groups. Furthermore, both for largest Lyapunov exponent and correlation dimension values, there is a statistically significant difference between energy drink and control group at 60 minutes of consumption (p = 0.014, p = 0.011). This is as the same for coffee group in correlation dimension parameter also (p = 0.023).

Table 2. Non-linear parameters of the groups before, after 30 minutes and 60 minutes of drink consumptions. P refers to the differences between the pre and 60th minute values of groups

Figure 1. The variation of largest Lyapunov exponent values of the groups during the experiment. Δ statistically significant versus pre-drink value of energy drink group (p = 0.016) and * statistically significant versus 60-minute value of control (p = 0.014).

Figure 2. The variation of correlation dimension values of the groups during the experiment. Δ statistically significant versus pre-drink value of energy drink group (p = 0.003) and * statistically significant versus 60-minute value of control (p = 0.011, p = 0.023, respectively).

Discussion

The acute effects of caffeinated beverages on complexity parameters of heart rate variability were investigated in this study. The principle finding is that energy drink consumption increases the complexity of the cardiovascular system in young adults significantly. Coffee and cola consumption have no significant effects on the non-linear parameters of heart rate variability.

In the last decades, especially young people consume caffeinated beverages more frequently for efficient alertness, physical strength, and mood enhancement. Insufficient information about the adverse effects of these beverages makes its consumption uncontrolled among young people. There are many case reports about their life-threatening adverse effects such as cardiac arrest, myocardial infarction, and atrial fibrillation.Reference Berger and Alford29Reference Ullah, Lakhani, Siddiq, Handa, Kahlon and Siddiqui31 There exist great health concerns worldwide because of the widespread use of these beverages. For this reason, its effects on young human must be investigated in more detail.

Heart rate variability analysis is a mostly preferred method since it divulges information about the influence of the autonomic nervous system activity on the heart and it is a reliable, informative, and applicable tool in daily practice.Reference Moraes, Silva and Massetti16 The most important step is the extraction of the features containing all the relevant information about the system dynamics.Reference Nayak, Bit, Dey, Mohapatra and Pal32 Time-domain, frequency-domain, and non-linear methods are generally used for this procedure. While time-domain and frequency-domain methods enable the quantification of heart rate variability on different timescales, non-linear methods provide additional information regarding the dynamics and structure of beat-to-beat time series.Reference Voss, Schulz, Schroeder, Baumert and Caminal9 Linear analysis methods have some limitations since it has sensitiveness about the band definition (low frequency, from 0.04 to 0.15 Hz and high frequency, from 0.15 to 0.40 Hz).Reference Neves, Takahashi and do Santos-Hiss13 Also, the physiological signals are not constant on time and do not show regular periodicity. Therefore, they cannot be explained accurately by linear analysis methods.Reference Lombardi33 Non-linear analysis outperforms time- and frequency-domain methods by providing a continuous analysis of the signal and detailed information about the system dynamics. These methods were applied in various studies in which the effects of different diseases on heart rate variability were investigated.Reference Pivatelli, dos Santos and Fernandes12,Reference Acharya, Faust and Sree21,Reference Rossi, Vanderlei and Bernardo34

In our study, we have shown that energy drink and cola consumption have risen the heart rate values when compared to pre-drink values. The increment for energy drink group was consistent with the studies in which reported a significant increase in heart rate after energy drink consumption.Reference Elitok, Öz and Panc5,Reference Grasser, Yepuri, Dulloo and Montani35 In addition, there are no studies that evaluate the acute effects of cola consumption on haemodynamic parameters in the literature. Therefore, there is no opportunity to compare the results we obtain. Furthermore, there was an interesting finding that in contradiction to other groups, there was no alteration in heart rate of coffee group. This is in line with the study of Hara et alReference Hara, Ohide and Miyagawa36 which reported no effect of coffee consumption in non-habitual coffee consumers. Also, Ammar et alReference Ammar, Song, Kluger and White6 reported that there was no alteration in heart rate with acute dosing of caffeine in healthy volunteers. The difference between energy drink–cola groups and coffee group might be due to the other ingredients of these beverages. One of the major contents of energy drink is taurine which is well known to rise heart rate.Reference Schaffer, Jong, Ramila and Azuma37 Also, glucose is one of the ingredients of both energy drink and cola and it has been shown that glucose ingestion or administration affects cardiac function and increases cardiac output.Reference Charrière, Loonam, Montani, Dulloo and Grasser38

On the other hand, we calculated largest Lyapunov exponent and correlation dimension parameters and found a statistically significant increase between the pre-drink and post-drink values of energy drink group. Also, these parameters increased for coffee and cola groups but these increments were trivial. These parameters give information about the system’s chaoticity and according to our results it can be deduced that caffeinated beverage consumption has a significant effect on the complexity of heart rate variability. In addition, there are no studies that evaluate the acute effects of caffeinated beverage consumption on heart rate variability by analysing these non-linear parameters in the literature. Therefore, there is no opportunity to compare the results we obtain.

Entropy is the measure of the degree to which the occurrence of a value depends on its predecessors in the input.Reference Koichubekov, Korshukov, Omarbekova, Riklefs, Sorokina and Mkhitaryan11 High values of the entropy indicate high irregularity and complexity in time series data.Reference Voss, Schulz, Schroeder, Baumert and Caminal9 In our study, we found an increment in approximate entropy values of the experimental groups when compared with pre-drink values for 30 and 60 minutes of caffeinated beverage consumption. This result is in agreement with Yeragani et al’s studyReference Yaregani, Krishnan, Engels and Gretebeck39 in which approximate entropy values are increased when compared to control group and concluded that caffeine consumption increased heart rate variability at rest. Furthermore, we investigated the detrended fluctuation analysis parameters which quantify the absence or presence of fractal correlation properties in RR intervals.Reference Moraes, Silva and Massetti16 In many studies, detrended fluctuation analysis parameters were investigated and the researchers suggested this analysis as an ideal way to have more accurate information about the complexity of heart rate variability.Reference Moraes, Silva and Massetti16,Reference Tsai, Lin and Ho25 However, there are very few studies related to the non-linear analysis of heart rate variability signals after caffeine consumption. Papaioannou et alReference Papaioannou, Vlachopoulos, Ioakeimidis, Alexopoulos and Stefanadis15 reported a significant increase in long-term scaling exponent α after caffeine ingestion, while the placebo induced no significant change. Our result is in agreement with this study. In the literature, it was reported that increased complexity and chaos were related to pathological phenomena.Reference Garfinkel, Spano, Ditto and Weiss40 Also, there exist various studies which observe a significant decrease in detrended fluctuation analysis of -α values of heart rate variability. It has shown that the complexity decreased in the presence of a systemic disease.Reference Pivatelli, dos Santos and Fernandes12,Reference Tsai, Ma and Lin14 The main reason of the difference may be the existence of a disease in the participants and the breakdown of the balance in the systems.

We acknowledge several limitations in this study. Relatively small number of the participants is one of the most important limitations of the present study. Also, we have evaluated only one dose of caffeinated beverage consumption which can be accepted as a low dose. Higher doses of these beverages may affect the results. We recorded electrocardiogram signals for 5 minutes which can be accepted as a short recording duration. A 24-hour recording can be more effective. The exact mechanism of heart rate variability complexity is remained unknown and further studies must be planned.

Conclusion

In conclusion, energy drink consumption increases the complexity of the cardiovascular system in young adults. Coffee and cola consumption have no significant effect on the non-linear parameters of heart rate variability. For investigating the effects of caffeinated beverage consumption more accurately, further studies which investigate the higher doses of caffeinated beverages are being planned.

Acknowledgements

We thank Research Assistant M. Alp Kılıç and the participants for their contribution to this paper.

Financial Support

This research received no specific grant from any funding agency, commercial or not-for-profit sectors.

Conflict of Interest

None.

Ethical Standards

The authors assert that all procedures contributing to this work comply with the Helsinki Declaration of 1975, as revised in 2008, and have been approved by the institutional committees (Institutional Medical Ethics Committee of Adnan Menderes University (Protocol number: 2017/1211)).

References

Gonzaga, LA, Vanderlei, LCM, Gomes, RL, Valenti, VE.Caffeine affects autonomic control of heart rate and blood pressure recovery after aerobic exercise in young adults: a crossover study. Sci Rep 2017; 7: 14091.10.1038/s41598-017-14540-4CrossRefGoogle ScholarPubMed
Monda, M, Viggiano, A, Vicidomini, , C, et al.Expresso coffee increases parasympathetic activity in young, healthy people. Nutr Neurosci 2009; 12: 4348.10.1179/147683009X388841CrossRefGoogle Scholar
Wolk, BJ, Ganetsky, M, Babu, KM.Toxicity of energy drinks. Curr Opin Pediatr 2012; 24: 243251.10.1097/MOP.0b013e3283506827CrossRefGoogle ScholarPubMed
Steinke, L, Lanfear, DE, Dhanapal, V, Kalus, JS.Effect of “energy drink” consumption on hemodynamic and electrocardiographic parameters in healthy young adults. Ann Pharmacother 2009; 43: 596602.10.1345/aph.1L614CrossRefGoogle ScholarPubMed
Elitok, A, Öz, F, Panc, C, et al.Acute effects of red bull energy drink on ventricular repolarization in healthy young volunteers: a prospective study. Anatol J Cardiol 2015; 15: 919922.10.5152/akd.2015.5791CrossRefGoogle ScholarPubMed
Ammar, R, Song, JC, Kluger, J, White, M.Evaluation of electrocardiographic and hemodynamic effects of caffeine with acute dosing in healthy volunteers. Pharmacotherapy 2001; 21: 437442.CrossRefGoogle ScholarPubMed
Hajsadeghi, S, Mohammadpour, F, Manteghi, MJ, et al.Effects of energy drinks on blood pressure, heart rate, and electrocardiographic parameters: an experimental study on healthy young adults. Anatol J Cardiol 2016; 16: 9499.Google Scholar
Peveler, WW, Sanders, GJ, Marczinski, CA, Holmer, B.Effects of energy drinks on economy and cardiovascular measures. J Strength Cond Res 2017; 31: 882887.10.1519/JSC.0000000000001553CrossRefGoogle ScholarPubMed
Voss, A, Schulz, S, Schroeder, R, Baumert, M, Caminal, P.Methods derived from nonlinear dynamics for analysing heart rate variability. Philos Trans Sries A Math Phys Eng Sci 2009; 367: 12231238.Google ScholarPubMed
Delliaux, S, Delaforge, A, Deharo, J-C, Chaumet, G.Mental workload alters heart rate variability, lowering non-linear dynamics. Front Physiol 2019; 10: 565. doi: 10.3389/fphys.2019.00565.CrossRefGoogle ScholarPubMed
Koichubekov, B, Korshukov, I, Omarbekova, N, Riklefs, V, Sorokina, M, Mkhitaryan, X.Computation of nonlinear parameters of heart rhythm using short time ECG segments. Comput Math Methods Med 2015; 2015: 983479.10.1155/2015/983479CrossRefGoogle ScholarPubMed
Pivatelli, FC, dos Santos, MA, Fernandes, GB, et al.Sensitivity, specificity and predictive values of linear and nonlinear indices of heart rate variability in stable angina patients. Int Arch Transl Med 2012; 5: 31.10.1186/1755-7682-5-31CrossRefGoogle ScholarPubMed
Neves, VR, Takahashi, ACM, do Santos-Hiss, MDB, et al.Linear and nonlinear analysis of heart rate variability in coronary disease. Clin Auton Res 2012; 22: 175183.10.1007/s10286-012-0160-zCrossRefGoogle ScholarPubMed
Tsai, CH, Ma, HP, Lin, YT, et al.Heart rhythm complexity impairment in patients with pulmonary hypertension. Sci Rep 2019; 9: 10710. doi: 10.1038/s41598-019-47144-1.CrossRefGoogle ScholarPubMed
Papaioannou, TG, Vlachopoulos, C, Ioakeimidis, N, Alexopoulos, N, Stefanadis, C.Nonlinear dynamics of blood pressure variability after caffeine consumption. Clin Med Res 2006; 4: 114118.10.3121/cmr.4.2.114CrossRefGoogle ScholarPubMed
Moraes, IAP, Silva, TD, Massetti, T, et al.Fractal correlations and linear analysis of heart rate variability in healthy young people with different levels of physical activity. Cardiol Young 2019; 29: 12361242. doi: 10.1017/S1047951119001793.CrossRefGoogle ScholarPubMed
Papakonstantinou, E, Kechribari, I, Sotirakoglou, K, et al.Acute effects of coffee consumption on self-reported gastrointestinal symptoms, blood pressure and stress indices in healthy individuals. Nutr J 2016; 15: 26.10.1186/s12937-016-0146-0CrossRefGoogle ScholarPubMed
Voskoboinik, A, Koh, Y, Kistler, PM.Cardiovascular effects of caffeinated beverages. Trends Cardiovas Med 2019; 29: 345350.10.1016/j.tcm.2018.09.019CrossRefGoogle ScholarPubMed
Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Eur Heart J 1996; 17: 354381.Google Scholar
Pan, J, Thompkins, WJ.A real-time QRS detection algorithm. IEEE Trans Biomed Eng 1985; 32: 230236.10.1109/TBME.1985.325532CrossRefGoogle ScholarPubMed
Acharya, UR, Faust, O, Sree, V, et al.Linear and nonlinear analysis of normal and CAD-affected heart rate signals. Comput Methods Programs Biomed 2014; 113: 5568.10.1016/j.cmpb.2013.08.017CrossRefGoogle ScholarPubMed
Cao, L.Practical method for determining the minimum embedding dimension of a scalar time series. Physica D 1997; 110: 4350.10.1016/S0167-2789(97)00118-8CrossRefGoogle Scholar
Fraser, AM, Swinney, HL.Independent coordinates for strange attractors from mutual information. Phys Rev A 1986; 33: 11341140.10.1103/PhysRevA.33.1134CrossRefGoogle ScholarPubMed
Caliskan, SG, Polatli, M, Bilgin, MD.Nonlinear analysis of heart rate variability of healthy subjects and patients with chronic obstructive pulmonary disease. J Med Eng Technol 2018; 42: 298305.10.1080/03091902.2018.1491650CrossRefGoogle ScholarPubMed
Tsai, CH, Lin, C, Ho, YH, et al.The association between heart rhythm complexity and the severity of abdominal aorta calcification in peritoneal dialysis patients. Sci Rep 2018; 8: 15627. doi: 10.1038/s41598-018-33789-x.CrossRefGoogle ScholarPubMed
Gronwald, T, Hoos, O.Correlation properties of heart rate variability during endurance exercise: a systematic review. Ann Noninvasive Electrocardiol 2019; 25: e12697. doi: 10.1111/anec.12697.Google ScholarPubMed
Peng, CK, Havlin, S, Stanley, HE, Goldberger, AL.Quantification of scaling exponents and crossover phenomena in nonstationary heart beat time series. Chaos 1995; 5: 8287.10.1063/1.166141CrossRefGoogle ScholarPubMed
Hardstone, R, Poil, SS, Schiavone, G, et al.Detrended fluctuation analysis: a scale-free view on neuronal oscillations. Front Physiol 2012; 30: 450. doi: 10.3389/fphys.2012.00450.Google Scholar
Berger, AJ, Alford, K.Cardiac arrest in a young man following excess consumption of caffeinated “energy drinks”. Med J Aust 2009; 190: 4143.10.5694/j.1326-5377.2009.tb02263.xCrossRefGoogle Scholar
Mattioli, AN, Pennella, S, Farinetti, A, Manenti, A.Energy drinks and atrial fibrillation in young adults. Clin Nutr 2018; 37: 10731074.10.1016/j.clnu.2017.05.002CrossRefGoogle ScholarPubMed
Ullah, MW, Lakhani, S, Siddiq, W, Handa, A, Kahlon, Y, Siddiqui, T.Energy drinks and myocardial infarction. Cureus 2018; 10: e2658.Google Scholar
Nayak, SK, Bit, A, Dey, A, Mohapatra, B, Pal, K.A review on the nonlinear dynamical system analysis of electrocardiogram signal. J Healthc Eng 2018; 2: 119.CrossRefGoogle Scholar
Lombardi, F.Chaos theory, heart rate variability, and arrhythmic mortality. Circulation 2000; 101: 810.10.1161/01.CIR.101.1.8CrossRefGoogle ScholarPubMed
Rossi, RC, Vanderlei, FM, Bernardo, AF, et al.Effect of pursed-lip breathing in patients with COPD: linear and nonlinear analysis of cardiac autonomic modulation. COPD 2016; 11: 3945.10.3109/15412555.2013.825593CrossRefGoogle Scholar
Grasser, EK, Yepuri, G, Dulloo, AG, Montani, JP.Cardio- and cerebrovascular responses to the energy drink Red Bull in young adults: a randomized cross-over study. Eur J Nutr 2014; 53: 15611571.10.1007/s00394-014-0661-8CrossRefGoogle ScholarPubMed
Hara, A, Ohide, H, Miyagawa, K, et al.Acute effects of caffeine on blood pressure and heart rate in habitual and non-habitual coffee consumers: a randomized, double blind, placebo-controlled study. Jpn J Pharm Health Care Sci 2014; 40: 383388.CrossRefGoogle Scholar
Schaffer, SW, Jong, CJ, Ramila, KC, Azuma, J.Physiological roles of taurine in heart and muscle. L Biomed Sci 2010; 14: S2.10.1186/1423-0127-17-S1-S2CrossRefGoogle Scholar
Charrière, N, Loonam, C, Montani, JP, Dulloo, AG, Grasser, EK.Cardiovascular responses to sugary drinks in humans: galactose present milder cardiac effects than glucose or fructose. Eur J Nutr 2017; 56: 21052113.10.1007/s00394-016-1250-9CrossRefGoogle ScholarPubMed
Yaregani, VK, Krishnan, S, Engels, HJ, Gretebeck, R.Effects of caffeine on linear and nonlinear measures of heart rate variability before and after exercise. Depress Anxiety 2005; 21: 130134.10.1002/da.20061CrossRefGoogle Scholar
Garfinkel, A, Spano, ML, Ditto, WL, Weiss, JN.Controlling cardiac chaos. Science 1992; 257: 12301235.10.1126/science.1519060CrossRefGoogle ScholarPubMed
Figure 0

Table 1. The characteristics of the participants

Figure 1

Table 2. Non-linear parameters of the groups before, after 30 minutes and 60 minutes of drink consumptions. P refers to the differences between the pre and 60th minute values of groups

Figure 2

Figure 1. The variation of largest Lyapunov exponent values of the groups during the experiment. Δ statistically significant versus pre-drink value of energy drink group (p = 0.016) and * statistically significant versus 60-minute value of control (p = 0.014).

Figure 3

Figure 2. The variation of correlation dimension values of the groups during the experiment. Δ statistically significant versus pre-drink value of energy drink group (p = 0.003) and * statistically significant versus 60-minute value of control (p = 0.011, p = 0.023, respectively).