Abstract
Background:
This research demonstrates how two distinct training modalities, high-intensity interval training (HIIT) and moderate-intensity continuous training (MICT), influence oxygen transport dynamics and microvascular remodeling.
Methods:
Twenty-five healthy sedentary men and women (median age 21 years) were randomly assigned to HIIT or MICT for 8 weeks. VO2max improvement was assessed in all participants. Non-invasive maximal cardiac output measurements (Qmax) were performed in 15 participants. Biopsies from vastus lateralis were obtained, cleared and immunolabeled for VE-cadherin and alpha-smooth muscle actin, in 10 subjects, to observe microvasculature architecture. A computational hemodynamic model was constructed to estimate muscle flow dynamics.
Results:
VO2max and Qmax increased significantly in both training groups, with a greater improvement for VO2max in HIIT that was accompanied by a significant increase in capillaries pericyte coverage. No formation of new capillaries nor anastomoses (angiogenesis) was detected in either group. Modelisation estimated higher shear stress during HIIT than MICT and pericyte recruitment was modelized to adapt to shear stress level limiting excessive capillary dilation.
Conclusion:
HIIT induces superior improvements in VO2max and distinct microvascular structural adaptations rather than angiogenesis. HIIT is thought to induce protective capillary adaptation, limiting dilation during maximal effort and improving oxygen diffusion.
Clinical trial registration:
https://clinicaltrials.gov/study/NCT07237854, identifier NCT07237854.
Introduction
Aerobic capacity, indexed by maximal oxygen consumption (VO2max), strongly predicts cardiovascular health and all-cause mortality (). Regular exercise provokes repeated perturbations of homeostasis that elicit coordinated central and peripheral adaptations, thereby increasing VO2max (). In this context, exercise training modalities such as moderate-intensity continuous training (MICT) and high-intensity interval training (HIIT) () have been successfully used for decades (). MICT is typically performed over a relatively extended period of time between 50% and 70% VO2max, corresponding to efforts below or near the first ventilatory threshold (VT1). VT1 reflects an individual’s aerobic capacity and marks a shift in substrate utilization, transitioning from predominantly aerobic pathways to increased engagement of anaerobic metabolism. HIIT, on the other hand, alternates between high-intensity exercise (above the VT1 threshold) and recovery phases. Ongoing debate, however, persists as to whether HIIT yields improvements in VO2max that are similar or superior to those of MICT and the exact mechanism of action remains elusive. Crucially, increases in exercise intensity produce larger hemodynamic forces within the microcirculation (arterioles, capillaries, and venules) of active muscles. Hemodynamic forces, represented by fluid shear stress (FSS) and vessel wall circumferential stretch (; ; ), are sensed by mechanoreceptors located between endothelial cells and vascular smooth muscle cells, triggering intracellular signalling cascades that may drive vascular remodeling (; ; ). Because direct in vivo quantification of microvascular hemodynamic forces within human skeletal muscle remains technically challenging, computational hemodynamic modeling may provide an important complementary framework for estimating how exercise-induced vascular remodeling influences flow distribution, shear stress, and impact the oxygen transport. The sequential nature and elevated intensities of HIIT may provoke greater, more dynamic, and more frequent hemodynamic stimuli, potentially leading to more pronounced vascular remodeling. According to the Fick principle, VO2 is the product of cardiac output and the arteriovenous oxygen difference; thus, improvements in aerobic capacity reflect both central enhancements (e.g., increased cardiac contractility and stroke volume) and peripheral changes (e.g., improved aerobic metabolism and enhanced muscle tissue perfusion) (; ). This study aims to elucidate how peripheral adaptations, particularly microvascular adaptive remodeling, contribute to improved cardiovascular efficiency and maximal oxygen uptake (VO2max) following different exercise regimens. By comparing HIIT and MICT, we sought to characterize intensity-dependent vascular adaptations to training and their potential role in optimizing tissue perfusion, oxygen delivery, and aerobic performance. To further strengthen this mechanistic framework, we incorporated computational hemodynamic modeling to evaluate the theoretical impact that structural vascular remodeling may have on muscle perfusion and oxygen transport. We hypothesized that: 1) both MICT and HIIT would improve VO2max and maximal cardiac output, with HIIT inducing greater improvements in VO2max; 2) HIIT would generate stronger hemodynamic stimuli, leading to more pronounced microvascular structural remodeling than MICT; and 3) these peripheral vascular adaptations would enhance muscle oxygen delivery and extraction during maximal exercise.
Methods
Study design
This prospective study was conducted between September 2022 and December 2024 within the Research Unit of Rehabilitation Sciences at Université libre de Bruxelles. All participants received detailed information about the study and provided written informed consent, which was approved by the Research Ethics Committee of the Brussels University Hospital (B4062021000227). Inclusion criteria required participants to be between 18 and 40 years old, to have a body mass index (BMI) below 30, and be classified as sedentary. Sedentary status was defined as self-reported physical activity levels below the World Health Organization guidelines, less than 150 minutes of moderate-intensity or 75 minutes of vigorous-intensity exercise per week (). Exclusion criteria included smoking, regular medication use, and any history of cardiovascular, respiratory, or metabolic diseases. The study protocol is illustrated in Figure 1. Twenty-five participants were recruited and underwent a cardiopulmonary exercise test (CPET) to determine their maximal oxygen uptake (VO2max). Of those 25 subjects, a subgroup of 15 participants completed a second exercise test after 24 hours and following a similar workload increment, with measurement of non-invasive cardiac output at rest and at maximal exercise workload (Qmax). VO2max and Qmax were subsequently used to derive the maximal arteriovenous difference (a-v O2 diff max) according to the Fick principle. Another subgroup of 10 participants underwent a resting muscle biopsy of the vastus lateralis (quadriceps). Muscle samples were collected, fixed and processed as detailed below (Muscle Biopsies section). All participants were then enrolled in an 8-week training program, 3 times a week, and randomly assigned to either a MICT group or HIIT group. At the end of the 8-week intervention, all baseline assessments were repeated. Cardiac output measurements and muscle biopsy data were used to construct a hemodynamic model (Hemodynamic model section) designed to link macro-level cardiovascular adaptations with micro-level vascular changes.
Figure 1
Training modalities
Cycling-based exercise training was performed under supervision three times per week. Throughout each session, heart rate was continuously monitored with visual feedback to facilitate adherence to prescribed intensity targets. The training workload was adjusted every five sessions, following individualized heart rate zones. Both exercise protocols were designed to elicit an equivalent hemodynamic load, defined as the product of training duration and targeted heart rate. All training sessions (HIIT and MICT) commenced with a standardized 3-minute warm-up at 50% HRmax. The 28-minute HIIT protocol consisted of seven 2-minute bouts at 90% HRmax, each separated by 2 minutes of moderate-intensity exercise at 65% HRmax. The 34-minute MICT protocol involved continuous exercise at 65% HRmax for 34 minutes. HRmax was determined from the initial cardiopulmonary exercise test (CPET).
Cardiopulmonary exercise test
CPET was conducted on a cycle ergometer (Ergoselect 100, Ergoline GmbH, Germany) using a one-minute workload increment protocol. Participants were instructed to maintain a pedaling cadence between 60 and 70 revolutions per minute throughout the test. Oxygen consumption (VO2), carbon dioxide production (VCO2), and ventilation (VE) were measured breath-by-breath via a face mask (COSMED, Rome, Italy). VO2max was defined as the highest VO2 value recorded over a 20-second interval at peak exercise. Additional methodological details, including test initiation procedures, criteria for maximal effort, and determination of the first ventilatory threshold, are provided in the Supplementary Methods (online appendix).
Non invasive measure of cardiac output
Cardiac output (Q) was measured non-invasively using the inert gas rebreathing method with the Innocor device (COSMED, Rome, Italy), which operates based on a single-alveolar lung model. Stroke volume was calculated by dividing the cardiac output (Q) by the heart rate (HR). Participants rebreathed a gas mixture containing 0.5% nitrous oxide (NO, blood-soluble) and 0.1% sulfur hexafluoride (SF6, blood-insoluble) diluted with ambient air. Pulmonary shunting is minimal in healthy individuals; pulmonary blood flow (PBF) was assumed to be equivalent to Q (PBF ≈ Q). Additional technical details are provided in the Supplementary Methods (online appendix).
Muscle biopsies
Resting muscle biopsies were obtained from the vastus lateralis of the quadriceps 14 days before the baseline CPET and 5 days after the post-training CPET to minimize any discomfort or pain that could interfere with exercise performance. To avoid sampling scar tissue and ensure anatomical consistency, the first biopsy was performed on the right leg and the second on the left. Muscle samples were processed using the PEGASUS protocol, a modified tissue-clearing technique optimized for immunolabeling and three-dimensional imaging. Vessel inner diameters (intima) were obtained from VE-cadherin staining of the endothelium, whereas external diameters (media) were determined using alpha-smooth muscle actin (alpha-SMA) staining of smooth muscle cells and pericytes. Confocal three dimensional imaging was performed using a Nikon AX R confocal microscope (Figure 2A). For each sample, z-stack acquisitions at 10× magnification were stitched together to generate wide-field composites that encompassed the full tissue volume. From these representative regions, areas containing arteriolar, capillary, and venular segments were selected for high-resolution imaging at 20× magnification. Vessel diameters were quantified transversely using NIS-Elements software, based on fluorescence-defined boundaries. Intima diameters were defined by VE-cadherin signal, while media diameters corresponded to α-SMA-positive structures. Measurements were performed independently on each fluorescence channel to ensure layer-specific morphometric accuracy. Capillary density was assessed in two anatomically distinct regions per sample using a volumetric approach that combined z-stack depth with a fixed surface area of 200 × 200 µm, with a mean imaging depth of 70 ± 6 µm at baseline and 83 ± 23 µm following the training intervention. A Full z-axis navigation allowed discrimination of individual vessels throughout tissue depth, minimizing errors related to vessel overlap or superimposition. Additional procedural details regarding biopsy collection, tissue preparation, immunostaining, and clearing steps are provided in the Supplementary Methods (online appendix).
Figure 2
Hemodynamic model
To improve our understanding of hemodynamic forces and how the microvascular system might structurally adapt to these forces, we developed a hemodynamic model representing the microvascular network. The hemodynamic model consisted of five dichotomous branches of the arteriole, resulting in 32 capillaries per arteriole. The venous side of the network consisted of five dual merges, giving one venule per arteriolar (Figure 2B).
(; ).
We distinguished three regimes of flows: Fåhræus-Lindqvist regime for vessel with a diameter above 8 µm, lubrication regime for vessels with a diameter below 7 µm and mixed regime for vessels with a diameter between 7 and 8 µm.
(; ; ; ; ).
The network was then deformed to account for the pressure of blood inside the vessels. For each generation, the final diameter was defined as the diameter providing a force balance between the pressure of the blood on the internal side and the elastic force of the walls and the intramuscular pressure on the external side.
(; ).
We assumed that the total flow in the muscles was 1 litre per minute at rest. For effort, we considered that the increase in cardiac output translated into an increased blood flow directed specifically toward the active muscles, taking into account that during intense exercise up to 80% of the cardiac output is redistributed to the working muscle (). The flow rate in the network was adjusted accordingly. Finally, to specifically assess the independent and combined effects of angiogenesis and vascular remodeling on network hemodynamics, simulations were performed under the post-training exercise flow conditions previously defined in the model (i.e., assuming that 80% of cardiac output was redistributed to active skeletal muscle). Three vascular architectures were compared. First, the baseline pre-training vascular network was simulated with the addition of 12.5% capillary angiogenesis to isolate the specific hemodynamic contribution of angiogenesis alone. This magnitude was selected based on its proximity to the ~15% capillary increase reported following aerobic exercise training (). Second, the experimentally derived post-training remodeled vascular network was simulated to evaluate the isolated effect of exercise-induced vascular remodeling. Third, the remodeled post-training network was combined with an additional 12.5% capillary angiogenesis to assess the cumulative effects of both remodeling and angiogenic expansion. This approach was designed to distinguish the relative contributions of angiogenesis and vascular remodeling to pressure-drop regulation under equivalent exercise flow demand. As a static experimental model, these simulations did not incorporate dynamic physiological mechanisms such as acute vasodilation, muscle contraction-induced vascular compression, enhanced cardiac pumping function, or skeletal muscle pump-mediated venous return. Across all simulations, vessel lengths were preserved so that pressure-drop differences specifically reflected network configuration and capillary expansion rather than additional geometric modifications. The length and diameter ratios between each vessel generation, the detailed descriptions of the three flow regimes, and the elastic vessel model are provided in the Supplementary Methods (online appendix).
Statistical analysis
All statistical analyses were performed using Jamovi (version 2.3.28). Data normality was assessed for each variable using the Shapiro-Wilk test. Results are expressed as mean ± standard deviation when normality was confirmed, or as median (interquartile range 25–75) when the distribution was non-normal. Paired Student’s t-tests were used to compare baseline and post-training (8 weeks) values. Two-way ANOVA was conducted with time and group as factors, and variance homogeneity was confirmed prior to analysis. Tukey’s post hoc test was applied when the ANOVA indicated a statistically significant effect. Statistical significance was accepted at p< 0.05.
Results
Twenty-five young, healthy, and sedentary individuals completed the 8-week training intervention (median age: 21 years; 15 women and 10 men), with 13 performing MICT and 12 HIIT. No significant differences were observed in baseline characteristics between the groups (Table 1). All participants who completed the intervention demonstrated full adherence to the prescribed training protocol. Aerobic capacity (VO2max) was similar between groups at baseline. As expected, VO2max increased in both groups following the training intervention, with a mean improvement of 4 ± 2 ml·kg-¹·min-¹ in the MICT group and 8 ± 5 ml·kg-¹·min-¹ in the HIIT group. The increase in VO2max was significantly greater in the HIIT group compared to MICT (p = 0.024). Consistently, maximal workload improved by 23 ± 13 watts in the MICT group and by 40 ± 18 watts in the HIIT group with a higher increase in HIIT (p=0.014). Both VO2 and workload at the first ventilatory threshold (VT1) showed significant improvements post-intervention in the MICT (p=0.032 and p=0.019, respectively) and HIIT groups (p<0.001 and p=0.007, respectively). These outcomes are summarized in Table 2.
Table 1
Characteristic | MICT (n = 13) | HIIT (n = 12) | p-value |
|---|---|---|---|
| Gender (F/M) | 9/4 | 6/6 | 0.347 |
| Age (years) | 21 (20 – 25) | 22 (20 – 27) | 0.722 |
| Weight (kg) | 70± 10 | 70 ± 11 | 0.896 |
| Height (cm) | 170 ± 12 | 171 ± 9 | 0.842 |
| BMI (kg/m²) | 24.5 ± 4.1 | 23.8 ± 2.2 | 0.595 |
Descriptive characteristics of participants.
BMI, body mass index. Statistics: Normality tests were performed by Shapiro-Wilk test. Student’s t-test was applied for normally distributed data (mean ± SD) and Mann-Whitney U test for non-parametric data (median (25th -75th percentile).
Table 2
| CPET values | MICT (N = 13) | HIIT (N = 12) | |||||
|---|---|---|---|---|---|---|---|
| Baseline | 8 weeks | Time | Baseline | 8 weeks | Time | Time* group | |
| Rest | |||||||
| VO2 (L. min-¹) | 0.42 ± 0.07 | 0.44 ± 0.12 | 0.582 | 0.39 ± 0.08 | 0.42 ± 0.07 | 0.091 | 0.278 |
| VO2 (mL·kg-¹·min-¹) | 6.1 ± 0.9 | 6.4 ± 1.4 | 0.532 | 5.6 ± 1.1 | 5.9 ± 0.8 | 0.145 | 0.167 |
| HR (beats·min-¹) | 93 ± 13 | 85 ± 13 | 0.082 | 84 ± 10 | 85 ± 13 | 0.743 | 0.113 |
| VT1 | |||||||
| VO2 (L. min-¹) | 1.75 ± 0.35 | 1.94 ± 0.51 | 0.052 | 1.82 ± 0.44 | 2.23 ± 0.42 | <0.001 | 0.038 |
| VO2 (mL·kg-¹·min-¹) | 25 ± 4 | 28 ± 6 | 0.032 | 27 ± 6 | 31 ± 4 | <0.001 | 0.538 |
| Workload (Watt) | 111 ± 27 | 129 ± 43 | 0.019 | 123 ± 43 | 149 ± 29 | 0.007 | 0.451 |
| HR (beats·min-¹) | 148 ± 19 | 150 ± 13 | 0.467 | 137 ± 24 | 144 ± 20 | 0.279 | 0.526 |
| Maximal | |||||||
| VO2 (L. min-¹) | 2.61 ± 0.64 | 2.92 ± 0.75 | 0.006 | 2.56 ± 0.45 | 3.18± 0.75 | <0.001 | 0.022 |
| VO2 (mL·kg-¹·min-¹) | 37 ± 7 | 41 ± 9 | 0.003 | 38 ± 6 | 46 ± 9 | <0.001 | 0.024 |
| Workload (Watt) | 187 ± 49 | 210 ± 57 | <0.001 | 198 ± 44 | 237 ± 51 | <0.001 | 0.014 |
| HR (beats·min-¹) | 184 ± 11 | 185 ± 11 | 0.567 | 177 ± 13 | 178 ± 16 | 0.585 | 0.978 |
| RER | 1.17 ± 0.07 | 1.16 ± 0.06 | 0.729 | 1.17 ± 0.08 | 1.16 ± 0.09 | 0.775 | 0.862 |
CPET data at before and following 8 weeks of training.
VO2, Oxygen uptake, HR,Heart rate, VT1, first ventilatory threshold, RER, Respiratory Exchange Ratio, Baseline: before 8-weeks training, 8 weeks: After 8-weeks training. Statistics: Normality was tested with a Shapiro-Wilk test. For pairwise comparisons, Student’s t-tests were used. ANOVA tests were performed to compare groups and Tukey post hoc tests were applied. Results are presented as mean ± SD.
Non-invasive determination of cardiac output
Out of the 25 participants, 15 performed non-invasive cardiac output (Q) measurements (8 assigned to MICT and 7 to HIIT). At peak exercise, Qmax increased similarly in both groups after 8 weeks of exercise training: from 16 ± 3 L/min to 19 ± 4 L/min (+ 13 ± 7%, p=0.002) in MICT and from 18 ± 3 L/min to 21 ± 2 L/min (+ 15 ± 9%, p=0.001) in HIIT. In parallel, stroke volume increased from 88 ± 13 to 101 ± 17mL (+ 15 ± 8%, p=0.01) in the MICT group and from 100 ± 15 to 119 ± 16mL (+ 19 ± 13%, p<0.001) in the HIIT group. No modification of maximal heart rate was observed. Within this subgroup, VO2max also increased, with a greater improvement in the HIIT group than in the MICT group (p = 0.049). Calculated peripheral O2 extraction remained unchanged in the MICT group (16 ± 3 to 16 ± 2 mlO2/100ml, p=0.984), while a tendency toward an increase was observed in the HIIT group (14 ± 2 to 16 ± 3 mlO2/100ml, p=0.066). However, a significant group-by-time interaction was detected (p = 0.041), suggesting a differential adaptation in peripheral oxygen extraction between training modalities. Qmax, VO2max and a-vO2 diff comparison are presented in Figure 3.
Figure 3
Resting muscle biopsies outcomes
In a subgroup of 10 participants equally divided between MICT (n=5) and HIIT (n=5), we performed vastus lateralis muscle biopsies before and after an 8-week training regimen. Arteriolar intima diameters exhibited non-significant enlargement from 20±3µm at baseline to 26 ± 4µm post-training in the MICT group (p=0.054), and from 22 ± 5µm to 27 ± 5µm in the HIIT group (p=0.171). Similarly, media diameters of arterioles showed a comparable dynamic, from 30 ± 4µm to 35 ± 7µm in MICT (p=0.181), and from 31 ± 9µm to 41 ± 12µm in HIIT (p=0.147). Capillary intima diameters increased following training, from 6.2 ± 0.6µm to 7.2 ± 0.5µm in MICT (p = 0.001) and from 6.3 ± 0.7µm to 7.4 ± 0.5µm in HIIT (p = 0.045). The media diameters of capillaries (representing pericytes coverage) showed no significant increase, from 7.1 ± 0.6µm to 7.7 ± 0.5µm in MICT (p = 0.06), and an eventual increase from 7.3 ± 0.6µm to 8.3 ± 0.5µm in HIIT (p = 0.047) (Figures 4E, J). No pro-angiogenic response was detected in both group, as three-dimensional analyses revealed no measurable expansion of the capillary network across conditions, including neither increased capillary number along vessel length nor increased anastomotic connections between capillaries. In the MICT group, mean capillary density was 3135 capillaries/mm³ at baseline and 2931 capillaries/mm³ post-training (p=0.231). Similarly, the HIIT group exhibited values of 3131 capillaries/mm³ at baseline and 3086 capillaries/mm³ following the intervention (p=0.805) (Figures 4D, I). Statistical analysis confirmed the absence of significant differences over time between groups (p = 0.499). Venular internal diameters demonstrated a pronounced increase, from 28 ± 8µm to 38 ± 10µm in MICT (p = 0.015), and from 23 ± 6µm to 36 ± 6µm in HIIT (p = 0.003). Correspondingly, venular external diameters expanded, from 36 ± 11µm to 45 ± 13µm in MICT (p = 0.009), and from 32 ± 10µm to 44 ± 7µm in HIIT (p = 0.021) (Figures 4G, L).
Figure 4
Figure 5
Hemodynamic model
Based on estimations from our vascular model capillary shear stress decreased in response to both training modalities after eight weeks. After training, and under training-intensity conditions (i.e., 65% HRmax for MICT and 90% HRmax for HIIT), capillary shear stress declined: from 7.45 Pa to 6.16 Pa in MICT and from 18.15 Pa to 10.47 Pa in HIIT. This indicates that, under the same conditions of exercise/blood flow intensity, shear stress before training was higher than after (Figure 5A). The relative reduction in shear stress during training yielded similar ratios with the increase in pericyte coverage for both modalities: 17.6 MPa·m-¹ for MICT and 18.35 MPa·m-¹ for HIIT (Figure 5B). Similarly to capillaries, venular shear stress decreased from baseline to post-training. Under the same exercise training intensity, venular shear stress dropped from 9.55 Pa to 4.69 Pa with MICT and from 26.90 Pa to 8.79 Pa with HIIT (Figure 5C). The ratios between changes in venular shear stress and corresponding alterations in venular diameter were consistent: 3.66 for MICT and 3.80 for HIIT (Figure 5D). These ratios reflect the relationship between the increase in venular diameter and the corresponding reduction in venular shear stress during training. Furthermore, we hypothesize that enhanced capillary pericyte coverage under HIIT conditions may contribute to mechanical stabilization of the capillary wall, potentially limiting deformation under elevated shear forces. This is illustrated in Figure 5E. Simulations of pressure-drop across distinct vascular architectures revealed marked differences in network hemodynamic efficiency depending on structural configuration. In MICT conditions (Figure 5F), the baseline vascular network supplemented with 12.5% angiogenesis exhibited a pressure drop of 140 mmHg, whereas the post-training remodeled vascular network demonstrated a substantially lower pressure drop of 74 mmHg. The addition of 12.5% angiogenesis to the remodeled network did not further alter this response, with pressure drop remaining at 74 mmHg. Similarly, under HIIT conditions (Figure 5G), the baseline vascular network supplemented with 12.5% angiogenesis exhibited a pressure drop of 252 mmHg, while the remodeled post-training vascular architecture reduced pressure drop to 98 mmHg. The subsequent addition of 12.5% angiogenesis to the remodeled network again produced negligible additional effect, with pressure drop measured at 99 mmHg.
Discussion
Our results demonstrate that both MICT and HIIT improved aerobic capacity, with HIIT eliciting greater increases in VO2max despite similar improvements in maximal cardiac output. At the microvascular level, both training modalities were associated with structural remodeling of pre-existing vessels, including capillary dilation and venular enlargement, without detectable capillary network expansion under our experimental conditions. Notably, only HIIT was associated with increased pericyte coverage around capillaries, whereas no such adaptation was observed following MICT, suggesting an intensity-dependent effect on capillary mural cell remodeling. Our hemodynamic model further suggests that these structural adaptations may enhance oxygen diffusion within muscle fibers. Collectively, these findings indicate that early exercise adaptation in young sedentary adults primarily involves remodeling of existing microvascular architecture rather than measurable angiogenesis.
Aerobic capacity
Our results demonstrate that both MICT and HIIT induce significant improvements in VO2max over an eight-week intervention in sedentary young adults. In parallel, the magnitude of improvement in VO2max was significantly greater following HIIT, reinforcing the growing body of evidence supporting its superior efficacy in enhancing aerobic capacity (; ). This differential effect appears to be particularly pronounced in healthy individuals aged 14 to 45 years () and influenced by baseline fitness level, while remaining largely independent of sex (). These effects were also observed by reporting greater improvement in VO2max following HIIT, compared to MICT under matched training volumes, suggesting that training intensity is a key driver of aerobic adaptations. It is noteworthy that the magnitude of VO2max improvement observed in our HIIT group exceeds that reported in the meta-analysis of that found a mean increase in VO2max of 4.9 ± 1.4 ml·kg-¹·min-¹ for MICT versus 5.5 ± 1.2 ml·kg-¹·min-¹ for HIIT what could partly be explained by the substantial heterogeneity in HIIT protocols across studies. Manipulation of the FITT principle which refers to frequency, intensity, time and type, has indeed emerged as a critical determinant of both acute and chronic, physiological and metabolic responses to interval training (). accordingly reported that protocols incorporating intervals of at least 2 minutes in duration, combined with a cumulative high-intensity workload of no less than 15 minutes per session over a period of 4 to 12 weeks, consistently produced greater improvements in aerobic fitness compared with MICT. Similarly, showed also that training with high intensity intervals lasting between 1 and 3 minutes, performed three times per week for a minimum of six weeks, optimize cardiorespiratory adaptations. The superior efficacy of longer HIIT bouts likely stems from sustained near-maximal intensity, enhancing central and peripheral cardiovascular stimuli (). The greater VO2max improvement observed in our study may reflect our protocol’s specific design of 14 minutes cumulative training of 2 minutes intervals at 90% of HRmax.
Peripheral oxygen extraction
In our subgroup of 15 participants VO2max improved across both training modalities with a parallel and similar increase in Qmax. These observations are in line with VO2max increase being primarily related to an increase in maximal oxygen transport (i.e. Qmax) to the exercising muscle, as described in several meta-analyses (; ). However, compared to the MICT group, VO2max increased slightly more in the HIIT group, suggesting a better oxygen extraction at maximal exercise. Although the arteriovenous oxygen difference (a-vO2 diff) remained generally stable among our participants, we observed a non-significant trend toward a greater increase in the HIIT group (14 ± 2 to 16 ± 3 mlO2/100ml, p=0.66) that was accompanied by a temporally divergent trajectory between training modalities (p=0.041). Oxygen extraction is considered to result from a combination of factors, including an intensified Bohr effect (), prolonged erythrocyte transit time, increased capillary density (), optimized vascular recruitment and flow distribution (), and greater oxidative capacity of muscle tissue (). In that context, greater exercise intensity during HIIT appears to amplify muscle intracellular signaling cascades, resulting in a superior increase in mitochondrial content compared to work-match MICT protocol (). Moreover, repeated exposure to elevated hemodynamic forces, which characterizes HIIT, appears to induce superior vascular adaptations compared to traditional endurance training (). A meta-analysis conducted by reported no significant improvement in a-vO2 diff following aerobic training (mixing MICT and HIIT protocols), although studies with longer duration or higher load of training induced a higher a-vO2 diff. However, as highlighted by Poole (), these reports should be interpreted with caution as of some of the studies included in this meta-analysis were subject to methodological limitations related to a limited use of catheter-based measurements and potential inaccuracies in VO2max and Qmax estimations. In contrast, a more recent meta-analysis by suggests that aerobic training, particularly when involving high-intensity intervals, can elicit modest but significant increases in a-vO2 diff. Finally, our results tend to confirm HIIT superiority in improving a-vO2 diff at maximal exercise that may stem from vascular remodeling processes, which could facilitate more efficient oxygen extraction as explained below.
Vascular remodeling
While macrovascular adaptations to exercise have been extensively documented, data regarding microcirculatory remodeling in humans remain limited. demonstrated a dose-dependent improvement in femoral artery flow-mediated dilation (FMD) and reduced arterial wall thickness following aerobic training. Meta-analyses further support a positive correlation between exercise intensity and endothelial function (), with HIIT yielding greater macrovascular benefits than MICT (). In contrast, little to no data exist describing structural adaptations at the microvascular level (arteriolar, capillary, and venular diameters), as much of the current understanding stems from animal models, and in vivo technical challenges restrict investigation within human quadriceps muscle. To our knowledge, our study is the first to assess microcirculation diameters and angiogenic adaptations in human skeletal muscle using advanced tissue clearing techniques combined with high-resolution confocal microscopy, which enables the three-dimensional visualization of entire human muscle samples, along with micron-level precision for measuring diameters from microcirculation and conducting longitudinal capillary quantification. Direct comparison with previous literature is therefore limited, as most studies on exercise-induced angiogenesis have relied on conventional two-dimensional histological sections with transverse slices stained with immunomarkers to estimate capillary density (capillaries per mm²) or capillary-to-fiber ratios. Other investigations have focused on the expression of angiogenesis-related biomarkers, especially the vascular endothelial growth factor (VEGF) (; ). From those studies, capillary density is estimated to range between ~300 and 600 capillaries per mm², depending on training status () with potential increase of 49.7 capillaries/mm² following endurance training (~15% improvement ()), and that exercise intensity is a key modulator of an angiogenetic response (). Interestingly, numerous studies () () () have identified shear stress as a key biomechanical stimulus for angiogenesis. While extensive literature supports exercise-induced angiogenic responses in human skeletal muscle, our study did not detect evidence of measurable capillary network expansion under our specific conditions. Using advanced tissue clearing, high-resolution confocal microscopy, and immunostaining, we achieved three-dimensional longitudinal visualization of vascular architecture with micron-level precision. This approach provides greater anatomical fidelity than conventional two-dimensional cross-sectional methods by enabling direct assessment of vessel continuity, branching, and anastomotic organization throughout tissue depth. In contrast, two-dimensional transverse sections may overestimate capillary number by counting vessel bifurcations or anastomoses as separate structures, potentially contributing to higher estimates of angiogenesis in previous studies. Moreover, we observed a general increase in capillary diameters after training, suggesting pre-existing capillaries dilatation rather than the formation of entirely new ones. The increased vascular caliber likely improves the detectability of pre-existing capillaries that existed before training but fell below the resolution threshold. Collectively, our findings indicate that, under our experimental conditions, microvascular remodeling was the predominant observable adaptation, without evidence of capillary angiogenesis after 8 weeks of training. Our participants were previously sedentary, and their vascular system was therefore likely more sensitive to exercise-induced hemodynamic stimuli, permitting relatively rapid early remodeling, leaving angiogenesis for a later adaptation. In contrast, trained individuals or athletes often already exhibit substantial cardiovascular and peripheral adaptations (). According to the shear stress set-point concept, repeated training may shift or expand the range of mechanical stimuli perceived as physiological, such that a habitual exercise stimulus may no longer be sufficient to induce further vascular remodeling (). Consequently, greater, more intense, or more specialized exercise stimuli may be required to exceed this expanded set-point and provoke additional vascular adaptations. Importantly, our computational hemodynamic model supports the hypothesis that remodeling of the baseline vascular network provides more effective pressure-drop regulation than physiologically relevant angiogenesis alone. Interrestingly, angiogenesis in combination with remodeling promotes effective pressure drop similar than when solely remodeling operates. However, these interpretation should be interpreted within the context of the model’s experimental and static design. Because dynamic physiological mechanisms could not be incorporated at the microvascular level, absolute pressure-drop values should be interpreted with caution. Specifically, this framework does not account for acute exercise-related factors such as muscle contraction-induced vascular compression, transient vasomotor responses, enhanced venous return through the skeletal muscle pump, or broader dynamic cardiovascular adjustments during exercise. Consequently, while the model provides valuable comparative insight into the relative structural hemodynamic effects of remodeling versus angiogenesis, it does not fully reproduce the complex dynamic physiology of exercising skeletal muscle.
Modelisation and interpretation
In humans at rest, typical shear stress values range between 0.5 and 5 Pa (). Numerous studies (; ; ) have introduced the concept of a “set point”, whereby endothelial cells drive vascular remodeling, adapting vessel diameter and luminal area to maintain shear stress within an optimal physiological range. A particularly compelling insight from the work of . is that the set point may differ across vessel types. In their study, comparison of endothelial cells from veins and lymphatic vessels highlighted the functional heterogeneity of endothelial responses. In our exercise-based model, structural adaptation following blood flow elevation and increased shear stress stimulated both capillaries and venules. Since the biopsies were obtained at rest, the adaptations observed in the capillary and venular segments represent permanent structural changes, in addition to functional vessel adaptation. However, a more pronounced arteriolar vasoreactivity, as evidenced by improved FMD in response to training (), might still exist without structural modifications. Supporting this interpretation, an experimental investigation of the rat mesenteric microvascular network reported that a ligature-induced rise in local pressure produced no diameter changes in small arterioles (10–20 µm), whereas larger arterioles (20–30 µm and >30 µm) exhibited clear structural adaptations. In contrast, venular segments displayed diameter alterations across all size classes (). Furthermore, our study demonstrated a significant increase in pericyte coverage at the capillary level following HIIT, whereas MICT did not produce a similar effect. This divergent pattern suggests that HIIT may trigger the activation of protective mechanisms to regulate microvascular architecture. Although both training modalities were matched for total exercise volume, exercise intensity appears to be a key determinant of the nature and magnitude of vascular remodeling. Within the shear stress set-point framework, both MICT and HIIT stimulated structural adaptation through elevations in hemodynamic forces (). However, the alternance of higher stimulation of HIIT may have more frequently and overall more intensively exceeded the physiological threshold promoting more frequent mechanical stimulation of greater magnitude than in MICT. This may therefore explain why both protocols induced remodeling of pre-existing vessels, while only HIIT promoted additional pericyte coverage. Furthermore, high-intensity intervals performed during HIIT are more representative of VO2max conditions and may promote microvascular adaptations optimized for maximal exercise demands. We hypothesize that the increased pericyte coverage observed under HIIT may help stabilize the microvessel structure in response to high mechanical stimuli, thereby preserving an optimal diameter for oxygen exchange (Figure 3F). From a rheological standpoint, vessel diameter has profound effects on blood viscosity and the spatial dynamics of red blood cells (RBCs) (). In vessels larger than 10 µm, multiple RBCs can travel side-by-side, interact and shift away from the vessel center, leading to peripheral displacement and increased apparent viscosity. In small vessels (diameter<10 µm), red blood cells (6–8 µm) deform and align in single file, moving centrally with a thin plasma layer separating them from the vessel wall, a behaviour described by the Fåhræus–Lindqvist effect. This microvascular flow configuration is especially favourable for oxygen diffusion, as it increases the surface area available for RBC gas exchange making microvascular oxygen transport tightly related to capillary geometry and RBC alignment. Studies have shown indeed that when RBCs travel centrally and in single file, oxygen extraction becomes more efficient due to minimized diffusion barriers (). We propose that the increased pericyte coverage observed in HIIT may help stabilize capillary diameter against pressure-induced deformation, thereby preserving this optimal flow configuration. This interpretation is further supported by our post-training analysis of the shear stress-to-pericyte coverage ratio yielded comparable ratios in both protocols inducing that pericyte coverage increases proportionally to the shear stress experienced during training. Such adaptation likely acts as a biomechanical safeguard, limiting excessive capillary deformation and helping maintain an optimal vessel diameter for oxygen exchange at specific training intensities. When extrapolated to maximal effort, this mechanism may carry important functional implications as when blood flow and shear stress surge, vessels conditioned through MICT, (i.e. lower shear stress) will exceed their tolerance, leading to over-dilation and reduced oxygen diffusing ability. In contrast, HIIT-trained vessels, conditioned to higher shear stress, are better equipped to withstand peak hemodynamic loads and preserve optimal flow geometry at maximal exercise intensity. Interestingly, veins adapted in accordance to the set point theory () and consistent with a homeostatic response to shear stress elevation, modulating luminal diameter to stabilize local hemodynamic stress. Indeed diameter increase following training was proportionate to hemodynamic forces increases in MICT or HIIT (Figures 3D, E). To note, mentioned shear stresses in veins ranging from 0.1 to 0.6 Pa in healthy individuals. In our cohort, resting shear stress levels decreased after training from 0.9 to 0.7 Pa in the MICT group and from 1.3 to 0.8 Pa in the HIIT group. These higher values are probably explained by the smaller diameters of venules compared to veins. Finally, no evidence of angiogenesis was observed following our exercise training. These findings may appear to contrast with the prevailing paradigm that exercise commonly induces angiogenesis. However, under our specific experimental conditions, our hemodynamic simulations suggest that the only remodeling of the pre-existing vascular architecture markedly enhanced hemodynamic efficiency. In both MICT and HIIT, the structurally remodeled network under post-training flow conditions substantially reduced pressure drop, while simulated physiologically relevant angiogenesis (12.5%) had negligible additional impact. Rather than excluding angiogenesis as a potential longer-term adaptation, these results suggest that, during early training in young sedentary adults, microvascular adaptation may primarily favor mechanotransduction-driven remodeling of existing vessels as a more efficient strategy to optimize perfusion and oxygen delivery. Our findings therefore support the concept that early exercise-induced vascular adaptation may be predominantly governed by structural optimization of the existing microvascular network before substantial angiogenic expansion becomes necessary.
Limitation
Vessel diameters were measured at rest, as it is impossible to fix perfused human biopsies. Second, blood samples were not collected, preventing direct assessment of hemoglobin concentration and limiting the precision of rheological estimations. Future studies could strengthen both functional and vascular remodeling interpretation by incorporating complementary approaches such as near-infrared spectroscopy (NIRS) to assess local muscle oxygenation and hemodynamic responses during exercise (). Third, cardiac output measurements and muscle biopsy data were collected from two distinct subject subgroups. Although both populations were matched for age, training status, and cardiovascular fitness, the absence of direct overlap introduces a methodological constraint. Consequently, while the findings offer valuable insight into microvessel behaviour, they may lack the granularity required to establish individual-level associations between flow dynamics and vascular remodeling.
In conclusion, our study demonstrates that high-intensity interval training elicits a more pronounced increase in VO2max compared to moderate-intensity continuous training, despite both protocols inducing similar enhancements in maximal cardiac output. We observed training intensity-dependent remodeling in pericyte capillaries, with no evidence of angiogenesis, to maintain optimal diameter and optimize muscle perfusion and oxygen extraction under the highest metabolic hemodynamic stress.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Research Ethics Committee of the Brussels University Hospital and had a clinical trial (identifier NCT07237854). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
EM: Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing, Investigation. CR: Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. CM: Investigation, Writing – review & editing. NB: Conceptualization, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing. GD: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Acknowledgments
We express our deep gratitude to all study participants for their time and commitment. The results of this study are presented clearly, honestly, and without fabrication, falsification, or inappropriate data manipulation. The results of the present study do not constitute endorsement by ACSM.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was used in the creation of this manuscript. Artificial intelligence–assisted tools were used during the preparation of this manuscript to improve the quality of the English language, particularly for refining scientific phrasing and enhancing clarity. All scientific content, interpretations, and conclusions were developed entirely by the authors.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2026.1840439/full#supplementary-material
References
1
Al-KhazrajiB. K.JacksonD. N.GoldmanD. (2016). A microvascular wall shear rate function derived from In Vivo hemodynamic and geometric parameters in continuously branching arterioles. Microcirculation23, 311–319. doi: 10.1111/micc.12279
2
AshorA. W.LaraJ.SiervoM.Celis-MoralesC.OggioniC.JakovljevicD. G.et al. (2015). Exercise modalities and endothelial function: a systematic review and dose–response meta-analysis of randomized controlled trials. Sports Med.45, 279–296. doi: 10.1007/s40279-014-0272-9
3
AstorinoT. A.CauserE.HazellT. J.ArhenB. B.GurdB. J. (2022). Change in central cardiovascular function in response to intense interval training: a systematic review and meta-analysis. Med. Sci. Sports Exerc54, 1991–2004. doi: 10.1249/mss.0000000000002993
4
BaeyensN.BandyopadhyayC.CoonB. G.YunS.SchwartzM. A. (2016). Endothelial fluid shear stress sensing in vascular health and disease. J. Clin. Invest.126 (3), 821–828. doi: 10.1172/JCI83083
5
BaeyensN.NicoliS.CoonB. G.RossT. D.Van Den DriesK.HanJ.et al. (2015). Vascular remodeling is governed by a VEGFR3-dependent fluid shear stress set point. eLife4, e04645. doi: 10.7554/eLife.04645
6
BaeyensN.SchwartzM. A. (2016). Biomechanics of vascular mechanosensation and remodeling. Mol. Biol. Cell. 27 (1), 7–11. doi: 10.1091/mbc.E14-11-1522
7
BalciG. A.AsH.OzkayaO.ColakogluM. (2022). Development potentials of commonly used high-intensity training strategies on central and peripheral components of maximal oxygen consumption. Respir. Physiol. Neurobiol.302, 103910. doi: 10.1016/j.resp.2022.103910
8
BartoliF.DebantM.Chuntharpursat-BonE.EvansE. L.MusialowskiK. E.ParsonageG.et al. (2022). Endothelial Piezo1 sustains muscle capillary density and contributes to physical activity. J. Clin. Invest.132, e141775. doi: 10.1172/JCI141775
9
BatemanR. M.SharpeM. D.EllisC. G. (2003). Bench-to-bedside review: Microvascular dysfunction in sepsis –hemodynamics, oxygen transport, and nitric oxide. Crit. Care7, 359–373. doi: 10.1186/cc2353. PubMed PMID: 12974969; PubMed Central PMCID: PMC270719.
10
BullF. C.Al-AnsariS. S.BiddleS.BorodulinK.BumanM. P.CardonG.et al. (2020). World Health Organization 2020 guidelines on physical activity and sedentary behaviour. Br. J. Sports Med.54, 1451–1462. doi: 10.1136/bjsports-2020-102955
11
ChebbiR. (2015). Dynamics of blood flow: modeling of the Fåhræus–Lindqvist effect. J. Biol. Phys.41, 313–326. doi: 10.1007/s10867-015-9376-1
12
CrowleyE.PowellC.CarsonB. P.DaviesW. (2022). The effect of exercise training intensity on VO2max in healthy adults: an overview of systematic reviews and meta-analyses. Transl. Sports Med.2022, 9310710. doi: 10.1155/2022/9310710
13
DengH.MinE.BaeyensN.CoonB. G.HuR.ZhuangZ. W.et al. (2021). Activation of Smad2/3 signaling by low fluid shear stress mediates artery inward remodeling. Proc. Natl. Acad. Sci.118, e2105339118. doi: 10.1073/pnas.2105339118
14
FiorenzaM.GliemannL.BrandtN.BangsboJ. (2020). Hormetic modulation of angiogenic factors by exercise-induced mechanical and metabolic stress in human skeletal muscle. 319 (4), 824–834. doi: 10.1152/ajpheart.00432.2020
15
FooteC. A.SoaresR. N.Ramirez‐PerezF. I.GhiaroneT.AroorA.Manrique‐AcevedoC.et al. (2022). “ Endothelial glycocalyx,” in Comprehensive physiology, 1st edn. Ed. PrakashY. S. (Hoboken, NJ, USA: Wiley), 3781–3811. doi: 10.1002/cphy.c210029
16
FungY. C.TsangW. C. O.PatitucciP. (1981). High-resolution data on the geometry of red blood cells. Biorheology18, 369–385. doi: 10.3233/BIR-1981-183-606
17
GreenD. J.HopmanM. T. E.PadillaJ.LaughlinM. H.ThijssenD. H. J. (2017). Vascular adaptation to exercise in humans: role of hemodynamic stimuli. Physiol. Rev.97, 495–528. doi: 10.1152/physrev.00014.2016
18
GuoZ.LiM.CaiJ.GongW.LiuY.LiuZ. (2023). Effect of high-intensity interval training vs. moderate-intensity continuous training on fat loss and cardiorespiratory fitness in the young and middle-aged a systematic review and meta-analysis. Int. J. Environ. Res. Public Health20, 4741. doi: 10.3390/ijerph20064741
19
HaasT. L.NwadoziE. (2015). Regulation of skeletal muscle capillary growth in exercise and disease. Appl. Physiol. Nutr. Metab.40, 1221–1232. doi: 10.1139/apnm-2015-0336
20
HellstenY.GliemannL. (2024). Peripheral limitations for performance: muscle capillarization. Scand. J. Med. Sci. Sports34, e14442. doi: 10.1111/sms.14442
21
HellstenY.NybergM. (2015). Cardiovascular adaptations to exercise training. Compr. Physiol., 1–32. doi: 10.1002/cphy.c140080
22
HendrickseP.DegensH. (2019). The role of the microcirculation in muscle function and plasticity. J. Muscle Res. Cell Motil.40, 127–140. doi: 10.1007/s10974-019-09520-2
23
HoierB.HellstenY. (2014). Exercise-induced capillary growth in human skeletal muscle and the dynamics of VEGF. Microcirculation21, 301–314. doi: 10.1111/micc.12117
24
HonigC. R.FeldsteinM. L.FriersonJ. L. (1977). Capillary lengths, anastomoses, and estimated capillary transit times in skeletal muscle. Am. J. Physiol. Heart Circ. Physiol.233, H122–H129. doi: 10.1152/ajpheart.1977.233.1.H122
25
HuoY.KassabG. S. (2012). Intraspecific scaling laws of vascular trees. J. R. Soc Interface9, 190–200. doi: 10.1098/rsif.2011.0270
26
JensenL.BangsboJ.HellstenY. (2004). Effect of high intensity training on capillarization and presence of angiogenic factors in human skeletal muscle. J. Physiol.557, 571–582. doi: 10.1113/jphysiol.2003.057711
27
JonesS.ChiesaS. T.ChaturvediN.HughesA. D. (2016). Recent developments in near-infrared spectroscopy (NIRS) for the assessment of local skeletal muscle microvascular function and capacity to utilise oxygen. Artery Res.16, 25. doi: 10.1016/j.artres.2016.09.001
28
JoynerM. J.CaseyD. P. (2015). Regulation of increased blood flow (hyperemia) to muscles during exercise: a hierarchy of competing physiological needs. Physiol. Rev.95, 549–601. doi: 10.1152/physrev.00035.2013
29
LiuY.ChristensenP. M.HellstenY.GliemannL. (2022a). Effects of exercise training intensity and duration on skeletal muscle capillarization in healthy subjects: a meta-analysis. Med. Sci. Sports Exerc54, 1714–1728. doi: 10.1249/MSS.0000000000002955
30
LiuY.ChristensenP. M.HellstenY.GliemannL. (2022b). Effects of exercise training intensity and duration on skeletal muscle capillarization in healthy subjects: a meta-analysis. Med. Sci. Sports Exerc54, 1714–1728. doi: 10.1249/MSS.0000000000002955
31
LockM.YousefI.McFaddenB.MansoorH.TownsendN. (2024). Cardiorespiratory fitness and performance adaptations to high-intensity interval training: are there differences between men and women? A systematic review with meta-analyses. Sports Med.54, 127–167. doi: 10.1007/s40279-023-01914-0
32
LoweG. D. O. (2019). Clinical blood rheology: volume 2 (Boca Raton: CRC Press), 254.
33
MacInnisM. J.GibalaM. J. (2017). Physiological adaptations to interval training and the role of exercise intensity. J. Physiol.595, 2915–2930. doi: 10.1113/JP273196
34
MairbäurlH. (2013). Red blood cells in sports: effects of exercise and training on oxygen supply by red blood cells. Front. Physiol.4, 332. doi: 10.3389/fphys.2013.00332
35
MaufroyE.BaeyensN.DeboeckG. (2026). Mechanoreceptor-mediated adaptations to physical exercise: from acute responses to long-term training effects. Cell. Mol. Life Sci.83, 180. doi: 10.1007/s00018-026-06140-1
36
MilanovićZ.SporišG.WestonM. (2015). Effectiveness of high-intensity interval training (HIT) and continuous endurance training for VO2max improvements: a systematic review and meta-analysis of controlled trials. Sports Med.45, 1469–1481. doi: 10.1007/s40279-015-0365-0
37
MølmenK. S.AlmquistN. W.Skattebo∈Ø. (2025). Effects of exercise training on mitochondrial and capillary growth in human skeletal muscle: a systematic review and meta-regression. Sports Med. Auckl Nz55, 115–144. doi: 10.1007/s40279-024-02120-2
38
MonteroD.Diaz-CañestroC.LundbyC. (2015a). Endurance training and V˙O2max: role of maximal cardiac output and oxygen extraction. Med. Sci. Sports Exerc47, 2024–2033. doi: 10.1249/MSS.0000000000000640
39
MonteroD.Diaz-CañestroC.LundbyC. (2015b). Endurance training and V˙O2max: role of maximal cardiac output and oxygen extraction. Med. Sci. Sports Exerc47, 2024–2033. doi: 10.1249/mss.0000000000000640
40
MyrkosA.SmiliosI.ZafeiridisA.KokkinouM. E.TzoumanisA.DoudaH. (2023). Aerobic adaptations following two iso-effort training programs: an intense continuous and a high-intensity interval. Appl. Physiol. Nutr. Metab.48, 583–594. doi: 10.1139/apnm-2022-0309
41
PooleD. C.KelleyG. A.MuschT. I. (2016). Training increases muscle O2 diffusing capacity intrinsic to the elevated V˙O2max. Med. Sci. Sports Exerc48, 762–763. doi: 10.1249/MSS.0000000000000853
42
PooleD. C.PittmanR. N.MuschT. I.ØstergaardL. (2020). August Krogh’s theory of muscle microvascular control and oxygen delivery: a paradigm shift based on new data. J. Physiol.598, 4473–4507. doi: 10.1113/JP279223
43
RamosJ. S.DalleckL. C.TjonnaA. E.BeethamK. S.CoombesJ. S. (2015). The impact of high-intensity interval training versus moderate-intensity continuous training on vascular function: a systematic review and meta-analysis. Sports Med.45, 679–692. doi: 10.1007/s40279-015-0321-z
44
RossM.KarglC. K.FergusonR.GavinT. P.HellstenY. (2023). Exercise-induced skeletal muscle angiogenesis: impact of age, sex, angiocrines and cellular mediators. Eur. J. Appl. Physiol.123, 1415–1432. doi: 10.1007/s00421-022-05128-6
45
ScribbansT. D.VecseyS.HankinsonP. B.FosterW. S.GurdB. J. (2016). The effect of training intensity on VO2max in young healthy adults: a meta-regression and meta-analysis. Int. J. Exerc Sci.9, 230–247. doi: 10.70252/HHBR9374
46
SecombT. W.SkalakR.ÖzkayaN.GrossJ. F. (1986). Flow of axisymmetric red blood cells in narrow capillaries. J. Fluid Mech.163, 405–423. doi: 10.1017/S0022112086002355
47
SkatteboØ.BjerringA. W.AuensenM.SarvariS. I.CummingK. T.CapelliC.et al. (2020). Blood volume expansion does not explain the increase in peak oxygen uptake induced by 10 weeks of endurance training. Eur. J. Appl. Physiol.120, 985–999. doi: 10.1007/s00421-020-04336-2
48
StavrinouP. S.AstorinoT. A.GiannakiC. D.AphamisG.BogdanisG. C. (2025). Customizing intense interval exercise training prescription using the “frequency, intensity, time, and type of exercise” (FITT) principle. Front. Physiol.16, 1553846. doi: 10.3389/fphys.2025.1553846
49
SunD.HuangA.KollerA.KaleyG. (1998). Adaptation of flow-induced dilation of arterioles to daily exercise. Microvasc. Res.56, 54–61. doi: 10.1006/mvre.1998.2083
50
ThijssenD. H. J.GreenD. J.HopmanM. T. E. (2011). Blood vessel remodeling and physical inactivity in humans. J. Appl. Physiol.111, 1836–1845. doi: 10.1152/japplphysiol.00394.2011
51
TiezziM.DengH.BaeyensN. (2022). Endothelial mechanosensing: A forgotten target to treat vascular remodeling in hypertension? Biochem. Pharmacol.206, 115290. doi: 10.1016/j.bcp.2022.115290
52
TiltonR. G.KiloC.WilliamsonJ. R. (1979). Pericyte-endothelial relationships in cardiac and skeletal muscle capillaries. Microvasc. Res.18 (3), 325–335. doi: 10.1016/0026-2862(79)90041-4
53
Van GiesonE. J.MurfeeW. L.SkalakT. C.PriceR. J. (2003). Enhanced smooth muscle cell coverage of microvessels exposed to increased hemodynamic stresses in vivo. Circ. Res. 92 (8), 929–936. doi: 10.1161/01.RES.0000068377.01063.79
54
WenD.UteschT.WuJ.RobertsonS.LiuJ.HuG.et al. (2019). Effects of different protocols of high intensity interval training for VO2max improvements in adults: a meta-analysis of randomised controlled trials. J. Sci. Med. Sport22, 941–947. doi: 10.1016/j.jsams.2019.01.013
55
ZhangX.GaoF. (2021). Exercise improves vascular health: role of mitochondria. Free Radic. Biol. Med.177, 347–359. doi: 10.1016/j.freeradbiomed.2021.11.002
Summary
Keywords
aerobic capacity, muscle, oxygen diffusion, training, vascularisation, VO2max
Citation
Maufroy E, Rigaut C, Maufroy C, Baeyens N and Deboeck G (2026) Vascular remodeling enhances high-flow muscle oxygen delivery following aerobic exercise training. Front. Physiol. 17:1840439. doi: 10.3389/fphys.2026.1840439
Received
27 March 2026
Revised
20 May 2026
Accepted
29 May 2026
Published
01 July 2026
Volume
17 - 2026
Edited by
Irena Levitan, University of Illinois Chicago, Chicago, United States
Reviewed by
Mark Daniel Ross, Heriot-Watt University, United KingdomMinyoung Kwak, University of Kentucky, Lexington, United States
Updates
Copyright
© 2026 Maufroy, Rigaut, Maufroy, Baeyens and Deboeck.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: G. Deboeck, gael.deboeck@ulb.be
†ORCID: E. Maufroy, orcid.org/0009-0006-4777-4144; C. Rigaut, orcid.org/0000-0003-4999-8601; C. Maufroy, orcid.org/0009-0009-0458-0480; N. Baeyens, orcid.org/0000-0001-9893-3041; G. Deboeck, orcid.org/0000-0002-4505-7815
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.