Oxygenation Threshold Derived from Near-Infrared Spectroscopy: Reliability and Its Relationship with the First Ventilatory Threshold

Abstract

Background

Near-infrared spectroscopy (NIRS) measurements of oxygenation reflect O₂ delivery and utilization in exercising muscle and may improve detection of a critical exercise threshold.

Purpose

First, to detect an oxygenation breakpoint (Δ[O₂HbMb-HHbMb]-BP) and compare this breakpoint to ventilatory thresholds during a maximal incremental test across sexes and training status. Second, to assess reproducibility of NIRS signals and exercise thresholds and investigate confounding effects of adipose tissue thickness on NIRS measurements.

Methods

Forty subjects (10 trained male cyclists, 10 trained female cyclists, 11 endurance trained males and 9 recreationally trained males) performed maximal incremental cycling exercise to determine Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds (VT1 and VT2). Muscle haemoglobin and myoglobin O₂ oxygenation ([HHbMb], [O₂HbMb], SmO₂) was determined in m. vastus lateralis. Δ[O₂HbMb-HHbMb]-BP was determined by double linear regression. Trained cyclists performed the maximal incremental test twice to assess reproducibility. Adipose tissue thickness (ATT) was determined by skinfold measurements.

Results

Δ[O₂HbMb-HHbMb]-BP was not different from VT1, but only moderately related (r = 0.58–0.63, p<0.001). VT1 was different across sexes and training status, whereas Δ[O₂HbMb-HHbMb]-BP differed only across sexes. Reproducibility was high for SmO₂ (ICC = 0.69–0.97), Δ[O₂HbMb-HHbMb]-BP (ICC = 0.80–0.88) and ventilatory thresholds (ICC = 0.96–0.99). SmO₂ at peak exercise and at occlusion were strongly related to adipose tissue thickness (r² = 0.81, p<0.001; r² = 0.79, p<0.001). Moreover, ATT was related to asymmetric changes in Δ[HHbMb] and Δ[O₂HbMb] during incremental exercise (r = -0.64, p<0.001) and during occlusion (r = -0.50, p<0.05).

Conclusion

Although the oxygenation threshold is reproducible and potentially a suitable exercise threshold, VT1 discriminates better across sexes and training status during maximal stepwise incremental exercise. Continuous-wave NIRS measurements are reproducible, but strongly affected by adipose tissue thickness.

About the Authors:

Stephan van der Zwaard

* E-mail: [email protected]

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Richard T. Jaspers

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Ilse J. Blokland

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Chantal Achterberg

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Jurrian M. Visser

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Anne R. den Uil

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Mathijs J. Hofmijster

Affiliations Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands, Faculty of Sports and Nutrition, Amsterdam University of Applied Sciences, Amsterdam, the Netherlands

Koen Levels

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Dionne A. Noordhof

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Arnold de Haan

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Jos J. de Koning

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Willem J. van der Laarse

Affiliation: Department of Physiology, VU Medical Center, Amsterdam, the Netherlands

Cornelis J. de Ruiter

Affiliation: Department of Human Movement Sciences, Vrije Universiteit Amsterdam, MOVE Research Institute Amsterdam, the Netherlands

Introduction

Performance benefits from adequate training. Even though training intensity can be quantified in many different ways, endurance training intensity is often quantified by the lactate thresholds (LT1 and LT2) obtained from blood sampling or the ventilatory thresholds (VT1 and VT2) obtained from gas exchange data [1]. Above LT1 and VT1, blood lactate concentration is elevated [1] and CO₂ production exceeds oxygen uptake due to anaerobic energy production in the muscle cells as bicarbonate buffers the produced lactic acid [2]. Above LT2 and VT2, blood lactate accumulates as lactate production exceeds removal rates [1] and respiratory compensation starts by an excessive increase in minute ventilation [3]. There is consensus that one can obtain ventilatory or lactate thresholds reproducibly from an incremental exercise test [1,4,5] and can use these thresholds for training purposes [1]. However, ventilatory and lactate thresholds are derived from measurements at the mouth or (closer to the muscle) in the blood. Therefore, the first ventilatory and lactate threshold rather indirectly and with a delay reflect an increase of anaerobic ATP resynthesis. Both exercise intensity and oxygen availability contribute to this transition. Increases in the former lead to recruitment of additional (larger) motor units with lower oxidative capacity, while the latter depends on O₂ consumption by the contracting muscle fibres and O₂ supply by the O₂ cascade from the air to the mitochondria [6]. Near-infrared spectroscopy (NIRS) provides non-invasive measures of the balance between oxygen delivery and consumption in the muscle [7], and thereby NIRS is potentially a very suitable technique to detect a critical exercise threshold directly in the exercising muscle.

In active muscle, the matching of oxygen supply () and oxygen consumption () is crucial for O₂ diffusion and metabolic control [8]. This matching can theoretically be assessed non-invasively by NIRS [8]. Hereto, concentration changes of deoxygenated haemoglobin and myoglobin ([HHbMb]), oxygenated haemoglobin and myoglobin ([O₂HbMb]) and tissue saturation (SmO₂) are measured. These NIRS signals present multiple breakpoints during incremental exercise. For instance, the [HHbMb] breakpoint during ramp exercise is suggested to resemble an upper limit of fractional oxygen extraction in the exercising leg muscles during incremental exercise [9–11] and approximates VT2 [10,12,11,13–15], maximal lactate steady state [14,16] and critical power [14]. At somewhat lower exercise intensities, a breakpoint in Hb difference ([O₂HbMb-HHbMb]) occurs during maximal stepwise incremental exercise and shows high correlations with VT1 and the lactate threshold (r>0.88) [17–19]. At this [O₂HbMb-HHbMb] breakpoint, the fractional oxygen extraction supposedly has not yet reached its upper limit in the exercising leg muscles during incremental exercise. Therefore, this breakpoint could reflect a transition in oxygen extraction that is likely attributed to capillary-venular PO₂ reaching its critical value (PO_2crit) and lactic acidosis inducing a rightward shift of the O₂Hb dissociation curve (i.e. Bohr effect) [17]. Hence, the [O₂HbMb-HHbMb] breakpoint may relate to the occurrence of oxygen limitation in the skeletal muscle fibres.

