 Research article
 Open Open Peer Review
 Published:
Quantification and physiological significance of the rightward shift of the Vslope during incremental cardiopulmonary exercise testing
BMC Sports Science, Medicine and Rehabilitationvolume 9, Article number: 9 (2017)
Abstract
Background
Ventilatory anaerobic threshold (VAT) is frequently used as a measure of exercise tolerance, with the Vslope method being the standard; however, this needs to be visually determined. Over the years, we have observed that the Vslope itself often appears to shift rightward before the appearance of the VAT (RtShift: rightward shift of Vslope). This phenomenon has long been known to occur during the first 1–2 min of steadystate exercise and disappears thereafter; it is attributed to CO_{2} storage, presumably in active muscle. However, during incremental exercise, we have observed that the RtShift persists; furthermore, it seems to be related to the level of VAT. Therefore, we attempted to objectively quantify the RtShift, and to confirm its relationship to an index of exercise tolerance (VAT).
Methods
This study was based on a retrospective analysis of data from 100 cardiopulmonary ramp exercise tests (submaximal) performed by patients with cardiac disease. VAT was determined with the visual Vslope method. The horizontal distances between the diagonal R = 1 line and each data point on the Vslope plot to the right of R = 1 were measured; the average of these measurements was used as an objectively determined estimate of RtShift.
Results
The predominant portion of RtShift occurred earlier than VAT. The mean RtShift was 33.9 ± 25.0 mL⋅min^{−1} VO_{2}, whereas the mean VAT was 635 ± 220 mL⋅min^{−1}. RtShift positively correlated with VAT (r = 718, p < 0.001), confirming previous visual observations. It also significantly correlated with ΔVO_{2}/Δwork rate, a marker of oxygen uptake efficiency (r = 0.531, p < 0.001).
Conclusions
We identified that among patients with cardiac disease, Vslope is shifted rightward to varying degrees. The objectively quantified rightward shift of Vslope is significantly correlated with an index of exercise tolerance (VAT). Furthermore, it appears to occur at even lower work rates. This may offer a new objective means of estimating exercise tolerance; however, its exact biological basis still needs to be elucidated.
Background
Ventilatory anaerobic threshold (VAT, or anaerobic threshold, AT) has been widely used as an index of exercise tolerance, primarily because it does not require maximal exercise [1, 2]. It is also recommended as an indicator of the optimal exercise training intensity during cardiac rehabilitation [3]. Among the methods for determining VAT, the Vslope method is considered to be the most basic; it directly assesses the relationship between VO_{2} and “excess CO_{2},” which is presumed to be derived from increased blood lactate levels [1]. It detects a breakpoint on the Vslope plotted on the x: VO_{2} versus y: VCO_{2} coordinates. The determination of the breakpoint (VAT), however, must be made visually, therefore making this parameter primarily a subjective measurement.
While using the Vslope method for determining VAT during routine cardiopulmonary incremental exercise tests (CPX) over a period of many years, we have found that the position of the Vslope itself is, from the early exercise stage, often shifted rightward to varying degrees from the reference diagonal line of the respiratory gas exchange ratio (R) of 1.0 in patients with cardiac disease as well as in normal subjects (Additional file 1: Figure S1). It also manifests itself as an initial drop in R. Since the 1960s, this phenomenon has been known to occur primarily in normal subjects during the first 1–2 min of steadystate exercise and disappears thereafter; it has been attributed to CO_{2} storage, presumably in active muscle [4–7]. However, we have noted that it also appears to occur during incremental exercise. We have also observed that the higher the VAT, the greater the rightward shift of the Vslope. We hypothesized that this rightward shift of the Vslope (RtShift) might be of clinical use as an index of exercise tolerance, if it could be quantified mathematically. This paper describes a method we have developed to mathematically derive RtShift and to elucidate whether this objective measure is in fact related to the level of VAT.
Methods
Patient characteristics
The CPX records of 100 patients with cardiac disease who underwent routine exercise testing and cardiac rehabilitation were retrospectively analyzed. There were 91 men and nine women, with a mean age of 63.8 ± 10.2 years. The underlying heart diseases were postacute myocardial infarction (n = 41), angina (n = 21), postcardiovascular surgery (n = 19), congestive heart failure (n = 14), and others. The characteristics of the study population are summarized in Table 1. New York Heart Association classification was not performed.
Data are presented as mean ± SD or number. BMI indicates body mass index; LVEF, left ventricular ejection fraction; LVDd, left ventricular diastolic dimension; Ca, calcium; ACE, angiotensinconverting enzyme; and ARB, Angiotensin II Receptor Blocker.
In a prospective substudy, the effect of different ramp exercise protocols on RtShift was assessed in 12 healthy young male students belonging to various college sports clubs; their mean age, body weight, and height were 20.8 ± 1.0 years, 66.0 ± 5.2 kg, and 172.6 ± 5.7 cm, respectively. This substudy was performed at a different institution (National Hospital Organization Hakodate Hospital, Hakodate, Japan).
The research plan was approved by the institutional review board of two institutions: Hokko Memorial Hospital and National Hospital Organization Hakodate Hospital. The study was conducted according to the Declaration of Helsinki.
Exercise test
Exercise tests were performed by using an upright bicycle with a breathbybreath gas analyzer (AE300S; Minato Ikagaku, Tokyo, Japan). The ramp protocol of 5–15 W⋅min^{−1} (10 W⋅min^{−1} in 84 cases, 15 W⋅min^{−1} in 13 cases, and 5 W⋅min^{−1} in three cases) was used, preceded by a 3min warmup. For the ramp protocol: 5 W⋅min^{−1}, 10 ⋅Wmin^{−1}, and 15 W⋅min^{−1}, the warmup load (W) was 0, 10 and 15, respectively; and the ramp start load (W) was 0, 10 and 15, respectively. The test was mostly terminated shortly after a VAT point was identified onscreen; therefore, the test protocol was submaximal, primarily performed for the purpose of identifying an initial exercise training intensity [3]: the work rate or heart rate at VAT. The Borg scale (ranging from 6 to 20) was used for the evaluation of perceived exertion. For the respiratory data analysis, a 10s average was used. The heart rate was monitored continuously and blood pressure was measured once every minute.
For the substudy to assess the effect of different ramp protocols, three ramps (15, 25, and 50 W⋅min^{−1}) were employed (symptomatic maximal exercise, on three different days 1 week apart).
Additionally, one normal volunteer performed an exercise protocol with three sixmin steady state steps; this was compared with a ramp protocol exercise to graphically demonstrate the effects of steady state exercise on CO_{2} storage in the framework of Vslope. In the past, the phenomenon has been repeatedly shown, but depicted graphically always with the elapsed time (s or min) on the xaxis, and VCO_{2} on the yaxis (1,7). Written informed consent was obtained before each exercise test.
Determination of VAT
VAT was primarily determined by using the Vslope method, as defined by Wasserman et al. [1]. However, the details of determination were based on those defined by Beaver et al. [8] and particularly by Sue et al. [9]. We drew the whole Vslope schematically as in Fig. 1. Two types of Vslope, types A and B, are shown. Type B is by far the most common type; the first segment of the slope starts from resting data points and moves both rightward and upward to varying degrees during the warmup and into the early ramp periods. This segment characteristically has an angle of <45°. We call this segment Str (slope transition). The Vslope then starts to ascend in parallel with the R = 1 line, at an angle of 45° (S1, preVAT slope). A break from S1 marks the appearance of VAT. Thereafter, the Vslope becomes steeper, with an angle of >45° (S2). In contrast, type A simply ascends with the R1 line before VAT. Further details of VAT determination and the Vslopes of all 100 cardiac subjects with VAT break points marked are included in Additional file 2: Figure S2. Overall, the VAT could not be determined in five cases.
The hypothesis that the S1 is parallel to the diagonal R = 1 was tested as described [10]. After all of the VAT points and S1 data points (those visually interpreted as being located on or around the hypothetical S1 line parallel to R = 1) had been finalized, the mean of the S1 slopes was tested against the slope of 1.0. Subsequently, the next data point on the Vslope was added and again tested against 1.0. The procedure was further repeated after sequentially adding the next data point.
The 95% limit of agreement for VAT determination was 170 mL⋅min^{−1} VO_{2} as per our experience [10]. This is based on the comparison of the averages of two readings of each assessor (i.e., A and B), therefore representing the interindividual agreement.
Rightward shift of the Vslope (RtShift)
Schematic representation of the concept of RtShift
The concept of RtShift is also schematically represented in Fig. 1. The S1 of Vslope A has no RtShift, being on the R = 1 line. The S1 of Vslope B is shifted rightward. The RtShift is defined as the horizontal distance between the two parallel lines B and A (b minus a), represented by the dotted line. However, the Vslope graph is conventionally drawn as a regular square with identical x and y scales; therefore, on the diagonal line, the xvalue is identical to the yvalue. The VO_{2} at “a” is equal to VCO_{2} at “a,” which then subsequently equals VCO_{2} at “b.” Therefore, the RtShift (b minus a) is simply calculated as a difference in mL⋅min^{−1} (VO_{2} minus VCO_{2}) at a single data point “b.” This diagram shows an interesting way in which RtShift affects a VAT value expressed as VO_{2} mL⋅min^{−1}; by merely shifting the whole S1 Vslope rightward, the VAT on B becomes greater than that on A. Therefore, it seems reasonable to assume that the level of VAT is variably augmented by the presence of RtShift.
Mathematical derivation of RtShift
Before the mathematical derivation of RtShift, we processed the data as follows. We employed only data points that are equal to or below R = 1 throughout the resting and exercise periods (warmup and ramp). Outlier data points were eliminated from further analysis by using the SmirnovGrubbs test; 24 data points were removed as a result. Next, we converted the conventional Vslope plot to the plot wherein RtShift is plotted against VO_{2} as in Fig. 2. On identifying a break point on the conventional Vslope, one is forced to detect a rising break point on an already rising S1 baseline. With the converted coordinate system, we have a normalized x–yaxis view of the relation between VO_{2} and the excess VCO_{2} (RtShift). A typical plot yielded a trapezoidal form; however, with a shorter S1, the plot became more like a triangle with a blunted apex.
We attempted to quantify RtShift by using following three methods.
The first method determined RtShift through the visual inspection of a Vslope. For this method we used the VO_{2} vs. RtShift graph. Additionally we also used the time (s) vs. RtShift graph, because this graph could avoid the overlapping of RtShift data points over a different time sequence, which is sometimes the problem on the VO_{2} vs. RtShift graph. The RtShift was visually taken as the highest parallel line to X = 0 (therefore, R = 1); occasionally, more than one horizontal line could be visualized before VAT. The second method involved the use of a curvefitting program. As the plot pattern on the converted Vslope graph mostly resembled a trapezoid or blunted triangle, we used quadratic regression to represent the curve (Fig. 2). Differentiation of the quadratic equation (ax^{2} + bx + c) yields (ax + b); solving (ax + b = 0) for x gives the tangent to the apex of the curve, RtShift. In 11 cases, the quadratic equation could not yield a regression curve with the convex within the data range (i.e., failure rate of 11%), because of a small data set and/or a large dispersion of data points, in which cases, a mathematical average of all data points equal to or below R = 1 was used instead. The third method was simply a mathematical average of all data points at or below R = 1, assuming no particular type of curve. This is equivalent to the integral of all RtShift values divided by the number of data points. In Additional file 2: Figure S2, each Vslope is drawn with the diagonal R = 1 line and a marked VAT point, so that the relation between RtShift and VAT may be individually ascertained.
Because RtShift defined this way included a shift at rest (calculated as an average resting VO_{2} minus an average resting CO_{2}), we may correct for this by subtracting this resting RtShift from the exercise RtShift. However, we decided that this would be unnecessary because we wanted to compare RtShift against VAT, which itself included a resting value; a resting VO_{2} is not subtracted from a VAT value.
The ΔVO_{2}/Δwork rate was calculated in the standard way [10]. In two cases, the ratio could not be calculated because of the small number and large dispersion of the data points.
Left ventricular function
Left ventricular dimension and function were assessed by echocardiography (Toshiba Aplio ™400, Toshiba, Tokyo, Japan). Left ventricular diastolic dimension (LVDd) was measured on the Mmode scan. Ejection fraction was calculated with the modified Simpson method.
Statistical analysis
Data are presented as mean ± SD. The distribution of RtShift was not Normal; therefore, the data for RtShift are presented both as mean ± SD and as median with interquartile range (IQR) and lower/upper quartile. Outlier detection was performed with the SmirnovGrubbs test. Correlation was assessed by using Pearson’s r. Comparison between two groups was performed with Student’s ttest. Comparisons between RtShift calculated through three different methods, and between RtShift in the three ramp protocols, were performed by using oneway repeatedmeasures analysis of variance (ANOVA) followed by the Bonferroni multiple comparison procedure. Multiple linear regression analysis was performed to assess whether VAT, endtidal CO_{2} concentration (as a sign of hyperventilation), resting RtShift, the use of betablocker, and factors such as age and weight independently contributed to RtShift. Because of the nonNormal nature of the data, we also applied nonparametric tests: namely, Spearman’s rank order correlation for correlation and the Friedman test as an alternative to repeatedmeasures ANOVA).
Results
The basic exercise data are summarized in Table 2.
Data are presented as mean ± SD. SBP indicates systolic blood pressure; VO_{2}, indicates minute oxygen uptake; RR, respiration rate; RER, respiratory exchange rate; VAT, ventilatory anaerobic threshold; VCO_{2}, minute carbon dioxide production; RtShift, rightward shift of Vslope; IQR, interquartile range.
First, typical examples of a case with very little rightward shift (type A) and a moderate shift (type B) of the Vslope are shown in Fig. 3(a, b). Next, the maintenance of RtShift during a ramp protocol (Fig. 4(a)) is contrasted with its lack during a steadystate protocol (Fig. 4(b)) consisting of a Vslope during three sixmin steps. The rightward shift is noted to occur only at the start of the step increase, ultimately reverting back later.
Among the visual method and the two mathematical methods for estimating RtShift, the first method (visual) yielded the highest mean RtShift, followed by the second (quadratic curve fitting) and then the third (simple averaging) method. The mean RtShift values were as follows: first, 50.8 ± 41.6; second, 42.0 ± 34.4; third, 33.9 ± 25.0 mL⋅min^{−1} VO_{2} (each different from the others, p < 0.001). However, the RtShift of each method correlated well with each other (first vs. second, r = 0.955; second vs. third, r = 0.989; first vs. third, r = 0.945). We decided to adopt the third method (simple averaging) to determine RtShift because there was no failure associated with the determination; however, with the second method, there was an 11% failure rate in determining RtShift. The third method also utilized all data points.
As mentioned in the Statistical analysis section, RtShift did not have a Normal distribution (Additional file 3: Figure S3(a)).
VAT (mL⋅min^{−1}) and RtShift (mL⋅min^{−1} VO_{2}) were significantly correlated (r = 0.718, p < 0.001; adjusted for weight, r = 0.543, p < 0.001); the greater the VAT, the greater the rightward shift (Fig. 5). VAT was also significantly correlated with RtShift after correction for the resting RtShift (r = 0.402, p < 0.001).
VAT, resting RtShift, and resting endtidal CO_{2} concentration were each significant (p < 0.001, p < 0.001, and p = 0.003; standardized correlation coefficient: 0.5699, 0.3966, and −0.2032, respectively). The univariate correlation between RtShift and resting endtidal CO_{2} concentration was, however, only 0.094 (p = 0.351). The use of a βblocker was not a significant variable in this analysis (p = 0.612).
Both VAT and RtShift were significantly correlated with ΔVO_{2}/Δwork rate, a marker of oxygen uptake efficiency (10 ⋅Wmin^{−1}, n = 83; r = 0.531, p < 0.001; adjusted for weight, r = 0.393, p < 0.001, respectively) (Additional file 3: Figure S3(b)).
The result of the hypothesis testing that S1 is parallel to the diagonal R = 1 was as follows: the first test of the mean S1 slope (1.013 ± 0.084) against the slope of 1.0 was found to be not significantly different from 1.0 (p = 0.132). Subsequently, the next data point on the Vslope was added and again tested against 1.0; the mean slope was 1.056 ± 0.090, significantly greater than 1.0 (p < 0.001). As the procedure was further repeated after sequentially adding the next data point, the mean S1 slope increased progressively. Furthermore, the individual linear regression line through the S1 data points varied depending on the number of data points included; a long S1 consisting of more data points appeared to converge to 1.0, whereas a short S1 with a few data points tended to diverge. However, the mean of all S1 linear regression lines (95 cases in total) was not significantly different from 1.0 (Additional file 3: Figure S3(c)).
The effect of the three different protocols (15, 25, and 50 W⋅min^{−1}) on RtShift (mL⋅min^{−1} VO_{2}) was not statistically significant (15 W⋅min^{−1}⋅, 132.2 ± 52.5; 25 W⋅min^{−1}, 137.8 ± 48.5; 50 W⋅min^{−1}, 136.1 ± 49.3; p = 0.946).
Nonparametric statistical tests yielded concordant results with the parametric tests. The detailed results are presented in the additional file (Additional file 4).
Discussion
During steadystate exercise, CO_{2} storage in the tissue occurs at the start and disappears after 1–2 min, when presumably the storage space has been filled. This has been observed since the 1960s [4–7]. In addition, the literature acknowledges that, even during incremental or ramp exercise testing, there is an artifact due to this phenomenon at the beginning of exercise and this segment of Vslope is to be excluded when VAT is determined [8, 11]. What has not been clearly pointed out is that during ramp exercise, this CO_{2} tissue storage effect never disappears and persists at least until VAT begins (Fig. 4). Because of its nature, ramp exercise never achieves a steady state. From the observed persistence of CO_{2} storage, it may be reasoned that the CO_{2} storage space also increases and its rate of increase remains constant as the work rate increases. We hypothesize that, with the constantly increasing work rate, new groups of muscle fibers are recruited, resulting in a steady increase in new CO_{2} storage space in active muscle.
CO_{2} storage was initially measured by using hyperventilation and/or rebreathing methods [4–6]. Later, CO_{2} storage was calculated as the difference between the measured CO_{2} and the predicted CO_{2}, based on the premise of a fixed respiratory quotient during sublactatethreshold exercise [12–14]. Estimates from the hyperventilation methods, with the apparent inclusion of an HCO_{3} ^{−} diffusing space in the whole body, were much larger than those from the rebreathing methods. The former measure of CO_{2} storage was presumed to result only from tissue with a high metabolic rate, such as active muscle [14]. There is a possibility that the phenomenon of RtShift may be chiefly metabolic in origin, such as a shift to fat utilization as energy source, instead of due to CO_{2} storage. However, we believe that this is highly unlikely because the RtShift occurs very early at the outset of exercise, whereas a major shift in the energy source to fat takes much longer to manifest, such as 20 min of steadystate exercise [15]. With ramp exercise too, there is no reason to believe that the energy source change to fat occurs so early, although no such basic study has ever been done during ramp protocols. Moreover and most important, the RtShift disappears (Vslope turns leftward) in 2–3 min while the work rate is being maintained at the same level (steady state). We assumed, on the basis of visual observations, that the slope of S1 of Vslope is approximately 1.0. However, the literature suggests that the RQ (R during steadystate exercise) increases by 0.05–0.1 from the resting value to moderate work rates [16]. This magnitude of RQ change is too small compared with R of either RtShift (far below 1.0) or S2 (far above 1.0). Therefore, we do not believe that the hypothesis of the slope of S1 being exactly 1.0 is an absolute requisite in the practical application of VAT and RtShift determination. Although there are no data on the energy utilization during ramp or rapid incremental testing, during work rate transition from rest to exercise, part of the energy fuel is estimated to come from phosphocreatine breakdown and glycolysis [17], both producing H^{+} and probably excess CO_{2}. Ramp exercise may be considered to consist of multiple work rate transitions; therefore, R may tend to increase slightly even during increments of light exercise intensity. Therefore, the observed slope of S1 of 1.0 (or its being parallel to R = 1 line) may be a result of a competing net effect of these factors (increasing slope of S1) and the RtShift (decreasing slope of S1). Wasserman states that the slope of S1 is slightly less than 1.0 [1].
The major finding of this research—that the size of the CO_{2} storage space (RtShift) is associated with increased exercise tolerance—may, at first, seem surprising. However, the above hypothesis suggests that a greater mass of aerobic muscle fibers recruited for a given workload may be a contributing mechanism. Additionally, a greater level of carbonic anhydrase (CA) activity in muscle may be a factor. As it is now known that there are subtypes of CA expressed in muscle [18], the hydration of CO_{2}, and therefore the production and retention of HCO_{3} ^{−} in muscle, may be facilitated by this mechanism. CA activity may increase as part of the overall increase in aerobic enzyme activity in muscle [19]. The only other study evaluating the relation between exercise tolerance and CO_{2} storage is that of Chuang et al. [20], in which the authors demonstrated a significant positive correlation between O_{2} deficit and CO_{2} store in 12 relatively young healthy subjects. However, in their study, the direct relation between indices of exercise tolerance (VO_{2}max or VAT) and the size of the CO_{2} store was not assessed.
The RtShift has several distinctive features that may be useful in the evaluation of exercise tolerance. Foremost among them is that RtShift does not require any assessor to determine its value. Second, its calculation is simple. Third, it has the additional advantage of an almost 100% determination rate. Fourth, RtShift, like VAT, does not require maximal exercise testing; it is mostly calculated from data points up to R = 1; therefore, it is unique in that it is primarily an aerobic phenomenon. Although it is significantly correlated with VO_{2} peak and VAT, the correlation coefficient is only moderate (~0.7). Therefore, although it may not be used as a substitute for VAT, it may indicate that a yetunidentified physiological factor that is absent in both VO_{2} peak and VAT is involved. That these two parameters do not always move in parallel is obvious from the only moderate correlation efficient between them and the inspection of each case of Vslope (Additional file 2: Figure S2).
Multiple regression analysis was used to assess whether factors other than VAT accounted for the association between RtShift and VAT. It has been reported that hyperventilation before exercise caused a rightward shift of the Vslope [21, 22]. However, that was not a major factor in this study.
As shown in this study, RtShift may be quantitatively defined by more than one method. In this study, we chose a method of simple averaging. However, any other method may be devised as long as it can provide a valid estimate of RtShift during incremental exercise testing.
The substudy was performed because we were not sure how a ramp change of exercise protocol affected the size of RtShift. If the RtShift is dependent on the recruited muscle mass, then the steeper ramp may result in a greater RtShift. As shown in the Results section, ramp change did not result in a significant change in RtShift.
The RtShift data distribution was clearly not Normal, with a deviation to the lower values of RtShift. Log transformation of the data was performed; however, the statistical results of various analyses did not materially differ from those obtained with raw data. We also ran nonparametric tests with the concordant results with parametric tests. Therefore, we chose to present only nontransformed data for analysis, which has the benefit of visually corresponding to the shift on the Vslope graph. In addition, we believe that because of the relatively large size of our data set (n = 100), the central limit theorem applies [23].
A limitation of this study is that the exact physiological mechanism through which RtShift is related to exercise tolerance is, at present, only speculative. The exact biological basis of RtShift must be elucidated through basic science research. Another limitation is that we do not have an age and sexmatched control population of sufficient size comparing RtShift between patients with cardiac disease and normal subjects. This issue must be systematically and prospectively addressed; in this study, we merely completed the process of quantifying RtShift.
Conclusions
In conclusion, during incremental CPX, there were varying degrees of rightward shift of the Vslope and they correlated significantly with VAT, suggesting that RtShift may be used as a completely objective measure of exercise tolerance. It is predominantly determined before the appearance of VAT. Its physiological meaning and clinical application requires further clarification.
Abbreviations
 ACE:

