Aerobic Interval Training Impacts Muscle and Brain Oxygenation Responses to Incremental Exercise.

The purpose of the present study was to assess the effects of aerobic interval training on muscle and brain oxygenation to incremental ramp exercise. Eleven physically active subjects performed a 6-week interval training period, proceeded and followed by an incremental ramp exercise to exhaustion (25 W min-1). Throughout the tests pulmonary gas exchange and muscle (Vastus Lateralis) and brain (prefrontal cortex) oxygenation [concentration of deoxygenated and oxygenated hemoglobin, HHb and O2Hb, and tissue oxygenation index (TOI)] were continuously recorded. Following the training intervention V . ⁢ O 2 peak had increased with 7.8 ± 5.0% (P < 0.001). The slope of the decrease in muscle TOI had decreased (P = 0.017) 16.6 ± 6.4% and the amplitude of muscle HHb and totHb had increased (P < 0.001) 40.4 ± 15.8 and 125.3 ± 43.1%, respectively. The amplitude of brain O2Hb and totHb had increased (P < 0.05) 40.1 ± 18.7 and 26.8 ± 13.6%, respectively. The training intervention shifted breakpoints in muscle HHb, totHb and TOI, and brain O2Hb, HHb, totHb and TOI to a higher absolute work rate and V . ⁢ O 2 (P < 0.05). The relative (in %) change in V . ⁢ O 2 peak was significantly correlated to relative (in %) change slope of muscle TOI (r = 0.69, P = 0.011) and amplitude of muscle HHb (r = 0.72, P = 0.003) and totHb (r = 0.52, P = 0.021), but not to changes in brain oxygenation. These results indicate that interval training affects both muscle and brain oxygenation, coinciding with an increase in aerobic fitness (i.e., V . ⁢ O 2 peak). The relation between the change in V . ⁢ O 2 peak and muscle but not brain oxygenation suggests that brain oxygenation per se is not a primary factor limiting exercise tolerance during incremental exercise.

Introduction

Incremental ramp exercise tests are commonly used in healthy and pathologic populations to determine physical fitness, to identify intensity zones for training/rehabilitation, and to assess the efficacy of training and/or treatment interventions. Next to the measurement of pulmonary gas exchange, which provides insight into whole-body physiological responses to increasing exercise intensity, also peripheral measures of regional oxygenation (predominantly at the level of the locomotor muscles) have been performed in the recent past. More specifically, the signals [the concentration of oxygenated (O2Hb) and deoxygenated hemoglobin (HHb)] derived from near-infrared spectroscopy reflect the relationship between O2 delivery and O2 utilization at the level of the microcirculation (Mancini et al., 1994). As such, more specific information can be obtained on the peripheral physiological responses and on the limiting factors of different populations to incremental exercise.

The NIRS responses at the level of the locomotor muscles have already been characterized extensively (for review see Boone et al., 2016) in healthy subjects, trained subjects, children, elderly and some patient populations. Traditionally, it is observed that muscle HHb, which is often considered as a reflection of fractional O2 extraction (Delorey et al., 2003; Grassi et al., 2003), increases following a sigmoid-like pattern, with a sluggish increase at the onset of the incremental exercise, followed by a linear increase to finally reach a plateau at high intensities (∼80–90%VO2max) (Ferreira et al., 2007; Spencer et al., 2012; Fontana et al., 2015). In this context, it has been shown from cross-sectional studies that trained subjects have a less steep increase but higher peak value in HHb, reflecting an improved matching between O2 supply and O2 demand and higher peak O2 extraction, respectively (Boone et al., 2009; Gifford et al., 2016; Okushima et al., 2016). The longitudinal effects of a training intervention on muscle oxygenation are scarcely documented. Prieur and Mucci (2013) found an increased amplitude of the HHb response (i.e., indicating an improved O2 extraction capacity) following 6 weeks of interval training at high intensities. Takagi et al. (2016) found a similar impact on the amplitude of HHb following an aerobic training program in post-myocardial infarction patients.

