Abstract
Introduction:
Skeletal muscles damaged by exercise exhibit disturbed energy metabolism and microvascular function for several days. However, it remains unclear whether these local changes might affect systemic cardiovascular responses to exercise. The present study aimed to investigate whether damaged muscles show changes in energy metabolism and oxygenation that influence systemic cardiovascular responses to exercise and post-exercise circulatory occlusion (PECO).
Methods:
A novel multi-parametric magnetic resonance imaging and spectroscopy approach was applied. Twelve healthy male participants completed assessments before and 48 h after 40 min of downhill running (20% decline). The assessments included muscle function, inflammation, and multi-parametric imaging at rest, exercise, and post-exercise occlusion using 31P spectroscopy, 1H- and muscle blood oxygen level-dependent imaging, and 23Na+ imaging to assess phosphate metabolism, oxygenation, and sodium disturbances. The mean arterial pressure (MAP) and heart rate (HR) were recorded throughout the MRI sequences.
Results:
Forty-eight hours after downhill running, muscle inflammation and Na+ disturbances were evident (both p < 0.05). Muscle oxygenation was lower and inorganic phosphates were higher during exercise and PECO than at baseline (both p < 0.05). However, MAP and HR during exercise and PECO remained unchanged at 48 h compared with baseline.
Conclusion:
Our multi-parametric MRI approach provides new insights into the local effects of muscle damage on energy metabolism, oxygenation, and Na+. Despite these local metabolic and microvascular disturbances, systemic cardiovascular responses, as indicated by MAP and HR, remained unchanged. These new findings suggest a dissociation between muscle metabolites, oxygenation, and the cardiovascular response to exercise and PECO 48 h after damaging exercise.
Introduction
Supporting exercising skeletal muscles by increasing blood flow while reducing metabolite accumulation is crucial for minimizing fatigue and preventing exercise intolerance (). This requires a redistribution of the cardiac output, primarily to the exercising muscle mass, which is regulated in part by muscle nerve afferents and the central nervous system (; ). In this context, the role of feedback is pivotal and driven by two types of skeletal muscle nerve afferents: myelinated group III (Aδ) and unmyelinated group IV (C-fiber). These afferents detect mechanical and metabolic stimuli, respectively, to match cardiovascular adjustments with local muscle metabolic demands during exercise (; ).
Afferent dysfunction has been reported in chronic inflammatory conditions, but it can also occur after acute, localized muscle damage (; ). However, the underlying mechanisms causing vascular dysfunction remain unclear. Previous studies have used exercise-induced muscle damage (EIMD) as a controlled inferential model to investigate the mechanisms underlying afferent-mediated changes in healthy populations (; ; ; ).
EIMD is characterized by localized inflammation (), muscular weakness, and soreness, and it usually peaks at 24 h–48 h post-exercise (). Additionally, EIMD is associated with increased phosphate muscle metabolism, ion perturbation (), and impaired microvascular function (; ; ; ). However, although these physiological alterations are well-known, the interaction between muscle afferents, metabolism, and cardiovascular regulation in response to exercise and post-exercise circulatory occlusion (PECO) post-EIMD remains poorly understood. EIMD alters local metabolic and inflammatory responses, which may, in turn, affect the muscle metaboreflex. This reflex helps match the cardiovascular output to metabolic demand. However, the extent to which EIMD disrupts this coupling remains unclear. Previous studies suggest that EIMD alters afferent sensitivity, particularly in mechano- and nociceptive receptors, which may lead to exaggerated cardiovascular responses to passive leg movements and reduced hemodynamic regulation (). However, some studies report intact blood pressure (BP) responses in the presence of central command activity (), while others demonstrate that the BP response during isometric exercise after eccentric exercise was enhanced, even when the workload on the damaged arm was reduced to match the perceived effort (i.e., central command) with the control arm (). On the other hand, an increase in the central command to the extent required to compensate for eccentric-induced muscle weakness does not affect the ventilatory response during dynamic exercises (). Other studies have also reported reduced or intact metaboreflex responses during PECO at 24 h and 48 h post-EIMD (; ). These discrepancies may arise due to methodological differences across the studies that cause different levels of muscle oxygenation and energy metabolic changes, both of which could significantly influence BP regulation during exercise and PECO (). Potential methodological differences across studies may be related to the exercise load used to assess cardiovascular responses during exercise and PECO, both before and after EIMD. Some studies have utilized the same absolute workload pre- and post-EIMD, while others have reduced the post-EIMD workload to account for increased central command because of muscle soreness or weakness. However, adjusting workload intensity can alter muscle energetics and the metabolic environment—particularly oxygenation and phosphate metabolism—both at baseline and following EIMD. These changes may vary between individuals and influence cardiovascular responses independently of afferent feedback.
To address this, we chose to maintain the same absolute workload across conditions while controlling for metabolic activation by matching the phosphocreatine (PCr) depletion levels at baseline. This approach allowed us to standardize exercise intensity and contextualize cardiovascular responses to exercise and PECO in relation to muscle metabolic changes induced by EIMD.
Therefore, the present study aimed to investigate how post-EIMD changes in muscle energy metabolism and oxygenation influence cardiovascular response during isometric exercise and PECO. We hypothesized that EIMD would increase muscle energy metabolism (inorganic phosphate) and reduce oxygenation, which would lead to a concomitant increase in the cardiovascular responses during exercise and PECO. By directly linking these variables to cardiovascular responses during sustained exercise in weakened muscles and the presence of inflammation, this study provides new insights into pathophysiological mechanisms leading to altered cardiovascular function.
Methods
Participants
The study received ethical approval from the Faculty of Science and Engineering Research Ethics and Governance Committee (reference number: 48097) and was conducted in accordance with the Declaration of Helsinki, except for pre-registration. Twelve healthy young men (age, 25.2 ± 5.0 years; mass, 76.3 ± 11.4 kg; and height, 178.0 ± 7.6 cm) completed the study. One participant dropped out of the study due to sickness, and only six participants completed the 31P assessment because of time constraints and scanner-associated costs.
All participants abstained from caffeine and food for 2 hours prior to participation and avoided alcohol intake and intense exercise for 2 days prior to participation in the study. The exclusion criteria included the use of non-steroidal anti-inflammatory medication and the presence of any injury or medical conditions that prevented downhill running exercise. Participants were also excluded from MRI if they had a cardiac pacemaker or any other standard contraindications. Due to the nature of the study design, participants were not randomized, and blinding was not implemented.
Experimental design
Participants attended the research laboratory for three visits. A familiarization session was conducted, during which a verbal explanation of the study procedures was provided, written informed consent was obtained, participants were familiarized with the MRI facility, and knee-extension assessments and appointments were scheduled for the first (pre-EIMD) and follow-up experimental sessions. Follow-up took place 48 h after baseline (post-EIMD). Pre- and post-EIMD assessments followed the same procedures. The EIMD exercise protocol was completed only once, at the end of the first session.
Questionnaires and soreness assessment
A physical activity readiness questionnaire (PAR-Q) was completed, and standing height and body mass were measured. Perceived muscle soreness of the knee extensors was measured using a visual analog scale as the participants held a squat with the knees bent at 90˚ (VASSQ) (). The assessment was performed by asking participants to mark an “X” on a 10-cm scale to indicate the level of soreness: with 0 representing no muscle soreness and 10 representing muscles so sore that movement is impossible (; ).
Maximal voluntary contraction assessment
The participants sat upright on a Kineo Multistation Dynamometer (Globus Italia, Treviso, Italy) with the hips and knees flexed at 90˚ and straps secured around the waist to minimize extraneous movements. Single-leg maximal voluntary knee extension isometric contraction (MVC) was tested on the dominant limb, and real-time feedback of force was recorded. The dominant leg was determined via self-report, based on the preferred leg used for kicking a ball. A short warm-up was permitted consisting of 10 brief isometric contractions, increasing progressively from approximately 50% to 80% maximal effort. After 1 min of rest, participants performed three MVCs, each separated by 1 min of rest. Each MVC trial lasted 3 s–5 s, with verbal encouragement provided to ensure the maximal effort. The highest external force value was accepted as each participant’s MVC (). This procedure was used to assess the initial load to be used for the exercise protocol that was later completed within the magnetic resonance (MR) scanner. For the MRI exercise load, a sandbag weighing 3 kg–5 kg was strapped firmly close to the ankle malleoli. The load was individualized depending on the initial MVC, allowing similar metabolism activation across participants, with a 20% ± 2% PCr reduction from baseline to the end of exercise. Standardization of PCr depletion between the participants was achieved by altering the ankle weight to a fraction of their MVC such that m = 0.131 MVC/(gּּּּL), where m is the ankle weight in kg, MVC is the torque in Nm, g is 9.81 ms−2, and L is the leg length in meters, which is defined as the distance from the lateral femoral condyle to the malleolus (). Specifically, participants with 177 Nm–204 Nm MVC used a 3-kg sandbag, participants with 204 Nm–235 Nm MVC used a 3.5-kg sandbag, participants with 235 Nm–267 Nm MVC used a 4-kg sandbag, participants with 267 Nm–298 Nm MVC used a 4.5-kg sandbag, and participants with 298 Nm–330 Nm MVC used a 5-kg sandbag. The mean ± SD sandbag load was 4 ± 1 kg. Pre- and post-EIMD loads were kept the same for each individual participant.
Quantitative multi-parametric MRI protocol
A Siemens MAGNETOM Vida 3 Tesla MRI system (Siemens Healthcare GmbH, Erlangen, Germany, housed in the university’s Wolfson ACTIVE Laboratory) was used for imaging/spectroscopy assessments of the thigh muscles. Participants were positioned feet-first supine, and all MRI sequences were centered on the mid-thigh, spanning the superior border of the patella to approximately 10 cm distal to the greater trochanter. The participants were provided with hearing protection, and the total time spent in the scanner, including the setup, was approximately 80 min. An integrated RF body coil was used for 1H signal transmission. A receive-only 16-channel 1H flex coil covered both thighs for structural imaging of the quadriceps. T1-weighted images (TR 800 ms, TE 12 ms) were acquired across 28 slices (0.5 × 0.5 × 7 mm voxels). T2 mapping utilized a 2D multi-spin-multi-echo (MSME) sequence with 17 equidistant echoes (TE1/ΔTE = 10 ms; TR = 3,000 ms; slice thickness = 5 mm; slice gap = 10 mm) (). Thereafter, the 1H flex coil was removed, and a dual-tuned single-channel 1H/31P surface coil (31P/1H Flex Coil 3T, RAPID Biomedical GmbH) was placed on the anterior thigh of the dominant leg. The knee was raised to rest on a firm bolster 14.5 cm high (thus, flexing the hip to approximately 20˚ and placing the RF coil at the MRI isocenter to maximize the signal-to-noise ratio) (). Sandbags (weight matching the relative load MVC assessment results) were strapped firmly around the distal tibia close to the malleoli (adapted from ). An 1H-based T2 localizer was acquired (3 slice packages; 15 slices; 256 mm2 field of view (FOV); 6 mm slice thickness; repetition time (TR) 606 ms; echo time (TE) 122 ms; number of averages (NA) 1; flip angle (FA) 150°; 60% phase resolution; 6/8 phase partial Fourier; 592 Hz/px; echo spacing 4.06 ms; turbo factor 154; RF pulse type, fast; and gradient mode, fast). This tri-planar localizer was used to position a cuboid adjustment window covering the entire thigh (the size and position covering the visible area of the individual thigh). Manual frequency, power, and 3D shimming adjustments were completed to achieve a 1H water peak with full width at half maximum (FWHM) of 19.9 ± 1.0 Hz and a 31P PCr peak with FWHM of 16.5 ± 0.97. 3D shimming utilized a multiple-gradient echo field mapping approach. The PCr resonance was set to 0 ppm. This was confirmed by a short 31P FID-based acquisition (TR 4,000 ms; NA 32; FA 90°; bandwidth (BW) 4,000 Hz; acquisition duration 512 ms; and spectral points 2,048). Scan parameters were as follows: NA = 32; TR = 12,000 ms; BW = 3,000 Hz; data points = 2,048; and acquisition duration = 692 ms. 31P spectra were acquired across the four workloads of the exercise protocol, including a 2-min of resting period, a 3-min submaximal sustained isometric contraction of the knee extensors, maintaining the shank with the ankle load parallel to the scanning bed, 2 min of PECO, and a 3-min resting recovery. At the beginning and end of the sustained knee extension contractions, participants were asked to provide their rating of perceived exertion (RPE) (). Following these procedures, the participants rested for 10 min–15 min outside the MRI scanner, while a flexible coil (23Na/1H Flex Coil 3T, RAPID Biomedical GmbH) was placed on the anterior thigh of the dominant leg. The participants then returned to the scanner to repeat the knee extension exercise. 2D radial FLASH-based 23Na images were captured pre- and post-exercise (TR 50 ms; TE 2.3 ms; FOV 250 mm2; matrix size 642; flip angle 90o; and thickness 100 mm). Echo planar-based blood-oxygen-level-dependent muscle oxygenation (mBOLD) acquisition was performed every 1 s (TR = 2 s) for a total of 10 min to capture rest, sustained isometric knee extension exercise, PECO, and recovery workloads (as described above).
Cardiovascular assessment and post-exercise cuff occlusion protocol
During the MRI exercise protocols, participants wore a beat-by-beat, non-invasive MRI-compatible blood pressure plethysmograph (Biopac System, United Kingdom) to monitor the mean arterial blood pressure (BP) and heart rate (HR). MAP and HR were collected continuously during the 31P MR spectroscopy and BOLD acquisition protocols. PECO used a rapid inflatable cuff (13 cm wide; E. Hokanson Inc., Bellevue, WA 98005, United States) around the proximal thigh. The cuff was inflated with compressed air to 250 mmHg–270 mmHg 15 s before the end of the 3-min sustained leg extension contraction. After 2 min, the cuff was instantly deflated. Arterial and venous occlusion was confirmed by the absence of a posterior tibial artery pulse ().
Exercise-induced muscle damage protocol
The EIMD protocol consisted of 40-min downhill running (DHR) at 20% decline on a treadmill (HP Cosmos Saturn 300/100, Nussdorf-Traunstein, Germany). The treadmill speed was set at 10 km/h, and a 5-min warm-up was provided as the treadmill gradient was gradually changed from 0% to 20% decline (). The treadmill speed, HR, mechanical power, and RPE were measured at 5 min intervals throughout the 40-min run. The treadmill speed was reduced to 8 Km/h if the participants reported RPE of 18 or asked for a slower running speed. A 2-min cool-down period at 0% decline at a self-selected speed was provided after the 40-min run.
Data analysis and MRI data processing
All MRI data (QCSA, T2 maps, muscle oxygenation, 31P and 23Na) were analyzed in MATLAB (2020b; MathWorks, Natick, MA, United States) using software routines developed in-house and available upon reasonable request. EIMD outcomes and cardiovascular responses to exercise were reported in Excel (Microsoft, Redmond, Washington, United States). Cardiovascular, muscle metabolism, and oxygenation metrics were analyzed for each given workload. Changes in sodium concentration were calculated by subtracting the resting values from post-exercise values for each time-point (pre-EIMD and post-EIMD).
For imaging data, a region of interest (ROI) was placed over the thickest portion of the quadriceps muscles (i.e., vastus lateralis, vastus intermedium, rectus femoris, and vastus medialis), avoiding subcutaneous fat and bone. The ROI mask was used to evaluate the quadriceps cross-sectional area (QCSA) on the T1-weighted image. Using the MSME images, quantitative muscle-water T2 maps were reconstructed based on a tri-exponential fitting procedure (; ). Muscle-water T2 values are reported as the mean from an ROI covering the interior of the vastus lateralis muscle. Muscle oxygenation (mBOLD) was analyzed from the imaging slice with the maximum QCSA. The vastus lateralis was used as the reference, motion artifacts were removed (using image registration), and the oxygenation signal was recorded and expressed as fold changes from rest. 23Na datasets were reconstructed offline using a non-uniform fast Fourier transform and calculated across. Before co-registration to the 1H images (tri-pilot images), a correction for gradient nonlinearities was performed. Then, a region-based partial volume correction using the tissue masks was applied to the 23Na images. Finally, sodium total concentration values were determined using the corrected signal intensities based on three external references (NaCl at 10/20/30 mmol). For 31P spectral data, raw frequency data underwent 5 Hz line broadening (Gaussian filter) to improve SNR. Spectra were manually phased (zeroth- and first-order correction). Following the Fourier transform, the baseline was fitted to a fourth-order polynomial and removed (elevated due to phospholipid contributions). The resultant spectral peaks were assigned to the following 31P resonances: phosphomonoesters (PMEs), inorganic phosphate (Pi; intracellular ∼4.9 ppm and extracellular ∼5.3 ppm), phosphodiesters (PDEs), phosphocreatine (PCr), and the three resonances of adenosine triphosphate (γ-ATP, α-ATP, and β-ATP). To quantify signal contributions from intracellular inorganic phosphate (Pi(i)), the spectra were further windowed to 4 ppm–6 ppm. Data were fitted using nonlinear least squares with a Levenberg–Marquardt algorithm to appropriate Gaussian/Lorentzian functions (to minimize the residual) for Pi(i) and Pi(e) with chemical shift, amplitude, and FWHM floating variables. Fitting parameters were used to isolate the signal contribution from intracellular Pi. The area under the curve (AUC) was extracted as a function of phase and time-point (averaged across the collected spectra). All data were normalized to the resting phase across the functional task.
Statistical analysis
Power analysis for sample size calculation was determined on the difference between pre- and post-changes in muscle oxygenation, as the primary outcome, reported by . Results from power calculation analysis using two-way repeated measures analysis of variance (ANOVA) for two groups were as follows: λ = 12.20, F = 4.9, and actual power = 0.87. This resulted in an estimated total sample size of n = 12 participants. While some secondary outcomes had smaller sample sizes (e.g., 31P MRS, n = 6), Cohen’s d effect sizes are reported to aid interpretation. Statistical analysis was carried out in GraphPad Prism (v.10.1 GraphPad software, San Diego, California United States). The normal distribution of the data was assessed using a Shapiro–Wilk test. If the sphericity assumption was violated, the Greenhouse–Geisser correction coefficient was reported. A paired t-test was performed for pre- and post-EIMD muscle soreness and neuromuscular outcomes (i.e., VASSQ, MVC, QCSA, T2 maps, resting [Na+], and delta [ΔNa+] concentrations). A two-way repeated measures ANOVA was performed for muscle oxygenation (mBOLD), Pi, Pi(i), Pi(e), MAP, and HR for workload (rest, contraction, PECO, and recovery) and time (pre-EIMD vs. post-EIMD). For all these variables, Cohen’s d effect sizes (ESs) were calculated as delta differences between delta changes from rest to exercise (Δexercise) and rest to PECO (ΔPECO) between pre- and post-EIMD measures by the pooled SD for each delta change. Cohen’s d effect size thresholds were considered at 0.2, 0.5, and 0.8 for small, moderate, and large effect sizes, respectively ().
Results
Neuromuscular function, soreness, and sodium concentrations following EIMD
Neuromuscular function, soreness, inflammation, and sodium concentrations before and after EIMD are summarized in Figure 1. The downhill running protocol significantly impaired neuromuscular function, with a decrease in MVC torque (pre-EIMD: 248.5 ± 51.2 Nm vs. post-EIMD: 216.1 ± 71.3 Nm; p < 0.01; ES = 0.52) and an increase in QCSA (pre-EIMD: 3518 ± 965 mm2 vs. post-EIMD: 3675 ± 970 mm2; p < 0.01; ES = 0.16) and VASSQ (pre-EIMD: 4.9 ± 4.9 mm vs. post-EIMD: 49.7 ± 23.3 mm; p < 0.01; ES = 2.6), indicating inflammation and tissue swelling. Consistent with these findings, muscle-water T2 mapping values increased significantly post-EIMD (pre-EIMD: 34.8 ± 0.9 ms vs. post-EIMD: 36.5 ± 0.8 ms; p < 0.01; ES = 1.99). The resting sodium concentration (Na+) did not change significantly (pre-EIMD: 29.8 ± 5.9 mmol vs. post-EIMD: 31.5 ± 6.2 mmol; p = 0.98; ES = 0.28). However, ΔNa+ (change in sodium from rest to post-exercise) was significantly increased post-EIMD (pre-EIMD: 1.7 ± 1.8 mmol vs. post-EIMD: 5.3 ± 3.6 mmol; p < 0.01; ES = 1.26), indicating a greater ionic imbalance during the sustained isometric contraction. Additionally, RPE during the exercise contraction inside the scanner was significantly higher post-EIMD (pre-EIMD: 14.0 ± 2.1 vs. post-EIMD: 18.3 ± 1.2; p < 0.01; ES = 2.5).
FIGURE 1
Muscle oxygenation and metabolism responses following EIMD
Muscle oxygenation and energy metabolism before and after EIMD are shown in Figures 2, 3. The reduction in muscle oxygenation (mBOLD) that occurred during exercise and PECO was greater for post-EIMD measurements than for pre-EIMD measurements [(workload effect: p < 0.01; F(1,20) = 39.41); time effect (p = 0.03; F(2,37) = 4.88); and time * workload interaction (p = 0.02; F(3,60) = 3.53)].
FIGURE 2
FIGURE 3
31P data were available from six participants. Pi at rest was significantly elevated post-EIMD compared with pre-EIMD but showed similar relative increases with exercise and PECO in the pre- and post-EIMD states (p < 0.01; workload effect: F(1,10) = 14.21; time effect: p = 0.02; F(2,20) = 6.63; and time * workload interaction: p = 0.17; F(3,30) = 1.78). Further analysis revealed that intracellular Pi concentration increased significantly with workload (p < 0.01; F(2,16) = 22.63), with a trend for the effect of time (p = 0.068; F(1,10) = 4.25). Extracellular Pi remained unchanged across workloads (p = 0.52; F(2,20) = 0.66) but increased significantly post-EIMD (p = 0.03; F(1,10) = 6.24). The Pi/PCr ratio showed a significant effect of workload (p < 0.01; F(2,17) = 30.92), time (p = 0.048; F(1,10) = 6.83), and time * workload interaction (p = 0.020; F(3,30) = 3.53). ATPγ, PCr, and pH values showed significant changes over phase (ATPγ: p < 0.01, F(2,23) = 39.50; PCr: p < 0.01, F(1,16) = 58.77; and pH: p < 0.01, F(1,13) = 15.83) with no changes in time (ATPγ: p = 0.75, F(1,10) = 0.10; PCr: p = 0.45, F(1,10) = 0.60; and pH: p = 0.12; F(1,10) = 2.73) or time * workload interaction (ATPγ: p = 0.77; F(3,30) = 0.37; PCr: p = 0.22, F(3,30) = 1.45; pH: p = 0.51, F(3,30) = 0.78). The corresponding effect sizes reporting changes in Δexercise and ΔPECO between pre- and post-EIMD values are reported in Figure 5.
Cardiovascular responses following EIMD
Cardiovascular responses (MAP and HR) to exercise are shown in Figure 4. Both MAP and HR increased significantly with workload (p < 0.001), but no significant main effect of time or time * workload interaction was observed, suggesting that the autonomic cardiovascular response remained unchanged post-EIMD despite altered muscle oxygenation and energy metabolism. The corresponding effect sizes reporting changes in Δexercise and ΔPECO between pre- and post-EIMD values are reported in Figure 5.
FIGURE 4
FIGURE 5
Discussion
This study provides the first accounts of the changes in muscle energy metabolism, oxygenation, and perceived exertion following EIMD and their relationship with the cardiovascular responses to exercise and PECO. This was achieved using a more comprehensive multi-parametric MRI approach. The primary novel findings were that damaged skeletal muscles exhibited lower oxygenation and higher Pi during exercise than in the non-damaged state, while MAP and HR were unchanged despite an increased RPE between pre- and post-EIMD. These results suggest that the systemic neuro-vascular regulation remains unaffected despite the substantial local muscle metabolic impairments and muscle weakness experienced by the participants that occurred following EIMD.
Skeletal muscle function following EIMD
The downhill running protocol caused muscular weakness, tenderness, and inflammation when measured 48 h later (Figure 1), which is in line with the anticipated levels of EIMD following downhill running (; ). In this regard, we observed 14% lower MVC at 48 h, which is similar to a past study using downhill running protocols in young healthy individuals (). We also observed an increase in muscle soreness () and inflammation (i.e., muscle-water T2 maps) () and QCSA, indicating localized swelling (; ). Inflammation accumulates extra fluid in the interstitial spaces, which increases the mobility of water molecules within the tissue, leading to a longer T2 relaxation time (). Muscle damage after downhill running occurs as sarcomeres are forcibly lengthened under tension (), affecting extracellular structures, the sarcolemma, and myofibrils (), leading to localized inflammation (), metabolic disruption (), and possible changes to muscle afferent sensitivity and activation ().
Sodium perturbation following EIMD
A novel finding of this study was the greater increase in muscle sodium concentration after a few minutes of damaged muscles, compared to the increase after exercising in non-damaged muscles than the increase after exercising non-damaged muscles. Na+ plays a critical role in maintaining action potential transmission during muscle contraction (). The failure to transport Na+ across the sarcolemma could contribute to excitation–contraction coupling failure, which may contribute to the loss of force production and increased muscle fatigue following eccentric contractions ().
Our findings build on the previous work that has shown increases in resting Na+ concentration 24 h following EIMD (; ), correlating muscle damage markers (i.e., creatine kinase) with increases in Na+ concentrations. Changes in Na+ concentration have been attributed to the disruption of the sarcolemma, involving increased concentration of Na+, calcium ions (), and phosphate (Pi) (). Na+ is expected to increase during exercise (2- to 2.5-fold during muscle contractions) () and remain elevated following fatiguing exercise, where increases in Δ[Na+] are reported to be approximately 30% at 15 min following exercise with slow recovery rates (approximately 60 min) (). Unfortunately, we cannot discriminate changes occurring between intra- and extracellular Na+ with our MRI sequence, which would have been helpful in explaining the mechanism underlying these changes. However, our findings point toward a link between increased ion perturbation and increased muscle fatigue and weakness in damaged muscles.
Muscle oxygenation following EIMD
By applying mBOLD sequences, we found a significant reduction in muscle oxygenation during exercise and PECO, which was greater for damaged muscles than for non-damaged muscles (Figure 2). These findings align with previous research in resting conditions (), where impaired muscle oxygenation was noted 48 h following EIMD (). However, our study provides new evidence of decreased oxygenation during exercise in the thigh muscles, which reflects a mismatch between local perfusion of oxygenated blood and local oxygen utilization. This may occur if damaged muscles utilize relatively more oxygen to hold the sustained isometric contraction without an increase in blood flow. Moreover, our results show that the relative increases in Pi and Pi/PCr from rest to exercise and PECO were greater in damaged than in non-damaged muscles, suggesting that energy turnover and metabolism differed between the two conditions, with a large effect size (Figures 2- 5). Another possible explanation is that the endothelial cells of the microvasculature were impaired due to changed circulation and impaired microvascular responses after muscle damage (). A link has been previously proposed between inflammation and endothelial cell impairments, which might be the underlying mechanism of an impaired microvascular response following EIMD ().
Interactions between muscle metabolism and the cardiovascular responses following EIMD
In line with our hypothesis, we observed significant changes in muscle energy metabolism (higher levels of Pi at rest and during exercise) and reduced oxygenation during exercise of damaged muscles. However, contrary to the second part of our hypothesis, we observed no change in cardiovascular responses (MAP and HR) during exercise and PECO following EIMD. This was surprising considering that an increase in phosphate metabolism and reduced pH has been linked to increased muscle sympathetic nerve activity and blood pressure response to exercise in non-damaged muscle (). In this regard, previous investigations examined the cardiovascular responses to exercise and the metaboreflex by applying PECO to assess the involvement of muscle nerve afferents following EIMD (; ). PECO avoids any involvement of the central command while maintaining activation of muscle afferent activity (). These recent studies used relative loads, which meant that the load at 48 h was lower than that used for baseline conditions. They found contrasting results in MAP responses following EIMD. found attenuated MAP responses during exercise and PECO at 24 h, while in our previous work, no differences in MAP were found at 48 h (). Therefore, it was suggested that these results might have been caused by the increased muscle weakness, the differences in workload used, or the changes in muscle energetics occurring following EIMD. However, in this current study, even after controlling for some of these variables, we found no differences in MAP and HR responses to exercise and PECO 48 h post-EIMD.
It is well-known that blood pressure increases with exercise intensity to help meet the higher demand for oxygen and nutrients by the muscles (; ). The exercise pressor reflex, involving muscle nerve afferents, plays a crucial role in this process (). Our results might suggest that localized damage changes the relationship between muscle mechanoreceptor and metaboreceptor afferent activity, which normally increases the heart rate and blood pressure in proportion to the increase in metabolites and muscle stretch (). The change may be linked to inflammation changing the afferent sensitivity and responsiveness, which can, in turn, affect cardiovascular responses (). Inorganic phosphate accumulation during isometric contractions contributes to increased blood pressure by activating metaboreceptors, which trigger sympathetic responses (). However, in muscle inflammation models, the expected increase in blood pressure might not occur despite Pi accumulation due to factors such as altered afferent signaling, changes in endothelial function, and modulation of the autonomic nervous system (). Our findings suggest a mismatch between cardiovascular responses and muscle Pi in damaged muscles as the increased Pi and Pi/PCr at rest and during exercise and PECO were not accompanied by corresponding increases in blood pressure (). This could be due to metabo-sensitive afferents becoming less responsive to metabolite accumulation, thus interfering with blood pressure regulation during exercise. Previous studies have found blunted metaboreflex activation in patients with chronic heart failure (), while others reported increases in mechanoreflex sensitization in response to cyclooxygenase (COX) but not lactic acid or adenosine produced during exercise (; ). EIMD is known to increase muscle inflammation and the production of prostaglandins and COX (; ), but at the same time, it may also produce higher amounts of lactate, which has been hypothesized to reduce metaboreceptor sensitivity (). Moreover, recent research has found an increased mechanoreflex activation in response to static stretching following delayed onset muscle soreness (DOMS) in the absence of central command, suggesting mechano-nociceptors sensitization following EIMD (), while other studies found no changes or blunted metaboreflex responses to EIMD (; ) when the central command was not activated.
Thus, following EIMD, there may be blunted metaboreflex activation with concomitant hyperactivation of the mechanoreflex, resulting in similar blood pressure responses during exercise and PECO. This mechanism may impact exercise performance as the increase in metabolic demand is not accompanied by a corresponding increase in blood pressure, oxygenation, or other factors contributing to the perception of effort and fatigue (; ; ). It is also important to clarify the relationship between the exercise pressor reflex (EPR)—a cardiovascular response mediated by group III/IV muscle afferents sensitive to mechanical and metabolic stimuli—and DOMS, which is typically associated with muscle damage and inflammation following eccentric exercise. Rather than continuous activation, eccentric exercise is thought to result in sensitization of these thin-fiber afferents, which may contribute to both DOMS and exaggerated reflex cardiovascular responses. Thus, the distinction between EPR and DOMS should not be viewed as a strict dichotomy as both phenomena involve overlapping afferent pathways. Prior work has questioned whether these represent two distinct classes of primary afferent neurons (). Our study did not directly assess nociceptive pathways, which may have influenced the cardiovascular responses observed, and this nuance should be considered when interpreting the findings.
Methodological consideration and study limitations
Unfortunately, due to time constraints and the type of study design we adopted, we could not control for the lack of time-matched controls, making it impossible to exclude time effects unrelated to EIMD. However, this study design has been utilized quite extensively in the literature, and the fact that we controlled for major confounding factors (exercise, medications, and diet) makes our inferences potentially more accurate. Furthermore, the subset of six participants involved in the 31P experiments makes it difficult to draw conclusive arguments about the presence of differences despite the large changes that were also found in the effect size analysis (). Nevertheless, our findings align with previous studies that have observed changes in muscle energetics following EIMD (; ). Additionally, our study acknowledges that some of the metabolites investigated, such as phosphate, might not contribute directly to the EPR, despite previous research in the field (; ). The increases in Pi observed might have been concomitant with exercise but not directly related to the EPR. Previous findings regarding the role of phosphate in the exercise pressor reflex might be over-interpreted, and the potential resetting of afferents at a higher operating point in response to acute alterations was not adequately addressed. However, we found a large effect size reduction in pH post-EIMD, which has been previously shown to be directly correlated with increased protons released and metaboreceptor activation (). Moreover, the observed increased in extracellular phosphate may represent a potential trigger for predominantly metabo-sensitive muscle nerve afferents (group IV) as the terminal endings of these afferents are located outside the cell (; ). However, future research needs to be conducted in this regard to confirm these assumptions. This study included only male participants to reduce variability related to hormonal fluctuations and sex-specific differences in muscle damage responses. However, this limits the generalizability of our findings, and future studies should include female participants to address potential sex-specific differences in EPR and muscle metabolism following EIMD.
Finally, the intensity of the exercise might not have been sufficiently intense to allow for a sufficient metabolite accumulation to trigger the metaboreflex. However, this intensity is still debatable, and no clear indication is present in the literature, making a priori assumptions in humans impossible. Based on our RPE data, we can assume that the intensity of the exercise, especially 48 h post EIMD, was increased, with a concomitant increase in muscle metabolism pointing toward the intense exercise domain.
Conclusion
Our findings provide new insights into the relationship between muscle oxygenation, energy metabolism, and the cardiovascular response to exercise and PECO following EIMD. The novel multi-parametric imaging approach revealed changes in the relationships when damaged skeletal muscles were exercised. We found that oxygenation decreased to lower levels and Pi levels were elevated during exercise with EIMD. This corresponded with elevated perceptions of effort but with unchanged cardiovascular response to exercise and PECO from the baseline condition. These results suggest that the neuro-vascular regulation and blood pressure response to exercise become uncoupled from the peripheral skeletal muscle energetic disturbance with EIMD.
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 the Faculty of Science and Engineering Research Ethics and Governance Committee (reference number: 48097). 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
FZ: Writing – original draft, Data curation, Methodology, Investigation, Visualization, Supervision, Formal analysis, Project administration, Conceptualization. J-CL: Project administration, Data curation, Methodology, Supervision, Writing – review and editing, Investigation. SP: Supervision, Investigation, Writing – review and editing, Project administration. JM: Investigation, Supervision, Methodology, Writing – review and editing. FS: Writing – review and editing, Data curation, Methodology. AK: Writing – original draft, Formal analysis, Writing – review and editing, Project administration, Supervision, Data curation, Methodology, Software, Conceptualization, Investigation. JM: Resources, Conceptualization, Validation, Funding acquisition, Writing – review and editing, Project administration, Investigation, Methodology, Supervision, Visualization, Data curation, Writing – original draft.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Conflict of interest
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.
References
1
AlamM.SmirkF. H. (1937). Observations in man upon a blood pressure raising reflex arising from the voluntary muscles. J. physiology89 (4), 372–383. 10.1113/jphysiol.1937.sp003485
2
AmannM.RunnelsS.MorganD. E.TrinityJ. D.FjeldstadA. S.WrayD. W.et al (2011). On the contribution of group III and IV muscle afferents to the circulatory response to rhythmic exercise in humans. J. Physiol.589 (Pt 15), 3855–3866. 10.1113/jphysiol.2011.209353
3
AmannM.WanH.-Y.ThurstonT. S.GeorgescuV. P.WeavilJ. C. (2020). On the influence of group III/IV muscle afferent feedback on endurance exercise performance. Exerc. Sport Sci. Rev.48 (4), 209–216. 10.1249/JES.0000000000000233
4
AzzabouN.Loureiro de SousaP.CaldasE.CarlierP. G. (2015). Validation of a generic approach to muscle water T2 determination at 3T in fat-infiltrated skeletal muscle. J. Magn. Reson Imaging41 (3), 645–653. 10.1002/jmri.24613
5
BansalN.SzczepaniakL.TernulloD.FleckensteinJ. L.MalloyC. R. (2000). Effect of exercise on 23Na MRI and relaxation characteristics of the human calf muscle. J. Magnetic Reson. Imaging11 (5), 532–538. 10.1002/(sici)1522-2586(200005)11:5<532::aid-jmri9>3.0.co;2-#
6
BontempsB.VercruyssenF.GruetM.LouisJ. (2020). Downhill running: what are the effects and how can we adapt? A narrative review. Sports Med.50 (12), 2083–2110. 10.1007/s40279-020-01355-z
7
BoushelM.NielsenQ. & S.NielsenH. B.QuistorffB.SecherN. H. (1998). Contribution of pH, diprotonated phosphate and potassium for the reflex increase in blood pressure during handgrip. Acta Physiol. Scand.164 (3), 269–275. 10.1046/j.1365-201X.1998.00429.x
8
BoushelR.PottF.MadsenP.RådegranG.NowakM.QuistorffB.et al (1998). Muscle metabolism from near infrared spectroscopy during rhythmic handgrip in humans. Eur. J. Appl. Physiology Occup. Physiology79 (1), 41–48. 10.1007/s004210050471
9
BurtD.LambK.NicholasC.TwistC. (2012). Effects of muscle-damaging exercise on physiological, metabolic, and perceptual responses during two modes of endurance exercise. J. Exerc. Sci. & Fit.10 (2), 70–77. 10.1016/j.jesf.2012.10.003
10
ByrneC.TwistC.EstonR. (2004). Neuromuscular function after exercise-induced muscle damage: theoretical and applied implications. Sports Med.34 (1), 49–69. 10.2165/00007256-200434010-00005
11
CaldwellJ. T.WardlowG. C.BranchP. A.RamosM.BlackC. D.AdeC. J. (2016). Effect of exercise-induced muscle damage on vascular function and skeletal muscle microvascular deoxygenation. Physiol. Rep.4 (22), e13032. 10.14814/phy2.13032
12
ChenT. C.NosakaK.LinM.-J.ChenH.-L.WuC.-J. (2009). Changes in running economy at different intensities following downhill running. J. Sports Sci.27 (11), 1137–1144. 10.1080/02640410903062027
13
DaviesR. C.EstonR. G.FulfordJ.RowlandsA. V.JonesA. M. (2011). Muscle damage alters the metabolic response to dynamic exercise in humans: a 31P-MRS study. J. Appl. Physiol.111 (3), 782–790. 10.1152/japplphysiol.01021.2010
14
DouglasJ.PearsonS.RossA.McGuiganM. (2017). Eccentric exercise: physiological characteristics and acute responses. Sports Med.47 (4), 663–675. 10.1007/s40279-016-0624-8
15
DucrocqG. P.KaufmanM. P. (2020). Inorganic phosphate and lactate potentiate the pressor response to acidic stimuli in rats. Exp. Physiol.105 (4), 613–621. 10.1113/EP088349
16
FouréA.WegrzykJ.Le FurY.MatteiJ. P.BoudinetH.VilmenC.et al (2015). Impaired mitochondrial function and reduced energy cost as a result of muscle damage. Med. Sci. Sports Exerc47 (6), 1135–1144. 10.1249/MSS.0000000000000523
17
FujiiY.OzakiN.TaguchiT.MizumuraK.FurukawaK.SugiuraY. (2008). TRP channels and ASICs mediate mechanical hyperalgesia in models of inflammatory muscle pain and delayed onset muscle soreness. Pain140 (2), 292–304. 10.1016/j.pain.2008.08.013
18
GastL. V.BaierL. M.MeixnerC. R.ChaudryO.EngelkeK.UderM.et al (2023). MRI of potassium and sodium enables comprehensive analysis of ion perturbations in skeletal muscle tissue after eccentric exercise. Invest Radiol.58 (4), 265–272. 10.1097/RLI.0000000000000931
19
GoodwinG. M.McCloskeyD. I.MitchellJ. H. (1972). Cardiovascular and respiratory responses to changes in central command during isometric exercise at constant muscle tension. J. Physiology226 (1), 173–190. 10.1113/jphysiol.1972.sp009979
20
HayashiK.LearyM. E.RoyS. J.LaosiripisanJ.PashaE. P.TanakaH. (2019). Recovery from strenuous downhill running in young and older physically active adults. Int. J. Sports Med.40 (11), 696–703. 10.1055/a-0951-0017
21
HigginsJ. P.GreenS. (2008). Cochrane handbook for systematic reviews of interventions.
22
HögerS. A.GastL. V.MartyB.HotfielT.BickelhauptS.UderM.et al (2023). Sodium and quantitative hydrogen parameter changes in muscle tissue after eccentric exercise and in delayed-onset muscle soreness assessed with magnetic resonance imaging. NMR Biomed.36 (2), e4840. 10.1002/nbm.4840
23
HolodovM.MarkusI.SolomonC.ShaharS.Blumenfeld-KatzirT.GepnerY.et al (2023). Probing muscle recovery following downhill running using precise mapping of MRI T(2) relaxation times. Magn. Reson Med.90 (5), 1990–2000. 10.1002/mrm.29765
24
HottaN.YamamotoK.OgataH.MaherP.OkumuraN.IshidaK. (2016). Does degree of alteration in effort sense caused by eccentric exercise significantly affect initial exercise hyperpnea in humans?J. Physiol. Anthropol.35 (1), 18. 10.1186/s40101-016-0107-5
25
HureauT. J.BroxtermanR. M.WeavilJ. C.LewisM. T.LayecG.AmannM. (2022). On the role of skeletal muscle acidosis and inorganic phosphates as determinants of central and peripheral fatigue: a (31) P-MRS study. J. Physiol.600 (13), 3069–3081. 10.1113/JP283036
26
KanoY.PadillaD. J.BehnkeB. J.HagemanK. S.MuschT. I.PooleD. C. (2005). Effects of eccentric exercise on microcirculation and microvascular oxygen pressures in rat spinotrapezius muscle. J. Appl. Physiol.99 (4), 1516–1522. 10.1152/japplphysiol.00069.2005
27
KaufmanM. P.HayesS. G. (2002). The exercise pressor reflex. Clin. Aut. Res.12 (6), 429–439. 10.1007/s10286-002-0059-1
28
KempG. J.BevingtonA. (1993). The regulation of intracellular orthophosphate concentration. J. Theor. Biol.161 (1), 77–94. 10.1006/jtbi.1993.1041
29
LarsenR. G.ThomsenJ. M.HirataR. P.SteffensenR.PoulsenE. R.FrøkjaerJ. B.et al (2019). Impaired microvascular reactivity after eccentric muscle contractions is not restored by acute ingestion of antioxidants or dietary nitrate. Physiol. Rep.7 (13), e14162. 10.14814/phy2.14162
30
LeeJ. B.KaterbergC.BommaritoJ.PowerG. A.MillarP. J. (2023). Blood pressure responses to postexercise circulatory occlusion are attenuated after exercise-induced muscle weakness. Med. Sci. Sports Exerc55, 1660–1671. 10.1249/MSS.0000000000003182
31
MacLeanD. A.ImadojemuV. A.SinowayL. I. (2000). Interstitial pH, K(+), lactate, and phosphate determined with MSNA during exercise in humans. Am. J. Physiol. Regul. Integr. Comp. Physiol.278 (3), R563–R571. 10.1152/ajpregu.2000.278.3.R563
32
MatsubaraT.HayashiK.WakatsukiK.AbeM.OzakiN.YamanakaA.et al (2019). Thin-fibre receptors expressing acid-sensing ion channel 3 contribute to muscular mechanical hypersensitivity after exercise. Eur. J. pain London, Engl.23 (10), 1801–1813. 10.1002/ejp.1454
33
McCloskeyD. I.MitchellJ. H. (1972). Reflex cardiovascular and respiratory responses originating in exercising muscle. J. Physiology224 (1), 173–186. 10.1113/jphysiol.1972.sp009887
34
McKennaM. J.BangsboJ.RenaudJ.-M. (2008). Muscle K+, Na+, and cl− disturbances and Na+-K+ pump inactivation: implications for fatigue. J. Appl. Physiology104 (1), 288–295. 10.1152/japplphysiol.01037.2007
35
McNeilP. L.KhakeeR. (1992). Disruptions of muscle fiber plasma membranes. Role in exercise-induced damage. Am. J. Pathol.140 (5), 1097–1109.
36
McpheeJ. S.Maden-WilkinsonT. M.NariciM. V.JonesD. A.DegensH. (2014). Knee extensor fatigue resistance of young and older men and women performing sustained and brief intermittent isometric contractions. Muscle & Nerve50 (3), 393–400. 10.1002/mus.24174
37
MenseS. (2010). “Functional anatomy of muscle: muscle, nociceptors and afferent fibers,” in Muscle pain: understanding the mechanisms. Editors MenseS.GerwinR. D. (Berlin, Heidelberg: Springer Berlin Heidelberg), 17–48.
38
MeyerspeerM.BoeschC.CameronD.DezortováM.ForbesS. C.HeerschapA.et al (2020). 31 P magnetic resonance spectroscopy in skeletal muscle: experts' consensus recommendations. NMR Biomed.34, e4246. 10.1002/nbm.4246
39
MiddlekauffH. R.ChiuJ. (2004). Cyclooxygenase products sensitize muscle mechanoreceptors in healthy humans. Am. J. Physiol. Heart Circ. Physiol.287 (5), H1944–H1949. 10.1152/ajpheart.00329.2004
40
MiddlekauffH. R.SinowayL. I. (2007). Increased mechanoreceptor stimulation explains the exaggerated exercise pressor reflex seen in heart failure. J. Appl. Physiol.102 (1), 492–496. 10.1152/japplphysiol.00994.2006
41
MilesM. P.LiY.RinardJ. P.ClarksonP. M.WilliamsonJ. W. (1997). Eccentric exercise augments the cardiovascular response to static exercise. Med. Sci. Sports Exerc29 (4), 457–466. 10.1097/00005768-199704000-00006
42
MurphyM. N.MizunoM.MitchellJ. H.SmithS. A. (2011). Cardiovascular regulation by skeletal muscle reflexes in health and disease. Am. J. Physiol. Heart Circ. Physiol.301 (4), H1191–H1204. 10.1152/ajpheart.00208.2011
43
PattenC.MeyerR. A.FleckensteinJ. L. (2003). T2 mapping of muscle. Semin. Musculoskelet. Radiol.7 (4), 297–305. 10.1055/s-2004-815677
44
PattersonS. D.HughesL.WarmingtonS.BurrJ.ScottB. R.OwensJ.et al (2019). Blood flow restriction exercise: considerations of methodology, application, and safety. Front. physiology10, 533. 10.3389/fphys.2019.00533
45
PaulsenG.EgnerI. M.DrangeM.LangbergH.BenestadH. B.FjeldJ. G.et al (2010). A COX-2 inhibitor reduces muscle soreness, but does not influence recovery and adaptation after eccentric exercise. Scand. J. Med. & Sci. Sports20 (1), e195–e207. 10.1111/j.1600-0838.2009.00947.x
46
PeakeJ. M.NeubauerO.Della GattaP. A.NosakaK. (2017). Muscle damage and inflammation during recovery from exercise. J. Appl. Physiol.122 (3), 559–570. 10.1152/japplphysiol.00971.2016
47
RayC. A.MahoneyE. T.HumeK. M. (1998). Exercise-induced muscle injury augments forearm vascular resistance during leg exercise. Am. J. Physiol.275 (2), H443–H447. 10.1152/ajpheart.1998.275.2.H443
48
RenaudJ. M.ØrtenbladN.McKennaM. J.OvergaardK. (2023). Exercise and fatigue: integrating the role of K(+), Na(+) and Cl(-) in the regulation of sarcolemmal excitability of skeletal muscle. Eur. J. Appl. Physiol.123, 2345–2378. 10.1007/s00421-023-05270-9
49
ReyngoudtH.SmithF. E.Caldas de Almeida AraújoE.WilsonI.Fernández-TorrónR.JamesM. K.et al (2022). Three-year quantitative magnetic resonance imaging and phosphorus magnetic resonance spectroscopy study in lower limb muscle in dysferlinopathy. J. Cachexia, Sarcopenia Muscle13 (3), 1850–1863. 10.1002/jcsm.12987
50
RottoD. M.SchultzH. D.LonghurstJ. C.KaufmanM. P. (1990). Sensitization of group III muscle afferents to static contraction by arachidonic acid. J. Appl. Physiol.68 (3), 861–867. 10.1152/jappl.1990.68.3.861
51
RowellL. B. (1992). Reflex control of the circulation during exercise. Int. J. Sports Med.13 (Suppl. 1), S25–S27. 10.1055/s-2007-1024583
52
SinowayL. I.SmithM. B.EndersB.LeuenbergerU.DzwonczykT.GrayK.et al (1994). Role of diprotonated phosphate in evoking muscle reflex responses in cats and humans. Am. J. Physiol.267 (2 Pt 2), H770–H778. 10.1152/ajpheart.1994.267.2.H770
53
SleighA.LupsonV.ThankamonyA.DungerD. B.SavageD. B.CarpenterT. A.et al (2016). Simple and effective exercise design for assessing in vivo mitochondrial function in clinical applications using 31P magnetic resonance spectroscopy. Sci. Rep.6 (1), 19057. 10.1038/srep19057
54
StaceyM. J. (1969). Free nerve endings in skeletal muscle of the cat. J. Anat.105 (Pt 2), 231–254.
55
StenvinkelP. (2001). Endothelial dysfunction and Inflammation—Is there a link?Nephrol. Dial. Transplant.16 (10), 1968–1971. 10.1093/ndt/16.10.1968
56
SternsD. A.EttingerS. M.GrayK. S.WhislerS. K.MosherT. J.SmithM. B.et al (1991). Skeletal muscle metaboreceptor exercise responses are attenuated in heart failure. Circulation84 (5), 2034–2039. 10.1161/01.cir.84.5.2034
57
StožerA.VodopivcP.Križančić BombekL. (2020). Pathophysiology of exercise-induced muscle damage and its structural, functional, metabolic, and clinical consequences. Physiol. Res.69 (4), 565–598. 10.33549/physiolres.934371
58
TaylorJ. L.AmannM.DuchateauJ.MeeusenR.RiceC. L. (2016). Neural contributions to muscle fatigue: from the brain to the muscle and back again. Med. & Sci. Sports & Exerc.48 (11), 2294–2306. 10.1249/MSS.0000000000000923
59
TeixeiraA. L.ViannaL. C. (2022). The exercise pressor reflex: an update. Clin. Auton. Res.32, 271–290. 10.1007/s10286-022-00872-3
60
TwistC.EstonR. G. (2009). The effect of exercise-induced muscle damage on perceived exertion and cycling endurance performance. Eur. J. Appl. Physiol.105 (4), 559–567. 10.1007/s00421-008-0935-z
61
WenkH. N.McCleskeyE. W. (2007). A novel mouse skeletal muscle-nerve preparation and in vitro model of ischemia. J. Neurosci. Methods159 (2), 244–251. 10.1016/j.jneumeth.2006.07.021
62
WesterbladH.AllenD. G.LännergrenJ. (2002). Muscle fatigue: lactic acid or inorganic phosphate the major cause?News Physiol. Sci.17, 17–21. 10.1152/physiologyonline.2002.17.1.17
63
WhiteheadN. P.AllenT. J.MorganD. L.ProskeU. (1998). Damage to human muscle from eccentric exercise after training with concentric exercise. J. Physiol.512, 615–620. 10.1111/j.1469-7793.1998.615be.x
64
WhiteheadN. P.WeerakkodyN. S.GregoryJ. E.MorganD. L.ProskeU. (2001). Changes in passive tension of muscle in humans and animals after eccentric exercise. J. Physiol.533 (Pt 2), 593–604. 10.1111/j.1469-7793.2001.0593a.x
65
WilliamsN. (2017). The borg rating of perceived exertion (RPE) scale. Occup. Medicine-Oxford67 (5), 404–405. 10.1093/occmed/kqx063
66
WilliamsS.CarsonR.TóthK. (2023). Moving beyond P values in the journal of physiology: a primer on the value of effect sizes and confidence intervals. J. Physiology601 (23), 5131–5133. 10.1113/JP285575
67
YeungE. W.BallardH. J.BourreauJ. P.AllenD. G. (2003). Intracellular sodium in mammalian muscle fibers after eccentric contractions. J. Appl. Physiol.94 (6), 2475–2482. 10.1152/japplphysiol.01128.2002
68
ZambolinF.GiuriatoG.LaginestraF. G.OttavianiM. M.FavarettoT.CalabriaE.et al (2022a). Effects of nociceptive and mechanosensitive afferents sensitization on central and peripheral hemodynamics following exercise-induced muscle damage. J. Appl. Physiology133 (4), 945–958. 10.1152/japplphysiol.00302.2022
69
ZambolinF.McPheeJ. S.Duro-OcanaP.GanseB.BagleyL.FaisalA. (2022b). The association of elevated blood pressure during ischaemic exercise with sport performance in master athletes with and without morbidity. Eur. J. Appl. Physiology122 (1), 211–221. 10.1007/s00421-021-04828-9
70
ZambolinF.PeçanhaT.PinnerS.VenturelliM.McPheeJ. S. (2023). Effects of exercise induced muscle damage on cardiovascular responses to isometric muscle contractions and post-exercise circulatory occlusion. Eur. J. Appl. Physiology123, 2747–2754. 10.1007/s00421-023-05255-8
71
ZambolinF.Duro OcanaP.GouldingR.SandersonA.VenturelliM.WoodG.et al (2024a). The corticomuscular response to experimental pain via blood flow occlusion when applied to the ipsilateral and contralateral leg during an isometric force task. Psychophysiology61 (3), e14466. 10.1111/psyp.14466
72
ZambolinF.LaginestraF. G.FavarettoT.GiuriatoG.OttavianiM. M.SchenaF.et al (2024b). Activation of skeletal muscle mechanoreceptors and nociceptors reduces the exercise performance of the contralateral homologous muscles. Am. J. Physiology-Regulatory, Integr. Comp. Physiology327 (4), R389–R399. 10.1152/ajpregu.00069.2024
Summary
Keywords
magnetic resonance spectroscopy, muscle oxygenation, muscle metabolism, group III–IV muscle afferents, exercise-induced muscle damage, muscle inflammation, delayed onset of muscle soreness
Citation
Zambolin F, Lagacé J-C, Pinner S, McStravick J, Smith FE, Kennerley AJ and McPhee JS (2025) Muscle energetics and the cardiovascular response to isometric exercise and post-exercise circulatory occlusion following exercise-induced muscle damage: insights from multiparametric MRI. Front. Physiol. 16:1636964. doi: 10.3389/fphys.2025.1636964
Received
28 May 2025
Accepted
26 September 2025
Published
29 October 2025
Volume
16 - 2025
Edited by
Anna Pedrinolla, University of Trento, Italy
Updates
Copyright
© 2025 Zambolin, Lagacé, Pinner, McStravick, Smith, Kennerley and McPhee.
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: Fabio Zambolin, f.zambolin@mmu.ac.uk
ORCID: Fabio Zambolin, orcid.org/0000-0002-4178-6965; James McStravick, orcid.org/0000-0001-6067-9893; Aneurin J. Kennerley, orcid.org/0000-0002-7599-7461
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.