Although reproducibility of ventilatory thresholds [4] and lactate threshold [5] are known to be high, data on reproducibility of NIRS signals [5,20–22] and NIRS derived breakpoints [23] is limited. Only in a single study with sedentary subjects, reproducibility of breakpoints in [O₂HbMb] has been determined (r = 0.67–0.85) [23]. However, it is yet unknown whether reproducibility of these [O₂HbMb] breakpoints also applies to trained subjects and to other NIRS derived breakpoints.

At present, the [O₂HbMb-HHbMb] breakpoint has been determined during cycling exercise in mountain climbers [17], fin swimmers [18] and college students [19]. However, the [O₂HbMb-HHbMb] breakpoint has not been assessed with subjects differing in training status and sex in one study. It is important to note that adipose tissue thickness (ATT) may vary widely between male and female subjects and across subjects with heterogeneous training status. Also, it is well-known that scattering and absorbance due to adipose tissue affect NIRS signals, reducing absorbance by underlying muscle tissue [24–26]. One may reduce these effects of ATT by normalizing oxygenation changes to the oxygenation at peak exercise [10,27], maximal voluntary contraction [26,28] or cuff occlusion [21,26], where the effect of ATT on the amplitude of separate [O₂HbMb] and [HHbMb] signals can be assessed in the absence of blood volume changes.

The main purpose of this study was to investigate how the Δ[O₂HbMb-HHbMb] breakpoint relates to ventilatory thresholds during incremental step exercise across sexes and training status. Secondary purposes of the study were to determine the reproducibility of NIRS signals and exercise thresholds and to assess the confounding effect of adipose tissue thickness on NIRS signals. We expected that exercise intensity is similar at the [O₂HbMb-HHbMb] breakpoint and VT1. In addition, we hypothesized that both the [O₂HbMb-HHbMb] breakpoint and VT1 differ across sexes and training status due to differences in oxygen supply [29–32] and oxygen utilization [29,30] within the exercising muscles.

Methods

Subjects

Forty subjects participated in the present study: 10 trained male cyclists (CM), 10 trained female cyclists (CF), 11 endurance trained males (EM) and 9 recreationally trained males (RM) as classified by their [33] (characteristics are summarized in Table 1). Female and male cyclists (CF, CM) were recruited from Dutch cycling teams, and RM and EM were recruited based on training volume (EM >5 hours endurance training per week and RM <3 hours training per week). Note that EM contained 7 cyclists and 4 non-cyclists, since many endurance athletes perform cycling exercise in the Netherlands. Reliability measurements were performed in CM and CF. Prior to participation, experimental procedures, risks and aims of the study were explained and all subjects provided written informed consent. The study was conducted according to the principles of the Declaration of Helsinki and was approved by the ethics committee of the Department of Human Movement Sciences, Vrije Universiteit, Amsterdam, the Netherlands.

[Figure omitted. See PDF.]

Table 1. Subject characteristics.

https://doi.org/10.1371/journal.pone.0162914.t001

Experimental design

Subjects performed a maximal incremental step exercise test to exhaustion on a friction-braked cycle ergometer (Monark Ergomedic 839E, Monark exercise AB, Sweden), pedalling at 90 rpm. After three minutes of seated rest, subjects started at a workload of 1.5W·kg^-1 (85–145W), which was increased with 0.5W·kg^-1 every three minutes (30–50W·3min^-1) until voluntary exhaustion. Trained cyclists performed two maximal incremental tests on separate days to assess reproducibility of NIRS signals and exercise thresholds. The test was terminated if cadence dropped below 80 rpm, despite verbal encouragement. Active cool-down consisted of five minutes cycling at 1 W·kg^-1 (60–100W). Thereafter, an arterial occlusion was performed to assess maximal desaturation of the m. vastus lateralis. Hereto, a pressure cuff was inflated to approximately ~300mmHg at the most proximal part of the thigh until maximal desaturation was reached (duration 314 ± 44 s). The cuff was secured by an extra strap that prevented it from unwrapping [34].

Prior to the maximal incremental test, subjects were instructed to avoid strenuous exercise and alcohol consumption within 24 hours before the test and to consume their last meal and caffeinated beverages at least three hours prior to testing. Subsequent measurements were separated by at least one day of rest. Saddle and handle bar position were adjusted to individual preferences and were recorded for subsequent measurements. Environmental circumstances were controlled in a climate-controlled environment (temperature 16.3±0.9°C, relative humidity 46.4±6.4%).

Data collection

Pulmonary measures.

Ventilation and gas exchange were analysed breath-by-breath using open circuit spirometry (Cosmed Quark CPET, Cosmed S. R. L., Rome, Italy). Prior to every test, the gas analyser and volume transducer were calibrated according to the manufacturer’s instructions.

Near-Infrared Spectroscopy.

Relative changes with respect to the start of maximal incremental test were measured for haemoglobin + myoglobin ([HHbMb], [O₂HbMb]) and muscle saturation (SmO₂) by a continuous-wavelength portable NIRS device using spatially resolved spectroscopy (Portamon, Artinis Medical Systems, Arnhem, the Netherlands). Concentration changes of the right leg were determined at 10Hz at the vastus lateralis (VL) muscle belly and calculated from light absorbance at 758nm and 847nm using the modified Lambert-Beer law. The NIRS device was positioned above the upper patella border at 1/3 distance of the patella and greater trochanter, parallel to the longitudinal femur axis. Light emitting optodes were positioned at 30, 35 and 40mm distal to the receiver, sending near-infrared light into the thigh with a corresponding penetration depth of 15–20mm [24]. How much light the muscle tissue can absorb depends on light absorbance and scattering by skin and subcutaneous adipose tissue (e.g. influenced by adipose tissue thickness, temperature, and pigment) and on secured contact of the NIRS device with the skin. Motion artefacts were minimized by fixating the optode position using adhesive tape (Fixomull®) and elastic bandages. Additionally, placement of a light absorbing cloth over the NIRS device minimized the influence from ambient light. Finally, probe position was marked for replication in subsequent measurements.

