Minute ventilation and heart rate relationship for estimation of the ventilatory compensation point at high altitude: a pilot study
Extreme Physiology & Medicine volume 2, Article number: 7 (2013)
The ventilatory compensation point (VCP) is an exercise threshold which has been used in the design of training programs in sports medicine and rehabilitation. We recently demonstrated that changes in the slope of the minute ventilation to heart rate relationship () can be utilized for estimation of the VCP during incremental exercise at sea level (SL). We hypothesized that in hypoxic conditions, such as high altitude (HA), VCP can be also reliably estimated by .
At SL and on immediate ascent to HA (5,050 m), six healthy subjects (42 ± 14 SD years) performed a maximal incremental exercise test on a cycle ergometer; O2 uptake (), CO2 output (), , and HR were measured breath-by-breath. The method for VCP estimation was compared to the standard method using the ventilatory equivalent for CO2 () and end-tidal PCO2 (PETCO2). The slope values below (S1) and above (S2) VCP were computed by linear regression analysis.
A significant difference between S1 and S2 was observed, at SL and HA, for both the and methods for VCP estimation. A good agreement between the two methods ( vs. ) was found for both environmental conditions; the mean difference ± 2 SD of at VCP (VCP-) was −22 ± 112 ml/min at SL and 39 ± 81 ml/min at HA. The VCP- was significantly lower at HA compared to SL; in addition, S1 and S2 mean values were significantly higher at HA compared to SL.
At HA, VCP may be reliably estimated by the method.
In normal individuals, minute ventilation () during exercise displays a more marked rapid increase when work rate exceeds the heavy-intensity domains, as a compensatory effect of the metabolic acidosis [1, 2], and this threshold is called the ventilatory compensation point (VCP) . The VCP has been utilized for the design of training programs in sports medicine and rehabilitation [4, 5]. Practical and, in some instances, simple approaches for estimation of the VCP have been evaluated in the past by different authors. Several strategies have been proposed, based either on  or heart rate (HR) responses [7, 8].
Since HR is usually reasonably linear with respect to oxygen uptake () over the entire work rate range , we have argued that the changes in the slope of increment in over HR response () can be utilized to estimate VCP during incremental exercise at sea level (SL) [10, 11]. The method has shown a good agreement with the standard method for VCP estimation (i.e., increased ventilatory equivalent for CO2) and decreased end-tidal PCO2 (PETCO2), independently of exercise modality  and protocol . The method therefore may be a potentially promising method to estimate VCP during field tests or in difficult environmental conditions.
At high altitude (HA), compared to SL, is greatly increased at any given work rate . This may have a detrimental effect on exercise tolerance because it may lead to a significant reduction in ventilatory reserve ; however, reduced regional blood flow and oxygen diffusion limitation in the lungs and exercising muscles are also believed to play an important role in limiting exercise tolerance at HA . The HR response at any given work rate is also increased at HA , although maximal HR appears to be reduced [13–16]. No data, however, are available on the slopes of the response to exercise at HA.
In the present study, we wanted to establish the reliability of the method for the estimation of VCP at HA, an environmental condition associated with an augmented ventilatory and cardiac response to exercise. We hypothesized that the method at HA would be as reliable as the standard method for VCP estimation. If so, this could lead to the development of less expensive exercise equipment, i.e., without gas analyzers, that could readily be easily used in sports medicine and in high altitude studies where calibaration of gas analyzers is more difficult.
Six healthy non-smoking SL resident subjects (two females, four males; Table 1) provided written informed consent to participate in the study. The procedures and protocols were approved by the National Research Council (CNR, Italy) as part of the Ev-K2-CNR research program and were conducted in accordance with the Declaration of Helsinki.
The SL phase was conducted in Rome, Italy, at approximately 60 m of altitude. On day 1, the subjects performed standard pulmonary function tests and got familiarized with the maximal incremental exercise test. On the subsequent day, the formal maximal incremental exercise test was conducted until exhaustion.
The HA phase of the study was conducted at the CNR Pyramid Laboratory, Lobuche, Khumbu, Nepal (5,050 m, barometric pressure approximately 410 Torr). After travelling by plane through Kathmandu (1,340 m) up to Lukla (2,860 m), all subjects performed a 10-day trek to the Pyramid Laboratory with the following intermediate stops for acclimatization: 3 days at Namche Bazaar (3,450 m) and then 2 days at Pheriche (4,252 m). Symptoms of acute mountain sickness and resting arterial O2 saturation (SpO2) were monitored. After 1 day of rest at the Pyramid Laboratory, the experimental phase started. On day 2, the subjects performed standard pulmonary function tests. Calibration of the spirometer (photoelectric digital turbine, diameter 28 mm, resolution 4 ml, Micro Kit, COSMED, Rome, Italy) was performed prior to each test, using a 3-l syringe.
Incremental exercise test
On day 3, each subject performed a maximal incremental exercise test on a cycle ergometer. The pedaling frequency was set at 60 ± 5 rpm by subjects following an electronic audio signal generated each second. The test protocol consisted (a) 2 min of rest, (b) 4 min of exercise at 20 W, (c) the incremental phase with work rate increments of 15 W/min at HA and 25 W/min at SL, (d) 6 min of recovery pedalling at 20 W. O2 uptake (, standard temperature and pressure, dry (STPD)), CO2 output (, STPD), minute ventilation (, body temperature and pressure saturated BTPS), and end-tidal partial pressures for O2 and CO2 (PETO2, PETCO2) were obtained breath-by-breath and edited to exclude occasional outlying breaths (>±4 SD of the local mean) as a result of coughs, swallows, sighing, or gasping. A 10-s average was utilized for subsequent analysis.
For VCP estimation by analysis of the relationship (), we utilized a least squares regression analysis. Using a dedicated software, we applied a ‘best fit’ line to the data which extended from end-exercise back to the sub-maximal point at which the linearity was lost (S2 region); this breakpoint was compared with the standard approach for VCP estimation () based on the relationship, where started to change out of proportion of and PETCO2 started to fall . From the aforementioned breakpoint, a second best fit line (S1) was applied to the data extending from the end of the warm-up phase upwards into the exercise data. If a second breakpoint was discernible, the S1 line was then divided in two different regions (S1” from the breakpoint to and S1’ from the end of warm-up up to the observed breakpoint), and the HR value at this breakpoint was then compared with the value observed at lactate threshold (θL, estimated by the ‘V-slope method’  and supported by standard ventilatory equivalent and end-tidal gas tension criteria ).
The incremental exercise test was performed on a mechanically braked cycle ergometer (828E, Monark Exercise AB, Varberg, Sweden). Prior to each test, the zero-load setting on the ergometer was checked, and a calibration was performed with a 4-kg weight. Ventilatory and pulmonary gas exchange variables were measured breath-by-breath in all tests using a portable system (K4b2, COSMED) which has been previously validated at HA ; the accuracy of the telemetric system has been previously established [18, 19]. The system comprised a face mask, analyzer unit (containing O2 and CO2 gas analyzers), heart rate monitor, and battery. The analyzer unit with battery pack, face mask, and tubing (weight 0.8 kg) was attached to the subject with a harness and connected to a personal computer by an Ethernet cable connection. The face mask contained a turbine for measurement of volume and flow; calibration was performed with a 3-l syringe (Hans Rudolph, Kansas City, MO, USA) over a range of different flow profiles. Respired gas, sampled continuously from a port within the turbine via a Nafion polymer capillary (PermaPure©, Toms River, NJ, USA), was analyzed at 100 Hz using rapid-response O2 (polarographic) and CO2 (infrared) analyzers (mean response time 120 ms) which were automatically thermostated and compensated for ambient variations in barometric pressure, humidity, and environmental temperature. Analyzers' calibration was performed using two precision-analyzed gas mixtures spanning the respired range. The volume and gas concentration signals were sampled and digitized every 10 ms, and time-aligned, i.e., correcting for the transport delay between the turbine and gas analyzers and for the analyzer rise time . HR was measured from a chest strip and recorded every breath. SpO2 was monitored non-invasively by finger pulse oximetry (Masimo Rad-5, Masimo Corporation, Irvine, CA, USA).
Differences among measured responses were determined by a Student's paired t test. Pearson's product–moment correlation coefficient (R2) was used to identify correlations between criterion variables. The level of statistical significance was set at P < 0.05. Group data are presented as mean ± SD. The limits of agreement between the and the methods for VCP estimation were evaluated by the Bland-Altman analysis , where the individual differences are plotted against their respective means. We proceeded with such a type of analysis if a significant linear correlation between methods was previously observed. The same statistical approach was also performed to compare the correspondence between θL and the breakpoint between S1” and S1’ in the relationship, when the latter was detectable.
The main results of the incremental tests are illustrated in Table 2. At HA, compared to SL, a significant reduction in , WRpeak, HRpeak, SpO2peak, and θL was observed; on the contrary, at HA was appreciably higher.
Figures 1 and 2 show , , PETO2, PETCO2 vs. , and relationship in two representative subjects at SL (upper panels) and HA (lower panels). As shown in Figure 3, in all subjects, a breakpoint in , which occurred at the VCP estimated by the method, was clearly discernible both at SL and HA. No significant differences were found in measured at VCP ( -VCP; Table 3) between methods utilized to identify the threshold (VCP- vs. VCP-).
The Bland-Altman analysis confirmed the agreement between the and methods in estimating VCP, both at SL and at HA (Figure 4, Table 3); the mean bias ± 95% confidence interval of the between-method differences were −22 ± 112 ml/min of for the SL protocol and 39 ± 81 ml/min of for the HA protocol (Figure 4).
Although in some instances a lower breakpoint in the S1 region was observed (four out of six subjects at SL and one out of six at HA), no significant agreement between such a breakpoint and the θL was obtained in terms of HR and values. -VCP was significantly lower at HA (approximately 500 ml of less) independently of the method utilized to estimate it, as shown in Table 2.
S1 and S2 at HA were significantly higher for both and methods (Table 3). The ratio S2/S1 tended to be higher at HA regardless of the estimation method utilized (Table 3), and this was particularly so for the method.
Discussion and conclusions
The main findings of this study are as follows: (1) VCP is clearly identifiable during incremental exercise at HA; (2) in conditions of persistent hypoxia, such as that elicited by HA exposure, VCP can be easily estimated by the analysis. In the present study on young healthy individuals, applying standard methods for the identification, VCP could be identified both at SL and HA. Our data are in accordance with the observation of Agostoni et al., who demonstrated that at HA the VCP was clearly identifiable .
Of note, we observed a 30% reduction in -VCP when moving from SL to HA. The lower VCP at HA likely reflected the shortening of the isocapnic buffering period because of the increased hypoxic sensitivity of the carotid chemoreflex . A further influence was presumably the premature development of metabolic acidosis during exercise at HA, as demonstrated by the lower θL likely due both to the reduction in arterial O2 content and impaired convective O2 transport to the exercising muscles .
Importantly, as demonstrated by the Bland-Altman analysis, we were able to demonstrate that the approach is a reliable method for VCP estimation during incremental exercise at HA. The limits of agreement and the coefficient of variation between the and methods for VCP estimation are quite narrow and likely to fall within the between-day intra-subject variability .
The results of the present study are in agreement with a previous work of our group that demonstrated the reliability of the breakpoint for VCP estimation at SL during incremental exercise [10, 11]. At SL, assuming a linear relationship between and HR (at least up to the heavy-intensity domain), the steepening of the slope is more evident above VCP where (1) a steeper increase in relationship is expected  and (2) the HR response could either maintain the same linearity as for more moderate exercise or in some instances (depending on exercise performance and protocol) reduce its rate of increment with respect to work rate . Interestingly, at HA, we observed higher S1 and S2 slopes. These findings are in accordance with previous studies demonstrating that at rest and at any given work rate during exercise is increased at HA because of the greater hypoxic drive [12–26] and also because of the reduction in the rate of increase in HR between rest and peak exercise that is commonly observed at HA [15, 16]. Moreover, compared to SL, we observed an increase in the S2/S1 ratio for at HA. This finding suggests that S2 is influenced both by hypoxic ventilatory drive and possibly by an augmented contribution to from the metabolic acidosis above the VCP . However, our reasoning remains speculative as, for technical reasons, we were unable to measure arterial (or arterialized) [lactate] or pH during the exercise.
Below VCP, an earlier (in the S1 region) but less evident breakpoint in was discernible in some instances (5 out of 12), but only in one out five cases such a breakpoint was coincident with the anaerobic threshold.
We are aware that our study has limitations, particularly with regard to the small sample size and the constrained characteristics of our study population, which prevented the evaluation of factors such as age, sex, and level of fitness, each of which is known to influence the variables of interest. Also, we are aware (1) that the choice of a rapid work rate incremental protocol may have influenced the results utilized, i.e., different results may be observed if a slower work rate increment protocol is used, and (2) of the lack of validation for pulmonary gas exchange measurement at HA. We took precautions to limit errors in the calibration procedures, with all equipment being calibrated before each test and all tests being performed indoors in the Pyramid Laboratory at an ambient temperature of approximately 20°C and the necessary corrections for barometric pressure, humidity, and environmental temperature being applied through the dedicated software. Thus, our investigation should be considered as a pilot study conducted at extreme altitude, with the conclusions only being applicable to healthy untrained young adults.
In conclusion, we were able to demonstrate the reliability of the method for the estimation of VCP at HA in a small group of sea-level residents, an environmental condition associated with different ventilatory and cardiac responses to exercise compared to SL. The attraction of the method is that it is a less expensive method for VCP estimation, compared to those that utilize expired gas measurements, and may therefore be utilized in sports medicine and in extreme conditions such as high altitude.
Wassermann K, Whipp BJ, Koyal SN, Beaver WL: Anaerobic threshold and respiratory gas exchange during exercise. J Appl Physiol. 1973, 35 (2): 236-243.
Meyer T, Faude O, Scharhag J, Urhausen A, Kindermann W: Is lactic acidosis a cause of exercise induced hyperventilation at the respiratory compensation point?. Br J Sports Med. 2004, 38 (5): 622-625. 10.1136/bjsm.2003.007815.
Beaver WL, Wasserman K, Whipp BJ: A new method for detecting anaerobic threshold by gas exchange. J Appl Physiol. 1986, 60: 2020-2027.
Azevedo LF, Perlingeiro PS, Brum PC, Braga AM, Negrão CE, de Matos LD: Exercise intensity optimization for men with high cardiorespiratory fitness. J Sports Sci. 2011, 29 (6): 555-561. 10.1080/02640414.2010.549613.
Carvalho VO, Mezzani A: Aerobic exercise training intensity in patients with chronic heart failure: principles of assessment and prescription. Eur J Cardiovasc Prev Rehabil. 2011, 18 (1): 5-14.
Neder JA, Stein R: A simplified strategy for the estimation of the exercise ventilatory thresholds. Med Sci Sports Exerc. 2006, 38 (5): 1007-1013. 10.1249/01.mss.0000218141.90442.6c.
Conconi F, Ferrari MP, Ziglio G, Droghetti P, Codeca L: Determination of the anaerobic threshold by a noninvasive field test in runners. J Appl Physiol. 1982, 52: 869-873.
Karapetian GK, Engels HJ, Gretebeck RJ: Use heart rate variability to estimate LT and VT. Int J Sports Med. 2008, 29 (8): 652-657. 10.1055/s-2007-989423.
Wassermann K, Hansen JE, Sue DY, Casaburi R, Whipp BJ: Principles of Exercise Testing and Interpretation. 2005, Baltimore: Lippincott Williams & Wilkins, 10-65. 4
Onorati P, Martolini D, Ora J, Valli G, Fedeli A, Palange P: Estimation of the exercise ventilatory compensation point by the analysis of the relationship between minute ventilation and heart rate. Eur J Appl Physiol. 2008, 104 (1): 3-17.
Onorati P, Martolini D, Valli G, Laveneziana P, Marinelli P, Angelici E, Palange P: A simplified approach for the estimation of the ventilatory compensation point. Med Sci Sports Exerc. 2012, 44 (4): 716-724. 10.1249/MSS.0b013e318237fb5f.
Sutton JR, Reeves JT, Wagner PD, Groves BM, Cymerman A, Malconian MK, Rock PB, Young PM, Walter SD, Huston CS: Operation Everest II: oxygen transport during exercise at extreme simulated altitude. J Appl Physiol. 1988, 64: 1309-1321.
Valli G, Cogo A, Passino C, Bonardi D, Morici G, Fasano V, Agnesi M, Bernardi L, Ferrazza AM, Ward SA, Palange P: Exercise intolerance at high altitude (5050m): critical power and W’. Respiratory Physiol & Neurobiol. 2011, 177 (3): 333-341. 10.1016/j.resp.2011.05.014.
Anholm JD, Foster GP: Con: Hypoxic pulmonary vasoconstriction is not a limiting factor of exercise at high altitude. High Alt Med Biol. 2011, 12 (4): 313-317. 10.1089/ham.2011.1059.
Vogel JA, Hartley LH, Cruz JC, Hogan RP: Cardiac output during exercise in sea level residents at sea level and high altitude. J Appl Physiol. 1974, 36: 169-172.
Reeves JT, Groves BM, Sutton JR, Wagner PD, Cymermann A, Malconian MK, Rock PB, Young PM, Houston CS: Operation Everest II: preservation of cardiac function at extreme altitude. J Appl Physiol. 1987, 63: 531-539.
American Thoracic Society/American College of Chest Physicians: ATS/ACCP Statement on cardiopulmonary exercise testing. Am J Respir Crit Care Med. 2003, 167: 211-277.
Pattinson K, Myers S, Gardner-Thorpe C: Apparatus. Problems with capnography at high altitude. Anesthesia. 2004, 59: 69-72. 10.1111/j.1365-2044.2004.03447.x.
Palange P, Forte S, Onorati P, Manfredi F, Carlone S, Serra P: Ventilatory and metabolic adaptation to walking and cycling in patients with COPD. J Appl Physiol. 2000, 88: 1715-1720.
Beaver WL, Wasserman K, Whipp BJ: On-line computer analysis and breath-by-breath graphical display of exercise function tests. J Appl Physiol. 1973, 34: 128-132.
Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1 (8476): 307-310.
Agostoni P, Valentini M, Magrì D, Revera M, Caldara G, Gregorini F, Bilo G, Styczkiewicz K, Savia G, Parati G: Disappearance of isocapnic buffering period during increasing work rate exercise at high altitude. Eur J Cardiovasc Prev Rehabil. 2008, 15: 354-358. 10.1097/HJR.0b013e3282f62982.
Lundby C, Sander M, van Hall G, Saltin B, Calbet JA: Maximal exercise and muscle oxygen extraction in acclimatizing lowlanders and high altitude natives. J Physiol. 2006, 573: 535-547. 10.1113/jphysiol.2006.106765.
Wergel-Kolmert U, Agehäll A, Rosenberg N, Wohlfart B: Day-to-day variation in oxygen consumption at submaximal loads during ergometer cycling by adolescents. Clin Physiol. 2001, 19 (2): 161-168.
Bodner ME, Rhodes EC: A review of the concept of the heart rate deflection point. Sport Med. 2000, 30 (1): 31-46. 10.2165/00007256-200030010-00004.
Richalet JP, Robach P, Jarrot S, Schneider JC, Mason NP, Cauchy E, Herry JP, Bienvenu A, Gardette B, Gortan C: Operation Everest III (COMEX’97): effects of prolonged and progressive hypoxia on humans during a simulated ascent to 8848m in a hypobaric chamber. Adv Exp Med Biol. 1988, 74: 323-333.
This study was carried out within the framework of the Ev-K2-CNR Project in collaboration with the Nepal Academy of Science and Technology, and the authors wish to thank the Italian National Research Council and the Italian Ministry of Foreign Affairs for their financial assistance.
The authors declare that they have no competing interests.
GV was the principal investigator, conceived the study, participated in all phases of the research, and drafted the manuscript. MI and AMF participated in the design of the study and in the data collection and performed the statistical analysis. PO helped conceive the study, participated in its design and coordination, and helped draft the manuscript. AC participated in the logistic organization, coordinated the research during the HA phase of the study that was conducted at the CNR Pyramid Laboratory, Lobuche, Khumbu, Nepal, and critically revised the final draft. PP helped conceive the study, supervised all phases of the research, helped in its design and coordination, and revised critically the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Valli, G., Internullo, M., Ferrazza, A.M. et al. Minute ventilation and heart rate relationship for estimation of the ventilatory compensation point at high altitude: a pilot study. Extrem Physiol Med 2, 7 (2013). https://doi.org/10.1186/2046-7648-2-7