Recently, also oxygenation responses to incremental exercise at the level of the brain have gained interest since it has been argued that the brain might be involved in the process of termination of maximal exercise (Robertson and Marino, 2016), especially since Nielsen et al. (Nielsen et al., 1999) found that O2 supplementation could maintain cerebral oxygenation at a higher level and as such increased performance. In healthy subjects cerebral oxygenation (cO2Hb) increases steadily during incremental exercise compared to baseline resting levels (Rooks et al., 2010). However, at high intensities (in close proximity to the respiratory compensation point) cO2Hb levels-off and even decreases (Bhambhani et al., 2007; Racinais et al., 2014; Oussaidene et al., 2015). In this context, it has even been shown that neural activity in the prefrontal cortex decreases at the respiratory compensation point (Robertson and Marino, 2015). Cross-sectional data show that trained subjects have a more pronounced increase in cO2Hb, indicating an improved cerebral oxygenation compared to less trained counterparts, which might add to the training-induced improvement in exercise tolerance (Rooks et al., 2010; Oussaidene et al., 2015). However, it is currently unclear whether cerebral oxygenation responses, as assessed in cO2Hb measured with NIRS, to incremental exercise are affected by short-term training interventions.

Therefore, the purpose of the present study was to assess the effects of a 6-week interval training program on muscle and cerebral oxygenation responses to incremental exercise. In line with cross-sectional studies we hypothesize first, that the amplitude of mHHb will have increased (as a reflection of an improved O2 extraction capacity) and that the slope of the increase in mHHb relative to work rate will be lower (as a reflection of improved matching between O2 supply and O2 demand) following the training intervention. Second, we hypothesized that the amplitude of cO2Hb will have increased (as a reflection of an improved cerebral oxygenation) and the leveling-off will occur at a higher absolute work rate, since it is has been proposed that the breakpoint in cO2Hb is mechanistically linked to the respiratory compensation point (Bhambhani et al., 2007; Racinais et al., 2014). Third, given the suggestion that cerebral oxygenation possibly affects exercise tolerance and the observation that cerebral oxygenation is improved in trained subjects, we will assess whether changes in muscle and/or cerebral oxygenation following the training intervention contribute to the improvement in VO2peak to obtain insight into the limiting factors of incremental ramp exercise.

Materials and Methods

Ethics Statement

This study was approved by the local ethical committee (Ghent University Hospital, Ghent, Belgium) with the code number EC/2015/1318 and followed the ethical recommendations for the study of humans as suggested by the Declaration of Helsinki. All participants give written informed consent prior to the start of the study.

Subjects

Eleven male physically active students (21.8 ± 1.2 year, 1.81 ± 0.08 m, 75.7 ± 4.0 kg) volunteered to take part in this study. Based on the effect (Cohen’s d: 0.8–2.0) of training interventions on the available NIRS responses (Prieur and Mucci, 2013; Takagi et al., 2016) a subject group of 5–11 subjects would be sufficient to obtain a power of 0.80 with an α of 0.05. All subjects participated in various recreational sports activities on a regular basis (1–3 times per week) and were habituated to maximal exercise efforts, although none of them had a history of cycling training. Prior to the study each participant underwent a medical examination. Each subject was declared to be in good health and no contra-indications for participation were detected.

Experimental Procedure

General Overview

Experimental testing (i.e., incremental ramp exercise) was conducted on an electromagnetically braked cycle ergometer (Lode Excalibur Sport, Groningen, Netherlands) and took place in the laboratory (Sport Science Laboratory – Jacques Rogge, Ghent University) on 2 different occasions prior to and following a 6-week training intervention. Subjects completed a maximal ramp incremental (RI) exercise test to assess their responses and adaptations in general cardiorespiratory fitness [gas exchange threshold (GET); peak oxygen uptake, peak] and in tissue oxygenation at the level of the M. Vastus Lateralis and prefrontal cortex. The study intervention itself consisted of 6-week cycling training with several bouts at a work rate corresponding to the Critical Power, as determined prior to the training period from 4 constant work rate trials (CWR trials) to exhaustion.

Experimental Testing

The incremental ramp exercise started from 3 min of baseline cycling at 50 W after which work rate (WR) increased continuously (25 W min–1). Participants were instructed to keep their cadence between 70 and 80 rpm and strong verbal encouragement was provided throughout the test to ensure maximum effort. The protocol was terminated at voluntary exhaustion or when the subjects’ cadence fell below 70 rpm for more than 5 consecutive seconds. Pulmonary gas exchange (,CO2) was measured breath-by-breath (Jaeger Oxycon Pro, Viasys Healthcare GmbH, Höchberg, Germany) and tissue oxygenation was registered using near-infrared spectroscopy (NIRO-200NX, Hamamatsu Phototonics, Hamamatsu, Japan) at a sampling rate of 0.5 Hz. This device records changes in oxygenated (O2Hb), deoxygenated (HHb), and total hemoglobin (totHb) from baseline values (i.e., seated rest) in μmol l–1 using the modified Beer Lambert law and TOI employing spatially resolved spectroscopy. Muscle oxygenation (mO2Hb, mHHb, mtotHb, and mTOI) was measured at the M. Vastus Lateralis of the right thigh. After shaving and cleaning, the probe was placed longitudinally on the distal section of the muscle belly. Simultaneously, cerebral oxygenation (cO2Hb, cHHb, ctotHb, and cTOI) was measured at the level of the right prefrontal cortex between Fp2 en F4 according to the modified international EEG 10-20 system (Robertson and Marino, 2015). Heart rate (HR) was monitored on a beat-by-beat basis (H7 Sensor, Polar, Kempele, Finland). During testing, environmental conditions were kept constant at a room air temperature of 19° and humidity of 50%.