Skinfold.

Skinfold thickness of the m. vastus lateralis at the site of the NIRS device was measured (median of three measurements) in seated position using Harpenden skinfold callipers (British Indicators Ltd, Burgess Hill, UK). Adipose tissue thickness was calculated by dividing skinfold thickness by two [25], resembling subcutaneous fat and skin. Subjects showed a wide range in adipose tissue thickness (2.0–11.7 mm), yet participants were excluded if adipose tissue thickness was above 12.5 mm.

Data analysis

Gas exchange, power output (PO) and NIRS data were converted to 1Hz samples using linear interpolation, filtered for extreme values (data points deviating > 2 SD from a local window of 6 data points were replaced by the local average) and time-assembled. Critical exercise thresholds and maximal physiological parameters were derived from 1Hz data using smoothed 30-second moving averages and peak PO was defined as average PO over the last 180 seconds prior to termination of the test. VT1 was derived from respiratory exchange measurements using the V-slope method and ventilatory equivalents [35] and VT2 was derived using the minute ventilation versus graph and ventilatory equivalents [3]. Both thresholds were detected by automatic regression and subsequently (re)evaluated by two independent observers. In addition, the difference between concentration changes in [HHbMb] and [O₂HbMb] were used to derive an oxygenation breakpoint (Δ[O₂HbMb-HHbMb]-BP) using a double linear regression method [17–19]. Both linear regressions were fitted through Δ[O₂HbMb-HHbMb] versus relative PO (%PO) plots. Δ[O₂HbMb-HHbMb]-BP was defined as the intercept of two congregating regression lines with the least combined residuals sum of squares (Eq 1). Warm-up stages were excluded and regression lengths were set at ≥10% of relative workload.

(1)

In addition, the confounding effect of ATT on NIRS signals was determined for SmO₂ at the end of maximal exercise and arterial occlusion. The amplitudes of [HHbMb] and [O₂HbMb] were assessed during rest and at the end of maximal exercise and arterial occlusion to determine the symmetry of these NIRS signals in relation to ATT. Moreover, next to the SmO₂ values provided by the NIRS apparatus, SmO₂ was also normalized to maximal changes during exercise (Eq 2). All data processing procedures were performed using custom written software in Matlab (Mathworks Co, Natick, Massachusetts, USA).

(2)

Statistics

All data are presented as individual values or as mean±SD, unless indicated otherwise. Differences in and PO between Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds were assessed by paired samples t-test and Pearson’s product-moment correlation analyses. A priori sample size calculations have been performed in GPower 3.1.2 (University Kiel, Germany) using the t-tests means (difference between two dependent means, matched pairs, a priori) with significance level 0.05 and power 80%. These calculations showed that differences in the Δ[O₂HbMb-HHbMb]-BP and VT1 could be detected with a sample size of 9 subjects for PO (W) and 4–8 subjects for (L/min) [18,19]. Differences between groups were assessed by one-way ANOVAs (between factor group: CM-CF-EM-RM). Moreover, reproducibility was determined by intra class-correlations (i.e. single measures ICC_3,1) and within-subject coefficient of variation (CV) after logarithmic transformation of the data. Qualification of correlation coefficients was performed according to Evans [36]. Relationships between ATT and SmO₂ were assessed by hyperbolic regressions and relationships between ATT and the symmetry of changes in [HHbMb] and [O₂HbMb] amplitude were assessed by linear regressions. Presented R-squared values were adjusted for the number of predictors (i.e. coefficients) in the regression model, being a more conservative measure of explained variance. Differences were considered to be significant if p < 0.05.

Results

Critical exercise threshold

Δ[O₂HbMb-HHbMb]-BP was determined with high regression model quality (r² = 0.97 ± 0.02, Fig 1). Including all subjects, Δ[O₂HbMb-HHbMb]-BP was not significantly different from VT1 (: -0.155±0.509 L/min; PO: -12.1±44.2 W), but significantly lower than VT2 (: -0.624±0.521 L/min; PO: -56.8±44.2W, p<0.001, Table 2). Subgroup analyses revealed that Δ[O₂HbMb-HHbMb]-BP was not significantly different from VT1 in CM, CF and RM, however, in EM the Δ[O₂HbMb-HHbMb]-BP was significantly smaller than VT1 (p<0.05). Δ[O₂HbMb-HHbMb]-BP was only moderately related to VT1 (: r = 0.63; PO: r = 0.58; p<0.001; Fig 2) and VT2 (: r = 0.68; PO: r = 0.62; p<0.001). For subgroups, Δ[O₂HbMb-HHbMb]-BP showed a better relationship with ventilatory thresholds in trained cyclists (r = 0.68–0.84, p<0.05) compared to endurance and recreationally trained males (r = 0.48–0.50, p<0.05). Although relative PO was not different between groups at Δ[O₂HbMb-HHbMb]-BP, relative PO at VT1 was significantly lower in recreational males compared to other groups (RM< CM~EM~CF, p<0.05). Therefore, ventilatory thresholds differed across sexes and training status, whereas Δ[O₂Hb-HHb]-BP differed only across sexes.

[Figure omitted. See PDF.]

Fig 1. Typical example of detection of Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds.

Ventilatory thresholds VT1 and VT2 (closed circles) are depicted within the ventilatory equivalent versus PO graph (upper graphs, a). Δ[O₂HbMb-HHbMb]-BP (open square) is derived from Δ([O₂HbMb]–[HHbMb]) versus PO graph by means of double linear regression analyses with the least combined residual sum of squares (lower graphs, b). Both Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds are determined for subsequent trials in a trained male cyclist for test (left graphs) and retest (right graphs).

https://doi.org/10.1371/journal.pone.0162914.g001

[Figure omitted. See PDF.]

Fig 2. Relationship between Δ[O₂HbMb-HHbMb]-BP and VT1.