Angiotensinconverting enzyme
 ARB:

Angiotensin II receptor blocker
 BMI:

Body mass index
 Ca:

Calcium
 CA:

Carbonic anhydrase
 CPX:

Cardiopulmonary exercise testing
 IQR:

Interquartile range
 LVDd:

Left ventricular diastolic dimension
 LVEF:

Left ventricular ejection fraction
 R:

Respiratory exchange rate
 RQ:

Respiratory quotient
 RR:

Respiration rate
 RtShift:

Rightward shift of Vslope
 S1:

preVAT slope of the Vslope or the S1 itself (if the slope is not mentioned)
 S2:

postVAT slope of the Vslope or the S2 itself (if the slope is not mentioned)
 SBP:

Systolic blood pressure
 VAT:

Ventilatory anaerobic threshold
 VCO_{2} :

Minute carbon dioxide production
 VO_{2} :

Minute oxygen uptake.
References
 1.
Wasserman K, Hansen JE, Sue DY, Stringer WW, Sietsema KE, Sun XG, et al. Principles of exercise testing and interpretation: including pathophysiology and clinical application. 5th ed. Philadelphia: Lippincott Williams & Wilkins; 2012. p. 30–4. 83–6.
 2.
Matsumura N, Nishijima H, Kojima S, Hashimoto F, Minami M, Yasuda H. Determination of anaerobic threshold for assessment of functional state in patients with chronic heart failure. Circulation. 1983;68:360–7.
 3.