Training Period

Six weeks of supervised cycling interval training on the ergometer was completed by all eleven subjects. Training sessions took place three times a week with a total of 18 visits. Each training lasted 49 min and included a 5 min warm-up and a 5 min cooling-down at a WR corresponding to the subject’s GET, as determined from the ramp incremental test. To determine the WR that would elicit a steady state corresponding to the GET, the linear /WR-relationship was shifted to the left to account for the mean response time (MRT; Fontana et al., 2015). The main part of each training was composed of six exercise bouts during which participants cycled at their Critical Power for 4 min, alternated with 3 min of active recovery at the level of the GET. Critical power of the subjects was determined prior to the training intervention from the relationship between work rate and time to exhaustion from four CWR trials to exhaustion at 75, 85, 95, and 105% Ppeak. During training sessions, subjects cycled at a self-selected cadence and HR sensors were worn.

Data Analysis

Peak power output (Ppeak) and peak heart rate (HRpeak) were defined as the highest values obtained during the RI test. Breath-by-breath data were transformed into 10s values for further analysis. peak was defined as the highest 30 s average achieved during the test. GET and RCP were determined by four independent researchers. GET was defined as (a) the point where CO2 increased disproportionate to O2, (b) the first departure from the linear increase in minute ventilation (E), and (c) an increase in E/O2 without a simultaneous increase in E/CO2 (Beaver et al., 1986). In case the physiologists encountered conflicting results, data were re-evaluated until mutual agreement was reached. RCP was defined as (a) the second departure from the linear increase in E, (b) an increase in E/CO2, and (c) a systemic fall in end-tidal PCO2 (Whipp et al., 1989).

The NIRS data (O2Hb, HHb, tHb, and TOI) were averaged into 10s-bins, expressed as function of power output and analyzed using a double-linear model (Osawa et al., 2011; Bellotti et al., 2013; Boone et al., 2015) (Sigmaplot 13, Systat Software Inc., San Jose, CA, United States). The data set used in this analysis was chosen on visual inspection by three independent researchers and included all data points between the middle portion of the RI test (i.e., the point where the signal showed a systematic linear pattern) and the end of the RI test. Piecewise linear regression analysis was applied and yielded two linear functions (expect for mO2Hb in which the signal did not show a clear breakpoint, see “Results” section):

where m represents the slope and b corresponds to the y-intercept. Subsequently, a breakpoint (BP) reflecting the intersection of these two linear functions could be determined. The WR at the time-point corresponding to the BPs was adjusted for the O2 MRT in each individual in order to account for the kinetics of O2 and the delay between muscles and lungs (Boone and Bourgois, 2012; Fontana et al., 2015). The MRT was defined as the time interval between the onset of the RI test and the intersection of the forward extrapolation of the baseline O2 and the backward extrapolation of the linear O2-time relationship below GET (Fontana et al., 2015). For each subject, the regression line (y = ax + b) of the O2/time relationship below GET was calculated. To make sure that the linear increase in O2 had already started, the first 2 min of the RI test were omitted from the analysis. Baseline O2 was defined as the mean O2 during the warm-up phase of the RI test, leaving out the first 90 s and the last 30 s. The individual MRT was then used to align the O2 data with the WR data in order to determine the O2 at which the BPs occurred. Additionally, the amplitude of mO2Hb, mHHb, mtotHb, mTOI, cO2Hb, cHHb, ctotHb, and cTOI was calculated as the largest change in the signal from the baseline values. To quantify the change in amplitude (Δ) following the training program the amplitude of the posttest was expressed relative (in %) to the amplitude of the pretest.

Statistical Analysis