Δ[O₂HbMb-HHbMb]-BP and VT1 are moderately related for (a) en PO (b). Individual values are shown for male cyclists (closed triangle), female cyclists (closed circle), endurance trained males (open square), recreationally trained males (open triangle). The thick solid line represents the best linear fit obtained from individual data at Δ[O₂HbMb-HHbMb]-BP and VT1 and the dotted line represents extrapolation of this best linear fit. The thin solid line is the line of identity.

https://doi.org/10.1371/journal.pone.0162914.g002

[Figure omitted. See PDF.]

Table 2. Group differences in SmO₂ values at rest, exercise and end arterial occlusion and in PO and values at the exercise thresholds.

https://doi.org/10.1371/journal.pone.0162914.t002

Reproducibility

Reproducibility was assessed in twenty trained cyclists (CM and CF). During the two maximal incremental tests, these cyclists displayed comparable values (mean difference: -0.012±0.179 L/min) and peak workload (mean difference: -1.3±8W). In rest, during exercise and during occlusion, SmO₂ measurements showed high reproducibility (mean differences between 0.5–1.8%, ICC = 0.69–0.97; Table 3). Moreover, reproducibility of and PO values was high at Δ[O₂HbMb-HHbMb]-BP (ICC = 0.80–0.88) and excellent at VT1 and VT2 (ICC = 0.96–0.99; Table 3).

[Figure omitted. See PDF.]

Table 3. Reproducibility of SmO₂ values at rest, exercise and arterial occlusion and reproducibility of PO and at the exercise thresholds among twenty trained cyclists.

https://doi.org/10.1371/journal.pone.0162914.t003

Adipose tissue

SmO₂ values during exercise and occlusion were higher in groups with larger adipose tissue thickness (Table 2). Accordingly, ATT explained ~80% of SmO₂ variance at peak exercise intensity and at the end of arterial occlusion (Fig 3). However, when SmO₂ was normalized to correct for ATT, at VT1 (75±15.3%) and VT2 (89±8.5%) group differences disappeared. During maximal incremental exercise and occlusion, we observed asymmetrical changes in the amplitude of [HHbMb] and [O₂HbMb] (Δ[HHbMb]/Δ[O₂HbMb] ≠ 1). For arterial occlusion as well as for maximal incremental exercise this asymmetry was significantly related to ATT (r = -0.50, p<0.05 and r = -0.64, p<0.01 respectively; Fig 3).

[Figure omitted. See PDF.]

Fig 3. ATT affects SmO₂ and symmetry of deoxygenated and oxygenated haemoglobin and myoglobin.

Adipose tissue thickness is hyperbolically related to SmO₂ at peak exercise (a) and at the end of arterial occlusion (b). Symmetry of (de)oxygenated haemoglobin (ΔHHbMb/ΔO₂HbMb) during maximal exercise (c) and from rest to the end of arterial occlusion (d) are linearly related to adipose tissue thickness. Individual values are shown for male cyclists (closed triangle), female cyclists (closed circle), endurance trained males (open square), recreationally trained males (open triangle). Note that occlusion data were successfully obtained in 36 out of 40 subjects.

https://doi.org/10.1371/journal.pone.0162914.g003

Discussion

The present findings show: 1) Δ[O₂HbMb-HHbMb]-BP was not different from VT1 in male and female cyclists and recreationally trained males, but was significantly smaller than VT1 in endurance trained males, and Δ[O₂HbMb-HHbMb]-BP was only moderately related to VT1 and VT2 in all groups, 2) Ventilatory thresholds showed differences across sexes and training status, whereas Δ[O₂HbMb-HHbMb]-BP differed only across sexes, 3) Reproducibility was high for SmO₂, Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds and 4) SmO₂ values were strongly affected by ATT and NIRS oxygenation present an asymmetry in [O₂HbMb] and [HHbMb] that was related to ATT.

Critical exercise threshold

The present study shows that Δ[O₂HbMb-HHbMb]-BP was not significantly different from VT1 in CF, CM and RM. However, Δ[O₂HbMb-HHbMb]-BP preceded VT1 in EM, which consisted of 7 male cyclists and 4 non-cyclists. Note that additional analysis revealed that in a combined group of these 7 EM cyclists and 10 CM cyclists, the Δ[O₂HbMb-HHbMb]-BP was not different from VT1 (i.e. average PO value of 17 male cyclists was 224 W for Δ[O₂Hb-HHb]-BP and 237 W for VT1, p>0.05). Therefore, it may be argued that in EM the Δ[O2HbMb-HHbMb]-BP was significantly smaller than VT1 due to the 4 non-cyclists (i.e. average PO of those 4 non-cyclists was 198 W for Δ[O₂HbMb-HHbMb]-BP and 259 W for VT1, p<0.01). Previous studies have shown that Δ[O₂HbMb-HHbMb]-BP preceded VT1 [18,19] and LT1 [18] or that Δ[O₂HbMb-HHbMb]-BP was not different from LT1 [17,19]. In addition, we showed that the correlation between Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds was rather moderate (r = 0.58–0.68). This deviates from previous findings in studies using similar methodology, demonstrating an excellent relationship between Δ[O₂HbMb-HHbMb]-BP and VT1 or lactate threshold during stepwise incremental cycling exercise (r>0.88) [17–19]. A possible explanation for these contrasting findings is that step increments (30–50W) in the present study were larger than increments in previous studies (30W) [17–19]. Even though the regression model quality is high (r² = 0.97), this could have led to more stepwise kinetics affecting Δ[O₂HbMb-HHbMb]-BP’s goodness of fit (illustrated in Fig 1). Moreover, measurement errors related to visual inspection of ventilatory thresholds and errors related to the detection of Δ[O₂HbMb-HHbMb]-BP (i.e. double linear regression model, motion artefacts, or adipose tissue (as discussed in detail below)) may have increased the unexplained variance in the relation between Δ[O₂HbMb-HHbMb]-BP and VT1. Note, however, that reproducibility of Δ[O₂HbMb-HHbMb]-BP and VT1 were high (ICC = 0.80–0.99). Between-subject variation, in Δ[O₂HbMb-HHbMb]-BP was higher than reported in previous studies [17–19], which may be due to a more heterogeneous subject group in the present study. This heterogeneous group could have resulted in larger inter-individual variation in NIRS patterns during exercise. Note that we observed better relationships between Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds in a more homogeneous group of trained cyclists (r = 0.68–0.84) as opposed to a more heterogeneous group of endurance and recreationally trained males (r = 0.48–0.50).

