Yuma Tamura1
Takafumi Kawashima2,3
Aoi Kodama4
Rui-Cheng Ji1,3,4
Yuta Itoh5
Nobuhide Agata6
Keisuke Kawakami1,3,4*- 1Physical Therapy Research Field, Graduate School of Medicine, Oita University, Yufu, Japan
- 2Department of Rehabilitation, Akeno-Central Hospital, Oita, Japan
- 3Faculty of Welfare and Health Science, Oita University, Oita, Japan
- 4Graduate School of Welfare and Health Science, Oita University, Oita, Japan
- 5Faculty of Rehabilitation Science, Nagoya Gakuin University, Nagoya, Japan
- 6Faculty of Health and Medical Sciences, Tokoha University, Hamamatsu, Japan
Introduction: Lymphatic vessels contribute to tissue homeostasis. Although the lymphatic vessels in skeletal muscle are known to undergo structural changes under certain conditions, such as atrophy and injury, effects of exercise on intramuscular lymphatic vessels remain unclear.
Methods: This study was aimed at investigating whether 8 weeks of voluntary wheel running (VWR) induces histological changes in lymphatic and blood capillaries, and whether these responses are related to age and myofiber type. Young (3-month-old) and aged (18-month-old) male C57BL/6 mice were assigned to sedentary or VWR groups. The soleus (SOL; slow-twitch) and plantaris (PLAN; fast-twitch) muscles were analyzed using immunohistochemistry and quantitative polymerase chain reaction.
Results: In young mice, VWR increased the quantity of type I myofibers and significantly enhanced the density of lymphatic vessels and blood capillaries in the SOL, besides upregulating the expression of vascular endothelial growth factors, VEGF-C and VEGF-D. These changes were not observed in aged mice or in the PLAN of mice in either age group.
Discussion: Although aged mice showed a similar increase in the quantity of type I myofibers, they did not exhibit corresponding vascular remodeling, which suggests that aging reduces responsiveness to exercise-induced angiogenic and lymphangiogenic signals. Overall, these findings indicate that VWR promotes lymphangiogenesis and angiogenesis in slow-twitch muscle in young mice, probably as an adaptive response to meet the increased oxygen demand. Exercise-induced vascular and lymphatic remodeling in skeletal muscle is significantly influenced by age and myofiber type, highlighting a reduced adaptive capacity of aged muscle that may impact strategies for promoting vascular health through physical activity.
1 Introduction
Lymphatic vessels are an integral component of the circulatory system, playing essential roles in lipid and immune cell transport as well as in the regulation of inflammatory response (Angeli and Lim, 2023; Hu et al., 2024). In peripheral tissues, lymphatic vessels are also responsible for collecting interstitial fluid that leaks from blood capillaries, and thereby contribute to fluid homeostasis (Oliver et al., 2020). Lymphatic vessels are also present within skeletal muscles. Kivelä et al. (2007a) were the first to report that small lymphatic vessels are located adjacent to blood capillaries between myofibers in human and murine skeletal muscles.
The distribution and morphology of intramuscular lymphatic vessels are altered under pathological conditions, such as muscle atrophy and injury (Kawashima et al., 2021; Tamura et al., 2024). Furthermore, recent report suggested that aging and myofiber type also influence the structural characteristics and distribution of lymphatic vessels (Taketa et al., 2024). These findings indicate that the architecture of lymphatic vessels within skeletal muscles is affected not only by pathological factors but also by physiological factors.
Changes in lymphatic vessels are primarily mediated by the expression of lymphangiogenic factors, such as vascular endothelial growth factors, VEGF-C and VEGF-D, and their receptors (VEGFR-3) (Tammela and Alitalo, 2010; Ji, 2012; Suarez et al., 2023). These molecules regulate the distribution and morphology of intramuscular lymphatic vessels in response to pathological and physiological stimuli (Ji, 2018; Kawashima et al., 2021; Pu et al., 2021; Tamura et al., 2024) and are critical for maintaining homeostasis in skeletal muscles.
The effects of exercise on intramuscular lymphatic vessels remain largely unclear. Although studies using treadmill-based exercise were previously conducted, no significant changes or decrease in intramuscular lymphatic vessels were observed after such interventions (Kivelä et al., 2007b; Greiwe et al., 2016). Voluntary wheel running (VWR) has recently gained attention as an exercise intervention that closely resembles natural locomotor patterns and does not disrupt the normal nocturnal circadian rhythm of animals (De Bono et al., 2006; Goh and Ladiges, 2015; Manzanares et al., 2018). VWR improves muscular endurance and enhances metabolic function in young mice (Pellegrino et al., 2005; LaBarge et al., 2016; Hedges et al., 2019). Furthermore, VWR prevents age-related decline in muscle function in aged mice (Yoshioka et al., 2019; Baek et al., 2022), underscoring its potential impact on skeletal muscles. VWR also increases blood capillary density and upregulates angiogenic factors such as VEGF-A in both slow- and fast-twitch skeletal muscles (Waters et al., 2004; Jackson et al., 2015; Hyatt et al., 2019; Zmudzka et al., 2025). However, the effects of VWR on intramuscular lymphatic vessels have not been examined. Investigating changes in intramuscular lymphatic vessels induced low-intensity and low-stress exercise like VWR may provide clues useful for understanding lymphedema rehabilitation and age-related muscle function decline.
Therefore, in this study, we investigated the histological changes in intramuscular lymphatic vessels and blood capillaries in young and aged mice subjected to VWR and evaluated the effects on the levels of lymphangiogenic and angiogenic factors in these animals.
2 Materials and methods
2.1 Animals
All experimental procedures were approved by the Animal Ethics Committee of Oita University (Approval No. 235701). Male C57BL/6 mice were purchased from The Jackson Laboratory Japan (Kanagawa, Japan). Three-month-old mice were assigned to the young group (n = 12), whereas 18-month-old mice were assigned to the aged group (n = 12). Each age group was further divided into two subgroups based on whether the mice engaged in spontaneous physical activity using a running wheel (VWR) or remained in a standard housing without access to the running wheel (SED). Accordingly, the following four groups were established: Young_SED, Young_VWR, Aged_SED, and Aged_VWR (n = 6 per group). All animals were housed under controlled environmental conditions (25 °C, 12 h light/dark cycle) with ad libitum access to standard laboratory chow and water. Only male mice were used in this study to avoid variability caused by the estrous cycle, which affects the voluntary running behavior. Female mice tend to run longer distances and exhibit greater day-to-day variability than males (Manzanares et al., 2018). The use of male mice helped minimize confounding factors and ensured consistency in activity levels.
2.2 Voluntary wheel running protocol
The VWR protocol was implemented by housing mice for 8 weeks in individual cages (RWC-15; Melquest, Toyama, Japan) equipped with a running wheel (diameter, 140 mm), allowing unrestricted access to voluntary running. Running distances were recorded every 24 h from the first to the last day of the protocol. The SED group mice were housed in standard cages without a running wheel for the same duration. The activity levels of mice in the SED groups were not recorded.
2.3 Muscle sampling
After the 8-week housing period, body weight was measured, and the mice were euthanized by cervical dislocation within 1 h of the final VRW. The soleus (SOL; predominantly slow-twitch) and plantaris (PLAN; predominantly fast-twitch) muscles were excised from the posterior region of both the hindlimbs, and their wet weight was measured. Muscles harvested from the left side were rapidly frozen in isopentane for subsequent histological and immunohistochemical analyses and stored at −80 °C until further analysis. Muscles from the right side were immediately frozen in liquid nitrogen for real-time polymerase chain reaction (RT-PCR) analysis and stored at −80 °C until further analysis.
2.4 Histological and immunohistochemical analyses
Muscle samples from the left hindlimb were used for histological and immunohistochemical analyses. Frozen transverse sections (8 μm thick) of the SOL and PLAN were prepared at −25 °C using a cryostat (CM 1860; Leica Biosystems, Nussloch, Germany). The sections were fixed in 4% paraformaldehyde (PFA; Fujifilm Wako Pure Chemical Co.) prepared in phosphate-buffered saline (PBS) for 11 min, rinsed with PBS, and stained with hematoxylin and eosin (H&E). The stained sections were mounted in 90% glycerol. Images were captured using a light microscope (IX70; Olympus, Tokyo, Japan) with a ×20 objective lens (UPLFLN; Olympus). The ImageJ software (ver. 1.53; NIH, Bethesda, MD, United States) was used to measure the whole muscle cross-sectional area (CSA). From the muscle belly, 150 myofibers were randomly selected for measuring the myofiber CSA.
For immunohistochemical analysis, the sections were fixed in 4% PFA in PBS for 11 min and blocked with 4% Block Ace powder (KAC, Kyoto, Japan) for 60 min at room temperature (RT). The sections were incubated overnight at 4 °C with the following primary antibodies: rabbit anti-lymphatic endothelial hyaluronan receptor 1 (LYVE-1; 1:1000; 103-PA50AG; Relia Tech GmbH, Wolfenbüttel, Germany; RRID: AB_2876870) and rat anti-CD31 (1:50; 550274; BD Biosciences, San Jose, US-CA, CA; RRID: AB_393571). The sections were washed with PBS and further incubated in the dark for 60 min at RT with Alexa Fluor 488-conjugated goat anti-rabbit IgG (1:400; A-11008; Thermo Fisher Scientific, Waltham, US-MA; RRID: AB_143165) and Alexa Fluor 568-conjugated goat anti-rat IgG (1:400; A-11077; Thermo Fisher Scientific; RRID: AB_2534121). The stained sections were mounted in 90% glycerol. Fluorescent images were acquired using a fluorescence microscope (IX70; Olympus) with a ×20 objective lens (UPLFLN; Olympus). The contrast of images was adjusted using Adobe Photoshop (Adobe, San Jose, CA, United States). The ImageJ software (ver. 1.53; NIH, Bethesda, MD, United States) was used to count LYVE-1-positive structures as lymphatic vessels and CD31-positive structures as blood capillaries (Kivelä et al., 2007a). Data were normalized by area (density per unit area) and number of myofibers (ratio to myofibers). Five random areas, measuring 300 μm × 300 μm, were selected from the immunohistochemical staining images and used for analysis. Please refer to Supplementary Table S1 for the mean ± standard deviation values of lymphatic vessels, blood capillaries, and myofibers analyzed in each group.
To identify myofiber types, the sections were incubated overnight at 4 °C with primary antibodies, including rabbit anti-slow skeletal myosin heavy chain (Type I myofiber; 1:200; ab234431; Abcam, Waltham, MA; RRID: AB_3076242) and rat anti-laminin alpha-2 (laminin; 1:100; sc-59854; Santa Cruz Biotechnology; Dallas, TX; RRID: AB_784266). The sections were washed with PBS and further incubated in the dark for 60 min at RT with Alexa Fluor 488-conjugated goat anti-rat IgG (1:400; A-11006; Thermo Fisher Scientific; RRID: AB_2534074) and Alexa Fluor 568-conjugated goat anti-rabbit IgG (1:400; A-11011; Thermo Fisher Scientific; RRID: AB_143157). The stained sections were mounted in 90% glycerol and imaged with a fluorescence microscope using a ×10 objective lens (UPLFLN; Olympus). Myofibers were defined as cells surrounded by laminin-positive boundaries, whereas those positive for slow myosin heavy chain were identified as Type I myofibers. The ImageJ software was used to quantify the total number of myofibers, and the quantity of Type I myofibers relative to the total myofiber count was calculated for each section of the SOL and PLAN.
2.5 Real-time reverse transcription PCR
Total RNA was isolated from the right SOL and PLAN using Isogen II (Nippon Gene, Tokyo, Japan) in accordance with the manufacturer’s instructions. RNA concentration was measured using a DS-11 spectrophotometer (DeNovix, Wilmington, DE, United States). Complementary DNA (cDNA) was synthesized from total RNA using SuperScript VILO Master Mix (Thermo Fisher Scientific). Real-time PCR was performed using Fast SYBR Green Master Mix (Thermo Fisher Scientific) on a StepOnePlus Real-Time PCR System (Applied Biosystems). The relative mRNA expression levels of VEGF-A, VEGF-C, VEGF-D, VEGFR-2, VEGFR-3, tumor necrosis factor (TNF)-α, and interleukin (IL)-1β were quantified using the comparative Ct method (ΔΔCt) and normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as an internal control. All primers were purchased from Greiner Bio-One (Kremsmünster, Austria). Primer sequences are listed in Supplementary Table S2.
2.6 Statistical analysis
All data are presented as mean ± standard deviation (SD). Two-sample t-tests were used to compare running distances for the Young_VWR and Aged_VWR groups. Two-way analysis of variance (ANOVA) was conducted to evaluate the effects of Age (Young vs. Aged), Run (SED vs. VWR), and their interaction. When a significant interaction was detected, one-way ANOVA followed by Tukey’s post hoc test was applied. For significant main effects without interaction, pairwise comparisons were conducted using the t-test. A p-value <0.05 was considered statistically significant. All statistical analyses were performed using Modified R Commander 4.0.2 (https://home.hirosaki-u.ac.jp/pteiki/r/).
3 Results
3.1 Running performance and physical parameters
Spontaneous running activity was compared between the Young_VWR and Aged_VWR groups. The total running distance over 8 weeks in the Young_VWR group was approximately twice as high (456.6 ± 93.1 km) as in the Aged_VWR group (221.0 ± 58.2 km) (p < 0.01).
At the end of the housing period, body mass was not significantly affected by VWR (p = 0.64); however, aged mice exhibited significantly higher body mass than young mice, regardless of the activity level. Both the Aged_SED (32.9 ± 1.2 g) and Aged_VWR (32.3 ± 1.3 g) groups exhibited significantly higher body mass than their respective Young groups (Young_SED: 26.7 ± 1.8 g; Young_VWR: 25.6 ± 0.8 g) (Table 1). Data on body mass changes over the 8-week intervention period in the VWR group are provided in Supplementary Table S3.
Table 1. Comparison of running performance and physical parameters in different groups.
The SOL mass was not significantly affected by VWR (p = 0.60). However, a significant main effect of age was detected (p < 0.05), with Aged_SED (11.8 ± 0.6 mg) and Aged_VWR (12.2 ± 1.0 mg) group mice showing greater muscle mass than the respective Young group mice (Young_SED: 10.6 ± 0.9 mg; Young_VWR: 11.4 ± 0.8 mg). For the PLAN mass, no significant interaction (p = 0.95) or main effect of Age (p = 0.19) and Run (p = 0.54) was detected. The values remained relatively comparable across the groups. Additionally, the mass of each muscle was normalized to the body weight (Supplementary Table S4).
3.2 Histological changes in muscle morphology
To assess morphological changes induced by VWR, we analyzed both whole muscle CSA and myofiber CSA. In the SOL (Figure 1A), whole muscle CSA showed no significant interaction between Age and Run (p = 0.09) and no main effects of Age (p = 0.77) and Run (p = 0.61) (Figure 1B). By contrast, myofiber CSA showed a significant interaction between Age and Run (p < 0.01); the Aged_VWR group (1572.0 ± 160.5 µm2) was significantly greater than the Young_VWR group (1242.8 ± 51.6 µm2), with no difference between Young_SED (1475.1 ± 161.9 µm2) and Aged_SED (1337.0 ± 229.1 µm2) (Figure 1C).
Figure 1. (A) Representative hematoxylin and eosin (H&E)-stained images of the soleus muscle. Enlarged views of the regions outlined by white boxes are also shown. (B) Quantification of the whole cross-sectional area of the soleus muscle. (C) Quantification of the myofiber cross-sectional area of the soleus muscle. (D) Representative H&E-stained images of the plantaris muscle. Enlarged views of the regions outlined by white boxes are also shown. (E) Quantification of the whole cross-sectional area of the plantaris muscle. (F) Quantification of the myofiber cross-sectional area of the plantaris muscle. Scale bars = 1 mm (full image) and 100 µm (enlarged image). Data are presented as means ± SD.
In the PLAN (Figure 1D), whole muscle CSA showed no significant Age × Run interaction (p = 0.17) but a significant main effect of Run (p < 0.05). Post hoc comparisons indicated that Old_VWR (1.89 ± 0.27 mm2) had greater whole muscle CSA than Old_SED (1.47 ± 0.19 mm2), whereas the Young groups did not differ (Figure 1E). For myofiber CSA, there was no significant Age × Run interaction (p = 0.26) and no main effects of Age (p = 0.39) or Run (p = 0.06) (Figure 1F).
3.3 Changes in myofiber type composition
The images of SOL and PLAN sections immunostained for Type I myofibers are presented in Figure 2A. The quantity of Type I myofibers in the SOL showed no significant interaction between Age and Run (p = 0.89); however, there was a significant main effect of Run (p < 0.01). Notably, both Young_VWR (52.8% ± 4.3%) and Aged_VWR (52.1% ± 4.0%) group mice showed significantly higher quantities of Type I myofibers in the SOL than Young_SED (41.0% ± 2.8%) and Aged_SED (40.9% ± 6.4%) group mice (Figure 2B). In the PLAN, the quantity of type I myofibers (two-factor model: Age, Run) showed no significant Age × Run interaction (p = 0.56). Using a two-sided α = 0.05 (significant if p < 0.05), the main effect of Age was borderline (p = 0.05) and was not considered significant; post hoc comparisons likewise showed no differences between young and aged mice within either activity condition (Young_SED vs. Aged_SED: p = 0.11; Young_VWR vs. Aged_VWR: p = 0.28) (Figure 2C). Data on the total number of myofibers and the number of Type I myofibers in each muscle are provided in Supplementary Table S5.
Figure 2. (A) Images of the soleus and plantaris muscle sections immunostained for laminin (green) and Type I myofibers (red). The quantity of Type I myofibers in the soleus (B) and plantaris muscle (C). Scale bars = 100 μm. Data are presented as means ± SD.
3.4 Changes in the lymphatic vessels and blood capillaries
To quantify lymphatic vessel density and the ratio of lymphatic vessels to myofibers in SOL, we immunostained the muscle sections for LYVE-1 and CD31 (Figure 3A). In the SOL, no significant interaction between Age and Run was observed for lymphatic vessel density (LVD) (p = 0.11), but a significant main effect of Run was detected in Young_SED vs. Young_VWR or Aged_SED vs. Aged_VWR (p < 0.01). The LVD in Young_VWR group mice (184.4 ± 26.1/mm2) was significantly higher than that in Young_SED group mice (133.0 ± 9.4/mm2), whereas no significant difference was found between the Aged groups (Aged_SED: 135.6 ± 17.0/mm2, Aged_VWR: 152.2 ± 33.1/mm2) (Figure 3B). The ratio of lymphatic vessels to myofibers showed no significant interaction between Age and Run (p = 0.77); however, a significant main effect of Run was detected in Young_SED vs. Young_VWR or Aged_SED vs. Aged_VWR (p < 0.05). The ratio of lymphatic vessels to myofibers in the Young_VWR group (0.29 ± 0.04) was significantly higher than that in the Young_SED group (0.24 ± 0.03), whereas no significant difference was observed between the aged groups (Aged_SED: 0.24 ± 0.04, Aged_VWR: 0.28 ± 0.06) (Figure 3C).
Figure 3. (A) Immunostaining images of the soleus muscles (gray: myofiber imaged via phase-contrast microscopy; green: LYVE-1-positive lymphatic vessels; red: CD31-positive blood capillaries). Quantification of lymphatic vessels density (B) and the ratio of lymphatic vessels to myofibers (C). Quantification of blood capillary density (D) and the ratio of blood capillaries to myofibers (E). Scale bars = 50 μm. Data are presented as means ± SD.
For blood capillary density in the SOL, a significant Age × Run interaction was detected (p < 0.05): mice in the Young_VWR group showed higher capillary density (1670.4 ± 121.6/mm2) than Young_SED (1351.5 ± 157.2/mm2), whereas no difference was observed between Aged_SED (1493.0 ± 169.3/mm2) and Aged_VWR (1448.5 ± 198.3/mm2) (Figure 3D). The ratio of blood capillaries to myofibers showed no significant Age × Run interaction (p = 0.27), and there were no main effects of Age (p = 0.12) or Run (p = 0.13) (Figure 3E).
To determine the density and myofiber-normalized ratios of lymphatic vessels and blood capillaries in the PLAN, we immunostained muscle sections for LYVE-1 and CD31 (Figure 4A). In the PLAN, LVD showed no significant Age × Run interaction (p = 0.29) and no main effects of Age (p = 0.22) or Run (p = 0.74) (Figure 4B). The ratio of lymphatic vessels to myofibers likewise showed no significant Age × Run interaction (p = 0.83) and no main effects of Age (p = 0.17) or Run (p = 0.11) (Figure 4C).
Figure 4. (A) Immunostaining images of the plantaris muscles (gray: myofiber imaged via phase-contrast microscopy; green: LYVE-1-positive lymphatic capillaries; red: CD31-positive blood capillaries). Quantification of lymphatic vessels density (B) and the ratio of lymphatic vessels to myofibers (C). Quantification of blood capillary density (D) and the ratio of blood capillaries to myofibers (E). Scale bars = 50 μm. Data are presented as means ± SD.
For blood capillary density in the PLAN, there was no significant Age × Run interaction (p = 0.15) and no main effects of Age (p = 0.76) or Run (p = 0.99) (Figure 4D). The ratio of blood capillaries to myofibers showed no significant Age × Run interaction (p = 0.41). There was a main effect of Age (p = 0.05); however, post hoc comparisons revealed no differences between SED and VWR within age (Young_SED vs. Young_VWR: p = 0.08; Aged_SED vs. Aged_VWR: p = 0.08) (Figure 4E).
3.5 Changes in mRNA expression of VEGF family members
In the SOL, VEGF-C expression showed a significant Age × Run interaction (p < 0.01): levels were higher in the Young_VWR and Aged_SED groups than in Young_SED, whereas the Aged_VWR group exhibited lower levels than both Young_VWR and Aged_SED (Figure 5A). In the PLAN, VEGF-C also showed a significant Age × Run interaction (p < 0.05), but post hoc comparisons did not reach significance (p = 0.12) (Figure 5B). The expression of VEGF-D in the SOL paralleled VEGF-C, except that no significant difference was found between Young_VWR and Aged_VWR (Figure 5C). In the PLAN, VEGF-D again showed a significant Age × Run interaction (p < 0.05), but group differences did not reach significance (p = 0.06) (Figure 5D). VEGFR-3 expression showed a significant Age × Run interaction in the SOL (p < 0.05), with Aged_VWR lower than the other groups (Figure 5E). In the PLAN, no Age × Run interaction was observed (p = 0.41), but a main effect of Run was detected (p < 0.05); post hoc comparisons indicated higher expression in Aged_VWR than Aged_SED (Figure 5F).
Figure 5. mRNA expression levels of VEGF family members in the soleus (A,C,E,G,I) and plantaris (B,D,F,H,J) muscles. mRNA levels of the lymphangiogenic factors VEGF-C and VEGF-D, and their receptor VEGFR-3 (A–F). mRNA levels of the angiogenic factor VEGF-A and its receptor VEGFR-2 (G–J). Data are presented as means ± SD.
VEGF-A expression in the SOL exhibited a significant Age × Run interaction (p < 0.01), with Aged_VWR lower than Young_SED and Aged_SED (Figure 5G). In the PLAN, an Age × Run interaction was also detected (p < 0.05), although pairwise differences were marginal and did not reach significance (p = 0.05) (Figure 5H). VEGFR-2 expression showed no Age × Run interaction and no main effects of Age or Run in either the SOL or PLAN (Figures 5I,J).
3.6 Changes in mRNA expression of inflammatory cytokines
In both SOL and PLAN muscles, TNF-α expression showed no significant Age × Run interaction and no main effects of Age or Run (Figures 6A,B). Similarly, IL-1β expression showed no Age × Run interaction and no main effects of Age or Run in either muscle (Figures 6C,D).
Figure 6. mRNA expression levels of TNF-α (A,B) and IL-1β (C,D) in the soleus and plantaris muscles. Data are presented as means ± SD.
4 Discussion
4.1 Morphological adaptations to VWR and aging
In this study, we investigated the effects of 8 weeks of VWR on histological changes in intramuscular lymphatic vessels and blood capillaries, using both young and aged mice. We focused on the SOL, a muscle rich in slow-twitch myofibers, and the PLAN, which is predominantly composed of fast-twitch myofibers (Bloemberg and Quadrilatero, 2012; Mobley et al., 2017).
No significant differences in whole muscle CSA or myofiber CSA were observed between the Young_SED and Aged_SED groups for both SOL and PLAN. However, in Aged_VWR, myofiber CSA was significantly larger for SOL than for the Young_VWR group. Previous study have shown that VWR reduces the myofiber CSA in young mice (Pellegrino et al., 2005). In contrast, VWR has been reported to increase myofiber CSA in aged mice (Baek et al., 2022). The results of this study indicate morphological adaptation due to changes in VWR responsiveness with aging. However, because muscle function, such as muscle strength or fatigue resistance, was not directly evaluated in this study, this interpretation should be considered preliminary and warrants further investigation. In young mice, VWR induced an increase in Type I myofibers in the SOL. VWR upregulates Type I myosin isoforms while downregulating Type IIa isoform in the SOL (Fuller et al., 2006; Hyatt et al., 2019); our findings are consistent with these observations. In contrast, no significant changes in myofiber-type composition were observed in the PLAN of young mice. VWR does not induce a shift from Type II to Type I myofibers in the PLAN (Waters et al., 2004), which is indicative of a limited response of the PLAN to endurance-based stimuli such as VWR.
4.2 Vascular remodeling in SOL induced by VWR
In the SOL of young mice, VWR significantly increased the LVD, the ratio of lymphatic vessels to myofibers. Lymphangiogenesis is primarily driven by the binding of VEGF-C and VEGF-D to their receptor VEGFR-3 (Tammela and Alitalo, 2010; Ji, 2012; Suarez et al., 2023). In our study, the expression of VEGF-C and VEGF-D was significantly increased by VWR, corresponding with the observed increase in LVD and the ratio of lymphatic vessels to myofibers. These results indicate that VEGF-C and VEGF-D may contribute to the VWR-induced increase in lymphangiogenesis. However, the VEGFR-3 mRNA levels did not change despite VWR. Kawashima et al. reported that the expression of VEGF-C and VEGF-D was decreased in muscle atrophy, whereas that of VEGFR-3 remained stable (Kawashima et al., 2021). Accordingly, the expression of these factors in skeletal muscle may not always be tightly correlated.
In contrast to our findings with VWR, previous treadmill-based exercise studies have not demonstrated lymphangiogenic responses in skeletal muscle (Kivelä et al., 2007b; Greiwe et al., 2016). Compared with VWR, the total daily running distance achieved with treadmill-based exercise is generally shorter. Moreover, treadmill-based exercise involves higher exercise intensity and exercise within a limited time frame. It is possible that differences in these exercise modalities may influence the development of intramuscular lymphatic vessels.
The increase in intramuscular blood capillary density in the SOL in response to VWR has also been documented (Hyatt et al., 2019), and our study corroborates these findings. While angiogenesis is mediated by the binding of VEGF-A to its receptor VEGFR-2 (Kudelka et al., 2013; Wang et al., 2020; Shah et al., 2025), no significant changes in VEGF-A or VEGFR-2 mRNA levels were observed in the SOL of young mice following VWR in this study. Prior research, including embryological studies and in vitro experiments, suggests that blood vessel formation precedes lymphangiogenesis (Martinez-Corral et al., 2015; Osaki et al., 2018; Hu et al., 2024; Schoofs and Mäkinen, 2024), and angiogenesis in skeletal muscles reportedly occurs earlier than lymphangiogenesis (Kawashima et al., 2021; Pu et al., 2021). Therefore, in our study, it is possible that at the time of muscle harvest, angiogenesis had already been completed and VEGF-A expression was no longer required. However, the reason for the unchanged VEGF-A expression remains unclear and warrants further investigation, including consideration of the timing of muscle sampling. Notably, VEGF-D has also been implicated in angiogenesis in skeletal muscles (Rissanen et al., 2003; Kärkkäinen et al., 2009). In our study, increased VEGF-D expression in the SOL of the Young_VWR group coincided with greater blood capillary density, suggesting a possible involvement of VEGF-D in capillary formation.
In contrast, no significant differences were observed in the blood capillary to myofiber ratio in the SOL among groups. This may be attributable to the evaluation being based on the numbers of myofibers, lymphatic vessels, and blood capillaries per unit area. Specifically, in the Young_VWR group, the slightly smaller myofiber CSA (Figure 1C) likely resulted in a greater number of myofibers within the analyzed regions, thereby influencing the ratio (Supplementary Table S1). Moreover, the quantity of type I myofibers increased in the VWR groups, suggesting that the observed changes in myofiber CSA may reflect physiological adaptations induced by exercise. However, further studies are needed to clarify these mechanisms in detail.
VWR increased the proportion of type I myofibers in the SOL, which are generally considered fatigue-resistant (Carrell et al., 2016). Type I myofibers possess high oxidative capacity (Xiao et al., 2021) and require a sustained supply of oxygen to maintain prolonged contractile activity. This continuous oxygen demand may have contributed to the observed increase in capillary density. The resulting increase in interstitial fluid leakage from blood capillaries could in turn enhance lymphatic vessel formation, which plays a role in fluid clearance (Angeli and Lim, 2023; Hu et al., 2024). In contrast, no changes in lymphangiogenesis or angiogenesis were observed in the PLAN. This muscle did not undergo a shift toward a slow-twitch phenotype following VWR, which provides a plausible explanation for the absence of vascular adaptations in this muscle.
It should be noted, however, that in mouse skeletal muscle, type IIa myofibers are known to exhibit the highest oxidative capacity (Gouspillou et al., 2014; Thomas et al., 2014), and this factor should be taken into account when interpreting the present findings. Therefore, a more comprehensive evaluation of myofiber type composition, including type IIa myofibers, together with functional assessments of oxidative capacity such as succinate dehydrogenase (SDH) activity, remains an important subject for future studies.
4.3 Age-related differences in vascular responsiveness
Although aging has been associated with a decline in the number of lymphatic vessels in dermal, heart, and rat skeletal muscles (Karaman et al., 2015; Taketa et al., 2024; Roh et al., 2025), we observed no significant differences in LVD between young and aged mice. This might be due to the lack of age-related changes in myofiber-type composition in our samples. Previous reports have indicated a decrease in the expression of lymphangiogenic factor with aging (Taketa et al., 2024); however, in our study, the expression levels of these factors were significantly increased in the SOL of aged mice. This may reflect a compensatory upregulation in response to an age-related decline in the lymphatic vasculature.
VEGF-C and VEGF-D are induced by proinflammatory cytokines, such as TNF-α and IL-1 (Zhang et al., 2007; Ji, 2012; Schwager and Detmar, 2019). Kawanishi and Machida (2021) reported elevated levels of these cytokines in the skeletal muscle of aged mice (27-month-old), which are near the end of their lifespan. In contrast, in the present study, mRNA levels of TNF-α and IL-1 did not differ significantly between young and aged mice, suggesting that the upregulation of VEGF-C and VEGF-D was not driven by proinflammatory cytokines. We selected 18-month-old mice, which are considered to retain relatively high activity levels (Shoji and Miyakawa, 2019). This selection was intended to evaluate whether VWR could help maintain physical activity and induce changes in intramuscular lymphatic vessels in aged mice. Notably, small differences in age among aged mice can lead to significant alterations in gene expression profiles (Kang et al., 2022). The discrepancy between our findings and those of Kawanishi and Machida (2021) may be attributable to differences in the age of the mice used.
Therapeutic strategies targeting VEGF-C or VEGF-D signaling may help overcome the impaired vascular responsiveness in aged skeletal muscle. Anisimov et al. (2009) reported that sustained overexpression of activated VEGF-C and VEGF-D in skeletal muscle using recombinant adeno-associated virus vectors significantly increased the densities of both blood and lymphatic vessels, improved tissue perfusion, and enhanced lymphatic drainage. These results implied that exogenous delivery of VEGF-C or VEGF-D could potentially restore lymphangiogenesis and might mitigate aging-associated muscle decline. Despite an increase in Type I myofibers following VWR, aged mice did not show corresponding increases in LVD and blood capillary density in the SOL. This discrepancy may be attributable to differences in the expression patterns of lymphangiogenic factors between young and aged mice. Blood capillary responsiveness to exercise is known to decline with age (Hendrickse et al., 2021), and a similar phenomenon may explain the lack of capillary adaptation observed in our aged cohort. Furthermore, age-related functional decline in lymphatic vessels has also been reported (Ji, 2018; Yang et al., 2023), which is suggestive of reduced responsiveness to exercise-induced stimuli. Because other molecules, such as angiopoietins and platelet-derived growth factors, are also critically involved in lymphangiogenesis (Cao et al., 2004; Korhonen et al., 2022; Hu et al., 2024), future studies should investigate their contributions and interactions in this context. In particular, the significantly shorter total running distance observed for aged mice compared with that for young mice over the 8-week VWR period may reflect an underlying decline in vascular and lymphatic function. Clarifying the relationship between vessel density and endurance performance will be essential to validate this hypothesis in future studies.
4.4 Methodological limitations and future studies
LYVE-1 is commonly used as a marker for lymphatic vessel structures; however, it can also be expressed by macrophages, particularly under inflammatory conditions (Hall et al., 2012; Zhang et al., 2021). On the other hand, voluntary wheel running (VWR) is reported to be a non-damaging form of exercise. For example, even in Duchenne muscular dystrophy model mice, VWR does not induce muscle injury (Monceau et al., 2022). In this study, no significant differences in the mRNA expression levels of the pro-inflammatory cytokines TNF-α and IL-1β were observed among groups, suggesting that inflammation was unlikely to have occurred. Nevertheless, protein-level analyses were not performed, and this remains an important limitation and direction for future studies. Moreover, the LYVE-1–positive cells identified in this study colocalized with VEGFR-3 (Supplementary Figure S1), another well-established marker of lymphatic vessels (Weber et al., 2018). Therefore, it is considered appropriate to regard LYVE-1–positive cells as lymphatic vessels in the present study.
In addition, our study was focused primarily on histological changes in lymphatic vessels and did not address their physiological functions, which represents a key limitation. Future studies using functional inhibition or stimulation approaches will be essential to elucidate the physiological roles of intramuscular lymphatic vessels and their contribution to exercise-induced muscle adaptation.
Furthermore, another significant limitation is that our analysis was confined to two-dimensional techniques such as immunohistochemical staining. With respect to blood capillaries in skeletal muscle, previous studies employing three-dimensional technologies such as confocal microscopy and two-photon laser scanning microscopy have demonstrated more comprehensive insights (Umek et al., 2019; Shimotsu et al., 2021). However, three-dimensional imaging of intramuscular lymphatic vessels has not yet been established and represents an important challenge for future studies. Once such techniques are developed, they may enable more precise analyses of exercise- or disease-induced alterations in the distribution and morphology of intramuscular lymphatic vessels. Moreover, they are expected to facilitate the evaluation of lymphatic fluid retention capacity using macromolecular tracer techniques, thereby contributing to more functional assessments.
As a side note, the lack of pre-intervention weight measurements in the SED group, coupled with the relatively small sample size (n = 6 per group), which should be considered when interpreting the results.
4.5 Conclusion
In conclusion, this study revealed that 8 weeks of VWR significantly increased LVD and blood capillary density in the SOL of young mice, possibly via VEGF-C and VEGF-D upregulation. These changes appear to represent an adaptive response to the increased oxygen demand associated with a VWR-induced shift toward a slow-twitch phenotype. In contrast, no changes in LVD and blood capillary density were observed in aged mice, possibly due to differences in their responsiveness to exercise stimulus.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by Animal Ethics Committee of Oita University (Approval No. 235701). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
YT: Data curation, Conceptualization, Writing – review and editing, Writing – original draft, Investigation. TK: Conceptualization, Investigation, Writing – review and editing, Data curation. AK: Investigation, Data curation, Conceptualization, Writing – review and editing. R-CJ: Writing – review and editing, Investigation, Funding acquisition, Conceptualization. YI: Writing – review and editing, Conceptualization, Funding acquisition. NA: Writing – review and editing, Conceptualization, Funding acquisition. KK: Conceptualization, Writing – review and editing, Funding acquisition, Investigation, Project administration.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Japan Society for the Promotion of Science under the grant numbers KAKENHI JP23K24698, 25K02980 (to KK), JP23K10429 (to R-CJ), JP24K14362 (to NA), and JP22K11351 and JP25K14520 (to YI).
Acknowledgments
We express our gratitude to Miyuki Hinago, Miwako Minami, and Kaho Tominaga from the Kawakami Lab for their helpful support. Additionally, we thank Editage (https://www.editage.jp/) for language editing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1654445/full#supplementary-material
References
Angeli V., Lim H. Y. (2023). Biomechanical control of lymphatic vessel physiology and functions. Cell. Mol. Immunol. 20, 1051–1062. doi:10.1038/s41423-023-01042-9
Anisimov A., Alitalo A., Korpisalo P., Soronen J., Kaijalainen S., Leppänen V. M., et al. (2009). Activated forms of VEGF-C and VEGF-D provide improved vascular function in skeletal muscle. Circ. Res. 104, 1302–1312. doi:10.1161/CIRCRESAHA.109.197830
Baek K. W., Kim S. J., Kim B. G., Jung Y. K., Hah Y. S., Moon H. Y., et al. (2022). Effects of lifelong spontaneous exercise on skeletal muscle and angiogenesis in super-aged mice. PLoS One 17, e0263457. doi:10.1371/journal.pone.0263457
Bloemberg D., Quadrilatero J. (2012). Rapid determination of myosin heavy chain expression in rat, mouse, and human skeletal muscle using multicolor immunofluorescence analysis. PLoS One 7, e35273. doi:10.1371/journal.pone.0035273
Cao R., Björndahl M. A., Religa P., Clasper S., Garvin S., Galter D., et al. (2004). PDGF-BB induces intratumoral lymphangiogenesis and promotes lymphatic metastasis. Cancer Cell 6, 333–345. doi:10.1016/j.ccr.2004.08.034
Carrell E. M., Coppola A. R., McBride H. J., Dirksen R. T. (2016). Orai1 enhances muscle endurance by promoting fatigue-resistant type I fiber content but not through acute store-o perated Ca2+ entry. FASEB J. 30, 4109–4119. doi:10.1096/fj.201600621R
De Bono J. P., Adlam D., Paterson D. J., Channon K. M. (2006). Novel quantitative phenotypes of exercise training in mouse models. Am. J. Physiol. Regul. Integr. Comp. Physiol. 290, R926–R934. doi:10.1152/ajpregu.00694.2005
Fuller P. M., Baldwin K. M., Fuller C. A. (2006). Parallel and divergent adaptations of rat soleus and plantaris to chronic exercise and hypergravity. Am. J. Physiol. Regul. Integr. Comp. Physiol. 290, R442–R448. doi:10.1152/ajpregu.00578.2005
Goh J., Ladiges W. (2015). Voluntary wheel running in mice. Curr. Protoc. Mouse Biol. 5, 283–290. doi:10.1002/9780470942390.mo140295
Gouspillou G., Sgarioto N., Norris B., Barbat-Artigas S., Aubertin-Leheudre M., Morais J. A., et al. (2014). The relationship between muscle fiber type-specific PGC-1α content and mitochondrial content varies between rodent models and humans. PLoS One 9, e103044. doi:10.1371/journal.pone.0103044
Greiwe L., Vinck M., Suhr F. (2016). The muscle contraction mode determines lymphangiogenesis differentially in rat skeletal and cardiac muscles by modifying local lymphatic extracellular matrix microenvironments. Acta Physiol. (Oxf.) 217, 61–79. doi:10.1111/apha.12633
Hall K. L., Volk-Draper L. D., Flister M. J., Ran S. (2012). New model of macrophage acquisition of the lymphatic endothelial phenotype. PLoS One 7, e31794. doi:10.1371/journal.pone.0031794
Hedges C. P., Bishop D. J., Hickey A. J. R. (2019). Voluntary wheel running prevents the acidosis-induced decrease in skeletal muscle mitochondrial reactive oxygen species emission. FASEB J. 33, 4996–5004. doi:10.1096/fj.201801870R
Hendrickse P. W., Krusnauskas R., Hodson-Tole E., Venckunas T., Degens H. (2021). Regular endurance exercise of overloaded muscle of young and old Male mice does not attenuate hypertrophy and improves fatigue resistance. GeroScience 43, 741–757. doi:10.1007/s11357-020-00224-x
Hu Z., Zhao X., Wu Z., Qu B., Yuan M., Xing Y., et al. (2024). Lymphatic vessel: origin, heterogeneity, biological functions, and therapeutic targets. Signal Transduct. Target. Ther. 9, 9. doi:10.1038/s41392-023-01723-x
Hyatt J. K., Brown E. A., Deacon H. M., McCall G. E. (2019). Muscle-specific sensitivity to voluntary physical activity and detraining. Front. Physiol. 10, 1328. doi:10.3389/fphys.2019.01328
Jackson J. R., Kirby T. J., Fry C. S., Cooper R. L., McCarthy J. J., Peterson C. A., et al. (2015). Reduced voluntary running performance is associated with impaired coordination as a result of muscle satellite cell depletion in adult mice. Skelet. Muscle. 5, 41. doi:10.1186/s13395-015-0065-3
Ji R. C. (2012). Macrophages are important mediators of either tumor- or inflammation-induced lymphangiogenesis. Cell. Mol. Life Sci. 69, 897–914. doi:10.1007/s00018-011-0848-6
Ji R. C. (2018). Recent advances and new insights into muscular lymphangiogenesis in health and disease. Life Sci. 211, 261–269. doi:10.1016/j.lfs.2018.09.043
Kang Y. K., Min B., Eom J., Park J. S. (2022). Different phases of aging in mouse old skeletal muscle. Aging (Albany, NY) 14, 143–160. doi:10.18632/aging.203812
Karaman S., Buschle D., Luciani P., Leroux J. C., Detmar M., Proulx S. T. (2015). Decline of lymphatic vessel density and function in murine skin during aging. Angiogenesis 18, 489–498. doi:10.1007/s10456-015-9479-0
Kärkkäinen A. M., Kotimaa A., Huusko J., Kholova I., Heinonen S. E., Stefanska A., et al. (2009). Vascular endothelial growth factor-D transgenic mice show enhanced blood capillary density, improved postischemic muscle regeneration, and increased susceptibility to tumor formation. Blood 113, 4468–4475. doi:10.1182/blood-2008-07-171108
Kawanishi N., Machida S. (2021). Alterations of macrophage and neutrophil content in skeletal muscle of aged versus young mice. Muscle Nerve 63, 600–607. doi:10.1002/mus.27158
Kawashima T., Ji R. C., Itoh Y., Agata N., Sasai N., Murakami T., et al. (2021). Morphological and biochemical changes of lymphatic vessels in the soleus muscle of mice after hindlimb unloading. Muscle Nerve 64, 620–628. doi:10.1002/mus.27402
Kivelä R., Havas E., Vihko V. (2007a). Localisation of lymphatic vessels and vascular endothelial growth factors-C and -D in human and mouse skeletal muscle with immunohistochemistry. Cell Biol. 127, 31–40. doi:10.1007/s00418-006-0226-x
Kivelä R., Silvennoinen M., Lehti M., Kainulainen H., Vihko V. (2007b). Effects of acute exercise, exercise training, and diabetes on the expression of lymphangiogenic growth factors and lymphatic vessels in skeletal muscle. Am. J. Physiol. Heart Circ. Physiol. 293, H2573–H2579. doi:10.1152/ajpheart.00126.2007
Korhonen E. A., Murtomäki A., Jha S. K., Anisimov A., Pink A., Zhang Y., et al. (2022). Lymphangiogenesis requires Ang2/Tie/PI3K signaling for VEGFR3 cell-surface expression. J. Clin. Invest. 132, e155478. doi:10.1172/JCI155478
Kudelka M. R., Grossniklaus H. E., Mandell K. J. (2013). Emergence of dual VEGF and PDGF antagonists in the treatment of exudative age-related macular degeneration. Expert Rev. Ophthalmol. 8, 475–484. doi:10.1586/17469899.2013.840095
LaBarge S. A., Migdal C. W., Buckner E. H., Okuno H., Gertsman I., Stocks B., et al. (2016). p300 is not required for metabolic adaptation to endurance exercise training. FASEB J. 30, 1623–1633. doi:10.1096/fj.15-281741
Manzanares G., Brito-da-Silva G., Gandra P. G. (2018). Voluntary wheel running: patterns and physiological effects in mice. Braz. J. Med. Biol. Res. 52, e7830. doi:10.1590/1414-431X20187830
Martinez-Corral I., Ulvmar M. H., Stanczuk L., Tatin F., Kizhatil K., John S. W. M., et al. (2015). Nonvenous origin of dermal lymphatic vasculature. Circ. Res. 116, 1649–1654. doi:10.1161/CIRCRESAHA.116.306170
Mobley C. B., Mumford P. W., Kephart W. C., Haun C. T., Holland A. M., Beck D. T., et al. (2017). Aging in rats differentially affects markers of transcriptional and translational capacity in soleus and plantaris muscle. Front. Physiol. 8, 518. doi:10.3389/fphys.2017.00518
Monceau A., Moutachi D., Lemaitre M., Garcia L., Trollet C., Furling D., et al. (2022). Dystrophin restoration after adeno-associated virus U7-mediated dmd exon skipping is modulated by muscular exercise in the severe D2-mdx Duchenne muscular dystrophy murine model. Am. J. Pathol. 192, 1604–1618. doi:10.1016/j.ajpath.2022.07.016
Oliver G., Kipnis J., Randolph G. J., Harvey N. L. (2020). The lymphatic vasculature in the 21st century: novel functional roles in homeostasis and disease. Cell 182, 270–296. doi:10.1016/j.cell.2020.06.039
Osaki T., Serrano J. C., Kamm R. D. (2018). Cooperative effects of vascular angiogenesis and lymphangiogenesis. Regen. Eng. Transl. Med. 4, 120–132. doi:10.1007/s40883-018-0054-2
Pellegrino M. A., Brocca L., Dioguardi F. S., Bottinelli R., D’Antona G. (2005). Effects of voluntary wheel running and amino acid supplementation on skeletal muscle of mice. Eur. J. Appl. Physiol. 93, 655–664. doi:10.1007/s00421-004-1237-8
Pu Z., Shimizu Y., Tsuzuki K., Suzuki J., Hayashida R., Kondo K., et al. (2021). Important role of concomitant lymphangiogenesis for reparative angiogenesis in hindlimb ischemia. Arterioscler. Thromb. Vasc. Biol. 41, 2006–2018. doi:10.1161/ATVBAHA.121.316191
Rissanen T. T., Markkanen J. E., Gruchala M., Heikura T., Puranen A., Kettunen M. I., et al. (2003). VEGF-D is the strongest angiogenic and lymphangiogenic effector among VEGFs delivered into skeletal muscle via adenoviruses. Circ. Res. 92, 1098–1106. doi:10.1161/01.RES.0000073584.46059.E3
Roh K., Li H., Freeman R. N., Zazzeron L., Lee A., Zhou C., et al. (2025). Exercise-induced cardiac lymphatic remodeling mitigates inflammation in the aging heart. Aging Cell 24, e70043. doi:10.1111/acel.70043
Schoofs H., Mäkinen T. (2024). Genetic targeting of lymphatic endothelial cells in mice: current strategies and future perspectives. Int. J. Dev. Biol. 68, 189–198. doi:10.1387/ijdb.230215tm
Schwager S., Detmar M. (2019). Inflammation and lymphatic function. Front. Immunol. 10, 308. doi:10.3389/fimmu.2019.00308
Shah F. H., Nam Y. S., Bang J. Y., Hwang I. S., Kim D. H., Ki M., et al. (2025). Targeting vascular endothelial growth receptor-2 (VEGFR-2): structural biology, functional insights, and therapeutic resistance. Arch. Pharm. Res. 48, 404–425. doi:10.1007/s12272-025-01545-1
Shimotsu R., Hotta K., Ikegami R., Asamura T., Tabuchi A., Masamoto K., et al. (2021). Vascular permeability of skeletal muscle microvessels in rat arterial ligation model: in vivo analysis using two-photon laser scanning microscopy. Am. J. Physiol. Regul. Integr. Comp. Physiol. 320, R972–R983. doi:10.1152/ajpregu.00135.2020
Shoji H., Miyakawa T. (2019). Age-related behavioral changes from young to old age in Male mice of a C57BL/6J strain maintained under a genetic stability program. Neuropsychopharmacol. Rep. 39, 100–118. doi:10.1002/npr2.12052
Suarez A. C., Hammel J. H., Munson J. M. (2023). Modeling lymphangiogenesis: pairing in vitro and in vivo metrics. Microcirculation 30, e12802. doi:10.1111/micc.12802
Taketa Y., Tamakoshi K., Hotta K., Maki S., Taguchi T., Takahashi H. (2024). Lymphatic capillarization in different fiber types of rat skeletal muscles with growth and age. Microcirculation 31, e12879. doi:10.1111/micc.12879
Tammela T., Alitalo K. (2010). Lymphangiogenesis: molecular mechanisms and future promise. Cell 140, 460–476. doi:10.1016/j.cell.2010.01.045
Tamura Y., Kawashima T., Ji R. C., Agata N., Itoh Y., Kawakami K. (2024). Histological and biochemical changes in lymphatic vessels after skeletal muscle injury induced by lengthening contraction in Male mice. Physiol. Rep. 12, e15950. doi:10.14814/phy2.15950
Thomas M. M., Trajcevski K. E., Coleman S. K., Jiang M., Di Michele J., O’Neill H. M., et al. (2014). Validity and reliability of a field technique for sweat na+ and K+ analysis during exercise in a hot-humid environment. Physiol. Rep. 2, e12007. doi:10.14814/phy2.12007
Umek N., Horvat S., Cvetko E., Kreft M., Janacek J., Kubinova L., et al. (2019). 3D analysis of capillary network in skeletal muscle of Obese insulin-resistant mice. Histochem. Cell Biol. 152, 323–331. doi:10.1007/s00418-019-01810-7
Wang X., Bove A. M., Simone G., Ma B. (2020). Molecular bases of VEGFR-2-mediated physiological function and pathological role. Front. Cell Dev. Biol. 8, 599281. doi:10.3389/fcell.2020.599281
Waters R. E., Rotevatn S., Li P., Annex B. H., Yan Z. (2004). Voluntary running induces fiber type-specific angiogenesis in mouse skeletal muscle. Am. J. Physiol. Cell Physiol. 287, C1342–C1348. doi:10.1152/ajpcell.00247.2004
Weber E., Sozio F., Borghini A., Sestini P., Renzoni E. (2018). Pulmonary lymphatic vessel morphology: a review. Ann. Anat. 218, 110–117. doi:10.1016/j.aanat.2018.02.011
Xiao L., Liu J., Sun Z., Yin Y., Mao Y., Xu D., et al. (2021). AMPK-Dependent and -independent coordination of mitochondrial function and muscle fiber type by FNIP1. PLoS Genet. 17, e1009488. doi:10.1371/journal.pgen.1009488
Yang Y., Wang X., Wang P. (2023). Signaling mechanisms underlying lymphatic vessel dysfunction in skin aging and possible anti-aging strategies. Biogerontology 24, 727–740. doi:10.1007/s10522-023-10016-3
Yoshioka K., Fujita R., Seko D., Suematsu T., Miura S., Ono Y. (2019). Distinct roles of Zmynd17 and PGC1α in mitochondrial quality control and biogenesis in skeletal muscle. Front. Cell Dev. Biol. 7, 330. doi:10.3389/fcell.2019.00330
Zhang Q., Lu Y., Proulx S. T., Guo R., Yao Z., Schwarz E. M., et al. (2007). Increased lymphangiogenesis in joints of mice with inflammatory arthritis. Arthritis Res. Ther. 9, R118. doi:10.1186/ar2326
Zhang Y., Zhang C., Li L., Liang X., Cheng P., Li Q., et al. (2021). Lymphangiogenesis in renal fibrosis arises from macrophages via VEGF-C/VEGFR3-dependent autophagy and polarization. Cell Death Dis. 12, 109. doi:10.1038/s41419-020-03385-x
Keywords: skeletal muscle, lymphatic vessels, blood capillary, voluntary wheel running, aged mice
Citation: Tamura Y, Kawashima T, Kodama A, Ji R-C, Itoh Y, Agata N and Kawakami K (2025) Voluntary wheel running promotes lymphangiogenesis in slow-twitch muscle in young mice. Front. Physiol. 16:1654445. doi: 10.3389/fphys.2025.1654445
Received: 26 June 2025; Accepted: 26 September 2025;
Published: 10 October 2025.
Edited by:
Irena Levitan, University of Illinois Chicago, United StatesReviewed by:
Nejc Umek, University of Ljubljana, SloveniaFan Li, Mayo Clinic Research in Arizona, United States
Sophie Møller, University of Copenhagen, Denmark
Copyright © 2025 Tamura, Kawashima, Kodama, Ji, Itoh, Agata and Kawakami. 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: Keisuke Kawakami, a2thd2FrYW1pQG9pdGEtdS5hYy5qcA==