All statistics were performed using SPSS Statistics 23 (IBM Corp., New York, NY, United States). Descriptive data are presented as mean values ± SD for n = 11 subjects. The Shapiro-Wilk test indicated that all variables were normally distributed. Paired samples t-tests were performed to detect training effects for Ppeak, O2peak, HRpeak, RERpeak. Also, the slopes, y-intercepts and amplitude of the NIRS-responses were compared using paired samples t-tests to identify pre–post-training differences. The breakpoints (BP) in muscle and cerebral oxygenation were compared by means of Repeated Measures ANOVA. Additionally, to assess the relationship between the change in aerobic fitness (i.e., relative increase from pre to post in O2peak) and changes in muscle and cerebral oxygenation (relative change from pre to post in amplitude and slope), a multiple linear regression analysis was performed with the inclusion of the NIRS variables that were affected by the training program. A stepwise backward elimination method was applied in which variables are eliminated that are redundant and/or do not contribute in predicting the outcome variable (i.e., O2peak) and thus, when their contribution to the coefficient of determination (R2) was not significant. The accuracy of the prediction across independent variables (i.e., contributing NIRS variables) was expressed as adjusted R2. Additionally, Pearson correlations were calculated between the significant variables of the multiple regression analysis and the relative increase in O2peak (in %). Statistical significance level was set at P < 0.05.

Results

General Cardiorespiratory Response

In Table 1 an overview is provided of the general cardiorespiratory response to the ramp incremental tests. The training intervention induced an increase in Ppeak of 7.9 ± 2.1% (P < 0.001), in O2peak of 7.8 ± 5.0% (P < 0.001). RERpeak and HRpeak did not differ significantly post-training compared to pre-training (P = 0.766 and P = 0.875, respectively).

TABLE 1

Overview of the cardiorespiratory responses to incremental ramp exercise pre- and post-training intervention.

	Pre	Post	Δ
Ppeak (Watt)	384 ± 36	415 ± 39	31 ± 8∗
V.O2peak (ml min–1)	3967 ± 308	4272 ± 357	305 ± 217∗
V.O2peak (ml min–1 kg–1)	52.4 ± 3.5	56.4 ± 3.8	4.0 ± 2.5∗
V.CO2peak (ml min–1)	4959 ± 429	5260 ± 471	301 ± 197∗
V.CO2peak (ml min–1 kg–1)	65.5 ± 5.3	69.4 ± 6.2	3.9 ± 2.1∗
RERpeak	1.25 ± 0.06	1.23 ± 0.06	0.01 ± 0.02
HRpeak (bts min–1)	190 ± 10	190 ± 9	0 ± 3

Muscle Oxygenation

In Figure 1 the pattern of mTOI, mHHb, mtotHb, and mO2Hb is presented as function of power output for a representative subject.

FIGURE 1

The pattern of mTOI, mHHb, mtotHb, and mO2Hb expressed as a function of power output for the pretest (white dots) and posttest (black dots) in a representative subject.

In Table 2 an overview is provided of the (double) linear fitting to mHHb, mtotHb, mO2Hb, and mTOI. For mTOI, m1 (−16.6 ± 6.4%, P = 0.017) was significantly lower and BP was significantly higher (P < 0.01) post-training compared to pre-training, whereas A of the response did not differ (P = 0.572). The mTOI at baseline cycling (70.1 ± 4.4% vs. 70.2 ± 4.1%, P = 0.881), at the BP (58.2 ± 3.4% vs. 58.4 ± 3.2%, P = 0.652) and at the end of the exercise test (57.7 ± 3.8% vs. 57.4 ± 3.6%, P = 0.857) did not differ between pre- and post-training. For mHHb, both m1 (61.5 ± 28.6%, P < 0.01) and A (+40.4 ± 15.8%, P < 0.01) were significantly higher post-training compared to pre-training. Also the BP occurred at a higher power output and O2 post-training (P < 0.001). For mtotHb, m1 (+118.8 ± 57.6%, P < 0.001) and A (+125.5 ± 43.1%, P < 0.001) were significantly higher post-training and also the BP occurred at a significantly (P = 0.02) higher power output and O2.

TABLE 2

The mean parameters of the (double) linear regression analysis of mHHb, mtotHb, mO2Hb, and mTOI to the incremental ramp exercise prior to (pre) and following (post) the training intervention.

	mHHb		mtotHb		mO2Hb			mTOI
	Pre	Post	Pre	Post	Pre	Post		Pre	Post