Physiological differences between sexes are reflected by differences in Δ[O₂HbMb-HHbMb]-BP and ventilatory thresholds. The absolute exercise intensity at VT1 and Δ[O₂HbMb-HHbMb]-BP (reflecting when aerobic ATP resynthesis can no longer match ATP use in the working muscles) was lower in females compared to males. The lower absolute exercise intensity in females may be explained by a lower oxygen carrying capacity and lower cardiac output [29,30,32] and possibly by smaller arteriovenous oxygen differences [29] in females. Expressed in percentage peak power output, Δ[O₂HbMb-HHbMb]-BP and VT1 were found to be similar for trained male and female cyclists, indicating that relative to PO_peak the matching and consequent occurrence of increased anaerobic energy production in the muscle are similar between sexes. Previously, VT1 (in %) has also shown to be similar between recreationally trained males and females [10]. However, at a higher intensity exercise, the Δ[HHbMb]-BP (in %) has demonstrated to be lower in females compared to males, but is likely explained by less sufficient blood distribution towards active muscle in females [10].

Physiological differences between males with different training status did become apparent by differences in ventilatory thresholds, but not by differences in Δ[O₂HbMb-HHbMb]-BP. In contrast, it was expected that exercise intensity (absolute and in %PO_peak) at both ventilatory thresholds and Δ[O₂HbMb-HHbMb]-BP would be higher in trained males, because of a higher gross efficiency in trained subjects [37] and since endurance training enhances oxidative metabolism (e.g. increases in mitochondrial density, oxidative enzymes and percentage slow-twitch myosin heavy chain types; [29,30]) and oxygen supply to the mitochondria (e.g. stroke volume, cardiac output, capillary density, haematocrit, Ca-vO₂ and O₂ diffusion; [29,30]). Also, subjects with a higher aerobic fitness level have been reported to demonstrate a higher percentage of slow-twitch fibres, which has been associated with more effective matching [31]. However, our results on Δ[O₂HbMb-HHbMb]-BP did not confirm these expectations. Possibly, in the present study, physiological differences in training status between male groups were too small and/or measurements of Δ[O₂HbMb-HHbMb]-BP were not sensitive enough to quantify differences in exercise intensity at Δ[O₂HbMb-HHbMb]-BP. Although the Δ[O₂HbMb-HHbMb] breakpoint is potentially a suitable exercise threshold revealing when anaerobic energy production starts to increase in the muscle, VT1 being a rather indirect measure of these changes in energy status of the muscle, discriminates better across sexes and training status in the present study.

Reproducibility

Saturation measured by our NIRS device was obtained in trained cyclists and showed high reproducibility in rest, during exercise and at the end of arterial occlusion. These findings are similar to previously reported reproducibility values of SmO₂ measurements during maximal incremental exercise (ICC = 0.81–0.95) [20], at (r = 0.99) [5] and following arterial arm or leg occlusion (ICC = 0.95–0.96) [21,22]. Thus, these results imply that in trained cyclists one is able to measure oxygenation reproducibly with a continuous-wave NIRS device.

Previous results [23] showed that in sedentary subjects, breakpoints in [O₂HbMb] are detected reproducibly during ramp exercise (r = 0.67–0.85, p<0.05). The present study shows that also in a group of trained cyclists the reproducibility was high for the Δ[O₂HbMb-HHbMb] breakpoint (ICC = 0.80–0.88). Note that in contrast to other studies that assessed an oxygenation threshold from NIRS signals [27,38–41], the present study did not exclude any of the participating subjects from analysis. In addition, our assessment of ventilatory thresholds demonstrates excellent reproducibility (ICC = 0.96–0.99). It is well-known that ventilatory thresholds can be obtained with good reproducibility (r = 0.91–0.98) and can accurately predict endurance time trial performance [4,5]. Our reproducibility results showed that systematic measurement errors related to visual inspection of ventilatory thresholds are smaller than the measurement errors related to detection of Δ[O₂HbMb-HHbMb]-BP, which may be due to effects of stepwise kinetics on the double linear regression model or effects of motion artefacts on the combined Δ[O₂HbMb-HHbMb] signal. Hence, our results indicated that measurement errors related to the reproducibility of NIRS and respiratory measurements only marginally account for the unexplained variance in the relation between Δ[O₂HbMb-HHbMb]-BP and VT1.

Adipose tissue thickness

ATT is known to affect the amplitude of NIRS measurements [7,24–26,42]. Essentially, the NIRS device assumes one homogeneous medium for extending the Lambert-Beer law of light attenuation [43]. This assumption is violated by light scattering and absorption by adipose tissue [7]. Even though the VL adipose tissue thickness was below 12.5mm in all subjects of the present study, ATT comprised 11.0–64.7% of the measurement depth [24,44]. Our results confirmed that SmO₂ values were substantially affected by ATT, explaining ~80% of the variance in SmO₂ at peak exercise and arterial occlusion. These results indicated that NIRS saturation measurements were largely determined by ATT. Although frequently assessed by linear regressions, the relationship between NIRS signals and ATT is better described hyperbolically [22,25,26]. Clearly, the hyperbola (Fig 3) showed that SmO₂ cannot be calculated for extremely low ATT values (i.e. below ~2 mm). This finding corresponds with our observation that the NIRS receiver detects no light in extremely lean subjects as the light is fully absorbed or scattered by underlying muscle tissue. Note that in extremely lean subjects the skinfold thickness predominantly consists of skin tissue. In addition to adipose tissue, cutaneous dilation in response to increased temperature will likely affect NIRS measurements as well [45]. Furthermore, lean subjects (trained males) showed SmO₂ values that were roughly in line with mixed venous saturation values reported in literature at rest (~70%), VT1 (~35%) and peak exercise (15–30%) [46–48]. Thus, NIRS saturation measurements (SmO₂) that are obtained from a combination of [HHbMb] and [O₂HbMb] signals are clearly affected by ATT.