JCS Joint Working Group.Guidelines for rehabilitation in patients with cardiac disease (JCS 2012). Circ J. 2014;78:2022–93.
 4.
Clode M, Clark TJ, Campbell EJ. The immediate CO_{2} storage capacity of the body during exercise. Clin Sci. 1967;32:161–5.
 5.
Fowle AS, Campbell EJ. The immediate carbon dioxide storage capacity of man. Clin Sci. 1964;27:41–9.
 6.
Jones NL, Jurkowski JE. Body carbon dioxide storage capacity in exercise. J Appl Physiol Respir Environ Exerc Physiol. 1979;46:811–5.
 7.
Yano T. Immediate CO_{2} storage capacity at the onset of exercise. Jpn J Physiol. 1986;36:1241–52.
 8.
Beaver WL, Wasserman K, Whipp BJ. A new method for detecting anaerobic threshold by gas exchange. J Appl Physiol. 1986;60:2020–7.
 9.
Sue DY, Wasserman K, Moricca RB, Casaburi R. Metabolic acidosis during exercise in patients with chronic obstructive pulmonary disease: use of the Vslope method for anaerobic threshold determination. Chest. 1988;94:931–8.
 10.
Kominami K, Nishijima H, Imahashi K, Katsuragawa T, Murakami M, Yonezawa K, et al. Very early lactate threshold in healthy young men as related to oxygen uptake kinetics. Medicine (Baltimore). 2015;94:e1559. doi:10.1097/MD.0000000000001559.
 11.