m1 (μmol l–1 Watt–1)	0.052 ± 0.029	0.084 ± 0.045∗	0.016 ± 0.012	0.035 ± 0.019∗	−0.022 ± 0.007	−0.014 ± 0.006∗	m1 (% Watt–1)	0.048 ± 0.012	0.040 ± 0.009∗
y-intercept1 (μmol l–1)	−5.1 ± 4.2	−15.1 ± 8.0∗	−0.6 ± 4.6	−3.8 ± 6.1	6.3 ± 3.8	3.8 ± 2.9	y-intercept1 (%)	73.2 ± 3.4	72.8 ± 2.8
m2 (μmol l–1 Watt–1)	−0.001 ± 0.004	−0.002 ± 0.003	−0.03 ± 0.02	−0.02 ± 0.03			m2 (% Watt–1)	0.009 ± 0.006	0.015 ± 0.006
y-intercept2 (μmol l–1)	8.8 ± 3.0	13.0 ± 3.4	11.7 ± 4.8	18.1 ± 6.9			y-intercept2 (%)	61.2 ± 4.3	64.1 ± 5.2
BP (Watt)	317 ± 20	342 ± 30∗	304 ± 25	330 ± 29∗			BP (Watt)	320 ± 23	347 ± 34∗
BP (ml min–1)	3820 ± 371	4051 ± 390∗	3674 ± 392	3969 ± 421∗			BP (ml min–1)	3847 ± 409	4094 ± 430∗
A (μmol l–1)	9.4 ± 4.8	13.2 ± 5.3∗	4.3 ± 4.1	9.7 ± 5.5∗	−6.0 ± 3.2	−4.7 ± 2.9∗	A (%)	−12.4 ± 3.1%	−12.8 ± 2.8%

The mO2Hb response did not show a uniform pattern for all subjects and based on visual inspection not all subjects showed a deflection of mO2Hb at high intensities (4 in the pretest and 5 in the posttest). Therefore, mO2Hb pattern was analyzed using a linear regression analysis to obtain insight into the slope of the decrease in mO2Hb. For mO2Hb the slope (m1) (−36.2 ± 17.5%, P = 0.017) and the total amplitude (A) (−22.7 ± 10.2%, P = 0.024) of the decrease was significantly lower post-training compared to pre-training.

Cerebral Oxygenation

In Figure 2 the pattern of cTOI, cHHb, ctotHb, and cO2Hb is presented as function of power output for a representative subject. In Table 3 an overview is provided of the double linear fitting to cHHb, ctotHb, cO2Hb, and cTOI. For cTOI, the BP occurred at a higher power output and O2 post-training (P < 0.001). Although m1 did not differ significantly (P = 0.081), cTOI was significantly higher at the power output corresponding to BP (72.6 ± 1.8% vs. 74.6 ± 1.6%, P = 0.039). The cTOI at baseline cycling (68.2 ± 3.6% vs. 68.5 ± 4.1%, P = 0.592) and the end of the exercise test did not differ (P = 0.534) between pre- and post-training (64.3 ± 3.9% vs. 64.8 ± 4.4%, P = 0.724). The cTOI at BP (70.8 ± 4.0% vs. 73.5 ± 3.3%, P = 0.019) was significantly higher post-training compared to pre-training. For cHHb, only the BP was significantly higher post-training compared to pre-training (P < 0.001). For ctotHb, m1 (50.1 ± 14.8%, P = 0.023) and A (+26.8 ± 13.6%, P = 0.031) were significantly higher post-training and also the BP occurred at a significantly higher power output (P = 0.037) and O2 (P = 0.029). For cO2Hb, m1 (+25.0 ± 11.5%, P = 0.034) and A (+40.1 ± 18.7%, P < 0.001) were significantly higher post-training compared to pre-training and also the BP occurred at a significantly higher (P < 0.001) power output and O2.

FIGURE 2

The pattern of cTOI, cHHb, ctotHb, and cO2Hb expressed as a function of power output for the pretest (white dots) and posttest (black dots) in a representative subject.

TABLE 3

The mean parameters of the double linear regression analysis of cHHb, ctotHb, cO2Hb and cTOI to the incremental ramp exercise prior to (pre) and following (post) the training intervention.

	cHHb		ctotHb		cO2Hb			cTOI
	Pre	Post	Pre	Post	Pre	Post		Pre	Post