One may reduce effects of ATT on the amplitude of [HHbMb] and [O₂HbMb] signals by normalizing oxygenation changes to the oxygenation at peak exercise [10,27], maximal voluntary contraction [26,28] or cuff occlusion [21,26]. Still, this correction does not account for the evident asymmetry in Δ[O₂HbMb] and Δ[HHbMb] amplitude. Asymmetrical changes during exercise may partially be explained by blood volume changes [26,49], increased haemoglobin concentration (i.e. haemoconcentration) [50] or scattering alterations [43]. For example, continuous-wave NIRS device incorrectly assumes fixed scattering coefficients and therefore can overestimate changes in [HHbMb] (i.e. short wavelengths <750nm) during exercise [43,51]. However, the asymmetry during arterial occlusion may not be explained by changes in blood volume and haemoglobin concentration (both assumed rather constant during occlusion) or scattering coefficients (contrastingly Δ[O₂HbMb] was generally larger than Δ[HHbMb]). Recent findings in canine muscle clearly showed that [O₂HbMb] and [HHbMb] signals change symmetrically when no skin and adipose tissue layer is present [52]. In the present study, the asymmetry during incremental exercise as well as occlusion was significantly related to ATT. Consequently, ATT may affect the amplitude of [O₂HbMb] and [HHbMb] signals to a different extent and therefore may affect kinetics when these are derived from combinations of [HHbMb] and [O₂HbMb] signals (i.e. tHb: [O₂HbMb+HHbMb], Hb difference: [O₂HbMb-HHbMb] and SmO₂). ATT may also affect determination of the Δ[O₂HbMb-HHbMb]-BP if breakpoints in the individual [O₂HbMb] and [HHbMb] signals do not occur at the same time. However, reproducibility of Δ[O₂HbMb-HHbMb]-BP and its relationship with VT1 did not improve when the Δ[O₂HbMb-HHbMb] signal was corrected for the observed asymmetry in [O₂HbMb] and [HHbMb] amplitude (results not presented). Hence, ATT scattering and absorbance likely affect [O₂HbMb] and [HHbMb] amplitude differently, but it remains to be established exactly how this occurs and whether this may confound the detection of Δ[O₂HbMb-HHbMb]-BP.

Conclusions

During maximal incremental step exercise, the Δ[O₂HbMb-HHbMb] breakpoint is reproducible and coincides with VT1, but is only moderately related to the ventilatory thresholds. Although the Δ[O₂HbMb-HHbMb] breakpoint is potentially a suitable exercise threshold revealing when anaerobic energy production starts to increase in the muscle, the VT1 (being a more indirect measure of these changes in muscle energy status) shows higher reproducibility and discriminates better across sexes and training status. Continuous-wave NIRS measurements are reproducible, but adipose tissue thickness strongly affects the amplitude of [O₂HbMb] and [HHbMb] signals to a different extent and thereby may also affect the kinetics of combined [O₂HbMb] and [HHbMb] signals. The large inter-individual variability observed at the Δ[O₂HbMb-HHbMb] breakpoint and its underlying nature is yet poorly understood, although some of this variability may be explained by differences in adipose tissue thickness. Correcting the Δ[O₂HbMb-HHbMb] breakpoint for ATT however is difficult, since it remains to be established how ATT affects the absorbance by O₂HbMb and HHbMb at the measured wavelengths. Therefore, the first ventilatory threshold is favoured over the Δ[O₂HbMb-HHbMb] breakpoint in context of determining an exercise threshold that reflects the metabolic demands of active muscle in maximal stepwise incremental exercise.

Supporting Information

[Figure omitted. See PDF.]

S1 File. Supporting excel-file with data of the present study.

https://doi.org/10.1371/journal.pone.0162914.s001

(XLSX)

Acknowledgments

The authors’ work was supported by a research grant from Technologiestichting STW, the Netherlands.

Author Contributions

1. Conceptualization: SvdZ RTJ MJH KL DAN AdH JJdK WvdL CJdR.

2. Data curation: SvdZ IJB CA JMV AdU CJdR.

3. Formal analysis: SvdZ IJB CA JMV AdU CJdR.

4. Funding acquisition: RTJ MJH DAN AdH JJdK WvdL CJdR.

5. Investigation: SvdZ IJB CA JMV AdU CJdR.

6. Methodology: SvdZ IJB CA JMV AdU RTJ MJH KL DAN AdH JJdK WvdL CJdR.

7. Project administration: SvdZ RTJ AdH JJdK CJdR.

8. Resources: SvdZ IJB CA JMV AdU CJdR.

9. Software: SvdZ, MJH CJdR.

10. Supervision: RTJ AdH JJdK CJdR.

11. Validation: SvdZ IJB CA JMV AdU CJdR.

12. Visualization: SvdZ.

13. Writing – original draft: SvdZ CJdR.

14. Writing – review & editing: SvdZ RTJ IJB CA JMV AdU MJH KL DAN AdH JJdK WvdL CJdR.

Citation: van der Zwaard S, Jaspers RT, Blokland IJ, Achterberg C, Visser JM, den Uil AR, et al. (2016) Oxygenation Threshold Derived from Near-Infrared Spectroscopy: Reliability and Its Relationship with the First Ventilatory Threshold. PLoS ONE 11(9): e0162914. https://doi.org/10.1371/journal.pone.0162914

References

1. Seiler KS, Kjerland GØ. Quantifying training intensity distribution in elite endurance athletes: is there evidence for an “optimal” distribution? Scand J Med Sci Sports. 2006;16: 49–56. pmid:16430681

2. Wasserman K. Coupling of external to cellular respiration during exercise: the wisdom of the body revisited. Am J Physiol. 1994;266: E519–539. pmid:8178973

3. Wasserman K. Determinants and detection of anaerobic threshold and consequences of exercise above it. Circulation. 1987;76: VI29–39. pmid:3315297