Cooper CB, Beaver WL, Cooper DM, Wasserman K. Factors affecting the components of the alveolar CO_{2} outputO_{2} uptake relationship during incremental exercise in man. Exp Physiol. 1992;77:51–64.
 12.
Hughson RL, Inman MD. Gas exchange analysis of immediate CO_{2} storage at onset of exercise. Respir Physiol. 1985;59:265–78.
 13.
Yano T. The differences in CO_{2} kinetics during incremental exercise among sprinters, middle, and long distance runners. Jpn J Physiol. 1987;37:369–78.
 14.
Barstow TJ, Landaw EM, Springer C, Cooper DM. Increase in bicarbonate stores with exercise. Respir Physiol. 1992;87:231–42.
 15.
Romijn JA, Coyle EF, Sidossis LS, Gastaldelli A, Horowitz JF, Endert E, et al. Regulation of endogenous fat and carbohydrate metabolism in relation to exercise intensity and duration. Am J Physiol. 1993;265(3 Pt 1):E380–91.
 16.
Goedecke JH, St Clair Gibson A, Grobler L, Collins M, Noakes TD, Lambert EV. Determinants of the variability in respiratory exchange ratio at rest and during exercise in trained athletes. Am J Physiol Endocrinol Metab. 2000;279(6):E1325–34.
 17.