m1 (μmol l–1 Watt–1)	0.02 ± 0.01	0.02 ± 0.01	0.12 ± 0.06	0.18 ± 0.05∗	0.16 ± 0.04	0.20 ± 0.07∗	m1 (% Watt–1)	0.025 ± 0.008	0.030 ± 0.011
y-intercept1 (μmol l–1)	−2.1 ± 1.9	−3.4 ± 2.2	−9.8 ± 5.2	−15.1 ± 6.2	−22.7 ± 11.4	−31.4 ± 11.7∗	y-intercept1 (%)	65.6 ± 4.4	66.2 ± 4.0
m2 (μmol l–1 Watt–1)	0.10 ± 0.03	0.10 ± 0.03	0.05 ± 0.02	0.05 ± 0.03	−0.05 ± 0.05	−0.07 ± 0.05	m2 (% Watt–1)	−0.067 ± 0.022	−0.069 ± 0.027
y-intercept2 (μmol l–1)	−15.5 ± 5.2	−17.9 ± 6.2	17.1 ± 7.8	27.8 ± 9.1	38.2 ± 10.7	41.4 ± 12.2	y-intercept2 (%)	92.6 ± 5.1	97.8 ± 6.1
BP (Watt)	251 ± 25	270 ± 33∗	312 ± 28	325 ± 32∗	289 ± 27	311 ± 25∗	BP (Watt)	258 ± 22	275 ± 27∗
BP (ml min–1)	3137 ± 385	3368 ± 402∗	3781 ± 441	3923 ± 381∗	3515 ± 379	3803 ± 337∗	BP (ml min–1)	3207 ± 380	3415 ± 419∗
A (μmol l–1)	16.2 ± 4.8	18.5 ± 5.3	36.5 ± 7.1	46.3 ± 6.8	22.0 ± 10.0	30.8 ± 13.1∗	A (%)	−4.0 ± 2.9%	−4.0 ± 3.1%

Breakpoints in Muscle and Cerebral Oxygenation

The statistical analysis revealed that there was no interaction effect (P = 0.744) indicating that the evolution of the BPs over the training period did not differ. However, there was a specific order in the occurrence of the BPs (expressed in ml min–1). BPs of cHHb and cTOI occurred at a significantly (P < 0.05) lower O2 compared to the other BPs, whereas BPs of mHHb and mTOI occurred at a significantly (P < 0.05) higher O2 compared to mtotHb, ctotHb and cO2Hb. Additionally, there was a significant main effect of training intervention (P < 0.001) indicating that all BPs shifted to a higher O2 post-training compared to pre-training.

O2peak vs. Muscle and Cerebral Oxygenation

The multiple linear regression analysis revealed that only the change in A of mHHb, the change in A of mtotHb and the change in m1 of mTOI contributed significantly to the relative increase in O2peak (Adjusted R2 = 0.68). Change in cerebal oxygenation parameters did not significantly add to the multiple regression analysis (P > 0.05). The relative increase in O2peak (in %) was correlated to the relative change (in %) in m1 of mTOI (Figure 3, upper panel) (r = 0.69, P = 0.011) and, in the amplitude of mHHb (Figure 3, lower left panel) (r = 0.73, P = 0.003) and mtotHb (Figure 3, lower right panel) (r = 0.52, P = 0.021). The improvement in O2peak was not correlated with any relative changes in m1 and A of cerebral oxygenation parameters.

FIGURE 3

Correlation between relative change in VO2peak and relative change in the slope (m1) of mTOI (upper panel), and amplitude of mHHb (lower left panel) and mtotHb (lower right panel).

Discussion

In the present study, the effects of a training intervention on muscle and brain oxygenation during incremental exercise were assessed. It was observed that aerobic interval training affects the patterns of muscle and cerebral oxygenation to incremental ramp exercise, in concert with an improvement with aerobic fitness indices (e.g., O2peak). More specifically, at the level of the muscle, it was found that mTOI decreased at a slower rate following the training intervention. The amplitude of mHHb and mtotHb increased from pre- to post-training, supporting an increased capacity for microvascular O2 extraction in combination with an improved O2 availability. Additionally, it was found that also cerebral O2Hb and totHb displayed a steeper increase at moderate intensities and increased to a higher amplitude, resulting in an improved cerebral oxygenation (cTOI) throughout the test. Finally, it was found that the increase in O2peak was predominantly related to changes in muscle oxygenation (i.e., decrease in the slope of mTOI, increase in amplitude of mHHb and mtotHb) and not to changes in cerebral oxygenation.

Effect of Training on Muscle Oxygenation