4. Amann M, Subudhi AW, Walker J, Eisenman P, Shultz B, Foster C. An evaluation of the predictive validity and reliability of ventilatory threshold. Med Sci Sports Exerc. 2004;36: 1716–1722. pmid:15595292

5. Austin KG, Daigle KA, Patterson P, Cowman J, Chelland S, Haymes EM. Reliability of near-infrared spectroscopy for determining muscle oxygen saturation during exercise. Res Q Exerc Sport. 2005;76: 440–449. pmid:16739682

6. Richardson RS, Noyszewski EA, Kendrick KF, Leigh JS, Wagner PD. Myoglobin O2 desaturation during exercise. Evidence of limited O2 transport. J Clin Invest. 1995;96: 1916–1926. pmid:7560083

7. Ferrari M, Muthalib M, Quaresima V. The use of near-infrared spectroscopy in understanding skeletal muscle physiology: recent developments. Philos Trans R Soc Math Phys Eng Sci. 2011;369: 4577–4590.

8. Koga S, Rossiter HB, Heinonen I, Musch TI, Poole DC. Dynamic Heterogeneity of Exercising Muscle Blood Flow and O2 Utilization. Med Sci Sports Exerc. 2014;46: 860–876. pmid:24091989

9. Spencer MD, Murias JM, Paterson DH. Characterizing the profile of muscle deoxygenation during ramp incremental exercise in young men. Eur J Appl Physiol. 2012;112: 3349–3360. pmid:22270488

10. Murias JM, Keir DA, Spencer MD, Paterson DH. Sex-related differences in muscle deoxygenation during ramp incremental exercise. Respir Physiol Neurobiol. 2013;189: 530–536. pmid:23994824

11. Boone J, Barstow TJ, Celie B, Prieur F, Bourgois J. The impact of pedal rate on muscle oxygenation, muscle activation and whole-body VO₂ during ramp exercise in healthy subjects. Eur J Appl Physiol. 2015;115: 57–70. pmid:25204279

12. Murias JM, Spencer MD, Keir DA, Paterson DH. Systemic and vastus lateralis muscle blood flow and O2 extraction during ramp incremental cycle exercise. Am J Physiol—Regul Integr Comp Physiol. 2013;304: R720–R725. pmid:23515617

13. Boone J, Barstow TJ, Celie B, Prieur F, Bourgois J. The interrelationship between muscle oxygenation, muscle activation and pulmonary vo2 to incremental ramp exercise: influence of aerobic fitness. Appl Physiol Nutr Metab. 2015;

14. Keir DA, Fontana FY, Robertson TC, Murias JM, Paterson DH, Kowalchuk JM, et al. Exercise Intensity Thresholds: Identifying the Boundaries of Sustainable Performance. Med Sci Sports Exerc. 2015;47: 1932–1940. pmid:25606817

15. Fontana FY, Keir DA, Bellotti C, De Roia GF, Murias JM, Pogliaghi S. Determination of respiratory point compensation in healthy adults: Can non-invasive near-infrared spectroscopy help? J Sci Med Sport Sports Med Aust. 2015;18: 590–595.

16. Bellotti C, Calabria E, Capelli C, Pogliaghi S. Determination of maximal lactate steady state in healthy adults: can NIRS help? Med Sci Sports Exerc. 2013;45: 1208–1216. pmid:23274611

17. Grassi B, Quaresima V, Marconi C, Ferrari M, Cerretelli P. Blood lactate accumulation and muscle deoxygenation during incremental exercise. J Appl Physiol Bethesda Md 1985. 1999;87: 348–355.

18. Wang B, Tian Q, Zhang Z, Gong H. Comparisons of local and systemic aerobic fitness parameters between finswimmers with different athlete grade levels. Eur J Appl Physiol. 2011;112: 567–578. pmid:21611824

19. Wang B, Xu G, Tian Q, Sun J, Sun B, Zhang L, et al. Differences between the Vastus Lateralis and Gastrocnemius Lateralis in the Assessment Ability of Breakpoints of Muscle Oxygenation for Aerobic Capacity Indices During an Incremental Cycling Exercise. J Sports Sci Med. 2012;11: 606–613. pmid:24150069

20. Thiel C, Vogt L, Himmelreich H, Hübscher M, Banzer W. Reproducibility of muscle oxygen saturation. Int J Sports Med. 2011;32: 277–280. pmid:21271493

21. Celie B, Boone J, Van Coster R, Bourgois J. Reliability of near infrared spectroscopy (NIRS) for measuring forearm oxygenation during incremental handgrip exercise. Eur J Appl Physiol. 2012;112: 2369–2374. pmid:21952981

22. Southern WM, Ryan TE, Reynolds MA, McCully K. Reproducibility of near-infrared spectroscopy measurements of oxidative function and postexercise recovery kinetics in the medial gastrocnemius muscle. Appl Physiol Nutr Metab. 2013;39: 521–529. pmid:24766233

23. Miura T, Takeuchi T, Sato H, Nishioka N, Terakado S, Fujieda Y, et al. Skeletal Muscle Deoxygenation During Exercise Assessed by Near-Infrared Spectroscopy and its Relation to Expired Gas Analysis Parameters. Jpn Circ J. 1998;62: 649–657. pmid:9766702

24. McCully KK, Hamaoka T. Near-infrared spectroscopy: what can it tell us about oxygen saturation in skeletal muscle? Exerc Sport Sci Rev. 2000;28: 123–127. pmid:10916704

25. van Beekvelt MC, Borghuis MS, van Engelen BG, Wevers RA, Colier WN. Adipose tissue thickness affects in vivo quantitative near-IR spectroscopy in human skeletal muscle. Clin Sci Lond Engl 1979. 2001;101: 21–28.

26. de Ruiter CJ, de Boer MD, Spanjaard M, de Haan A. Knee angle-dependent oxygen consumption during isometric contractions of the knee extensors determined with near-infrared spectroscopy. J Appl Physiol Bethesda Md 1985. 2005;99: 579–586.

27. Chin LMK, Kowalchuk JM, Barstow TJ, Kondo N, Amano T, Shiojiri T, et al. The relationship between muscle deoxygenation and activation in different muscles of the quadriceps during cycle ramp exercise. J Appl Physiol Bethesda Md 1985. 2011;111: 1259–1265.

