- Department of Neuromedicine and Movement Science, Norwegian University of Science and Technology, Trondheim, Norway
A decline in skeletal muscle function is a key factor contributing to reduced functional capacity and quality of life with aging. While the mechanisms are multifactorial, impairments in mitochondrial and microvascular function are believed to play a central role. Near-infrared spectroscopy (NIRS) has emerged as a promising non-invasive tool for identifying aging biomarkers and evaluating interventions to preserve muscle health. However, aging often coincides with reduced physical fitness due to lifestyle changes, while age-related changes in body composition further complicate the identification of biomarkers using NIRS. This study investigated aging effects on microvascular and mitochondrial function in healthy adults and explored how physical performance and body composition influence these functions. Eighteen healthy young adults (25.8 ± 2.3 years) and eighteen healthy older adults (69.8 ± 6.0 years) participated. Mitochondrial and microvascular function in arm and leg muscles were simultaneously measured using NIRS during and following vascular occlusion. Physical fitness was evaluated through graded cycling and handgrip tests, and body composition by skinfold thickness and bioimpedance analysis. Data were analyzed using repeated measures ANOVA and multiple regression analysis. Although baseline muscle oxygen saturation (SmO2) was significantly lower in older adults (FDS: 62.8% ± 1.5% vs. 66.6% ± 1.5%; VL: 73.9% ± 2.9% vs. 79.5% ± 2.2%; both P < 0.001), no significant group differences were found in muscle oxygen consumption (mVO2), reperfusion rate (RR), or desaturation/resaturation responses. A significant interaction between muscle and age suggested muscle-specific age-related differences. Multiple regression analysis using forward selection revealed age as a moderate predictor of forearm mVO2, explaining 20.3% of the variance, while adipose tissue thickness (ATT) was the strongest predictor of leg mVO2, explaining 53.9% of the variance. Although the regression analyses should be considered exploratory and interpreted with caution, given the limited sample size and use of forward selection, these findings highlight the complexity of investigating aging effects on muscle function and underscore the importance of accounting for individual fitness and body composition, particularly local ATT, when interpreting NIRS data. To accurately assess age-related differences, future studies should report detailed information on both factors. Addressing these confounders is essential for understanding the physiological impact of aging on muscle function measured by NIRS.
1 Introduction
Ageing is commonly associated with an increased susceptibility to disease and a decline in physical function (Burton and Sumukadas, 2010). One of the most evident changes that happen with advancing age is the gradual loss of skeletal muscle strength and mass (Mitchell et al., 2012), as well as a progressive decline in aerobic capacity (Tanaka and Seals, 2008). The association between aging and a decline in aerobic exercise capacity, commonly expressed as maximal oxygen uptake (VO2max), is well-known and VO2max is thought to decline at a rate of about 10% per decade in sedentary people after the age of 25 years (Betik and Hepple, 2008). Although mostly presented as a linear decline, recent evidence indicates a nonlinear age-related decline in VO2max with accelerated rates above 70 years of age (Hawkins and Wiswell, 2003). The reduction in aerobic capacity is of relevance as it is considered to be an independent predictor of cardiovascular morbidity and mortality (Lang et al., 2024), regardless of age.
Much less clear are the declines in muscle mass and muscle strength. Estimates of the age-related decline in muscle mass vary greatly and have shown to be affected by gender and ethnicity (Gallagher et al., 1997; Silva et al., 2010). Median values of 4.7% loss of peak muscle mass per decade for men and 3.7% for women have been reported across studies (Mitchell et al., 2012). Moreover, while skeletal muscle mass is associated with muscle strength, the loss of strength seems to proceed at a faster rate than the loss of mass (Goodpaster et al., 2006; Delmonico et al., 2009), suggesting that other factors beyond muscle mass contribute to strength loss in older adults (Clark and Manini, 2008). While eventually all aging individuals will lose muscle mass, strength and physical function, and this loss can thus be considered inherent to the typically aging process, the extent and speed of the decline can vary greatly (Mitchell et al., 2012). Over time, time, severe loss of muscle strength, mass and function may result in sarcopenia (Cruz-Jentoft et al., 2010; Cruz-Jentoft et al., 2018), which has recently been formally recognized as a muscle disease (Cruz-Jentoft et al., 2018; Vellas et al., 2018). Sarcopenia is an independent predictor of adverse health outcomes (Landi et al., 2016), further underscoring the health implications of aging. Although operational definitions for sarcopenia enable diagnosis and tracking, they reflect a prolonged trajectory of decline. Identifying physiological markers that detect changes earlier, such as those potentially measurable by NIRS, could help mitigate or prevent functional deterioration.
A prominent challenge in aging research is the differentiation between changes attributable to aging and those resulting from reduced physical activity as both contribute to declines in muscle mass, strength and aerobic capacity. Some authors even propose that physical inactivity is the main contributor to these changes (Tieland et al., 2018) since substantial evidence shows that physical activity can preserve muscle mass, strength and function in older adults (Liu et al., 2014; Grevendonk et al., 2021; Hortobágyi et al., 2023; Choi et al., 2024) and may even benefit sarcopenic individuals (Lozano-Montoya et al., 2017). However, despite the positive effects of physical activity, structural changes occur with aging that appear to be independent of activity level. These include degradation of elastic fibers, increased stiffness in tissues such as the lungs and vasculature (Bloom et al., 2023), which may lead to higher strain during exercise. Reduced vessel compliance increases systolic blood pressure and compromises peripheral blood flow while increasing cardiac workload. Endothelial function becomes less responsive to physiological stimuli due to cellular aging, impairing vasodilation and regulation of local blood flow (Scioli et al., 2014; Han and Kim, 2023; Ajoolabady et al., 2024). At the muscle level, aging is associated with reduced capillary density and diminished mitochondrial function, contributing to lower ATP production and increased fatigue (Grevendonk et al., 2021; Distefano and Goodpaster, 2018; Gonzalez-Freire et al., 2015; Conley et al., 2013; Short et al., 2005). These microvascular and mitochondrial changes are central to age-related declines in muscle performance and metabolic health, and may persist even in physically active older adults (Grevendonk et al., 2021). Understanding the physiological changes that occur during typically aging is essential, as this represents the baseline trajectory against which pathological deviations can be identified and interpreted. However, a major challenge is disentangling intrinsic aging effects from those that are influenced by physical fitness as both contribute to reductions in muscle mass, strength, and aerobic capacity.
To better interpret age-related changes in muscle function and distinguish between intrinsic and lifestyle-related effects, it is important to clarify key concepts related to physical capacity. In the context of aging and muscle function, it is important to distinguish between physical activity, physical fitness, and physical performance. Physical activity refers to any bodily movement produced by skeletal muscles that requires energy and is typically assessed through self-reported measures such as frequency, duration, and intensity. In this study, physical fitness refers specifically to cardiorespiratory fitness, assessed through objective measures like VO2max and peak work rate. Physical performance refers to the ability to carry out specific tasks and was evaluated using incremental cycling and handgrip tests. Although these definitions are distinct, physical activity and fitness are related, as increased activity generally improves physical fitness. Physical activity was assessed via validated questionnaires, while fitness and performance were evaluated through graded exercise tests to provide a comprehensive view of participants’ functional status.
Near-infrared spectroscopy (NIRS) has emerged as a promising non-invasive tool to assess muscle oxygen consumption and vascular function in vivo. By measuring tissue saturation dynamics during and after vascular occlusion, NIRS provides insight into both microvascular and mitochondrial function. However, interpretation of NIRS data can be complicated by individual differences in adipose tissue thickness (ATT), as ATT is known to confound the NIRS measurements. While the general effect of ATT on NIRS is well established, aging is also associated with an increase in intramuscular fat deposits (Farsijani et al., 2020; Wang et al., 2014). Thus, although NIRS holds potential for studying age-related changes in muscle function, factors such as physical fitness and body composition may mask or exaggerate these changes, making it difficult to isolate the effects of aging per se.
Therefore, the main purpose of this study was to investigate possible aging effects on microvascular and mitochondrial function in typically aging older adults. In addition, we investigated the relationship between microvascular and mitochondrial function with performance variables, as an indication of physical fitness, and with body composition variables, as possible confounding variables for the continuous-wave NIRS method. By combining group comparisons with regression analyses, we sought to disentangle the effects of age from those of body composition and training status, and to contribute to a more nuanced understanding of muscle function in typically aging populations. We hypothesized that aging would be associated with reduced microvascular and mitochondrial function, and that both physical fitness and ATT would significantly influence NIRS-derived outcomes.
2 Materials and methods
2.1 Participants
Nineteen young (25.9 ± 2.3 years) and twenty older adults (69.7 ± 5.8 years) participants volunteered to participate in this study. All were healthy, non-obese, non-smoking and recreational active participants without a history of cardiovascular, pulmonary or metabolic diseases. None of the participants were using medication that might affect the acute hemodynamic responses during exercise. All but two of the participants were right-handed. Activity level was self-reported, and participants were asked how often (frequency), how long (duration) and how hard (intensity) they performed physical activity on a weekly basis. The questions for training frequency, training duration and training intensity were previously used in the Nord-Trøndelag Health study and have recently been validated (Kurtze et al., 2007; Vie et al., 2023). Training frequency, duration and intensity were recoded as described by Vie et al. (2023), and physical activity (PA) level was subsequently quantified by multiplying training frequency, duration and intensity and reported in minutes per week (Vie et al., 2023). The study was approved by the Norwegian Regional Ethics Committee, and all participants gave their written, informed consent prior to data collection.
2.2 Experimental design
All participants visited the lab on 2 separate occasions, consisting of a short visit for the measurement of body composition using bioelectrical impedance analysis (BIA) and one test day. BIA measurements were done during a separate visit in the early morning and after an overnight fast, typically, after the test day. The measurements on the test day started with a vascular occlusion test, followed by an incremental handgrip test and an incremental cycling test, to obtain both resting as well as peak exercise values. Participants were instructed to refrain from alcohol intake and vigorous exercise 24 h prior to testing, limit their caffeine-intake 6 h prior to testing, and to avoid intake of any large meal less then 3 h prior to testing. Participants were encouraged to be well rested and hydrated when reporting to the lab for exercise testing.
2.3 Experimental procedures
The measurements on the test day started with a vascular occlusion test, followed by an incremental handgrip test and an incremental cycling test. NIRS equipment was attached prior to the 3 tests and recordings were continued during the complete session. Markers were set in the NIRS software (Oxysoft, Artinis Medical Systems, BV, Netherlands) to identify the start and end of each test and used to divide the data file in separate tests during the analysis process. Start and end of the various events were likewise marked. Prior to the tests, blood pressure was measured using an automatic portable blood pressure monitor (OSZ 5 easy, Welch Allyn Jungingen, Germany) that measured systolic blood pressure (SBP) and diastolic blood pressure (DBP). Mean arterial pressure (MAP) was calculated as MAP = DBP +1/3(SBP - DBP).
2.3.1 Vascular occlusion test (VOT)
The vascular occlusion test (VOT) was used to measure oxygenation and saturation characteristics of the muscles in resting state. The participants were placed in a comfortable semi-supine position, supporting arms and knees to prevent overextension and to optimize relaxation. Vascular occlusion of arm and leg was applied using two pneumatic cuffs (Hokanson SC5 and Hokanson SC10, Marcom Medical ApS, Denmark) that were attached proximally around the right upper arm and the right thigh. To ensure rapid inflation and deflation (<0.5 s), an automatic cuff inflator (Hokanson E20 Rapid Cuff Inflator, Marcom Medical ApS, Denmark) was used in combination with an external pressurized air source (8 bar, NTNU). Cuff pressure was set to 300 mmHg to ensure blockage of both inflow and outflow to the limbs and both cuffs were inflated and deflated simultaneously. NIRS signals were measured continuously during the whole test and simultaneously in both muscles. The vascular occlusion test started with a 7.5 -minute period to establish baseline saturation levels, followed by a 1-min occlusion for familiarization purposes. After cuff release and a 3-min recovery period, a 10-min occlusion was applied followed by at least 5 min of recovery. The duration of the occlusion was based on laboratory experience and observations indicating that approximately 30% of participants do not reach a desaturation plateau within 5 min (unpublished data). The longer duration ensured consistent and reliable deoxygenation across all individuals. The participants were instructed to refrain from any movements during the entire test, with special emphasis on the initial phase of the occlusion and the initial phase of recovery. The duration for measurement of the recovery phase was guided by the occurrence of the hyperemic response and lasted either 5 min or as long as it took to monitor the signals past a clear maximum response after the release of the cuffs.
2.3.2 Incremental handgrip test (IHT)
The vascular occlusion test was followed by the incremental handgrip test without a change in position. Prior to the incremental handgrip test, handgrip strength (HGS) of the right hand was measured using a handheld handgrip dynamometer (Lafayette Instruments, Model 5030L1, Indiana, United States). Testing was done in seated position with the elbow in 90° flexion and the shoulder and forearm in neutral position. HGS was defined as the absolute highest handgrip force value of three 3-s trials separated by 1 min of recovery. After a short break, forearm muscle endurance was measured with an incremental handgrip test (IHT) that consisted of dynamic handgrip (HG) exercise using a custom-made handgrip dynamometer (Figure 1). Participants were in the same comfortable semi-supine position as during the vascular occlusion test, and at least 15 min separated both tests. The right hand rested on the custom-made handgrip dynamometer with the upper arm at heart level and the forearm in an upward angle of 30° to avoid venous pooling of the blood (Figure 1). The arm was supported at the wrist and elbow to avoid contact between forearm and dynamometer, ensuring unrestricted circulation of the forearm. Participants raised and lowered a bucket through a distance of approximately 5 cm by squeezing a handgrip device at a contraction rate of 30 min–1 (duty cycle 50%) and guided by a metronome, resulting in dynamic exercise without a static component. The opened handgrip position was determined by an ∼90° angle between the second and third phalanx of the index finger and marked with blue tape on the dynamometer as a visual guidance marker. The closed handgrip position was defined by the handle touching the handlebar. After a minimum of 3 min with baseline measurements, the incremental handgrip test started with 1 min dynamic HG exercise at an initial load of 2.5 kg followed by a gradually increasing load with 250 g increments every 15 s until volitional exhaustion. The test was terminated when contraction rate could no longer be maintained without compensatory movements or lack of contact between handle and handlebar. Handgrip force was measured using a high precision load cell (Revere Transducers Europe B.V., Model 9363-D3-50kg-20T1) mounted in series with the lifting cable of the dynamometer. Calibration of the force cell was done twice, prior to and following the complete measurement of all participants. Force was measured continuously during the handgrip test and recorded as an analog channel within the NIRS acquisition software, ensuring simultaneous sampling and automatic synchronization with the NIRS data. Performance variables derived from the incremental handgrip test were handgrip strength (HGS) and peak dynamic handgrip load (DHGpeak) where HGS was defined as the absolute highest handgrip strength and DHGpeak as the final workload completed for a full 15-s period prior to termination of the incremental handgrip test.
Figure 1. Schematic representation of the custom-made handgrip dynamometer used during the incremental handgrip test (IHT). The right hand rested on the handgrip dynamometer with the upper arm at heart level and the forearm in an upward angle of 30° to avoid venous pooling of the blood. The arm was supported at the wrist and elbow to avoid contact between forearm and dynamometer, ensuring unrestricted circulation of the forearm. The NIRS device was placed on top of the flexor digitorum muscle and was not in contact with the dynamometer. The dynamometer cable was connected with a bucket that was gradually increased in load (see text for details). Participants raised and lowered the bucket through a distance of ca. 5 cm by squeezing the handgrip device at a contraction rate of 30 min–1 and guided by a metronome, resulting in dynamic exercise without a static component. Handgrip force was measured using a high precision load cell (Revere Transducers Europe B.V., Model 9363-D3-50kg-20T1) mounted in series with the lifting cable of the dynamometer (not visible in the drawing).
2.3.3 Incremental cycling test (ICT)
Following the handgrip test, and without disconnecting the equipment, participants switched to a cycle ergometer (Lode Excalibur Sport, Lode B. V., Groningen, Netherlands) to perform an incremental cycling test (ICT). A warm-up period of 5–10 min preceded the test, with external work rates ranging from 50 to 100 W depending on age, sex, and cycling experience. Participants were instructed to cycle at a self-selected cadence between 80 and 100 RPM throughout the tests. To determine the work rate (WR) at the onset of blood lactate accumulation (WROBLA; 4 mmol/L), participants cycled in 4-min stages using age- and sex-specific starting work rates and increments as previously described by Skovereng et al. (2016). Blood lactate was measured at the end of each stage, and WROBLA was defined as the last completed stage before exceeding 4 mmol/L. After 5–10 min of active recovery (50–100 W), the incremental test began at a WR that was one level below WROBLA. WR increments were 25 W/min for young men, 20 W/min for young women and older men, and 15 W/min for older women and were continued until volitional exhaustion or cadence dropped below 60 RPM. Strong verbal encouragement was given throughout the test to promote maximal effort. Immediately following termination of the test, blood lactate concentration was measured, and the participants were asked to rate their perceived exertion (RPE) using the modified Borg scale [6–20 Borg scale (Borg, 1970)]. Blood samples (0.3 µL) for lactate measurements were taken before warm-up, during the lactate phase, and after the incremental test using a handheld analyzer (Lactate Pro2 LT-1730, Arkray KDK, Kyoto, Japan). Muscle oxygenation, heart rate, pulmonary oxygen uptake, and work rate were measured continuously during the whole test. Performance variables derived from the incremental cycling test were peak oxygen uptake (VO2peak), peak work rate (WRpeak) and peak heart rate (HRpeak). VO2peak was defined as the highest 30-s average value, WRpeak as the highest WR sustained for 1 min and HRpeak as the highest 5-s average value. WRpeak and VO2peak were used to construct the model for multiple linear regression.
2.4 Measurements
2.4.1 Near-infrared spectroscopy (NIRS)
A continuous-wave near-infrared spectrophotometer (Oxymon MKIII, Artinis Medical Systems, Netherlands), and two portable NIRS devices (Portamon, Artinis Medical Systems, Netherlands) were used to measure relative changes in optical density at rest as well as during handgrip and cycling exercise. The portamon devices were placed in a consisted order on top of the flexor digitorum superficialis (FDS) and the vastus lateralis (VL) muscles as the main active muscles in handgrip and cycling exercise, respectively. Simultaneous measurement of both muscles was done through a single data acquisition system (Oxysoft, Artinis Medical Systems, BV, Netherlands). The Oxymon system was dormant but used to enable simultaneous data acquisition of handgrip force and cuff pressure within Oxysoft. Excessive body hair at the measurement site, judged visually and based on prior experience, was removed to ensure optimal light propagation. The portamon devices were placed directly on the skin on top of the muscle belly, fixed with adhesive tape to prevent movement between skin and portamon device, and then covered with a lightproof cloth and elastic tape. Both channels consisted of three transmitters and one receiver. The LED transmitters of the portamon devices generated light at 844.5 ± 3.2 nm and 761.7 ± 1.0 nm, and the fixed source-detector distances were 30, 35 and 40 mm (transmitter separation of 5 mm). These CW devices estimate relative changes in tissue oxygenation and do not provide absolute chromophore concentrations. Spatially resolved spectroscopy (SRS) was used to continuously measure muscle oxygen saturation (SmO2) (Patterson et al., 1989) during vascular occlusion as well as during the handgrip and cycling tests. Data were sampled at 10 Hz, displayed real-time and stored on disk for off-line analysis. All NIRS measurements were done on the right side of the body.
After completing the tests, NIRS optodes were detached and skinfold thickness was measured on top of both muscles and in between source and detector using a skinfold caliper (Holtain Ltd, Crymmych, United Kingdom). All skinfold measurements were done twice and in series where the mean of the two measurements per muscle was used as an estimate of adipose tissue thickness (ATT = fat + skin layer = skinfold/2). In addition, limb circumference was measured at the same site as that used for skinfold thickness measurements. Skinfold thickness was measured 10–15 min after the cycling test. Given the short duration of the test (≤20 min), we consider the impact of acute exercise on compressibility and ATT values to be negligible. Antropometric measurements were done by the same skilled investigator.
2.4.2 Cardiorespiratory assessments
Heart rate (HR) was measured continuously during all tests using a heart rate monitor (Polar RS800, Polar Electro OY, Kempele, Finland) with a sensor belt (Wearlink, Polar Electro OY, Kempele, Finland) and providing average HR samples over 5 s. Pulmonary gas exchange was measured during the cycling test by open-circuit indirect calorimetry and using a mixing chamber (Oxycon Pro, Jaeger GmbH, Hoechberg, Germany). Data were sampled every 10 s. Every test day, the equipment was calibrated using a 3-liter calibration syringe (Hans Rudolph Inc, Kansas City, MO, United States) to calibrate the flow turbine and a high-precision gas mixture (16.0% ± 0.04% O2 and 5.0% ± 0.1% CO2; Riessner-Gase GmbH & Co., Lichtenfels, Germany) to calibrate the gas analyzers. Both HR and gas exchange recordings were started simultaneously, and a marker was manually set in both the Jaeger and Oxymon software to synchronize with the NIRS measurements.
2.4.3 Body composition assessments
Body composition and weight, recorded to the nearest 0.1 kg, were measured using bioelectrical impedance analysis (BIA) (Inbody 770, Biospace, Seoul, Korea). Body height was measured during the same occasion using a wall-mounted stadiometer (Seca 222, Seca GmbH & Co, Hamburg, Germany) and recorded to the nearest 0.5 cm. Measurements were performed in fasted state in the early morning, wearing light sports clothes and no shoes. BIA measurements were done in accordance with the manufacturer’s guidelines. Participants were informed about the procedure and asked to urinate and/or defecate before the measurements. BIA measurements of interest were body weight and percent body fat (PBF) as well as segmental estimates of fat free mass (SFFM), and body fat mass (SBFM). SFFM and SBFM for the right arm and leg were used to calculate limb mass (LM = SFFM + SBFM), and percent body fat of the limb (SPBF (%) = SBFM/LM). This resulted in an estimation of PBF for the right arm (PBF-RA) and for the right leg (PBF-RL). PBF-RA and PBF-RL were used to construct the model for multiple linear regression as a peripheral factor possibly affecting the NIRS measurements. Although BIA accuracy can be influenced by age, sex, and hydration status, the standardized protocol used in this study was intended to minimize variability across participants.
2.5 Data analysis
2.5.1 Peripheral measurements (NIRS)
NIRS signals were imported in matlab for further analysis and calculation of peripheral variables. Figure 2 shows a schematic representation of the SmO2 signal during the vascular occlusion test as well as the various variables that were derived from this test. All variables were calculated for both muscles using the same procedures. Baseline saturation (BL) was calculated as the mean saturation level (in %) over the last 30-s period prior to the start of the vascular occlusion or exercise period. Muscle oxygen consumption (mVO2) during rest was calculated from the rate of decrease in SmO2 (in %/s) during the initial phase of the occlusion. To omit possible artefacts due to inflation, the first 10 s following cuff inflation were omitted, and linear regression was done over the subsequent 120-s period (1200 samples). The maximum desaturation level (VOT-min) during vascular occlusion was calculated as the mean saturation level (in %) over the last 30-s period of the 10-min occlusion period and expressed as the relative change (Δ%) from baseline saturation levels. Reperfusion rate (RR) following cuff release was calculated from the rate of increase in SmO2 (in %/s) during the initial period of recovery. To omit possible artefacts due to the cuff deflation, the first 0.5 s following the release of the cuff was omitted. Linear regression was done over the subsequent 10-s period following the release of the cuff. Based on the short period that is used for this analysis, it is thought that the rapid increase in saturation primarily depends on endothelial function (Creteur et al., 2007; Doerschug et al., 2007; Skarda et al., 2007; Nanas et al., 2008; Gerovasili et al., 2010; Barstow, 2019) and, therefore, provides information about microvascular function. The maximum resaturation level (HP-max) was calculated as the highest peak value after cessation of the cuff (in %) and relative (Δ%) to the start of recovery and was used as a measure to quantify the hyperemic response. Calculation of HP-max was performed after smoothing the recovery signal using a simple moving average with a window size of 100 samples, corresponding to approximately 10 s. Group responses in SmO2 responses during and following vascular occlusion were calculated as mean and SD for 20 intervals between start and end of the test, as well as 5 pre-occlusion intervals and 10 post-occlusion intervals, all with 30 s duration.
Figure 2. Schematic representation of the NIRS signals during the vascular occlusion test as well as the variables that were derived from this test. 1) Baseline saturation levels (BL) were calculated as the mean saturation level in percentage (%) over the last 30-s period prior to the start of the vascular occlusion. 2) Oxygen consumption (mVO2) during rest was calculated from the rate of decrease in SmO2 during the initial phase the occlusion, omitting the first 2 s to omit possible inflation artefacts. Linear regression was done over the subsequent 45-s period and expressed in %/s Maximum desaturation levels (VOT-min) during vascular occlusion were calculated as the mean desaturation over the last 30-s period of the 10-min occlusion and expressed as the relative change from baseline saturation levels. 4) Reperfusion rate (RR) following cuff release was calculated from the rate of increase in SmO2 during the initial period of recovery, omitting the first 0.5 s to omit possible deflation artefacts. Linear regression was done over the subsequent 10-s period and expressed in %/s Maximum resaturation (HP-max) was calculated as the highest peak after cessation of the cuff and relative to the start of recovery after smoothing the recovery signal using a simple moving average with a window size of 100 samples, corresponding to approximately 10 .
2.6 Statistical analysis
To determine the required sample size for this study, we used published data from the literature for the major outcome variables: muscle oxygen consumption (Rosenberry et al., 2018; De Oliveira et al., 2019b; Landers-Ramos et al., 2024), reperfusion rate (George et al., 2018; Rosenberry et al., 2018; Iannetta et al., 2019; De Oliveira et al., 2019b; Rosenberry et al., 2019), VOT-min (Rosenberry et al., 2018; Iannetta et al., 2019; De Oliveira et al., 2019b; Landers-Ramos et al., 2024), and HP max (Costes et al., 1999; Rosenberry et al., 2018; Iannetta et al., 2019; De Oliveira et al., 2019b; Landers-Ramos et al., 2024). Effect sizes were calculated based on these data, and using an alpha level of 0.05 and a power of 0.90, the estimated required sample size ranged from 6 to 14 participants per group, assuming independent samples and equal group sizes.
All variables were tested for normal distribution using the Shapiro-Wilk test. Apart from some variables related to the leg (ATT, mVO2 and RR), all variables were normally distributed. When tested per group, all VL variables apart from mVO2 in HO and RR in HY were normally distributed. To compare group mean values for physical characteristics and performance variables, an independent samples t-test was used for all normally distributed variables, while a non-parametric Mann-Whitney U test was used for variables that were not normally distributed. Self-reported training frequency, duration and intensity, as well as computed PA minutes were compared between the two groups using Student’s t-test and cross-tabulation chi-square tests. For (simple) regression analysis, Pearson correlation was used for normally distributed variables, while Spearman correlation was used when normality was not met. Investigating the effect of aging on mVO2 and RR was done using a combination of ANOVA and multiple linear regression analysis to provide complementary insights and to strengthen the robustness of our analysis. The regression analysis aimed to identify individual-level predictors of mVO2 and RR, and to assess the relative contributions of age, body composition, and physical performance. This approach complements group comparisons by revealing relationships that may be masked in categorical analyses.
To identify significant differences between our predefined age groups, NIRS data were analyzed using a 2 × 2 repeated measures ANOVA with group (HY and HO) as the between-subject factor and muscle (FDS and VL) as the within-subject factor. Due to the uneven distribution of sex across both age groups, sex was included as an additional between-subjects factor in the repeated measures ANOVA to account for potential sex-related differences in muscle oxygenation. When assumption of sphericity was violated according to Mauchly’s test of sphericity, the Greenhouse-Geisser correction was used. When the interactions were significant, post hoc analyses (Tukey) were applied to test between-group effects, whereas Student’s t-test for paired samples was used to test changes between the two muscles. To evaluate the practical significance of the results, effect sizes (ES) for t-tests and ANOVA’s were calculated using Hedges’s g. Interpretation of the absolute Hedges’s g values was as follows: trivial >0.1, small ES = 0.2, medium ES = 0.5, large ES = 0.8, very large ES = 1.3.
Multiple linear regression with forward selection was used to investigate to what degree mVO2 and RR were affected by age, body composition (BC), and physical performance (PP). Due to the limited sample size, only 3 predictor variables were considered for each model. Apart from age, one variable for BC and one variable for PP were included as predictors. The NIRS variables (mVO2 and RR) were defined as dependent variables. Candidate variables for BC were whole-body and segmental percentages of body fat (PBF, PBF-RA and PBF-RL), as well as ATT (ATT-FDS and ATT-VL). Candidate variables for PP were handgrip strength (HGS) and peak dynamic handgrip load (DHGpeak) for handgrip performance, and peak WR (WRpeak) and peak VO2 (VO2peak) for cycling performance. Simple linear regression analysis using either Pearson or Spearman correlation was used to select the BC and PP variables that had the strongest correlation with the dependent variables. These BC and PP variables were subsequently entered into the multiple linear regression model using a forward selection method, allowing only the most significant predictor(s) to be included in each model. Model variance was reported using adjusted R2 values from the multiple linear regression analyses, reflecting the proportion of variance explained by the final model. Assumptions of linearity, homoscedasticity, normality of residuals, independence of errors, and absence of multicollinearity were evaluated using residual plots, Q-Q plots, Durbin-Watson statistics, and variance inflation factors (VIF).
All values are presented as mean ± standard deviation (SD) and the level of significance was set to P < 0.05. The analyses were done in matlab (R2022b) and SPSS [29.0.2.0 (20)].
3 Results
Of the 39 participants that were included in the study, two of the older participants dropped out after day 1 due to health issues and were therefore excluded from the data set. The data of one younger participant was excluded from data analysis due to poor quality in all NIRS signals. In addition, NIRS data for FDS was partly excluded in 5 participants (3 in HY, 2 in HO), due to poor quality of the NIRS signals. All participants tolerated the 10-min vascular occlusion without interruption or premature termination. The physical characteristics of the 36 participants that fulfilled all tests are shown in Table 1. As expected, the two groups varied significantly in terms of age, covering an age span of about 44 years. No differences with respect to height, weight, PBF, and segmental PBF were found. BMI was slightly higher in the older participants, while ATT for FDS and VL were significantly lower in the older group as compared to the younger group. The range of ATT values for the arm (FDS) were 1.2–4.1 mm in HO and 1.8–5.0 mm in HY, while for the leg (VL), ATT varied between 1.5–9.1 mm in HO and 3.1–17.0 mm in HY. In addition, ATT was significantly lower in male participants compared to female participants in both muscles (P < 0.001 for both), with values of 2.3 ± 0.9 mm vs. 3.4 ± 0.7 mm in FDS and 4.4 ± 1.4 mm vs. 8.7 ± 2.9 mm in VL, respectively. All participants were considered to be active based on their self-reported physical activity level. There was a medium to large group effect for activity level (PA), where HO reported significantly more minutes of weekly physical activity than HY (Table 1). Although there were no differences in self-reported frequency, intensity and duration of the exercise sessions (all P (χ2) > 0.18), HO tended to exercise slightly more often and slightly longer than HY, whereas exercise intensity was similar in both groups. Consequently, this contributed to a higher PA in HO since PA was calculated as a multiplication of exercise frequency, intensity and duration (Vie et al., 2023). Resting blood pressure (SBP, DBP and MAP) was significant higher in HO compared to HY (Table 1) with MAP being 103.5 ± 11.3 mmHg in HO and 87.8 ± 11.3 mmHg in HY (P < 0.001, Hedges’ g = −1.56). Performance variables from the incremental handgrip and cycling tests are presented in Table 2. No significant differences in handgrip performance were found as both handgrip strength (HGS) and peak dynamic handgrip force (DHGpeak) were similar in both groups. However, cycling performance was different between the two groups as peak WR and peak VO2 were both lower in HO as compared to HY. Peak HR was, as expected, also lower in HO.
3.1 Group differences in NIRS variables
Repeated measures ANOVA was used to identify possible differences in the NIRS variables between our two predefined age groups. Baseline SmO2 values showed a significant main effect of group [F (1,32) = 26.6, P < 0.001, ES = 0.45] and muscle [F (1,32) = 177.3, P < 0.001, ES = 0.85], with no interaction effect [F (1,32) = 0.00, P = 0.95, ES = 0.00], indicating that baseline values for SmO2 prior to VOT10 were significantly lower in HO compared to HY and significantly lower in FDS as compared to VL. SmO2 at rest was 62.4% ± 1.6% and 66.8% ± 1.5% for FDS in HO and HY, while 74.6% ± 2.0% and 78.9% ± 2.0% for VL in HO and HY, respectively. No main effect of sex was found for baseline SmO2 [F (1,32) = 3.5, P = 0.072, ES = 0.10], but a significant interaction of muscle x sex was found [F (1,32) = 14.3, P < 0.001, ES = 0.31], indicating that baseline SmO2 values were similar for male and female participants in FDS (P = 0.087) while higher in females for VL (P = 0.001).
The rescaled NIRS responses during and after the 10-min arterial occlusion (VOT10), are shown in Figure 3 while the quantitative values for the variables obtained from the signals during VOT10 are shown in Figure 4. Figures 3A,B show that the mean responses during VOT10 for each muscle were reasonable similar in both groups where desaturation and resaturation responses were more blunted in the leg (VL) compared to the arm (FDS). This was confirmed by statistics done on the variables extracted from the vascular occlusion test. No significant between-group differences were found for any of the variables (VOT-min: F (1,27) = 3.2, P = 0.087, ES = 0.11; HP-max; F (1,27) = 2.2, P = 0.150, ES = 0.08; mVO2; F (1,27) = 2.4, P = 0.135, ES = 0.08; RR; F (1,27) = 1.6, P = 0.216, ES = 0.06, indicating that the overall NIRS values derived from the vascular occlusion test were similar for young and older participants (Figure 4).
Figure 3. Group responses for simultaneously measured muscle oxygen saturation (SmO2) during and following 10 min of arterial occlusion in (A) flexor digitorum superficialis (FDS) and (B) vastus lateralis (VL) for healthy young (●; HY) and healthy older (◯; HO) participants.
Figure 4. Group mean values for (A,E) muscle oxygen consumption (mVO2) calculated from the initial decrease in SmO2 during occlusion, (B,F) maximum desaturation levels at the end of occlusion (VOT-min), (C,G) reperfusion rate (RR) calculated from the initial increase in SmO2 following cuff release, and (D,H) maximum resaturation levels during hyperemia following cuff release (HP-max), in flexor digitorum superficialis (FDS) and vastus lateralis (VL).
Apart from RR (F (1,27) = 2.7, P = 0.109, ES = 0.09), a main effect of muscle was found for VOT-min (F (1,27) = 97.2, P < 0.001, ES = 0.78), HP-max (F (1,27) = 76.7, P < 0.001, ES = 0.74) and mVO2 (F (1,27) = 60.2, P < 0.001, ES = 0.69), indicating a more pronounced desaturation, resaturation and a higher mVO2 in FDS compared to VL. The interaction effect of muscle x group was significant for VOT-min (F (1,27) = 4.8, P = 0.037, ES = 0.15), but not for HP-max (F (1,27) = 3.3, P = 0.082, ES = 0.11), mVO2 (F (1,27) = 3.7, P = 0.065, ES = 0.12), or RR (F (1,27) = 0.01, P = 0.923, ES = 0.00), suggesting that the found main effect of muscle was similar for both age groups, apart from VOT-min that showed a more pronounced difference (P = 0.017) between FDS and VL in the younger participants compared to the older participants (Figure 4).
Apart from RR (F (1,27) = 1.2, P = 0.283, ES = 0.04), a main effect of sex was found for VOT-min (F (1,27) = 24.3, P < 0.001, ES = 0.47), HP-max (F (1,27) = 31.7, P < 0.001, ES = 0.54) and mVO2 (F (1,27) = 16.1, P < 0.001, ES = 0.37), indicating a more pronounced desaturation, resaturation and higher mVO2 in male compared to female participants. A significant interaction effect of muscle x sex was found for all variables (VOT-min: F (1,27) = 5.5, P = 0.027, ES = 0.17; HP-max: F (1,27) = 4.5, P = 0.043, ES = 0.14; mVO2: F (1,27) = 8.2, P = 0.008, ES = 0.23; RR: F (1,27) = 13.2, P = 0.001, ES = 0.33), indicating that sex influenced the muscle-specific responses during and after vascular occlusion.
No significant three-way interaction (muscle × group × sex) was found for any of the variables. The regression fit (R2) for calculation of mVO2 over 120 s was 0.994 ± 0.01 and 0.992 ± 0.01 in FDS and VL, respectively, while for RR calculated over 10 s, it was 0.974 ± 0.04 and 0.981 ± 0.02, respectively.
3.2 Influence of age, body composition and physical fitness on mVO2 and RR
To investigate to what degree mVO2 and RR were affected by age, body composition (BC) and physical performance (PP), multiple linear regression was used. Prior to building these models, simple linear regression was performed to identify the BC and PP variables that most strongly correlated with each outcome (Table 3). All models met the assumptions for multiple linear regression, as indicated by residual diagnostics, and no collinearity was found between the independent predictor variables (VIF range: 1.0–2.7). Of the BC variables for the arm, the percentage of body fat in the arm (PBF-RA) showed the best correlation with both dependent variables (mVO2 and RR), and handgrip strength (HGS) for the PP variables. For the leg, ATT at VL showed the strongest correlation of the BC variables, and peak cycling WR (WRpeak) for the PP variables (Table 3).
For the arm, age, PBF-RA and HGS were entered into the regression model as candidate predictor variables, while mVO2 in FDS, as measured by NIRS SmO2, was entered as the outcome variable. A similar model was used for RR in FDS. A summary of the model fit statistics is shown in Table 4. Using the forward selection method, age emerged as the only significant predictor of mVO2 in FDS (Table 5), indicating a moderate and meaningful association and accounting for 20.3% of the variance in mVO2 in FDS (Table 4; Figure 5A). Neither PBF-RA nor HGS were included in the final model (Table 5). For RR in FDS, PBF-RA was selected as the only significant predictor and emerged as a moderate and meaningful predictor, accounting for 13.3% of the variance in RR in FDS, while age and HGS were excluded (Table 5).
Table 5. Regression coefficients for predictors tested and retained in multiple linear regression models.
Figure 5. Estimated marginal means from multiple linear regression. (A) Expected decline in mVO2 in FDS with increasing age. (B) Expected decline in in mVO2 in VL with increasing ATT. mVO2 is presented as the absolute slope in SmO2 during the initial 120 s of vascular occlusion.
For the leg, age, ATT-VL and WRpeak were entered into the regression model as candidate predictor variables, while mVO2 in VL was entered as the outcome variable. A similar model was used for RR in VL. The model fit statistics for these analyses are summarized in Table 4. Using the forward selection method, ATT-VL emerged as the only significant predictor of mVO2 in VL (Table 5), indicating a strong association and accounting for 53.9% of the variance in mVO2 in VL (Table 4; Figure 5B). Neither age nor peak cycling work rate were included in the final model (Table 5). ATT-VL also emerged as the only significant predictor of RR in VL, accounting for 9.9% of the variance, while age and WRpeak were excluded (Table 5). Substituting PBF-RA with ATT-FDS in the arm model and WRpeak with VO2peak in the leg model did not alter the results (data not shown).
3.3 Self-reported physical activity vs. lab-based performance variables
As might be expected, strong correlations were found between peak work rate and peak VO2 (rho = 0.84, P < 0.001) as well as between handgrip strength and maximal dynamic handgrip load (rho = 0.73, P < 0.001). However, none of these objectively measured performance variables correlated with the PA index quantified from self-reported exercise frequency, intensity and duration according to Vie et al. (2023), i.e., between PA index and peak cycling work rate (rho = 0.14, P = 0.423), peak oxygen uptake (rho = 0.11, P = 0.518), handgrip strength (rho = 0.03, P = 0.850) or maximum handgrip load (rho = 0.29, P = 0.088).
4 Discussion
The main purpose of the study was to investigate possible aging effects on microvascular and mitochondrial function in healthy typically-aging older adults. The key findings were: (1) baseline muscle oxygen saturation was lower in older adults compared to younger adults; (2) no significant group differences were observed in muscle oxygen consumption (mVO2), reperfusion rate (RR), or desaturation/resaturation responses; and (3) regression analysis revealed that age was a moderate predictor of forearm mVO2, while ATT was the strongest predictor of leg mVO2. These results suggest that age-related changes in mVO2 may be muscle-specific and influenced by local adiposity.
While no direct effect of age was found on the desaturation and resaturation responses, the interaction effects indicated that the relationship between age and several NIRS-derived variables (mVO2, VOT-min, HP-max) differed between forearm (FDS) and leg (VL) muscles. The multiple linear regression analysis further supported this muscle-specific difference. In the arm, age emerged as the only significant predictor of mVO2 in FDS, suggesting a moderate and meaningful age-related effect when ATT was relatively low. In contrast, in the leg, ATT was the strongest predictor of mVO2, explaining over half of the variance, while age was not included in the model. This suggests that in regions with greater variation in adiposity, such as the leg, ATT may mask age-related effects, whereas in regions with less adiposity, such as the arm, these effects may be more detectable.
4.1 Agreement with previous NIRS-based findings
Over the past few decades, near-infrared spectroscopy (NIRS) has been explored as a potential tool for investigating age-related changes in muscle function. A number of quantitative variables derived during and/or after vascular occlusion have been reported for young and older adults, starting in the ´90s with the studies of McCully et al. (McCully et al., 1991), and Costes et al. (Costes et al., 1999), reporting age-related alterations in the recovery phase following vascular occlusion. These initial findings laid the groundwork for the use of NIRS to assess changes in muscle physiology associated with aging and especially in the last 10 years, interest in this approach has grown. To date, approximately thirteen studies have used NIRS to compare young and older adults, ten of which incorporated vascular occlusion protocols and reported quantitative results obtained during and/or after occlusion. Of these ten studies, approximately half used continuous-wave NIRS devices and half used frequency-domain NIRS. In this study, we used continuous-wave NIRS devices, which estimate relative changes in tissue oxygenation. Although this study used continuous-wave NIRS, both ATT and physical fitness have been shown to influence NIRS-derived variables across modalities. Cross-modality comparisons, therefore, remain relevant for interpreting muscle function. Among the most commonly reported outcomes are recovery-related parameters, such as the recovery rate and the maximum level of resaturation or reoxygenation observed during the hyperemic response following occlusion release. In contrast to most studies that use a 5-min occlusion period, we applied a 10-min occlusion to ensure a distinct plateau in the NIRS signal. Although we did not systematically compare different durations, previous literature suggests that longer occlusions may enhance the magnitude of desaturation and hyperemic responses. This may have contributed to more pronounced signal characteristics observed in our data and should be considered when comparing results across studies.
4.1.1 Baseline SmO2 values
Baseline SmO2 (BL) values found in this study were comparable with the values of around 70% that were reported in other studies investigating NIRS derived skeletal muscle saturation in healthy young and older adults (Rosenberry et al., 2018; Rosenberry et al., 2019; Iannetta et al., 2019; De Oliveira et al., 2019a) and in older adults that were at high risk for cardiovascular disease cardiovascular disease (De Oliveira et al., 2019b; De Oliveira et al., 2019a). We found that baseline saturation was slightly lower (4.7%) in our older participants as compared to our younger participants. A similar tendency, though not significant, was seen in other studies as well (Hart et al., 2014; Iannetta et al., 2019; Rosenberry et al., 2019). We also found that baseline saturation levels differed between the two muscles, with lower values in the forearm (FDS) as compared to the thigh (VL). While the moderate group difference in (ES = 0.49) indicated a meaningful difference in baseline saturation between young and older adults, a strong difference (ES = 0.80) was found between arm and leg muscles, indicating that substantial differences in baseline saturation can be expected when comparing different muscles.
4.1.2 Muscle oxygen consumption
Muscle oxygen consumption (mVO2) is typically defined as the initial linear decline in oxygenation or saturation during arterial occlusion. Whereas in earlier days mVO2 was estimated from changes in oxyhemoglobin (+myoglobin), deoxyhemoglobin (+myoglobin) or the difference between the two (HbDiff), recent studies predominantly report the slope of the initial decrease in tissue/muscle saturation (in %/s). Of the ten studies employing vascular occlusion to assess age-related changes in NIRS-derived muscle function, only four reported resting values for mVO2, either as deoxygenation/desaturation slopes or as slope-derived estimates. The magnitude of mVO2 values that we found during rest were comparable to the values that were reported in the literature. In agreement with other studies (De Oliveira et al., 2019a; Landers-Ramos et al., 2024), we did not find differences in mVO2 between our young and older healthy adults. This in comparison to the lower mVO2 values found in older adults by Rosenberry et al. (2018) and in older adults with high risk for cardiovascular disease (De Oliveira et al., 2019b; De Oliveira et al., 2019a). Unfortunately, two other studies that investigated mVO2 recovery kinetics in healthy young and older adults (Chung et al., 2018; Lagerwaard et al., 2020), did not report resting values for mVO2 for the various muscles that were measured.
Independent of age, we found a difference in mVO2 between the two muscles with higher values in FDS as compared to VL. This difference might be due to intrinsic properties of the two muscles, as they typically vary in muscle mass, fiber type composition, capillary density, and mitochondrial content, and these muscle-specific characteristics likely affect their metabolic turnover. However, the overall higher ATT in the VL (6.5 ± 2.7 mm) as compared to FDS (2.8 ± 0.9 mm) cannot be ruled out as another likely contributor to the differences in mVO2 between the two muscles as ATT has been a known confounder for mVO2 measured by NIRS (Beekvelt et al., 2001).
4.1.3 Maximum desaturation
Maximum desaturation (VOT-min) during the 10-min vascular occlusion was similar in both groups for both FDS and VL. This is in agreement with the studies of Iannetta et al. (2019) and De Oliveira et al. (2019a) where maximum desaturation was neither found to be different between young and healthy older adults, though desaturation was less in older adults at risk for cardiovascular disease (De Oliveira et al., 2019a). A less pronounced desaturation in older adults was also found by Rosenberry et al. (2019). The maximum desaturation values that we found were higher than those reported in other studies, but this is likely attributed to the occlusion duration of 10 min that we used, in comparison to the 5-min occlusion that is most often used. This is supported by Iannetta et al. (2019) who found an effect of occlusion duration on the magnitude of desaturation, independent of age. The longer occlusion duration that we used was done to ensure a distinct plateau in NIRS signals and based on our experience that 5 min often is not enough. The progressive increase in desaturation that was shown by Iannetta et al. (2019) occurred at least up to an occlusion duration of 5 min. To what extend a longer occlusion will lead to even more desaturation remains unclear, but our data suggests that maximal desaturation might be even more pronounced after 10 min of occlusion. Although we did not find any differences between groups, we did find that maximum desaturation differed between the two muscles with more pronounced desaturation in the arm (FDS) as compared to the leg (VL). To what extend these differences are due to muscle characteristics or rather to the differences in ATT remains unclear.
4.1.4 Reperfusion rate
Reperfusion rate (RR), sometimes simply called “slope 2”, is thought to reflect microvascular function and is typically calculated as the resaturation or reoxygenation slope during the first 10 seconds after the release of the vascular occlusion. RR has often been used in NIRS studies with focus on microvascular function, e.g., in septic patients (Creteur et al., 2007; Doerschug et al., 2007; Skarda et al., 2007; Nanas et al., 2008; 2009) and has shown strong reliability (McLay et al., 2016b; Iannetta et al., 2019; McLay et al., 2016a; Iannetta et al., 2019). RR has also shown a moderate correlation with flow-mediated dilation (FMD) (McLay et al., 2016a; Soares et al., 2019), a noninvasive measure predicative for cardiovascular disease. The methods are relatively similar in that both RR and FMD measure the vascular response following a cuff-induced transient period of ischemia. However, FMD is measured using ultrasound instead of NIRS and thereby reflects macrovascular function instead of microvascular function.
Typical values for RR in arm or leg range between 0.75%/s and 2.5%/s, though values seem to vary depending on the muscle of interest. No differences in RR following 10 min of arterial occlusion were found between our young and older adults. Our results are in agreement with other studies that did not find differences in RR between healthy young and older adults (George et al., 2018; De Oliveira et al., 2019a; Rosenberry et al., 2019). However, slower forearm recovery rates in healthy older adults have been reported immediately following low intensity dynamic handgrip exercise (Kutsuzawa et al., 2001) and following arterial occlusion (Rosenberry et al., 2018). In the two studies of De Oliveira et al., (De Oliveira et al., 2019b; De Oliveira et al., 2019a), slower recovery rates were found for older adults with high risk for cardiovascular disease, though no differences were found between young and healthy older adults (De Oliveira et al., 2019a). Overall, these results imply that healthy older adults show a tendency towards slower recovery rates, though not enough to reach significance levels. A sound conclusion is further compromised by other factors that might affect the hyperemic response, such as occlusion duration, fitness level of the participants, and skinfold thickness, as will be discussed later.
4.1.5 Maximum resaturation
Similar values for maximum resaturation during post-occlusion reactive hyperemia (HP-max) were found for both groups in this study, for forearm (FDS) and thigh (VL). Though, as for the other variables, HP-max was more pronounced in the arm (FDS) as compared to the leg (VL). The lack of a difference between healthy young and healthy older adults was in agreement with Oliveira et al. (De Oliveira et al., 2019a), and Landers-Ramos et al., (Landers-Ramos et al., 2024). However, Rosenberry et al. (Rosenberry et al., 2019), did find a difference between young and older adults. Although less pronounced than for maximum desaturation and reperfusion rate, HP-max showed a dependency with occlusion duration as well (Iannetta et al., 2019). This was true for healthy young and healthy older adults, and was again, independent of age (Iannetta et al., 2019).
4.2 Confounding effect of adipose tissue thickness
The thickness of the subcutaneous adipose layer is a known confounding factor for NIRS measurements (McCully and Hamaoka, 2000; Beekvelt et al., 2001; Barstow, 2019; Hendrick et al., 2024), as NIR light is assumed to penetrate to a depth of approximately half the source-detector distance (Cui et al., 1991; Homma et al., 1996). With a 35 mm source-detector separation, as used in this study, and an adipose tissue thickness (ATT) of 17.5 mm or more, none of the NIR light is thought to reach the muscle but will instead travel through the adipose tissue. Given the lower metabolic activity of adipose tissue compared to muscle, this leads to attenuated oxygenation responses and underestimation of quantitative NIRS values (Beekvelt et al., 2001; Bopp et al., 2011; Bowen et al., 2013).
In this study, ATT values for the arm (FDS) ranged between 1.2–5.0 mm, while ATT for the leg (VL) ranged between 1.5–17.0 mm. Our older participants had significantly less subcutaneous adipose tissue than our younger participants, and this was most pronounced in the leg. Although we recognize that it would be best to keep ATT below 25% of the source-detector distance according to the recommendations of Barstow (Barstow, 2019), this might be difficult to achieve in a non-athletic general population, especially for leg muscles. Instead of the recommended 25% that would have resulted in an ATT threshold of 8.8 mm and consequently excluding seven, mainly young and all female participants, we chose to rely on signal quality. This resulted in the exclusion of 1 younger participant, as described in the methods, while all other participants were well within the 95% CI for the NIRS variables and including or excluding these participants did not affect the results.
The saturation responses during and following vascular occlusion for both muscles and both groups are shown in Figure 3. Compared to the younger group, our older group showed a less pronounced response in the arm (FDS) while the assumed propagation of light traveling through the muscle was roughly similar (∼81% and ∼85% in HY and HO, respectively) (Figure 3A). In the leg (VL), however, group responses of older and younger adults were roughly similar but with assumably substantially less light traveling through the muscle in the younger participants (∼54% and ∼72% in HY and HO, respectively) due to the higher ATT values in this group (Figure 3B). It is, therefore, likely that the higher ATT values in the younger group masks to some extent the true effects of age between these groups, at least in the leg (VL).
Multiple linear regression showed that ATT on top of the VL was a strong predictor for mVO2 in VL, but not for RR. Using simple linear regression, strong significant correlations were found between ATT on top of the VL and mVO2 (rho = 0.80, P < 0.001), changes in desaturation at the end of the occlusion (VOT-min; rho = 0.83, P < 0.001) and resaturation during the hyperemic phase (HP-max; rho = −0.84, P < 0.001), as well as a moderate significant correlation between ATT-VL and RR (rho = −0.63, P < 0.001). These correlations indicate slower rates for mVO2 and RR and less pronounced changes in desaturation (VOT-min) and resaturation (HP-max) with higher ATT. Thus, ATT did affect the NIRS measurements in the leg and may have attenuated the age-related changes in mitochondrial and microvascular function in the leg. This was not the case for the measurements in the arm, most likely due to a narrower range of ATT values and the smaller differences in ATT between both groups (Δ 0.7 mm in FDS vs. Δ 3.1 mm in VL), thereby not leading to any significant correlations between ATT and measured NIRS variables in FDS. Moreover, it enabled age to emerge as a moderate predictor of mVO2 in the arm independent of the other predictors in the multiple linear regression model.
Whereas ATT was lower in older adults in both arm and leg, we did not find any differences in the estimated segmental body fat percentages (Table 1: PBF-RA and PBF-RL). With similar values for segmental fat, but lower values for ATT, one could speculate that the ratio of fat within the tissue might have been higher in the older group compared to the younger group. An age-related body fat redistribution from subcutaneous sites to other tissues such as, e.g., muscle tissue, has been reported previously (Farsijani et al., 2020) and is thought to often occur in humans as well as in aged animals (Wang et al., 2014). This difference in fat distribution with more subcutaneous fat in our younger adults and possibly more intramuscular fat in our older adults might have affected the NIRS measurements differently for the two age groups. Although it is tempting to assume similar confounding effects of inter/intramuscular fat deposits on the NIRS measurements as we know for subcutaneous fat, the effect is, as far as we know, unknown.
In addition to differences in ATT between the age groups, ATT also differed significantly between male and female participants. Due to the uneven sex distribution between the young and older groups in our study, it remains unclear whether the observed sex effects in NIRS-derived variables were primarily driven by differences in ATT or by other biological differences between sexes. However, no significant three-way interaction (muscle x group x sex) was found, suggesting that sex did not differentially affect age-related responses across muscles.
Although the effect of ATT on NIRS measurements has been clearly emphasized in review papers (Perrey and Ferrari, 2018), and despite the relative ease of the measurement, ATT is, unfortunately, still frequently omitted in published NIRS studies, particularly in the context of aging and muscle function. Of the nine studies that used vascular occlusion for comparison between healthy young and healthy older adults, only three studies reported muscle specific ATT values (De Oliveira et al., 2019a; Lagerwaard et al., 2020; Landers-Ramos et al., 2024). These studies reported that apart from a group difference in ATT for the gastrocnemius muscle (Lagerwaard et al., 2020), no differences in ATT were found for the other leg muscles (Lagerwaard et al., 2020; Landers-Ramos et al., 2024) or arm muscle (De Oliveira et al., 2019a). The study of De Oliveira et al. (De Oliveira et al., 2019a), reported the same NIRS variables as the current study, and no differences in any of the variables were found between young and older healthy adults. Landers-Ramos et al. (Landers-Ramos et al., 2024), neither found differences between young and older healthy adults for mVO2, VOT-min and HP-max, but a significantly faster half-recovery time in older adults. An opposite slower half-recovery time was found by Lagerwaard et al. (Lagerwaard et al., 2020), in the tibialis and vastus lateralis muscles, while half-recovery time was similar in the gastrocnemius. In both latter studies vascular occlusion was preceded by voluntary exercise or electrostimulation.
4.3 Influence of physical fitness level
One factor that complicates the investigation of age-related effects on muscle function is the coinciding tendency for physical activity to decline with age. This decline can lead to reductions in physical fitness, which in turn affects physical performance. Although these concepts are interrelated, they are not synonymous. In this study, physical activity was assessed via validated self-report questionnaires, while physical fitness and physical performance were evaluated through graded cycling and handgrip tests. This distinction is important for interpreting the results, as it allows us to consider both habitual activity patterns and objective measures of functional capacity.
Self-reported PA measures, as an indirect indicator for fitness level, are typically used in large population studies where it is a cost-effective tool (Sun et al., 2013) to assess activity patterns in large groups. However, it is also known to be questionably in terms of its accuracy on individual level and might be less feasible in older adults due to factors like fluctuations in health status, cognitive problems, and the idea that older adults tend to engage most frequently in light to moderate intensity activities (Garatachea et al., 2009; Kowalski et al., 2012). Self-reported physical activity in this study was merely used in combination with physical performance measurements to obtain a better view on the overall training status of our population. PA was quantified according to the method of Vie et al. (2023) and was based on three questions related to frequency, intensity and duration of leisure time physical activity. According to the definition of Vie et al. (2023), all participants were physically active, scoring a PA index ≥1, corresponding to 1 h or more of PA per week. With respect to activity minutes per week, our older adults were more active (Table 1) than our younger participants and trained at similar self-reported intensity.
Other indicators for fitness level were measured in the lab during physical tests. Physical performance variables were measured by means of objective values for peak work rate and peak oxygen consumption from the incremental cycling test as well as handgrip strength and maximum dynamic handgrip load from the incremental handgrip test. As expected, we found strong correlations between peak work rate and peak VO2 (rho = 0.84, P < 0.001) and between handgrip strength and maximum dynamic handgrip load (rho = 0.73, P < 0.001). However, no correlations were found between the quantified PA index obtained from the self-reported training habits and any of the measured performance values, i.e., peak cycling work rate (rho = 0.14, P = 0.423), peak oxygen uptake (rho = 0.11, P = 0.518), handgrip strength (rho = 0.03, P = 0.850) or maximum dynamic handgrip load (rho = 0.29, P = 0.088).
Concerning the overall physical fitness level, we can conclude that our older population was healthy and fit, and engaged in more activity minutes on a weekly basis as compared to our younger group. Handgrip strength (HGS) was similar in both age-groups with mean HGS values for younger (38.0 ± 6.2 kg) and older women (28.1 ± 3.3 kg) that were well above the centile reference values reported by Dodds et al. (2014). Mean HGS for our older men (49.1 ± 4.5 kg) was also well above their age-comparable reference values, while HGS for our younger men (51.4 ± 6.2 kg) was roughly similar to that of the age-comparable reference values (Dodds et al., 2014). Contrary to the performance variables for arm exercise (i.e., handgrip strength and incremental handgrip load), physical performance during whole body exercise (cycling) was less for our older participants as compared to our younger group, indicated by the lower peak work rate and VO2max. These findings underscore the challenge of isolating age-related effects from physical fitness, given the multidimensional nature of fitness and its limited correspondence with laboratory-measured physical variables.
With respect to the use of NIRS investigating age-related effects on muscle function, studies addressing the effect of physical fitness in addition to that of age are scarce. George et al. (2018) showed that the magnitude of the hyperemic response following vascular occlusion was positively related to training level in both young and older adults, where endurance trained men had significantly faster recovery rates than both inactive and recreational active men regardless of age. This is supported by the results of Hart et al. (2014) who neither found an effect of age between young and older participants that were matched for physical activity. Although Hart et al. (2014) reported the time constant (Tc) of recovery instead of RR, this is an established method to characterize exercise recovery, especially in studies using magnetic resonance spectroscopy. Both phosphocreatine and NIRS deoxyhemoglobin recovery were similar in their young and older adults (Hart et al., 2014). An age-related prolongation of half-recovery time following an incremental cycling test was reported by Ichimura et al. (2006). However, in line with the above findings, this age-associated slowing of recovery was attenuated in physically active middle-aged and older women (40–80 years) compared to their sedentary counterparts, suggesting that sedentary behavior may be associated with less efficient recovery from exercise.
Two other studies (Kutsuzawa et al., 2001; Rosenberry et al., 2018) reported slower recovery rates in older adults following vascular occlusion, although the extent to which physical fitness influenced these outcomes remains unclear. In the study of Kutsuzawa et al. (2001), both young and older participants were self-reported sedentary, though older adults had ∼16% lower grip strength. Similarly, Rosenberry et al. (2018) reported ∼47% lower grip strength and 33%–57% fewer weekly physical activity minutes for their older participants, indicating differences in physical fitness between groups. De Oliveira et al. (2019a) found a non-significant moderate age effect on recovery rate in sedentary adults. However, data on physical fitness or handgrip strength is lacking. In our study, young and older participants had comparable grip strength, dynamic handgrip performance, and recovery rates in the arm. However, our older adults showed lower peak cycling WR and VO2, indicating reduced overall fitness, yet recovery rates in the leg did not differ by age. To date, no additional data are available, and the influence of physical fitness on other NIRS-derived parameters remains unclear.
4.4 Complementary insights from group comparisons and regression analysis
To gain a comprehensive understanding of how aging affects muscle function, this study combined group comparisons with multiple linear regression analysis. While group comparisons revealed no significant differences in vascular occlusion responses between young and older adults, regression analysis revealed nuanced relationships between age and mVO2 that varied by muscle group. Whereas age emerged as a moderate predictor of mVO2 in the forearm (FDS), explaining 20.3% of the variance, ATT was the strongest predictor in the leg (VL), accounting for 53.9% of the variance. This dual approach allowed us to identify broad group effects while also detecting individual predictors that might not emerge from group comparisons alone. These findings suggest that in muscle regions with relatively low adiposity and minimal variation, such as the arm, age-related effects on mVO2 may be more readily detectable in comparison to regions with higher adiposity, such as the leg. Although this study focused on comparing age groups, the broader goal is to understand how age and muscle function relate across the full age spectrum. To support this aim, we complemented group comparisons with multiple linear regression analysis. Given the relatively small sample size, the use of forward selection, and the categorical nature of the age groups, the regression analyses should be considered exploratory and interpreted with caution.
4.5 Limitations
This study focused on healthy older adults to investigate physiological changes associated with normal aging. While this provides valuable insight into typical aging trajectories, future research should include a broader range of health statuses and a more continuous age spectrum to better capture the full complexity of aging. Although the confounding effect of ATT on NIRS measurements is well-established, many studies still fail to report ATT, limiting interpretability. In our study, ATT strongly affected NIRS measurements, particularly in the leg where ATT values were higher and showed greater variability than in the arm. Our findings reinforce the importance of routinely measuring and reporting ATT, not only in the field of muscle aging, but in all NIRS studies. Moreover, while subcutaneous fat is easily measured and known to confound NIRS measurements, aging is also known to coincide with increased intramuscular fat infiltration. The influence of intramuscular fat infiltration on NIRS signals remains unclear, and represents an important area for future investigation. Additionally, all skinfold measurements were taken after the final test was completed. This enabled precise measurements at the site of the NIRS measurements, but also meant that skinfold measurements were taken after the exercise tests, and although the impact of acute exercise on compressibility is considered negligible, hydration differences between groups cannot be ruled out. NIRS devices were not randomly assigned across limbs, and all tests were performed in a fixed order, which may have introduced systematic effects. Finally, while this study employed continuous-wave NIRS and findings therefore apply most directly to this modality, comparisons across NIRS types remain relevant for interpreting muscle function.
The sample size provided sufficient power for detecting group differences using ANOVA but was limited for muscle × group interactions and regression models, which were therefore exploratory. This constraint also restricted inclusion of multiple covariates such as sex. Although sex was included as a factor in ANOVA and did not affect the main outcomes, it was not included in regression models because of limited power and risk of overfitting. Although assumptions of linear regression were met, the use of forward selection in a small sample may increase the risk of overfitting, and results should be interpreted with caution. Moreover, forward selection may limit the inclusion of predictor combinations and is sensitive to variable entry order, which reinforces the exploratory nature of the regression analysis and the need for cautious interpretation. To overcome these constraints, future studies should include larger, more balanced samples, enabling the incorporation of sex and other relevant covariates in multivariable models.
4.6 Practical implications
Understanding how aging affects muscle oxygenation and vascular function requires careful consideration of both physiological and methodological factors. In the forearm, where ATT was low and relatively uniform, age emerged as a significant predictor of mVO2. In contrast, in the leg, where ATT varied more widely, ATT was the dominant predictor, possibly masking the age-related effects. This, together with the typical decline in physical activity with age, underscores the need for careful study design when investigating aging effects on muscle function using NIRS. Moreover, regardless of study design or population, ATT should be routinely measured and reported in all studies to improve data interpretation and comparability. While muscle-specific differences in NIRS responses are well-known, they are thought to be particularly relevant in aging research, as the aging process itself likely varies between muscles. More research is needed to explore these muscle-specific aging patterns. Complementary methods such as ultrasound and MRI that can assess muscle architecture, muscle adiposity, and quantify intramuscular fat and perfusion, may enhance interpretation and provide further insight into muscle aging. Moreover, given the known influence of sex on body composition and aerobic capacity, future studies should aim for balanced sex distribution or include sex as a covariate to account for potential confounding effects. Finally, more knowledge of the effects of occlusion duration and its impact on NIRS-derived variables can broaden our understanding of the physiological responses to vascular occlusion in general, and muscle aging in particular.
5 Conclusion
Investigating the effect of aging on muscle oxygen consumption (mVO2) and reoxygenation rate (RR) using both ANOVA and multiple linear regression analysis provided complementary insights and improved the interpretability of our findings. ANOVA was used to identify potential differences between the two predefined age groups and showed that, apart from baseline SmO2 values, no differences were found between young and older participants. In contrast, multiple linear regression analysis revealed that age was a moderate predictor of mVO2 in the forearm, explaining 20.3% of its variance, particularly in a context where body composition and performance were relatively comparable between the two age groups. In the leg, however, ATT emerged as the strongest predictor of mVO2, explaining 53.9% of its variance. Although RR did not differ significantly between age groups, it showed moderate associations with ATT in the leg and percent body fat in the arm, suggesting a potential influence of adiposity on microvascular function. These findings underscore the importance of considering individual-level factors such as ATT and physical fitness when assessing age-related changes in muscle function using continuous-wave NIRS, as these factors may confound or interact with the effects of aging, making them more difficult to detect.
Data availability statement
The datasets generated and analyzed for this study are not publicly available due to NTNU regulations and GDPR compliance. Requests to access the datasets should be directed to bWlyZWlsbGUudmFuLmJlZWt2ZWx0QG50bnUubm8=.
Ethics statement
The studies involving humans were approved by NTNU/REK midt, Faculty of Medicine and Health Sciences, PO Box 8905, 7491 Trondheim. The studies were conducted in accordance with the local legislation and institutional requirements and adhered to the principles of the Declaration of Helsinki. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from all participants for participation in the study.
Author contributions
MvB: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. MW: Conceptualization, Data curation, Methodology, Writing – review and editing. AP: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review and editing.
Funding
The authors declare that no financial support was received for the research and/or publication of this article.
Acknowledgements
The authors would like to thank Erwin Marsi for his valuable guidance in the early stages of MATLAB programming, which formed the foundation for the data analysis presented in this study, and Xiang-Chun Tan for her assistance in the laboratory. Parts of this study appear in the master’s thesis of the last author (Andreas Parviz Gaarden, 2021).
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 authors 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
Ajoolabady A., Pratico D., Ren J. (2024). Endothelial dysfunction: mechanisms and contribution to diseases. Acta Pharmacol. Sin. 45, 2023–2031. doi:10.1038/s41401-024-01295-8
Andreas Parviz Gaarden (2021). Age-related differences in muscle microvascular response and muscle oxygen uptake: a near infrared spectroscopy approach. Physical Activity and Health - Exercise Physiology. (Master’s thesis) Available online at: https://ntnuopen.ntnu.no/ntnu-xmlui/bitstream/handle/11250/2783259/no.ntnu:inspera:82525395:9098907.pdf?sequence=1
Barstow T. J. (2019). Understanding near infrared spectroscopy and its application to skeletal muscle research. J. Appl. Physiol. 126, 1360–1376. doi:10.1152/japplphysiol.00166.2018
Beekvelt M. van, Borghuis M., Engelen B. van, Wevers R., Colier W. (2001). Adipose tissue thickness affects in vivo quantitative near-IR spectroscopy in human skeletal muscle. Clin. Sci. 101, 21–28. doi:10.1042/cs20000247
Betik A. C., Hepple R. T. (2008). Determinants of VO2 max decline with aging: an integrated perspective. Appl. Physiol. Nutr. Metab. 33, 130–140. doi:10.1139/h07-174
Bloom S. I., Islam M. T., Lesniewski L. A., Donato A. J. (2023). Mechanisms and consequences of endothelial cell senescence. Nat. Rev. Cardiol. 20, 38–51. doi:10.1038/s41569-022-00739-0
Bopp C. M., Townsend D. K., Barstow T. J. (2011). Characterizing near-infrared spectroscopy responses to forearm post-occlusive reactive hyperemia in healthy subjects. Eur. J. Appl. Physiol. 111, 2753–2761. doi:10.1007/s00421-011-1898-z
Borg G. (1970). Perceived exertion as an indicator of somatic stress. Scand. J. Rehabil. Med. 2, 92–98. doi:10.2340/1650197719702239298
Bowen T. S., Rossiter H. B., Benson A. P., Amano T., Kondo N., Kowalchuk J. M., et al. (2013). Slowed oxygen uptake kinetics in hypoxia correlate with the transient peak and reduced spatial distribution of absolute skeletal muscle deoxygenation. Exp. Physiol. 98, 1585–1596. doi:10.1113/expphysiol.2013.073270
Burton L. A., Sumukadas D. (2010). Optimal management of sarcopenia. Clin. Interv. Aging 5, 217–228. doi:10.2147/cia.s11473
Choi Y., Kim D., Kim S. K. (2024). Effects of physical activity on body composition, muscle strength, and physical function in old age: bibliometric and meta-analyses. Healthcare 12, 197. doi:10.3390/healthcare12020197
Chung S., Rosenberry R., Ryan T. E., Munson M., Dombrowsky T., Park S., et al. (2018). Near-infrared spectroscopy detects age-related differences in skeletal muscle oxidative function: promising implications for geroscience. Physiol. Rep. 6, e13588. doi:10.14814/phy2.13588
Clark B. C., Manini T. M. (2008). Sarcopenia =/= dynapenia. J. Gerontol. A Biol. Sci. Med. Sci. 63, 829–834. doi:10.1093/gerona/63.8.829
Conley K. E., Jubrias S. A., Cress M. E., Esselman P. (2013). Exercise efficiency is reduced by mitochondrial uncoupling in the elderly. Exp. Physiol. 98, 768–777. doi:10.1113/expphysiol.2012.067314
Costes F., Denis C., Roche F., Prieur F., Enjolras F., Barthélémy J.-C. (1999). Age-associated alteration of muscle oxygenation measured by near infrared spectroscopy during exercise. Arch. Physiol. Biochem. 107, 159–167. doi:10.1076/apab.107.2.159.4343
Creteur J., Carollo T., Soldati G., Buchele G., Backer D. D., Vincent J.-L. (2007). The prognostic value of muscle StO2 in septic patients. Intensive Care Med. 33, 1549–1556. doi:10.1007/s00134-007-0739-3
Cruz-Jentoft A. J., Baeyens J. P., Bauer J. M., Boirie Y., Cederholm T., Landi F., et al. (2010). Sarcopenia: european consensus on definition and diagnosis: report of the european working group on sarcopenia in older people. Age Ageing 39, 412–423. doi:10.1093/ageing/afq034
Cruz-Jentoft A. J., Bahat G., Bauer J., Boirie Y., Bruyère O., Cederholm T., et al. (2018). Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing 48, 16–31. doi:10.1093/ageing/afy169
Cui W., Kumar C., Chance B. (1991). Experimental study of migration depth for the photons measured at sample surface. Time-Resolv. Spectrosc. Imaging Tissues 1431, 180–191. doi:10.1117/12.44189
De Oliveira G. D., Soares R., Volino-Souza M., Leitão R., Murias J., Alvares T. (2019a). The effects of aging and cardiovascular risk factors on microvascular function assessed by near-infrared spectroscopy. Microvasc. Res. 126, 103911. doi:10.1016/j.mvr.2019.103911
De Oliveira G. D., Soares R., Volino-Souza M., Murias J., Alvares T. (2019b). The association between near-infrared spectroscopy assessment of microvascular reactivity and flow-mediated dilation is disrupted in individuals at high risk for cardiovascular disease. Microcirculation 26, e12556. doi:10.1111/micc.12556
Delmonico M. J., Harris T. B., Visser M., Park S. W., Conroy M. B., Velasquez-Mieyer P., et al. (2009). Longitudinal study of muscle strength, quality, and adipose tissue infiltration. Am. J. Clin. Nutr. 90, 1579–1585. doi:10.3945/ajcn.2009.28047
Distefano G., Goodpaster B. H. (2018). Effects of exercise and aging on skeletal muscle. Cold Spring Harb. Perspect. Med. 8, a029785. doi:10.1101/cshperspect.a029785
Dodds R. M., Syddall H. E., Cooper R., Benzeval M., Deary I. J., Dennison E. M., et al. (2014). Grip strength across the life course: normative data from twelve British studies. PLoS ONE 9, e113637. doi:10.1371/journal.pone.0113637
Doerschug K. C., Delsing A. S., Schmidt G. A., Haynes W. G. (2007). Impairments in microvascular reactivity are related to organ failure in human sepsis. Am. J. Physiol. Heart Circ. Physiol. 293, H1065–H1071. doi:10.1152/ajpheart.01237.2006
Farsijani S., Santanasto A. J., Miljkovic I., Boudreau R. M., Goodpaster B. H., Kritchevsky S. B., et al. (2020). The relationship between intermuscular fat and physical performance is moderated by muscle area in older adults. J. Gerontol. Ser. A 76, 115–122. doi:10.1093/gerona/glaa161
Gallagher D., Visser M., Meersman R. E. D., Sepúlveda D., Baumgartner R. N., Pierson R. N., et al. (1997). Appendicular skeletal muscle mass: effects of age, gender, and ethnicity. J. Appl. Physiol. 83, 229–239. doi:10.1152/jappl.1997.83.1.229
Garatachea N., Luque G. T., Gallego J. G. (2009). Physical activity and energy expenditure measurements using accelerometers in older adults. Nutr. Hosp. 25, 224–230.
George M. A., McLay K. M., Doyle-Baker P. K., Reimer R. A., Murias J. M. (2018). Fitness level and not aging per se,, determines the oxygen uptake kinetics response. Front. Physiol. 9, 277. doi:10.3389/fphys.2018.00277
Gerovasili V., Dimopoulos S., Tzanis G., Anastasiou-Nana M., Nanas S. (2010). Utilizing the vascular occlusion technique with NIRS technology. Int. J. Ind. Ergon. 40, 218–222. doi:10.1016/j.ergon.2009.02.004
Gonzalez-Freire M., de Cabo R., Bernier M., Sollott S. J., Fabbri E., Navas P., et al. (2015). Reconsidering the role of mitochondria in aging. J. Gerontol. A Biol. Sci. Med. Sci. 70, 1334–1342. doi:10.1093/gerona/glv070
Goodpaster B. H., Park S. W., Harris T. B., Kritchevsky S. B., Nevitt M., Schwartz A. V., et al. (2006). The loss of skeletal muscle strength, mass, and quality in older adults: the health, aging and body composition study. J. Gerontol. A Biol. Sci. Med. Sci. 61, 1059–1064. doi:10.1093/gerona/61.10.1059
Grevendonk L., Connell N. J., McCrum C., Fealy C. E., Bilet L., Bruls Y. M. H., et al. (2021). Impact of aging and exercise on skeletal muscle mitochondrial capacity, energy metabolism, and physical function. Nat. Commun. 12, 4773. doi:10.1038/s41467-021-24956-2
Han Y., Kim S. Y. (2023). Endothelial senescence in vascular diseases: current understanding and future opportunities in senotherapeutics. Exp. Mol. Med. 55, 1–12. doi:10.1038/s12276-022-00906-w
Hart C. R., Layec G., Trinity J. D., Liu X., Kim S.-E., Groot H. J., et al. (2014). Evidence of preserved oxidative capacity and oxygen delivery in the plantar flexor muscles with age. J. Gerontol. A 70, 1067–1076. doi:10.1093/gerona/glu139
Hawkins S. A., Wiswell R. A. (2003). Rate and mechanism of maximal oxygen consumption decline with aging: implications for exercise training. Sports Med. 33, 877–888. doi:10.2165/00007256-200333120-00002
Hendrick E., Jamieson A., Chiesa S. T., Hughes A. D., Jones S. (2024). A short review of application of near-infrared spectroscopy (NIRS) for the assessment of microvascular post-occlusive reactive hyperaemia (PORH) in skeletal muscle. Front. Physiol. 15, 1480720. doi:10.3389/fphys.2024.1480720
Homma S., Eda H., Ogasawara S., Kagaya A. (1996). Near-infrared estimation of O2 supply and consumption in forearm muscles working at varying intensity. J. Appl. Physiol. 80, 1279–1284. doi:10.1152/jappl.1996.80.4.1279
Hortobágyi T., Vetrovsky T., Brach J. S., Haren M. van, Volesky K., Radaelli R., et al. (2023). Effects of exercise training on muscle quality in older individuals: a systematic scoping review with meta-analyses. Sports Med. Open 9, 41. doi:10.1186/s40798-023-00585-5
Iannetta D., Inglis E. C., Soares R. N., McLay K. M., Pogliaghi S., Murias J. M., et al. (2019). Reliability of microvascular responsiveness measures derived from near-infrared spectroscopy across a variety of ischemic periods in young and older individuals. Microvasc. Res. 122, 117–124. doi:10.1016/j.mvr.2018.10.001
Ichimura S., Murase N., Osada T., Kime R., Homma T., Ueda C., et al. (2006). Age and activity status affect muscle reoxygenation time after maximal cycling exercise. Med. Sci. Sports Exerc. 38, 1277–1281. doi:10.1249/01.mss.0000227312.08599.f1
Kowalski K., Rhodes R., Naylor P.-J., Tuokko H., MacDonald S. (2012). Direct and indirect measurement of physical activity in older adults: a systematic review of the literature. Int. J. Behav. Nutr. Phys. Act. 9, 148. doi:10.1186/1479-5868-9-148
Kurtze N., Rangul V., Hustvedt B.-E., Flanders W. D. (2007). Reliability and validity of self-reported physical activity in the nord-Trøndelag health study (HUNT 2). Eur. J. Epidemiol. 22, 379–387. doi:10.1007/s10654-007-9110-9
Kutsuzawa T., Shioya S., Kurita D., Haida M., Yamabayashi H. (2001). Effects of age on muscle energy metabolism and oxygenation in the forearm muscles. Med. Sci. Sports Exerc. 33, 901–906. doi:10.1097/00005768-200106000-00008
Lagerwaard B., Nieuwenhuizen A. G., Boer V. C. J. de, Keijer J. (2020). In vivo assessment of mitochondrial capacity using NIRS in locomotor muscles of young and elderly males with similar physical activity levels. Geroscience 42, 299–310. doi:10.1007/s11357-019-00145-4
Landers-Ramos R. Q., Merriman K., Chartier G. L., Parks J., Park H., Knuth N. D. (2024). Comparison of skeletal muscle oxidative function in young, middle-aged, and older males and females with similar physical activity levels. J. Appl. Physiol. 136, 618–629. doi:10.1152/japplphysiol.00545.2023
Landi F., Calvani R., Tosato M., Martone A. M., Bernabei R., Onder G., et al. (2016). Impact of physical function impairment and multimorbidity on mortality among community-living older persons with sarcopaenia: results from the ilSIRENTE prospective cohort study. BMJ Open 6, e008281. doi:10.1136/bmjopen-2015-008281
Lang J. J., Prince S. A., Merucci K., Cadenas-Sanchez C., Chaput J.-P., Fraser B. J., et al. (2024). Cardiorespiratory fitness is a strong and consistent predictor of morbidity and mortality among adults: an overview of meta-analyses representing over 20.9 million observations from 199 unique cohort studies. Br. J. Sports Med. 58, 556–566. doi:10.1136/bjsports-2023-107849
Liu C., Shiroy D. M., Jones L. Y., Clark D. O. (2014). Systematic review of functional training on muscle strength, physical functioning, and activities of daily living in older adults. Eur. Rev. Aging Phys. Act. 11, 95–106. doi:10.1007/s11556-014-0144-1
Lozano-Montoya I., Correa-Pérez A., Abraha I., Soiza R. L., Cherubini A., O’Mahony D., et al. (2017). Nonpharmacological interventions to treat physical frailty and sarcopenia in older patients: a systematic overview – the SENATOR project ONTOP Series. Clin. Interv. Aging 12, 721–740. doi:10.2147/cia.s132496
McCully K. K., Hamaoka T. (2000). Near-infrared spectroscopy: what can it tell Us about oxygen saturation in skeletal muscle? Exerc. Sport Sci. Rev. 28, 123–127.
McCully K. K., Kakihira H., Vandenborne K., Kent-Braun J. (1991). Noninvasive measurements of activity-induced changes in muscle metabolism. J. Biomech. 24, 153–161. doi:10.1016/0021-9290(91)90385-z
McLay K. M., Fontana F. Y., Nederveen J. P., Guida F. F., Paterson D. H., Pogliaghi S., et al. (2016a). Vascular responsiveness determined by near-infrared spectroscopy measures of oxygen saturation. Exp. Physiol. 101, 34–40. doi:10.1113/ep085406
McLay K. M., Nederveen J. P., Pogliaghi S., Paterson D. H., Murias J. M. (2016b). Repeatability of vascular responsiveness measures derived from near-infrared spectroscopy. Physiol. Rep. 4, e12772. doi:10.14814/phy2.12772
Mitchell W. K., Williams J., Atherton P., Larvin M., Lund J., Narici M. (2012). Sarcopenia, dynapenia, and the impact of advancing age on human skeletal muscle size and strength; a quantitative review. Front. Physiol. 3, 260. doi:10.3389/fphys.2012.00260
Nanas S., Gerovasili V., Dimopoulos S., Pierrakos C., Kourtidou S., Kaldara E., et al. (2008). Inotropic agents improve the peripheral microcirculation of patients with end-stage chronic heart failure. J. Card. Fail. 14, 400–406. doi:10.1016/j.cardfail.2008.02.001
Nanas S., Gerovasili V., Renieris P., Angelopoulos E., Poriazi M., Kritikos K., et al. (2009). Non-invasive assessment of the microcirculation in critically ill patients. Anaesth. Intensive Care 37, 733–739. doi:10.1177/0310057X0903700516
Patterson M. S., Chance B., Wilson B. C. (1989). Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties. Appl. Opt. 28, 2331–2336. doi:10.1364/AO.28.002331
Perrey S., Ferrari M. (2018). Muscle oximetry in sports science: a systematic review. Sports Med. 48, 597–616. doi:10.1007/s40279-017-0820-1
Rosenberry R., Munson M., Chung S., Samuel T. J., Patik J., Tucker W. J., et al. (2018). Age-related microvascular dysfunction: novel insight from near-infrared spectroscopy. Exp. Physiol. 103, 190–200. doi:10.1113/ep086639
Rosenberry R., Trojacek D., Chung S., Cipher D. J., Nelson M. D. (2019). Interindividual differences in the ischemic stimulus and other technical considerations when assessing reactive hyperemia. Am. J. Physiol. Regul. Integr. Comp. Physiol. 317, R530–R538. doi:10.1152/ajpregu.00157.2019
Scioli M. G., Bielli A., Arcuri G., Ferlosio A., Orlandi A. (2014). Ageing and microvasculature. Vasc. Cell 6, 19. doi:10.1186/2045-824X-6-19
Short K. R., Bigelow M. L., Kahl J., Singh R., Coenen-Schimke J., Raghavakaimal S., et al. (2005). Decline in skeletal muscle mitochondrial function with aging in humans. Proc. Natl. Acad. Sci. 102, 5618–5623. doi:10.1073/pnas.0501559102
Silva A. M., Shen W., Heo M., Gallagher D., Wang Z., Sardinha L. B., et al. (2010). Ethnicity-related skeletal muscle differences across the lifespan. Am. J. Hum. Biol. 22, 76–82. doi:10.1002/ajhb.20956
Skarda D. E., Mulier K. E., Myers D. E., Taylor J. H., Beilman G. J. (2007). Dynamic near-infrared spectroscopy measurements in patients with severe sepsis. Shock 27, 348–353. doi:10.1097/01.shk.0000239779.25775.e4
Skovereng K., Ettema G., Van Beekvelt M. (2016). Local muscle oxygen consumption related to external and joint specific power. Hum. Mov. Sci. 45, 161–171. doi:10.1016/j.humov.2015.11.009
Soares R. N., Somani Y. B., Proctor D. N., Murias J. M. (2019). The association between near-infrared spectroscopy-derived and flow-mediated dilation assessment of vascular responsiveness in the arm. Microvasc. Res. 122, 41–44. doi:10.1016/j.mvr.2018.11.005
Sun F., Norman I. J., While A. E. (2013). Physical activity in older people: a systematic review. BMC Public Health 13, 449. doi:10.1186/1471-2458-13-449
Tanaka H., Seals D. R. (2008). Endurance exercise performance in masters athletes: age-associated changes and underlying physiological mechanisms. J. Physiol. 586, 55–63. doi:10.1113/jphysiol.2007.141879
Tieland M., Trouwborst I., Clark B. C. (2018). Skeletal muscle performance and ageing. J. Cachexia Sarcopenia Muscle 9, 3–19. doi:10.1002/jcsm.12238
Vellas B., Fielding R. A., Bens C., Bernabei R., Cawthon P. M., Cederholm T., et al. (2018). Implications of ICD-10 for sarcopenia clinical practice and clinical trials: report by the international conference on frailty and sarcopenia research task force. J. Frailty Aging 7, 2–9. doi:10.14283/jfa.2017.30
Vie T. L., Hufthammer K. O., Rangul V., Andersen J. R., Meland E., Breidablik H. J. (2023). Patterns of physical activity over 34 years in a large sample of adults: the HUNT study, Norway. Scand. J. Public Health. 52, 669–677. doi:10.1177/14034948231174947
Keywords: near-infrared spectroscopy, aging, muscle oxygen consumption, microvascular function, adipose tissue thickness, physical fitness, vascular occlusion, post-occlusion reactive hyperemia
Citation: van Beekvelt M, Wilson M and Parviz Gaarden A (2025) Muscle oxygen consumption and microvascular function in healthy aging assessed by near-infrared spectroscopy: effect of adipose tissue thickness and physical fitness. Front. Physiol. 16:1659756. doi: 10.3389/fphys.2025.1659756
Received: 04 July 2025; Accepted: 11 November 2025;
Published: 05 December 2025.
Edited by:
Re Rebecca, Polytechnic University of Milan, ItalyReviewed by:
Stephane Perrey, Université de Montpellier, FranceAlessandro L. Colosio, Université Jean Monnet, France
Copyright © 2025 Van Beekvelt, Wilson and Parviz Gaarden. 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: Mireille van Beekvelt, TWlyZWlsbGUudmFuLmJlZWt2ZWx0QG50bnUubm8=
†Present addresses: Marte Wilson, University Hospital of North Norway, Tromsø, Norway; Andreas Parviz Gaarden, Trondheim Municipality, Trondheim, Norway