The effects of a training intervention on muscle oxygenation responses, assessed with NIRS, to incremental exercise are scarcely documented. In the present study an aerobic interval training program, consisting of work bouts at the level of the critical power (i.e., the boundary between the heavy and the severe intensity domain), resulted in an improved muscle oxygenation throughout the incremental exercise test, as can be deducted from the slower decrease in mTOI post-training. As mtotHb reflects the amount of Hb under the NIRS probe, the increase in the slope and amplitude of mtotHb indicates that the O2 availability at the muscle level throughout the incremental exercise has improved. Since the change with the interval training in mHHb, which can be considered to reflect microvascular O2 extraction, is less pronounced than the change in mtotHb, this will result in a slower decrease in overall tissue oxygenation (i.e., mTOI). The increase in O2 availability (i.e., higher slope and amplitude of mtotHb) has likely resulted from training-induced adaptations in O2 diffusive capacity and/or convective O2 supply. The O2 diffusing capacity is dictated predominantly by capillary hematocrit and the volume density of red blood cell-flowing capillaries (Davis and Barstow, 2013). The strong increase in mtotHb following interval training is thus, likely in part the result of an improvement in capillary density (and thus capillary-to-muscle fiber ratio) and adaptations in vascular control inducing a local redistribution of blood flow to strongly recruited areas (Laughlin and Roseguini, 2008). Next to these adjustments at the level of the microcirculation, affecting the O2 diffusing capacity, also functional and structural adaptations in the conduit arteries (Dinenno et al., 2001; Rowley et al., 2011; Green et al., 2017) and with regards to cardiac output (MacInnis and Gibala, 2017), can improve the convective O2 supply. It has been shown, that there is a certain time course in the adaptations in the different vascular properties, with rapid initial improvements in vascular dilation function to normalize shear stress during exercise bouts. Afterward, a more prolonged training program will induce more permanent structural enlargements (Green et al., 2017). Also, at the level of the microcirculation the adaptations have been reported to occur at an early stage during a training program (Duscha et al., 2011). In this context, it would have been interesting to observe the changes in NIRS responses also in the middle of the training program (i.e., after 3 weeks), next to those at the end (i.e., after 6 weeks) to obtain information on the time course of the adaptations in O2 diffusive capacity and convective O2 supply.

Additionally, the amplitude of mHHb, reflecting microvascular O2 extraction, strongly increased following the training intervention. This response indicates that structural and/or functional adaptations have been induced at the level of the skeletal muscles. In this context Jacobs et al. (2013) have shown that the increase in exercise performance with high-intensity interval training is strongly linked to an increase in skeletal muscle mitochondrial content and function, which has been documented as a general effect of high-intensity interval training (Dinenno et al., 2001). This will in turn affect the driving O2 pressure gradient between capillaries and muscle with a consequent increase in O2 extraction capacity. Previous studies in healthy (Murias et al., 2010; Jacobs et al., 2013; Prieur and Mucci, 2013) and patient populations (Fu et al., 2016; Takagi et al., 2016) have found a similar effect on mHHb following training, whereas cross-sectional studies have also shown a positive relationship between O2peak and the amplitude of the mHHb response during exercise (Okushima et al., 2016).

Interestingly, the change in O2peak showed to be related to the changes in muscle oxygenation parameters. Individuals that could maintain mTOI at a higher level compared to pre-training and that could enhance mHHb (microvascular O2 extraction) to a higher extent showed the strongest improvements in O2peak. The observation that mTOI at the breakpoint, where a leveling-off in mTOI occurs, and at the end of the incremental exercise were similar pre-training compared to post-training suggests that a critical level of mTOI dictates the termination of the incremental exercise test or at least triggers a cascade that leads to this termination.

Effect of Training on Cerebral Oxygenation