28. de Ruiter CJ, Maas EA, Wesseling MGH, de Haan A. Knee extensor fatigue threshold is related to whole-body VO2max. Med Sci Sports Exerc. 2012;44: 1366–1374. pmid:22246220

29. Mitchell JH, Tate C, Raven P, Cobb F, Kraus W, Moreadith R, et al. Acute response and chronic adaptation to exercise in women. Med Sci Sports Exerc. 1992;24: S258–265. pmid:1625551

30. Ogawa T, Spina RJ, Martin WH, Kohrt WM, Schechtman KB, Holloszy JO, et al. Effects of aging, sex, and physical training on cardiovascular responses to exercise. Circulation. 1992;86: 494–503. pmid:1638717

31. McDonough P, Behnke BJ, Padilla DJ, Musch TI, Poole DC. Control of microvascular oxygen pressures in rat muscles comprised of different fibre types. J Physiol. 2005;563: 903–913. pmid:15637098

32. Peltonen JE, Hägglund H, Koskela-Koivisto T, Koponen AS, Aho JM, Rissanen A-PE, et al. Alveolar gas exchange, oxygen delivery and tissue deoxygenation in men and women during incremental exercise. Respir Physiol Neurobiol. 2013;188: 102–112. pmid:23707876

33. De Pauw K, Roelands B, Cheung SS, de Geus B, Rietjens G, Meeusen R. Guidelines to classify subject groups in sport-science research. Int J Sports Physiol Perform. 2013;8: 111–122. pmid:23428482

34. de Ruiter CJ, Goudsmit JFA, Van Tricht JA, de Haan A. The isometric torque at which knee-extensor muscle reoxygenation stops. Med Sci Sports Exerc. 2007;39: 443–453. pmid:17473770

35. Beaver WL, Wasserman K, Whipp BJ. A new method for detecting anaerobic threshold by gas exchange. J Appl Physiol. 1986;60: 2020–2027. pmid:3087938

36. Evans JD. Straightforward statistics for the behavioral sciences. Pacific Grove: Brooks/Cole Pub. Co.; 1996.

37. Hopker JG, Coleman DA, Wiles JD. Differences in efficiency between trained and recreational cyclists. Appl Physiol Nutr Metab Physiol Appliquée Nutr Métabolisme. 2007;32: 1036–1042.

38. Racinais S, Buchheit M, Girard O. Breakpoints in ventilation, cerebral and muscle oxygenation, and muscle activity during an incremental cycling exercise. Front Physiol. 2014;5.

39. Osawa T, Kime R, Hamaoka T, Katsumura T, Yamamoto M. Attenuation of Muscle Deoxygenation Precedes EMG Threshold in Normoxia and Hypoxia. Med Sci Sports Exerc. 2011;43: 1406–1413. pmid:21266933

40. Wang L, Yoshikawa T, Hara T, Nakao H, Suzuki T, Fujimoto S. Which common NIRS variable reflects muscle estimated lactate threshold most closely? Appl Physiol Nutr Metab Physiol Appliquée Nutr Métabolisme. 2006;31: 612–620.

41. Rissanen A-PE, Tikkanen HO, Koponen AS, Aho JM, Hägglund H, Lindholm H, et al. Alveolar gas exchange and tissue oxygenation during incremental treadmill exercise, and their associations with blood O2 carrying capacity. Front Physiol. 2012;3: 265. pmid:22934021

42. Koga S, Poole DC, Fukuoka Y, Ferreira LF, Kondo N, Ohmae E, et al. Methodological validation of the dynamic heterogeneity of muscle deoxygenation within the quadriceps during cycle exercise. Am J Physiol Regul Integr Comp Physiol. 2011;301: R534–541. pmid:21632845

43. Ferreira LF, Hueber DM, Barstow TJ. Effects of assuming constant optical scattering on measurements of muscle oxygenation by near-infrared spectroscopy during exercise. J Appl Physiol Bethesda Md 1985. 2007;102: 358–367.

44. Ferrari M, Quaresima V. A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. NeuroImage. 2012;63: 921–935. pmid:22510258

45. Tew GA, Ruddock AD, Saxton JM. Skin blood flow differentially affects near-infrared spectroscopy-derived measures of muscle oxygen saturation and blood volume at rest and during dynamic leg exercise. Eur J Appl Physiol. 2010;110: 1083–1089. pmid:20700602

46. Stringer W, Wasserman K, Casaburi R, Pórszász J, Maehara K, French W. Lactic acidosis as a facilitator of oxyhemoglobin dissociation during exercise. J Appl Physiol Bethesda Md 1985. 1994;76: 1462–1467.

47. Mortensen SP, Damsgaard R, Dawson EA, Secher NH, González-Alonso J. Restrictions in systemic and locomotor skeletal muscle perfusion, oxygen supply and VO2 during high-intensity whole-body exercise in humans. J Physiol. 2008;586: 2621–2635. pmid:18372307

48. McArdle WD, Katch VL. Exercise Physiology: Energy, Nutrition, and Human Performance. 6th Revised edition edition. Philadelphia: Lippincott Williams and Wilkins; 2006.

49. Bhambhani YN. Muscle Oxygenation Trends During Dynamic Exercise Measured by Near Infrared Spectroscopy. Can J Appl Physiol. 2004;29: 504–523. pmid:15328597

50. Davis ML, Barstow TJ. Estimated contribution of hemoglobin and myoglobin to near infrared spectroscopy. Respir Physiol Neurobiol. 2013;186: 180–187. pmid:23357615

51. Mancini DM, Bolinger L, Li H, Kendrick K, Chance B, Wilson JR. Validation of near-infrared spectroscopy in humans. J Appl Physiol. 1994;77: 2740–2747. pmid:7896615

52. Goodwin ML, Hernández A, Lai N, Cabrera ME, Gladden LB. VO2 on-kinetics in isolated canine muscle in situ during slowed convective O2 delivery. J Appl Physiol Bethesda Md 1985. 2012;112: 9–19.