ORIGINAL RESEARCH article
Sec. Autonomic Neuroscience
Variations of Time Irreversibility of Heart Rate Variability Under Normobaric Hypoxic Exposure
- 1School of Instrument Science and Engineering, Southeast University, Nanjing, China
- 2Air Force Medical Center, Beijing, China
In the field of biomedicine, time irreversibility is used to describe how imbalanced and asymmetric biological signals are. As an important feature of signals, the direction of time is always ignored. To find out the variation regularity of time irreversibility of heart rate variability (HRV) in the initial stage of hypoxic exposure, the present study implemented 2 h acute normobaric hypoxic exposure on six young subjects who have no plateau or hypoxia experiences; oxygen concentration was set as 12.9%. Electrocardiogram (ECG) signals were recorded in the whole process and RR interval sequences were extracted. Mathematical operations were executed to transform the difference of adjacent RR intervals into proportion and distance with delay time to conduct time irreversibility analysis of HRV. The same calculating method was implemented on six items randomly picked out from the MIT-BIH normal sinus rhythm database as a control group. Results show that variation of time irreversibility of HRV in a hypoxic environment is different from that in a normoxic environment, time irreversibility indices of a hypoxic group decreases continually at a delay time of 1 and 2, and indices curves of time irreversibility gradually tend to be steady and gather with each other at a delay time of 3 or 4. The control group shows no consistent tendency no matter what the delay time is in the range of 1–4. Our study indicates that in short-time hypoxic exposure, as hypoxic time goes by, regulation of the cardiovascular autonomic nervous system weakens; regulation times and intensity of sympathetic and parasympathetic nerves tend to be equal.
Time irreversibility is one of the basic features of nonlinear systems. It describes the process of how the physical state and statistical characteristics of a system depend on the direction of time (Yao et al., 2019). If a sequence has the same joint probability distribution with its reversal sequence, it is said to be time invertible. Invertible processes include Gaussian linear processes and isentropic processes. On the contrary, time irreversibility of series indicates the existence of nonlinear processes in the underlying dynamics, including non-Gaussian stochastic processes and dissipative chaos (Lacasa et al., 2012). The heartbeat of normal individuals is jointly controlled by multiple factors of sympathetic and parasympathetic nerves, and thus the time series of cardiac intervals shows characteristics of chaos, nonlinearity, and complexity (Rajbhandari Panday and Panday, 2018). Costa et al. (2005) pointed out that time irreversibility of cardiac intervals is influenced by age and pathological factors; healthy young people have significantly higher time irreversibility of cardiac intervals than that of healthy elderly people, and healthy aged people have distinctly higher time irreversibility of cardiac intervals than patients with congestive heart failure (CHF) or atrial fibrillation. This conclusion was also supported by Hou et al. (2011) with multi-scale time irreversibility method, and she also indicated that there are significant differences of time irreversibility in healthy young people between day and night (Hou et al., 2012).
Many methods have been proposed for the measurement of time irreversibility of cardiac interval sequences. Differentiate operations between adjacent RR intervals is used by Porta et al. (2006) who took the occurrence frequency of adjacent intervals’ difference greater than or less than zero as the time irreversibility index. And Costa et al. (2008) introduced the coarse-grained multi-scale method on Porta’s work. Guzik et al. (2006) performed scaling operations on the square of RR intervals. Lacasa et al. (2012) traced the amplitude of each point in the sequence with a graph, calculated the number of points that each point could reach without crossing with other points, and constructed a new sequence as the calculating basis of time irreversibility. Ordinal mode has been applied to calculate the emergence frequency of various patterns (Graff et al., 2013). Hou et al. (2011) defined high-dimensional irreversibility calculation index on the work of Porta (Graff et al., 2013) and Guzik (Guzik et al., 2006) by introducing different delay times for irreversibility analysis, and the mean value of all calculation results is the value of high-dimensional irreversibility. It reflects the comprehensive influence of different levels’ time delay of nerves and fluid in the human body on time irreversibility of heart rate variability (HRV). And research shows that both age and exercise have a great effect on high-dimensional irreversibility (Hou et al., 2011, 2013).
Hypoxia is one of the stimulating factors that make HRV change. When entering a hypoxic environment, blood oxygen saturation (SpO2) of the human body decreases, and regulation of the autonomic nervous system changes, which is usually manifested as an accelerated heart rate and shortened breath (Xue et al., 2010; Liu et al., 2019). Lower oxygen concentration may even result in altitude sickness, with symptoms such as headache and nausea (Zhang et al., 2014). Different methods have been applied to assess how HRV changes when people are exposed to a hypoxic environment. Most of the research results found that the standard deviation (SDNN) and root mean square (RMSSD) between adjacent heartbeat intervals and low frequency (LF) and high frequency (HF) in the spectrum decrease while the ratio of LF to HF increases (Krejčí et al., 2018) and sample entropy declines (Zhang et al., 2014) at the same time. However, other researchers gave different results (Iwasaki et al., 2006; Vigo et al., 2010). Normobaric hypoxic chambers are often used in laboratories to simulate the plateau hypoxic environment. It creates hypoxic effects by filling the chamber with excess nitrogen to dilute the air, which is safe, economic, and easy to use.
Existing studies mostly compare the physiological state under normoxia and hypoxia and find out how physiological indicators change after people are exposed to a hypoxic environment. For example, variation of electrocardiogram (ECG) in hypoxic exposure was studied in Coustet et al. (2015); variation tendencies of HRV and SPO2 within 10 min of acute hypoxic exposure were revealed in Krejčí et al. (2018). Some other studies explored people’s hypoxic acclimatization speed or physical training methods through experiments (Beidleman et al., 2008; Zeng et al., 2013, 2019; Sinex and Chapman, 2015). As one of the methods to evaluate the health of the physiological system, time irreversibility analysis method has been widely used in medical research. It can reflect the pathological information of the cardiovascular system and aging status, however, this method has never been used in the study of human physiological variation in a hypoxic environment. Some symptoms that occur in hypoxia are similar to other diseases and there must be connections between time irreversibility and those symptoms. We need to conduct research over a continuous period of time at the beginning of hypoxia to study the subtle variation trends of physiological indicators. This research is of great importance for the prevention of acute altitude sickness, the physical training of pilots and astronauts in response to accidents of loss of pressure, and the revelation of dynamics and autonomic rhythm changes of the human cardiovascular nervous system in the early stage of hypoxic exposure. In the process of cardiac feedback control, delay happens at all levels of the physiological system due to the reaction speed of chemical transmitters (Casali et al., 2008; Alvarez-Ramirez et al., 2009), and some useful information that the data can express may be lost if time irreversibility is calculated only based on the difference between adjacent RR intervals. In order to explore variation laws of time irreversibility of HRV in the early stage of hypoxic exposure, experiments will be carried out in a normobaric hypoxic chamber. By introducing different delay times, time irreversibility variation rules of HRV will be studied in an acute 2 h hypoxic exposure. Few reports have been posted on the research of time irreversibility of HRV in a hypoxic environment.
Materials and Methods
Six young and healthy male volunteers were recruited for our experiment. Over 24 h stay at an altitude >3,000 m in the previous 6 months, born at an altitude >1,000 m, being a smoker, or having a history of severe respiratory or cardiopulmonary diseases are the exclusion criteria. All subjects were forbidden to consume tea, coffee, alcohol, medicine that could get the neural system excited, or do intense physical exercise 24 h prior to the experiment. Characteristics of the subjects are presented in Table 1. The study was approved by the Ethics Committee of Southeast University, China. All the subjects were well informed of the aim and risks of the experiment and they all signed their informed consent before the study.
Experimental Design and Data Acquisition
This experiment was implemented in a normobaric hypoxic chamber. Before the experiment, nitrogen was filled into the chamber to reduce the oxygen concentration to 12.9% (equivalent to the altitude of 3,600 m), then subjects entered the chamber. During the experiment, the nitrogen flow was controlled by an automatic valve (Yu et al., 2011), so that the oxygen concentration in the chamber was maintained at 12.9% ± 0.2%. The experiment lasted for 2 h, and subjects were asked to keep quiet and lie in their bunks. They could watch cellphones, read books, and listen to gentle music. ECG signals were recorded by tiny equipment inlaid in a chest belt for each subject (SensEcho-5A, Health Regulation Co. Ltd., Beijing, China). The sampling rate was 200 Hz, and digitized ECG signals were stored in a built-in hard disk, which is capable of storing 24 h of data. The temperature and humidity in the chamber were maintained at 22 ± 1°C and 25 ± 1%, respectively.
Data Processing and Calculating Method
First, we identified each QRS complex automatically by P&T methods (Pan and Tompkins, 1985) for its extensively tested accuracy and efficiency. False and missing detection were calibrated artificially. Then we extracted R-R intervals series by corrected R wave positions.
Porta et al. (2006) counted the number of times that the difference of one RR interval with the next was less than zero and when the number was unequal to zero; the two numbers’ ratio was defined as time irreversibility index P%. Guzik et al. (2006) performed arithmetical operations on the square of the difference between adjacent RR intervals and defined time irreversibility index G%. We denote any element in RR intervals sequence as RRi(1≤i≤n), where n is length of RR intervals sequence. The difference between adjacent RR intervals is shown as ΔRR = RRi−RRi−1. ΔRR is also expressed as ΔRR+ when ΔRR > 0, and ΔRR− when ΔRR < 0. N(ΔRR+) is the number of ΔRR+ in the entire ΔRR sequence. Delay time τ is introduced in this study. For τ = 1,2,…,n−1, ΔRRτ = xi + τ−xi(1≤i≤n−τ), time irreversibility indices P%(τ) and G%(τ) for each τ are given by
Distance between P%(τ), G%(τ) and 50 are expressed as
In order to make comparisons with the hypoxic group, we randomly selected six items of 2 h from the MIT-BIH Normal Sinus Rhythm Database (Goldberger et al., 2000) to perform the same time irreversibility analysis on ECG data. We will take QPτ and QGτ as the measurement of time irreversibility.
Since QPτ and QGτ are all statistics independent of the sequence length, we calculated each index value from the beginning of the experiment and moved forward by 1 min each time. The first computation scope is [0,5 min]. That is, we calculated the indexes in the scope of [0,5 min], [0,6 min], [0,7 min], …, [0,120 min], respectively, and drew line graphs. In this way, data volume for calculation is gradually increased and the curve is relatively flat. Interference for observation and analysis because of a short or sharp change in a certain data section can be avoided.
Variation of QP1, QG1, QP2, and QG2
If delay time was set to be 1 or 2 (as shown in Figure 1), indexes of all subjects in the hypoxic group fluctuated greatly in the first 15 min after the experiment began, which is possibly due to the data volume for calculation being small meaning the autonomic nervous system of the human body must undergo the process of adjustment and adaptation when first entering the hypoxic environment. In the entire experiment, the trends of QP1, QG1, QP2, and QG2 are generally toward zero, which indicates time irreversibility of HRV of the hypoxic group tends to decline.
Figure 1. (A–D) show the corresponding results of QP1, QP2, QG1, and QG2. Variation tendency of QP1, QP2, QG1, and QG2 within 2 h of hypoxic group (each color represents a subject, the same below).
Different from the hypoxic group, which showed a significant downward trend of QP1, QG1, QP2, and QG2, indicators of the MIT normoxic group did not show an obvious consistent changing direction (as shown in Figure 2).
Figure 2. (A–D) show the corresponding results of QP1, QP2, QG1, and QG2. Variation tendency of QP1, QP2, QG1, and QG2 within 2 h of MIT normoxic group.
If we calculate indicators every 15 min and take both τ = 1 and τ = 2 into consideration, we can draw an error bar chart with the average value of QP1 and QP2 (as shown in Figure 3), which clearly shows the average time irreversibility of all subjects with a delay time of 1 and 2 in different groups. As in the charts above, indicators of the hypoxic group decreases and the MIT normoxic group shows no clear differences.
Figure 3. Average value error bar chart of QP1 and QP2 of two groups. (A) Error bar chart of hypoxic group. (B) Error bar chart of MIT normoxic group.
Variation of QP3, QG3, QP4, and QG4
If delay time was set to be 3 or 4 (as shown in Figure 4), a great difference in variation trend happened in QP and QG from that when delay time is 1 or 2. Indicators no longer show a slow downward trend toward zero with time going by, but gradually tend to be stable after the drastic transition period of the first 30 min. And most of the subjects’ indicator curves gradually converge to a narrow range. In contrast, there was no significant trend change in the MIT normoxic group (as shown in Figure 5).
Figure 4. (A–D) show the corresponding results of QP3, QP4, QG3, and QG4. Variation tendency of QP3, QP4, QG3, and QG4 within 2 h of hypoxic group.
Figure 5. (A–D) show the corresponding results of QP3, QP4, QG3, and QG4. Variation tendency of QP3, QP4, QG3, and QG4 within 2 h of MIT normoxic group.
Research on time irreversibility of HRV in a hypoxic environment has not been carried out before. Previous studies have shown that age and pathology will lead to a decrease of time irreversibility of HRV (Costa et al., 2005) and, compared to the daytime, time irreversibility of both healthy people and patients with congestive heart failure (CHF) declines at night (Porta et al., 2009). Our study in a hypoxic environment reveals a similar rule: prolonged hypoxic exposure leads to a sustained decline of time irreversible of HRV within 2 h when the delay time is 1 or 2.
Hypoxia is one of induction factors that can make regulatory mechanisms of cardiovascular nervous system change: parasympathetic activity was significantly reduced and autonomic nervous system regulation was greatly inhibited (Zhu et al., 2010). Since healthy physiological systems have the strongest nonlinearity and complexity, the decrease of time irreversibility in an anoxic state indicates a depression of nonlinearity and complexity.
In RR interval series, the positive difference between adjacent RR intervals represents a sudden increase of heart rate due to sympathetic modulation, whereas the negative difference represents a decrease in instantaneous heart rate due to parasympathetic modulation. The indicators variation of an MIT normoxic group reveals that there is always a certain difference of regulation times between sympathetic and parasympathetic nerves in normal individuals under normal oxygen concentration, while the difference tends to diminish in hypoxic environments. By definition, QP1 reflects the absolute difference of regulation times between sympathetic and parasympathetic nerves, and QG1 is a measurement of regulation intensity by calculating the square of RR intervals difference. The fact that both of them are approaching zero indicates that regulation times and intensity of sympathetic and parasympathetic nerves tend to be close to each other.
The LancetIn the present study, we observed a decrease in time irreversibility of the hypoxic group when the delay time was 1 and 2, but the variation trend changed when the delay time was increased to 3 or 4. It can be seen that the trend of QP3, QG3, QP4, and QG4 in the hypoxic group gradually becomes stable and converges with each other, while there is no significant difference whatever the delay time is of the MIT normoxic group, which indicates that time irreversibility of HRV is closely related to the delay time in a hypoxic environment. The results of mutual information analysis shows that the delay time of RR interval sequence in the elderly was significantly longer than that in the young (Ma and Zhang, 2010). Further research is needed to find out what generates the difference of delay time between hypoxic and normoxic environments. Hou et al. (2011) defined the average value at different delay times of QP and QG as the high dimensional time irreversibility index Pm, Gm, and the square root of their sum of squares was defined as their comprehensive measurement index Dm. Her study pointed out that the relationship between age, CHF, and time irreversibility of HRV was affected by the embedded dimension. When the value range of the embedded dimension changed, the variation rule revealed by time irreversibility indicators also changed. If the embedding dimension is m, the essential meaning of Pm and Gm is the average value of QP and QG when the delay time is taken from 1 to m. In the multi-level control loop and feedback network system of the human body, they can reflect the average effect of signals derived from different levels and different delays, but researchers cannot make observations on the specific value of time irreversibility index under a certain delay time. When the value of delay time increases continuously, variation of Pm and Gm must become gentler and have the characteristics of lag. It was found that time irreversibility of HRV during exercise was significantly increased when taking multiple dimensions into consideration (Hou et al., 2013). Time irreversibility of HRV in a short period of time at the beginning of hypoxic exposure was investigated in this paper. Studies have shown that intermittent hypoxic training over several days significantly improves SDNN, RMSSD, LF HF (Xu et al., 2004), and SPO2 (Fulco et al., 2013), however, the impact on time irreversibility needs further research.
Variation laws of time irreversibility of HRV in a normobaric hypoxic environment were researched in the present study. We found that in 2 h hypoxic exposure, time irreversibility indicators QP1, QG1, QP2, and QG2 decrease continuously, while QP3, QG3, QP4, and QG4 gradually tend to be stable and converge to a narrow range, which means there is a relationship between time irreversibility of HRV and delay time, and that proper selection of delay time is important for the observation of the relationship. In contrast, normoxic environments do not have any influence on time irreversibility of HRV whatever the delay time is. Results when delay time is 1 or 2 suggests a decline in automatic nervous system activities in a hypoxic environment, and regulation times and intensity of sympathetic and parasympathetic nerves get close to each other.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving human participants were reviewed and approved by the Ethics Committee of Southeast University, Nanjing, China.
YL did all the analysis work and wrote the article. JqL checked research protocol. JL performed the data calibration. YX conducted the experiment and collected data. ZC provided the experimental site and equipment. CL guided the research method, wrote the article and examined the research results. All authors contributed to the article and approved the submitted version.
This work was supported by the National Key Research and Development Program of China (2019YFE0113800), theNational Natural Science Foundation of China (81871444 and 62071241), and the Distinguished Young Scholars of Jiangsu Province (BK20190014).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful for the experimental site and equipment provided by the Aeromedicine Engineering Center in Air Force Medical Center.
Beidleman, B. A., Muza, S. R., Fulco, C. S., Cymerman, A., Sawka, M. N., Lewis, S. F., et al. (2008). Seven intermittent exposures to altitude improves exercise performance at 4300 m. Med. Sci. Sports Exerc. 40, 141–148. doi: 10.1249/mss.0b013e31815a519b
Casali, K. R., Casali, A. G., Montano, N., Irigoyen, M. C., Macagnan, F., Guzzetti, S., et al. (2008). Multiple testing strategy for the detection of temporal irreversibility in stationary time series. Phys Rev. E Stat. Nonlin. Soft. Matter. Phys. 77(6 Pt 2):66204. doi: 10.1103/PhysRevE.77.066204
Costa, M. D., Goldberger, A. L., and Peng, C. K. (2005). Broken asymmetry of the human heartbeat: loss of time irreversibility in aging and disease. J. Electrocardiol. 38(Suppl. 4):141. doi: 10.1016/j.jelectrocard.2005.06.076
Coustet, B., Lhuissier, F. J., Vincent, R., and Richalet, J. (2015). Electrocardiographic changes during exercise in acute hypoxia and susceptibility to severe high-altitude illnesses. Circulation 131, 786–794. doi: 10.1161/CIRCULATIONAHA.114.013144
Goldberger, A., Amaral, L., Glass, L., Hausdorff, J., Ivanov, P. C., Mark, R., et al. (2000). PhysioBank, physiotoolkit, and physionet: components of a new research resource for complex physiologic signals. Circulation [Online] 101, 215–220.
Graff, G., Graff, B., Kaczkowska, A., Makowiec, D., Amigó, J. M., Piskorski, J., et al. (2013). Ordinal pattern statistics for the assessment of heart rate variability. Eur. Phys. J. Spec. Top. 222, 525–534. doi: 10.1140/epjst/e2013-01857-4
Hou, F. Z., Huang, X. L., Zhuang, J. J., Huo, C. Y., and Ning, X. B. (2012). Multi-scale strategy and data surrogating test: two elements for the detection of time irreversibility in heart rate variability. Acta Phys. Sin. 61:220507. doi: 10.7498/aps.61.220507
Guzik, P., Piskorski, J., Krauze, T., Wykretowicz, A., and Wysocki, H. (2006). Heart rate asymmetry by Poincaré plots of RR intervals. Biomed. Tech. Biomed. Eng. 51, 272–275. doi: 10.1515/BMT.2006.054
Hou, F. Z., Ning, X. B., Zhuang, J. J., Huang, X. L., Fu, M. J., and Bian, C. H. (2011). High-dimensional time irreversibility analysis of human interbeat intervals. Med. Eng. Phys. 33, 633–637. doi: 10.1016/j.medengphy.2011.01.002
Iwasaki, K., Ogawa, Y., Aoki, K., Saitoh, T., Otsubo, A., and Shibata, S. (2006). Cardiovascular regulation response to hypoxia during stepwise decreases from 21% to 15% inhaled oxygen. Aviat. Space Environ. Med. 77, 1015–1019.
Krejčí, J., Botek, M., and McKune, A. J. (2018). Dynamics of the heart rate variability and oxygen saturation response to acute normobaric hypoxia within the first 10 min of exposure. Clin. Physiol. Funct. Imaging 38, 56–62. doi: 10.1111/cpf.12381
Porta, A., D’Addio, G., Bassani, T., Maestri, R., and Pinna, G. D. (2009). Assessment of cardiovascular regulation through irreversibility analysis of heart period variability: a 24 hours Holter study in healthy and chronic heart failure populations. Philos. Trans. A Math. Phys. Eng. Sci. 367, 1359–1375. doi: 10.1098/rsta.2008.0265
Porta, A., Guzzetti, S., Montano, N., Gnecchi-Ruscone, T., Furlan, R., and Malliani, A. (2006). “Time reversibility in short-term heart period variability,” in Proceedings of the 2006 Computers in Cardiology (Valencia: IEEE), 77–80.
Vigo, D. E., Perez, L. S., Videla, A. J., Perez, C. D., Hunicken, H. M., Mercuri, J., et al. (2010). Heart rate nonlinear dynamics during sudden hypoxia at 8230 m simulated altitude. Wilderness Environ. Med. 21, 4–10. doi: 10.1016/j.wem.2009.12.022
Xue, Y., Yang, J., Feng, Y. T., Zhou, Y. B., Qin, Y. F., Li, Y., et al. (2010). Effects of mindful breathing on rapid hypoxia preacclimatization training. J. Med. Imaging Health Inform. 10, 718–723. doi: 10.1166/jmihi.2020.2923
Yao, W., Yao, W., and Wang, J. (2019). Equal heartbeat intervals and their effects on the nonlinearity of permutation-based time irreversibility in heart rate. Phys. Lett. A 383, 1764–1771. doi: 10.1016/j.physleta.2019.03.002
Yu, M. S., Yang, J., Wang, H. T., Miao, Y., Wang, N. Z., and Li, Y. (2011). Normobaric Hypoxic Cabin for Intermittent Hypoxic Training and Method for Manipulate Oxygen Concentration. CN Patent No 102198321 B. Beijing: State Intellectual Property Office of the People’s Republic of China.
Zeng, N., Wang, Z., Li, Y., and Du, M. (2013). Time series modeling of nano-gold immunochromatographic assay via expectation maximization algorithm. IEEE Trans. Biomed. Eng. 60, 3418–3424. doi: 10.1109/TBME.2013.2260160
Zeng, N., Wang, Z., Zhang, H., Kim, K.-E., Li, Y., and Liu, X. (2019). An improved particle filter with a novel hybrid proposal distribution for quantitative analysis of gold immunochromatographic strips. IEEE Trans. Nanotechnol. 18, 819–829. doi: 10.1109/TNANO.2019.2932271
Zhang, D., She, J., Zhang, Z., and Yu, M. (2014). Effects of acute hypoxia on heart rate variability, sample entropy and cardiorespiratory phase synchronization. Biomed. Eng. Online 13:73. doi: 10.1186/1475-925X-13-73
Keywords: time irreversibility, hypoxia, heart rate variability, delay time, autonomic nervous system
Citation: Li Y, Li J, Liu J, Xue Y, Cao Z and Liu C (2021) Variations of Time Irreversibility of Heart Rate Variability Under Normobaric Hypoxic Exposure. Front. Physiol. 12:607356. doi: 10.3389/fphys.2021.607356
Received: 17 September 2020; Accepted: 15 February 2021;
Published: 05 March 2021.
Edited by:Nianyin Zeng, Xiamen University, China
Copyright © 2021 Li, Li, Liu, Xue, Cao and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.