Baker JS, McCormick MC, Robergs RA. Interaction among skeletal muscle metabolic energy systems during intense exercise. J Nutr Metab. 2010;2010:905612. doi:10.1155/2010/905612. Epub 2010 Dec 6.
 18.
Geers C, Gros G. Carbon dioxide transport and carbonic anhydrase in blood and muscle. Physiol Rev. 2000;80:681–715.
 19.
Zoll J, Ponsot E, Dufour S, Doutreleau S, VenturaClapier R, Vogt M, et al. Exercise training in normobaric hypoxia in endurance runners. III. Muscular adjustments of selected gene transcripts. J Appl Physiol. 2006;100:1258–66. doi:10.1152/japplphysiol.00359.2005.
 20.
Chuang ML, Ting H, Otsuka T, Sun XG, Chiu FY, Beaver WL, et al. Aerobically generated CO(2) stored during early exercise. J Appl Physiol. 1999;87:1048–58.
 21.
Ozcelik O, Ward SA, Whipp BJ. Effect of altered body CO_{2} stores on pulmonary gas exchange dynamics during incremental exercise in humans. Exp Physiol. 1999;84:999–1011.
 22.
Ward SA, Whipp BJ, Koyal S, Wasserman K. Influence of body CO_{2} stores on ventilatory dynamics during exercise. J Appl Physiol Respir Environ Exerc Physiol. 1983;55:742–9.
 23.