Also at cerebral level the training-induced adaptations support an enhanced oxygenation level. The slope (m1) of cTOI is slightly but not significantly steeper, cTOI at the breakpoint is higher, and also the observation that the slope and amplitude of cO2Hb and ctotHb are higher support an improved O2 availability throughout the incremental exercise. It has been found (e.g., Peltonen et al., 2009; Rooks et al., 2010; Oussaidene et al., 2013), also in the present study, that cO2Hb and ctotHb increase up to a high intensity (∼75% Ppeak), where the response first levels off and then starts to decrease. cTOI follows a similar pattern albeit that the breakpoint occurs at a slightly lower intensity. These responses indicate that at high intensities there is a mismatch between O2 supply and O2 demand at the level of the prefrontal cortex. The occurrence of the breakpoint in cO2Hb and ctotHb, following the increase at low to moderate intensities, can be related to the pattern of cerebral blood flow during incremental exercise which depends on the exercise intensity (Jorgensen et al., 1992; Gonzalez-Alonso et al., 2004; Querido and Sheel, 2007). Cerebral blood flow is tightly regulated by the cerebrovascular responsiveness to alterations in partial pressures of O2 (PaO2) and CO2 (PaCO2) (Paulsen et al., 1990). At moderate intensities cerebral blood flow increases in relation to the increase in cardiac output. Once the GET is exceeded, PaCO2 increases and induces a vasodilation resulting in a more pronounced blood flow to the brain. It appears from the higher slope and amplitude of cO2Hb and ctotHb that cerebral blood flow is enhanced following interval training, dictated by a commonly reported increase in cardiac output following interval training (MacInnis and Gibala, 2017) and/or redistribution of the blood flow favoring the prefrontal cortex. These results correspond to the study of Oussaidene et al. (2015) were similar responses were found in trained vs. untrained subjects and indicate that the O2 supply exceeded the demand more following the training intervention. However, at the respiratory compensation point the PaCO2 is reduced (i.e., hypocapnia) due to the ventilatory response to the metabolic acidosis and this will reduce cerebral blood flow. This mechanistic link between RCP and the BP in cO2Hb also explain the tight relationship between the two parameters observed in the present study.

Given this decrease in cerebral oxygenation at a given intensity during incremental exercise, the role of the brain (i.e., cerebral oxygenation) in the termination of incremental exercise has been put forward (e.g., Subudhi et al., 2008). The decrease in cO2Hb, ctotHb, and cTOI following the breakpoint (i.e., at high intensities) in cerebral oxygenation might have a negative impact on central motor command (Ide and Secher, 2000; Subudhi et al., 2007; Rupp and Perrey, 2008) and as such set a limitation to the maximal exercise performance. The present study showed that the changes in cerebral oxygenation parameters (i.e., higher cTOI at BP, higher slope and amplitude of cO2Hb and ctotHb) following the training intervention were not related to the changes in O2peak. This indicates that cerebral oxygenation per se does not dictate the termination of the incremental exercise test, in contrast to the suggestion of Robertson and Marino (2016).

The results of the present study, however, suggest that a tight coupling exists between the muscle and brain during incremental exercise. The breakpoints in the muscle and cerebral oxygenation parameters occur in a specific order which is unchanged by the interval training program. First cHHb and cTOI start to level off probably as a consequence of the higher cerebral metabolic load to increase motor command to fast twitch fibers in response to the increase in power output and to account for fatigue in earlier recruited muscle fibers. The progressive recruitment of these less oxidative fibers enhances the occurrence of a metabolic acidosis, which has an impact on the ventilatory responses (RCP). As indicated above, this will reduce PaCO2 which will in turn affect cerebral blood flow (Paulsen et al., 1990) and thus, cO2Hb and ctotHb. The origin of the plateau in mHHb (and here also mTOI) and possible link with the cerebral oxygenation responses is currently unclear. Although originally the plateau in mHHb was considered as a leveling-off in O2 extraction, recent studies have shown an O2 extraction reserve following incremental ramp exercise when a blood flow occlusion is applied (Inglis et al., 2017; Iannetta et al., 2018). It is argued that the leveling off in mHHb near the end of exercise is related to locally released vasoactive compounds (H+, ATP, lactate, K+) triggering a local vasodilation and thus a redistribution of the blood flow. In order to establish a mechanistic link between muscle and cerebral oxygenation it would be interesting to directly measure the effect of training on blood gas (arterial O2 and CO2 pressures) and metabolite (e.g., lactate, H+) parameters.

Conclusion

The present study shows that aerobic interval training impacts both muscle and brain oxygenation, in concert with an increase in aerobic fitness. More specifically, the training program slowed the decrease in mTOI and increased the amplitude of mtotHb and mHHb, pointing at an improved O2 availability and O2 extraction capacity, respectively. Also the amplitude of cerebral O2Hb and totHb was increased, suggesting an enhanced brain perfusion at high intensities. The observation that the increase in O2peak correlated with the changes in the slope of mTOI and, amplitudes of both muscle HHb and totHb, but not to those of cerebral oxygenation, indicate that brain oxygenation per se will not be a primary limiting factor to incremental exercise tests.

Data Availability

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

KC, KV, JGB, and JB: study design. KC, KV, and JB: data collection. KC, KV, AM, VN, SP, and JB: data analysis. KC, KV, AM, VN, SP, JGB, and JB: data interpretation. KC, VN, SP, and JB: writing manuscript. AM, VN, SP, and JGB: revising manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.