Altman DG. Practical statistics for medical research. London: Chapman & Hall/CRC; 1991. p. 154.
Acknowledgments
We offer our special thanks to Sapporo Medical University Scholarly Communication Center for acquiring the needed literature.
Funding
None.
Availability of data and materials
The data set supporting the conclusions is provided in Additional file 5 (Table S1).
Authors’ contributions
Conceptualization and methodology: HN, KK, KY, MS. Formal analysis: HN, KK, HH. Investigation: HN, KK, KY, HH. Writing the paper: HN, KK, KY, HH, MS. Supervision: HN, KK, MS. All authors critically revised the manuscript. All gave final approval and agreed to be accountable for all aspects of this work ensuring integrity and accuracy.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
This research was approved by the institutional review board of two institutions: Hokko Memorial Hospital and National Hospital Organization Hakodate Hospital. The informed consent for the 100 patients who underwent cardiopulmonary exercise testing was judged not necessary by the institutional review board because the study was observational based on the data collected in the past. For the substudy of 12 normal subjects, a written informed consent was obtained from each participant.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Additional files
Additional file 1: Figure S1.
Our past experience on the relation between ventilatory anaerobic threshold (VAT) and RtShift. (PPTX 138 kb)
Additional file 2: Figure S2.
Individual Vslope plots (n = 100). (PPTX 3053 kb)
Additional file 3: Figure S3.
(a) Histogram of RtShift data distribution. (b) RtShift versus ΔVO_{2}/Δwork rate. (c) Test of S1 being parallel to R = 1. (PPTX 108 kb)
Additional file 4:
Result: Results of nonparametric statistical analyses. (DOCX 13 kb)
Additional file 5: Table S1.
Complete data set. (XLSX 26 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Exercise tolerance
 Ventilatory anaerobic threshold
 CO_{2} storage