Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 06 October 2021
Sec. Exercise Physiology

Acute Effects of Foam Rolling on Hamstrings After Half-Marathon: A Muscle Functional Magnetic Resonance Imaging Study

\r\nDingbo Shu,,Dingbo Shu1,2,3Chuan ZhangChuan Zhang4Siyu Dai,Siyu Dai2,3Shubo WangShubo Wang5Jie Liu,Jie Liu1,3Jianping Ding,*Jianping Ding1,3*
  • 1Department of Radiology, Affiliated Hospital of Hangzhou Normal University, Hangzhou, China
  • 2School of Clinical Medicine, Hangzhou Normal University, Hangzhou, China
  • 3Institute of Sport Medicine, Hangzhou Normal University, Hangzhou, China
  • 4School of Physical Education and Sports, Central China Normal University, Wuhan, China
  • 5Department of Mechanical Engineering, University of Delaware, Newark, DE, United States

Purpose: Foam rolling (FR) is widely used for post-exercise muscle recovery; yet, the effects of FR on skeletal muscle inflammation and microvascular perfusion following prolonged exercise are poorly understood. We aim to address the gap in knowledge by using magnetic resonance imaging (MRI) T2 mapping and intravoxel incoherent motion (IVIM) sequences to study the acute effects of FR on hamstrings following half-marathon running in recreational runners.

Methods: Sixteen healthy recreational marathon runners were recruited. After half-marathon running, FR was performed on the hamstrings on the dominant side, while the other limb served as a control. MRI T2 and IVIM scans were performed bilaterally at baseline (pre-run), 2–3 h after running (post-run), immediately after FR (post-FR0), 30 min after FR (post-FR30) and 60 min after FR (post-FR60). T2, a marker for inflammatory edema, as well as IVIM microvascular perfusion fraction index f for biceps femoris long head (BFL), semitendinosus (ST) and semimembranosus (SM) were determined. Total Quality Recovery (TQR) scale score was also collected.

Results: Both T2 and f were higher at post-run compared to pre-run in all hamstrings on both sides (all p < 0.05; all d > 1.0). For the FR side, T2 decreased, and f increased significantly at post-FR0 and post-FR30 compared to post-run in all muscles (p < 0.05; all d > 0.4) except for f at BFL and SM at post-FR30 (both p > 0.05), though f at BFL was still marginally elevated at post-FR30 (p = 0.074, d = 0.91). Both parameters for all muscles returned to post-run level at post-FR60 (all p > 0.05; all d < 0.4) except for T2 at SM (p = 0.037). In contrast, most MRI parameters were not changed at post-FR0, post-FR30 and post-FR60 compared to post-run for the control side (p < 0.05; d < 0.2). TQR scores were elevated at post-FR0 and post-FR30 compared to post-run (both p < 0.05; both d > 1.0), and returned to the post-run level at post-FR60 (p > 0.99; d = 0.09). Changes in TQR scores compared to post-run at any time points after FR were correlated to T2 for ST at post-FR30 (r = 0.50, p = 0.047) but not T2 for other muscles and any changes in f values.

Conclusions: Hamstrings inflammatory edema and microvascular perfusion were elevated following half-marathon running, which were detectable with MRI T2 mapping and IVIM sequences. FR resulted in acute alleviation in inflammation and greater microvascular perfusion; however, the effects seemed to last only for a short period of time (30–60 min). FR can provide short-term benefits to skeletal muscle after prolonged running.

Introduction

Running is a popular sport activity which can provide physical and psychological benefits to the runners. However, long-distance running induced muscular changes may also contribute to lower extremity musculoskeletal injury risks. Prolonged running has been reported to induce muscular strength decrease, plasma creatine kinase elevation, as well as delayed onset muscle soreness (DOMS) and feeling of fatigue (Hill et al., 2014; Clifford et al., 2016; Takayama et al., 2017), and may cause acute and chronic imbalance on multiple aspects of the human body (Junker and Stöggl, 2019). Repetitive overuse of muscles could lead to microdamage (Hooijmans et al., 2020), which triggers inflammatory responses, and the accumulative effects may ultimately result in muscle injury and dysfunction (Curran et al., 2008; Wan and Shan, 2016).

Greater microvascular perfusion in the hamstrings than the quadriceps was observed following running (Jungmann et al., 2019), suggesting the great involvement of the hamstrings for running activities. Hamstrings are consisted of biarticular muscles biceps femoris long head (BFL), semitendinosus (ST), and semimembranosus (SM), and play an important role in both knee and hip joints movements. Repeated eccentric contractions associated with running could lead to the accumulation of hamstrings microdamage, which puts hamstrings at risk for strain injury (Fyfe et al., 2013). Interestingly, Higashihara et al. (2020) showed that compared to the proximal region, the distal and middle portion of the hamstrings are more vulnerable to marathon-induced injuries. Therefore, post-running recovery methods that are effective in reducing local muscle damage at the middle and distal portion of the hamstrings could be promising remedies to prevent running related hamstrings injuries.

Foam rolling (FR) is a simple, time-efficient and effective post-exercise recovery tool that is able to reduce the risk of sports injury by shortening the recovery period with less negative impact on performance (Škarabot et al., 2015; Su et al., 2016; Jo et al., 2018; Behm and Wilke, 2019; Hodgson et al., 2019). This technique requires individuals to put their own body weight on a foam roller to exert pressure on the soft tissues, and roll back and forth starting at the proximal portion of the muscle (MacDonald et al., 2013). Evidence suggests that FR has a number of physiological benefits, including but not limited to increase blood flow, decrease tissue stiffness, promote muscle recovery, and reduce pain level (Aboodarda et al., 2015; Cavanaugh et al., 2016; Pablos et al., 2020; Phillips et al., 2021). However, whether FR can be used to reduce muscle damage induced by prolonged exercise, such as long-distance running, and therefore lower potential injury risks associated with such exercise modality remains to be addressed.

Therefore, the purpose of this study was to evaluate the acute effects of FR on hamstrings after a half-marathon in recreational runners using muscle functional magnetic resonance imaging (mfMRI). Specifically, we employed T2 sequence and intravoxel incoherent motion (IVIM) to assess inflammatory edema and microvascular perfusion respectively. T2 value reflects the level of muscle inflammation and is related to the severity of muscle injury (Maillard et al., 2004), and has been applied to detect muscle changes induced by exercise and monitor the recovery from muscle injury (Maillard et al., 2004; Higashihara et al., 2020; Hooijmans et al., 2020). MRI-IVIM model is a technique sensitive to the diffusion characteristics of water molecules, which can be used to evaluate the microvascular perfusion in human tissues, including the skeletal muscle (Iima and Le Bihan, 2016; Ogura et al., 2020). We hypothesized that the inflammation and microvascular perfusion as assessed using mfMRI at the mid portion of the hamstrings will both be elevated following half-marathon running, and that FR will result in decreased inflammation and increased microvascular perfusion.

Materials and Methods

Participants

A total of 16 recreational marathon runners (12 males) were recruited from the Zhejiang University Entrepreneurs Outdoor Association. This study received ethical approval (approval number 2021 (E2) - KS - 036) from the Affiliated Hospital of Hangzhou Normal University Research Ethics Committee, and informed consent was obtained from all participants prior to any data collection. Anthropometric characteristics for all participants are summarized in Table 1. Demographic and running related information were obtained from a questionnaire filled by the participant, including name, age, sex, monthly running distance and the pace during marathon. Recreational marathon runners were defined as runners who had not participated in formal distance running training and whose occupation was not running related. The main inclusion criteria were as follows: (1) history of regular long distance running for more than 1 year, (2) three or more regular weekly running sessions, (3) run for longer than 10 km at least once every week, (4) monthly running distance no less than 100 km, (5) having participated in a formal half-marathon or longer race at least once in the past year, and (6) no contraindications for MRI examinations. Those who have had musculoskeletal injuries at the lower extremities were excluded from the study.

TABLE 1
www.frontiersin.org

Table 1. Participant and race characteristics of the 16 half-marathon runners.

Experimental Design

MRI data collections were performed five times for all participants, including baseline (pre-run), 2–3 h after half-marathon running (post-run), immediately after FR (post-FR0), 30 min after FR (post-FR30), and 60 min after FR (post-FR60). Participants lay on the MRI examination table for 30 min before the pre-run and post-run examinations, and were kept on the table throughout the rest of the experiment except for FR intervention. Half-marathon (21.0975 km) running was completed on the flat roads near our laboratory, and running related information, including average speed, time and distance were recorded using a phone app (Ledong Information Technology Co., Ltd., Chengdu, Sichuan). Subjects were asked to keep their RPE scale level below 15 throughout the half-marathon running. FR was performed immediately after the post-run MRI examination was finished. All participants were prohibited from running for long distance for at least 7 days, and refrained from performing any form of exercise for at least 5 days prior to the baseline MRI examination. Participants were asked to complete a Total Quality Recovery (TQR) scale before each MRI examination (except for baseline) to assess their overall sense of recovery (Kenttä, 1996; Kenttä and Hassmén, 1998), and one score was obtained at each time point to indicate the general sense of recovery for the participant. This is a scale ranging from 6 to 20 that was designed to assess the degree of fatigue recovery after exercise, with higher values indicating better sense of recovery. A full experimental protocol is shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Example of the experimental protocol. FR, foam rolling; pre-run, baseline; post-run, 2–3 h after half-marathon running; post-FR0, immediately following FR; post-FR30, 30 min after FR; post-FR60, 30 min after FR; TQR, Total Quality Recovery scale; MRI, magnetic resonance imaging examination.

Foam Rolling Protocol

An experienced researcher (DS) showed all participants about the FR protocol prior to the first MRI data collection, and actual FR was performed after the post-run MRI examination was finished. A foam roller (LiNing, Shanghai, China; size 33 × 14 cm) constructed of a hollow polyvinylchloride pipe and surrounded by ethylene vinyl acetate copolymer was used for all FR massage in the current study. FR was performed on the dominant side for all participants, and the other side served as a control. To roll the hamstrings, the foam roller was placed on the ground, and the participant sat on the foam roller with hands placed on the ground to support their body weight. Participant was instructed to place as much of his/her weight as possible onto the foam roller and to roll gradually and continuously down from the hip to the knee. Each roll was completed in 5 s and the speed was controlled using a metronome. Once the foam roller reached the distal thigh, the participants were instructed to return the foam roller to the starting position (hip). The FR massage protocol was 60 s each with 5 repetitions. A 30 s rest interval was given between two consecutive FR massages, and pressure was controlled at 6–7 on visual analog scale (VAS) score (Aboodarda et al., 2015; Behm et al., 2020). VAS is a numeric scale ranging from 0 to 10 that measures the perceived pain rating, with 0 being no pain at all, and 10 being the worst pain imaginable. During the massage process, the same researcher (DS) assisted the participants to complete the FR as needed. All FR activities were performed at the MRI examination suite.

Muscle Functional Magnetic Resonance Imaging

T2 mapping and IVIM sequences were employed to study the inflammatory edema and microvascular perfusion changes of the hamstrings in responses to half-marathon running and FR. The main perfusion indexes from IVIM sequences are perfusion fraction “f” and pseudo diffusion coefficient “D*.” However, D* value is highly influenced by noise, and accurately obtaining D* value requires high-quality data and long scanning time (Mastropietro et al., 2018), which may not be possible in post-exercise recovery time course studies that only permits short scanning times between consecutive scan sessions. In contrast, f value can be obtained accurately with relatively low-quality images and shorter scanning time. Therefore, f value was chosen to indicate microvascular perfusion from IVIM sequence in the current study. A 1.5 T MRI scanner (Magnetom Avanto, Siemens Healthcare, Germany) with an 18-channel body coil was used for all MRI data collection. Participants lay supine with legs fully extended within the coils. Sponges were added to the gaps in the coil in order to prevent movement artifact. T1 and T2 weighted image sequences were collected bilaterally from the greater trochanter to the medial condyle of the femur. T1-weighted anatomical scans in axial plane (TR = 850 ms, TE = 12 ms, FOV = 400 mm, thickness = 3.5 mm, average = 2, Bandwidth = 182 Hz/Px; scan time 2 min 34 s) were first collected. T1 sequence was collected to provide anatomical reference for IVIM images. T2 mapping sequence (TR = 2,000 ms, TE = 25, 50, 75, and 100 ms, FOV = 400 mm; thickness = 10 mm; average = 1; Bandwidth = 160 Hz/Px; scan time 3 min) was also collected. IVIM sequence (TR = 5,300 ms, TE = 78 ms, FOV = 380 mm, thickness = 3.5 mm, total slice number = 25, Bandwidth = 1,044 Hz/Px; b values [average = 0 (1), 10 (1), 20 (1), 40 (1), 80 (1), 110 (1), 140 (1), 170 (2), 200 (2), 300 (2), 400 (2), 500 (2), 800 (3) s/mm2; scan time 5 min and 20 s] in 3 orthogonal directions were subsequently collected in axial plane. Note that the slice thickness for T1 and IVIM sequences was set at 3.5 mm and was 10 mm for T2 mapping. The selection of different slice thicknesses was done to ensure image quality while keep the scan time short.

Image Analysis

For all participants, the mid-point line for the thigh region covered by the scans (from the greater trochanter to the medial condyle of the femur) was first identified and used as the anatomical reference (Higashihara et al., 2020). Four most adjacent slices to the mid-point line (2 above and 2 below) from the T2 sequence and 12 most adjacent slices to the line (6 above and 6 below) from the IVIM sequence were further processed. Therefore, the processed T2 (total thickness: 40 mm) and IVIM images (total thickness: 42 mm) represent approximately the same region at the mid-thigh. Example for IVIM images is shown in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Example of intravoxel incoherent motion images at different b values.

For T2 images processing, hamstrings, including BFL, ST, and SM were manually segmented on the selected T2-weighted images using syngo MR Workplace (syngo MR E11, Siemens, Germany). Care was taken to avoid the inclusion of any non-contractile tissues. The obtained T2 relaxation times were averaged over the 4 selected slices for each muscle and reported as the T2 value for each participant. Example for image segmentation (A) and color-coded T2 maps at different time points (C) are shown in Figure 3.

FIGURE 3
www.frontiersin.org

Figure 3. Example magnetic resonance images for raw T2 (A) and raw intravoxel incoherent motion (IVIM) b0 images (B) images with region of interest outlined for biceps femoris long head (BFL), semitendinosus (ST) and semimembranosus (SM). Panels (C,D) are color-coded maps for T2 mapping and IVIM f at different time points, respectively, with images on the left side being the limb that received foam rolling intervention. Light colors on the color-coded maps indicate lower values. FR, foam rolling; pre-run, baseline; post-run, 2–3 h after half-marathon running; post-FR0, immediately following FR intervention; post-FR30, 30 min after FR; post-FR60, 60 min after FR.

For IVIM images processing, a free open-source software platform Start MITK Diffusion (Division of Medical Image Computing, German Cancer Research Center) was used. Regions of interest corresponding to BFL, ST, and SM were manually segmented on b0 images, with the visual aid using T1-weighted anatomical images as reference when needed, and care was taken to avoid the inclusion of any non-contractile tissues. Pixels in which the fitting algorithm failed (f lower than 10–4) were excluded from the analysis (Mastropietro et al., 2018). The obtained fitted f values in percentage were averaged over the 12 selected slices reported as the f value for each muscle. Example for image segmentation (B) and color-coded IVIM f map at different time points (D) are shown in Figure 3.

Statistics

Statistical analyses were conducted using SPSS version 20.0 (IBM Corp., Armonk, NY). MRI parameters, including T2 and f values were analyzed using mixed-ANOVA, with sides (FR or control) as the between groups factor and time as the within groups factor. Bonferroni adjustment was applied for post hoc analysis as needed. Importantly, although we made all the pairwise comparisons, our primary interests for statistical analyses were (1) to compare the MRI parameters between the time points after half-marathon running vs. baseline, and (2) to compare the MRI parameters after FR vs. post-run, therefore the results were reported accordingly for the ease of interpretation. Results for the comparisons between each single time point vs. pre-run within each group, as well as between each single time point vs. post-run within each group are presented in Figure 4. TQR scores were also compared using one-way repeated ANOVA with time as the within subject factor, and Bonferroni post hoc adjustment was applied as needed. The relationship between changes in TQR scores and changes in MRI parameters (compared to post-run) at different post FR time points was examined using Pearson correlation. Statistical significance was set at p < 0.05. Mean ± SD was used for descriptive statistics. Cohen’s d (d) was calculated to indicate the effect size whenever applicable, with 0.2, 0.5, and 0.8 represent small, medium, or large effect size, respectively (Cohen, 2013). Post-examination power analysis was performed using G power and effect sizes were calculated for all muscles for T2 and f between post-run and post-FR0.

FIGURE 4
www.frontiersin.org

Figure 4. Bar graphs showing muscle functional magnetic resonance imaging results for foam rolling (FR) side and control (CON) side for the hamstring muscles at five different time points. (A–C) The T2 value changes from T2 mapping sequence. Panels (D–F) show f value changes from intravoxel incoherent motion sequence. BFL, biceps femoris long head; ST, semitendinosus; SM, semimembranosus; pre-run, baseline; post-run, 2–3 h after half-marathon running; post-FR0, immediately following FR intervention; post-FR30, 30 min after FR; post-FR60, 60 min after FR. aTime effect for each single time point vs. pre-run within each group; btime effect for each single time point vs. post-run within each group. Data are presented as Mean ± SE.

Results

Participants

All participants (12 men and 4 women) completed the half-marathon and all TQR and MRI examinations. Physical characteristics and running volume and speed for all participants are summarized in Table 1.

Quantitative Muscle Functional Magnetic Resonance Imaging Parameters

Significant interactions were observed for T2 values and f values for all muscles (all p < 0.001).

Post hoc analysis revealed that T2 values were higher at post-run compared to pre-run for BFL, ST, and SM (all p < 0.001; all d > 1.0; Figure 4) for both FR side and control side. For the FR side, T2 values were lower at post-FR0 and post-FR30 compared to post-run for BFL, ST, and SM (all p < 0.05; all d > 0.5); T2 value was still significantly different at post-FR60 compared to post-run for SM (p = 0.037; d = 0.32); however, T2 value returned to post-run level at post-FR60 for BFL and ST (all p > 0.05; all d < 0.3). For the control side, T2 value was not significantly different at post-FR0, post-FR30 or post-FR60 compared to post-run for BFL, ST, and SM (all p > 0.05; all d < 0.1).

Post hoc analysis suggested that f values were higher at post-run compared to pre-run for BFL, ST, and SM (all p < 0.05; all d > 1.0; Figure 4) for both FR side and control side. For the FR side, f values were higher at post-FR0 compared to post-run for BFL, ST and SM (all p < 0.05; all d > 1.0), and was higher at post-FR30 compared to post-run for ST (p = 0.009; d = 0.77). f value was also higher at post-FR30 compared to post-run for BFL, although the result was only marginally significant (p = 0.074; d = 0.91). f value was not significantly different at post-FR30 compared to post-run for SM (p = 0.485; d = 0.4). f value returned to post-run level at post-FR60 in all three muscles (all p > 0.05; all d < 0.2). For the control side, f values were not significantly different at post-FR0, post-FR30 and post-FR60 compared to post-run for all three muscles (all p > 0.05; all d < 0.2), except for post-FR30 for BFL, which was slightly lower than post-run (p = 0.003; d = 0.35).

TQR scores were higher at post-FR0 and post-FR30 compared to post-run (15.06 ± 0.44 and 13.38 ± 0.50 vs. 12.63 ± 0.72, both p < 0.001, both d > 1.0). TQR score was not significantly different at post-FR60 compared to post-run (12.69 ± 0.6 vs. 12.63 ± 0.72, p > 0.99, d = 0.09). Changes in TQR scores compared to post-run at any time points after FR were not related to changes in f values or T2 values for any muscles (all p > 0.05), except for T2 for ST at post-FR30 (r = 0.50, p = 0.047). Post-examination power analysis suggested that we had high statistical power for this study (>0.81 for both T2 and f for all three muscles).

Discussion

To the best of our knowledge, this is the first study that employed mfMRI to evaluate the efficacy of recovery strategies on skeletal muscle recovery following prolonged running. In this study, MRI T2 mapping and IVIM sequences were used to non-invasively evaluate inflammatory edema and microvascular perfusion in hamstrings before and after FR. The microscopic physiological changes of BFL, ST and SM at five time points of pre-run, post-run, post-FR0, post-FR30, and post-FR60 were demonstrated bilaterally in recreational runners. After half-marathon, T2 and f values for all hamstrings were increased compared to pre-run on both sides. After FR (post-FR0), T2 values decreased while f values increased compared to post-run for BFL, ST muscle, and SM muscle on the FR side, with minimal changes observed on the control side. TQR score was significantly increased, indicating improved sense of recovery following FR. All changes on f values on the FR side and TQR returned to post-run level at post-FR60, indicating that the effect of FR is acute and lasts approximately 30–60 min following one session of FR. Our data suggest that FR massage protocol employed in the current study acutely reduced hamstrings inflammation and increased microvascular perfusion after the running.

After Half-Marathon

We examined the participants shortly after (2–3 h) the half-marathon running and found that T2 values for all three muscles on both sides were significantly higher than pre-run, which is consistent with findings from previous studies that investigated T2 changes in skeletal muscle following running (Higashihara et al., 2020; Hooijmans et al., 2020). In addition, f values for all hamstrings were higher compared to pre-run, which is consistent with a previous study that found increased f value in the skeletal muscle after running (Jungmann et al., 2019). These changes in MRI signals may be caused by the loss of calcium homeostasis, the pro-inflammatory response, as well as the increased oxidative stress after long-distance running (Siegel et al., 2001; McAnulty et al., 2003; Overgaard et al., 2004; Hooijmans et al., 2020). Therefore, our results suggest that MRI T2 and IVIM sequences can be used to monitor the microscopic changes in the muscle, particularly to evaluate the presence of inflammatory edema and microdamage after exercise.

After Foam Rolling

We found that T2 values for BFL, ST, and SM muscles were significantly lower compared to post-run immediately following acute FR intervention for the FR side. In contrast, no significant changes were observed for the control side for any muscles. The likely explanation for the decreased post-FR0 T2 value is the alleviation of inflammatory edema. Myofascia, including blood and lymph, has close interaction with surrounding tissues (Barnes, 1997). An animal study suggested that the water content of myofascia has a high adaptability to mechanical stimulation (Schleip et al., 2011). It is likely that the mechanical force exerted on the muscle during FR intervention caused a water content redistribution between muscle and the surrounding tissues, leading to decreased water volume retained within the muscle tissue and manifest as decreased T2 values. Our results indicate that FR can effectively reduce the level of edema in muscle tissue after long-distance running, which could be beneficial for muscle recovery.

We found that the f values for BFL, ST, and SM were significantly higher than post-run following acute FR intervention. In contrast, there were no significant changes in f for any of the hamstrings in the control group. The results indicate that the microvascular perfusion for BFL, ST, and SM was significantly increased after FR. This is consistent with results from a previous study which suggested that blood flow in the lateral femoral posterior artery was increased after FR (Hotfiel et al., 2017). Similarly, using near-infrared spectroscopy, Soares et al. (2020) found that FR could improve skeletal muscle oxygenation and microvascular reactivity. One possible explanation for the increased microvascular perfusion is that FR changes the viscoelasticity and thixotropy of the fascia, and these changes combined with the heat production and friction induced by FR, would result in increased blood flow in the muscle (Bradbury-Squires et al., 2015; Cheatham et al., 2015). In addition, the increased shear stress on the vascular wall caused by external compression through FR is a powerful stimulus for nitric oxide (NO) production (Okamoto et al., 2014), which is an vasoactive substance known to modulate microvascular function and therefore could lead to microvascular perfusion changes (Kurose et al., 1994; Paniagua et al., 2001). The enhanced microvascular blood flow, as well as greater oxygen delivery as induced by acute FR intervention, could stimulate the mitochondria to accelerate adenosine triphosphate and phosphocreatine filling, a process that contributes to the elimination of metabolic waste products during recovery (Romero-Moraleda et al., 2019). Therefore, the increased microvascular perfusion as a result of FR may play an important role in promoting the recovery of muscle microdamage after exercise.

Results of TQR scores suggested that the participants’ overall sense of recovery was improved immediately following FR. The improvement is likely due to the combination of both physiological and psychological factors. Physiologically, the increased microvascular perfusion after FR could accelerate the elimination of metabolic wastes (Cheatham et al., 2015; Romero-Moraleda et al., 2019). In addition, massage may facilitate the reconstitution of exercise induced Ca2+ -ATPase activity decrease (Ogai et al., 2008) and therefore reduce muscle stiffness and provide better sense of recovery. Psychologically, FR intervention modifies the perception of stretch (Krause et al., 2019), increases pressure pain threshold (Cheatham and Kolber, 2018) and improve mood (Arroyo-Morales et al., 2008), which could all contribute to the increased TQR scores in response to FR.

Thirty and Sixty Minutes After Foam Rolling

At post-FR30, T2 values for BFL, ST, and SM for the FR side were still lower compared to post-run. f values for SM and BFL were still elevated at post-FR30 compared to post-run, but f value for SM returned to post-run level. In contrast, both T2 and f values remained unchanged at post-FR30 compared to post-run and post-FR0 for the control side for all muscles, except for BFL. Sixty minutes after FR, T2 value for BFL returned to post-run level, while remained elevated for ST and SM compared to post-run. f values for BFL, ST, and SM all returned to post-run level at post-FR60. TQR scores were still elevated at post-FR30 and returned to post-run level at post-FR60. Our data, combined with data at post-FR0, indicate that (1) mfMRI is capable of detecting the physiological changes induced by FR, and (2) FR could provide short-term benefits to muscle recovery that generally lasts 30–60 min.

As previously stated, the increase in T2 value may be due to the increase of microvascular perfusion or muscle inflammatory edema after exercise (Esposito et al., 2013; Sharafi et al., 2017; Messer et al., 2018). In this study the lowest T2 for the FR side were observed at post-FR30 for all muscles, while the peak values for f were presented at post-FR0 for all muscles. Combing the pattern changes for both T2 and f, we postulate that FR is capable to both reduce inflammatory edema as well as improve microvascular circulation, with the later effect being more acute while persistently reduce edema for a relatively longer period. It is possible that the improvement in inflammatory edema and microcirculation could facilitate the microdamage recovery experienced by the hamstrings following half-marathon running, and therefore reduce future injury risks and shorten the recovery time needed to prepare for the next race; however, more studies with longer follow-up times are needed in order to confirm this inference. At post-FR60, T2 value remained lower than post-run for ST and SM, indicating that the effects of FR to reduce inflammatory edema could last for at least an hour for certain muscles. Studies that investigate the time course for FR to reduce muscle inflammatory edema are needed to determine if such effect could persist even longer than an hour. In addition, it is possible that more frequently applied FR during recovery could shorten the inflammation process through the removal of interstitial fluid, the dissolved proteins and lymphocytes via the lymphatic system, although clearly more research is needed to confirm this notion (Crane et al., 2012; Macdonald et al., 2014; Pablos et al., 2020).

Our results that microvascular perfusion was still elevated at post-FR30 for the FR side is consistent with a previous study that measured blood flow response to FR. Using Doppler ultrasound, Hotfiel et al. (2017) measured arterial perfusion after FR intervention at two time points (post-FR0 and post-FR30), and found that the increase in blood flow was more prominent at post-FR0 than post-FR30, but remains higher at post-FR30 compared to baseline. In addition, in our study we found that all f values return to post-run level at post-FR60. Therefore, it can be concluded that effects of FR to increase blood flow is somewhat transient (for about 30 min), and future studies that manipulate different variables of the FR protocols, such as rest interval or the number of repetitions, to determine whether more significant blood flow effects can be achieved are warranted.

Overall, changes in TQR scores were not significantly related to either MRI parameters for any of the three measured muscles. The lack of association between changes in TQR scores and MRI parameters may be attributed to several possibilities. First, since this is not the primary goal of this study and we only enrolled 16 participants, the relatively small sample size may not provide enough power to detect the true relationship between the two measurements. Second, underlying factors driving the improved TQR scores following FR may be multifaceted and cannot be specifically attributed to the improved microvascular perfusion or alleviated inflammatory edema as measured in the current study. It was previously reported that FR can induce endorphin release (Weerapong et al., 2005) and increase pain threshold (Cheatham and Kolber, 2018), which are probably important factors that could modulate participant’s overall sense of recovery. It should be noted that while the TQR scale used in this study only evaluates the overall sense of recovery for the participants, studies have demonstrated that unilateral FR can elicit crossover effect on the contralateral limb to alleviate pain (Aboodarda et al., 2015; Cavanaugh et al., 2017; Behm and Wilke, 2019). Thus future studies with larger sample size and sham group (light rolling) are needed to determine whether improvement in microvascular perfusion and inflammation are contributing factors to the better sense of recovery following FR.

To sum up, the use of FR after half-marathon could provide both physiological benefits as well as overall sense of recovery at least for a short period of time. In this study we only evaluated the effects of FR on hamstrings, but it can be inferred that FR may also be used to facilitate the acute recovery of other muscles following strenuous exercise. Importantly, it should be pointed out that FR should be avoided in case of acute muscle injury. External loading exerted through FR intervention may further stretch the injured muscle tissue.

Limitation

This study has several limitations. First, as previously discussed, high quality images and longer scan time are required to accurately determine D*; however, due to the time constraint of our study protocol, obtaining high-quality images to calculate D* may not be feasible. Future studies with better image quality are needed to determine the D* changes following FR intervention. Second, TQR scale used in the current study is designed to assess the overall sense of recovery. Due to a lack of control group, it is possible that the increased TQR score following FR was partly due to passive recovery. However, this is unlikely to be the case, as TQR score returned to post-run level 60 min after FR. Third, we monitored the participants up to 60 min following FR, therefore it is unclear how the microscopic physiological changes observed in the current study can be translated into better longitudinal recovery for recreational runners. Future studies with longer follow-up periods are needed to determine whether FR provides long-term benefits for post-exercise recovery and injury prevention. Last but not the least, the relatively small number of participants precludes the possibility of gaining conclusive evidence regarding whether improvement in microvascular perfusion and inflammation are driving factors for the alleviated TQR scores. Studies with larger sample size and sham group (light rolling) are therefore needed before definitive conclusion can be drawn.

Conclusion

Hamstrings inflammatory edema and microvascular perfusion were elevated following half-marathon running, which were detectable with mfMRI T2 mapping and IVIM sequences. FR resulted in acute alleviation in inflammation and greater microvascular perfusion, as well as improvement in overall sense of recovery; however, the effects only seemed to last for a short period of time (30–60 min). FR can provide short-term benefits to skeletal muscle after prolonged running. More research is needed to explore whether the short-term benefits seen in the current study can be translated into reduced injury risks and faster recovery.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Affiliated Hospital of Hangzhou Normal University Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

DS designed the study, performed data collection and data analysis, and drafted the initial manuscript. CZ designed the study and performed data analysis and interpretation. SD designed the study and assisted with data collection and data analysis. SW and JL assisted with data collection, analysis, and interpretation. JD designed the study, performed data interpretation and analysis, and contributed to the structure, planning, and execution for this project. All authors approved the final manuscript as submitted.

Funding

This study was funded by the Medical and Health Science & Technology Project of Zhejiang Province (Grant No. 2020KY711) and the Key Medical Disciplines of Hangzhou (Grant No. YDYX).

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.

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

Aboodarda, S. J., Spence, A. J., and Button, D. C. (2015). Pain pressure threshold of a muscle tender spot increases following local and non-local rolling massage. BMC Musculoskelet Disord. 16:265. doi: 10.1186/s12891-015-0729-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Arroyo-Morales, M., Olea, N., Martínez, M. M., Hidalgo-Lozano, A., Ruiz-Rodríguez, C., and Díaz-Rodríguez, L. (2008). Psychophysiological effects of massage-myofascial release after exercise: a randomized sham-control study. J. Altern. Complement. Med. 14, 1223–1229. doi: 10.1089/acm.2008.0253

PubMed Abstract | CrossRef Full Text | Google Scholar

Barnes, M. F. (1997). The basic science of myofascial release: morphologic change in connective tissue. J. Bodyw. Mov. Ther. 1, 231–238.

Google Scholar

Behm, D. G., Alizadeh, S., Hadjizadeh Anvar, S., Mahmoud, M. M. I., Ramsay, E., Hanlon, C., et al. (2020). Foam rolling prescription: a clinical commentary. J. Strength Cond. Res. 34, 3301–3308. doi: 10.1519/JSC.0000000000003765

PubMed Abstract | CrossRef Full Text | Google Scholar

Behm, D. G., and Wilke, J. (2019). Do self-myofascial release devices release myofascia? rolling mechanisms: a narrative review. Sports Med. 49, 1173–1181. doi: 10.1007/s40279-019-01149-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradbury-Squires, D. J., Noftall, J. C., Sullivan, K. M., Behm, D. G., Power, K. E., and Button, D. C. (2015). Roller-massager application to the quadriceps and knee-joint range of motion and neuromuscular efficiency during a lunge. J. Athl. Train. 50, 133–140. doi: 10.4085/1062-6050-49.5.03

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavanaugh, M. T., Döweling, A., Young, J. D., Quigley, P. J., Hodgson, D. D., Whitten, J. H., et al. (2016). An acute session of roller massage prolongs voluntary torque development and diminishes evoked pain. Eur. J. Appl. Physiol. 117, 109–117. doi: 10.1007/s00421-016-3503-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavanaugh, M. T., Döweling, A., Young, J. D., Quigley, P. J., Hodgson, D. D., Whitten, J. H., et al. (2017). An acute session of roller massage prolongs voluntary torque development and diminishes evoked pain. Eur. J. Appl. Physiol. 117, 109–117.

Google Scholar

Cheatham, S. W., and Kolber, M. J. (2018). Does roller massage with a foam roll change pressure pain threshold of the ipsilateral lower extremity antagonist and contralateral muscle groups? An exploratory study. J. Sport Rehabil. 27, 165–169. doi: 10.1123/jsr.2016-0196

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheatham, S. W., Kolber, M. J., Cain, M., and Lee, M. (2015). The effects of self-myofacial release using a foam roll or roller massager on joint range of motion, muscle revovery, and performance: a systematic review. Int. J. Sports Phys. Ther. 10, 827–838.

Google Scholar

Clifford, T., Allerton, D. M., Brown, M. A., Harper, L., Horsburgh, S., Keane, K. M., et al. (2016). Minimal muscle damage after a marathon and no influence of beetroot juice on inflammation and recovery. Appl. Physiol. Nutr. Metab. 42, 263–270. doi: 10.1139/apnm-2016-0525

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, J. (2013). Statistical Power Analysis For The Behavioral Sciences. Cambridge, MA: Academic press.

Google Scholar

Crane, J. D., Ogborn, D. I., Cupido, C., Melov, S., Hubbard, A., Bourgeois, J. M., et al. (2012). Massage therapy attenuates inflammatory signaling after exercise-induced muscle damage. Sci. Transl. Med. 4:119ra13. doi: 10.1126/scitranslmed.3002882

PubMed Abstract | CrossRef Full Text | Google Scholar

Curran, P. F., Fiore, R. D., and Crisco, J. J. (2008). A comparison of the pressure exerted on soft tissue by 2 myofascial rollers. J. Sport Rehabil. 17, 432–442. doi: 10.1123/jsr.17.4.432

PubMed Abstract | CrossRef Full Text | Google Scholar

Esposito, A., Campana, L., Palmisano, A., De Cobelli, F., Canu, T., Santarella, F., et al. (2013). Magnetic resonance imaging at 7T reveals common events in age-related sarcopenia and in the homeostatic response to muscle sterile injury. PLoS One 8:e59308. doi: 10.1371/journal.pone.0059

CrossRef Full Text | Google Scholar

Fyfe, J. J., Opar, D. A., Williams, M. D., and Shield, A. J. (2013). The role of neuromuscular inhibition in hamstring strain injury recurrence. J. Electromyogr. Kinesiol. 23, 523–530. doi: 10.1016/j.jelekin.2012.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Higashihara, A., Nakagawa, K., Inami, T., Fukano, M., Iizuka, S., Maemichi, T., et al. (2020). Regional differences in hamstring muscle damage after a marathon. PLoS One 15:e0234401. doi: 10.1371/journal.pone.0234401

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, J. A., Howatson, G., van Someren, K. A., Walshe, I., and Pedlar, C. R. (2014). Influence of compression garments on recovery after marathon running. J. Strength Cond. Res. 28, 2228–2235. doi: 10.1519/JSC.0000000000000469

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodgson, D. D., Quigley, P. J., Whitten, J. H. D., Reid, J. C., and Behm, D. G. (2019). Impact of 10-minute interval roller massage on performance and active range of motion. J. Strength Cond. Res. 33, 1512–1523. doi: 10.1519/JSC.0000000000002271

PubMed Abstract | CrossRef Full Text | Google Scholar

Hooijmans, M. T., Monte, J. R. C., Froeling, M., van den Berg-Faay, S., Aengevaeren, V. L., Hemke, R., et al. (2020). Quantitative MRI reveals microstructural changes in the upper leg muscles after running a marathon. J. Magn. Reson. Imaging 52, 407–417. doi: 10.1002/jmri.27106

PubMed Abstract | CrossRef Full Text | Google Scholar

Hotfiel, T., Swoboda, B., Krinner, S., Grim, C., Engelhardt, M., Uder, M., et al. (2017). Acute effects of lateral thigh foam rolling on arterial tissue perfusion determined by spectral doppler and power doppler ultrasound. J. Strength Cond. Res. 31, 893–900. doi: 10.1519/JSC.0000000000001641

PubMed Abstract | CrossRef Full Text | Google Scholar

Iima, M., and Le Bihan, D. (2016). Clinical intravoxel incoherent motion and diffusion mr imaging: past. present, and future. Radiology 278, 13–32. doi: 10.1148/radiol.2015150244

PubMed Abstract | CrossRef Full Text | Google Scholar

Jo, E., Juache, G. A., Saralegui, D. E., Weng, D., and Falatoonzadeh, S. (2018). The acute effects of foam rolling on fatigue-related impairments of muscular performance. Sports 6:112. doi: 10.3390/sports6040112

PubMed Abstract | CrossRef Full Text | Google Scholar

Jungmann, P. M., Pfirrmann, C., and Federau, C. (2019). Characterization of lower limb muscle activation patterns during walking and running with Intravoxel Incoherent Motion (IVIM) MR perfusion imaging. Magn. Reson. Imaging 63, 12–20. doi: 10.1016/j.mri.2019.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Junker, D., and Stöggl, T. (2019). The training effects of foam rolling on core strength endurance. balance, muscle performance and range of motion: a randomized controlled trial. J. Sports Sci. Med. 18, 229–238.

Google Scholar

Kenttä, G. (1996). Överträningssyndrom: En Psykofysiologisk Process (Overtraining: A Psychophysiological Process). Luleå: Högskolan i Luleå.

Google Scholar

Kenttä, G., and Hassmén, P. (1998). Overtraining and recovery. a conceptual model. Sports Med. 26, 1–16. doi: 10.2165/00007256-199826010-00001

PubMed Abstract | CrossRef Full Text | Google Scholar

Krause, F., Wilke, J., Niederer, D., Vogt, L., and Banzer, W. (2019). Acute effects of foam rolling on passive stiffness, stretch sensation and fascial sliding: a randomized controlled trial. Hum. Mov. Sci. 67:102514. doi: 10.1016/j.humov.2019.102514

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurose, I., Wolf, R., Grisham, M. B., and Granger, D. N. (1994). Modulation of ischemia/reperfusion-induced microvascular dysfunction by nitric oxide. Circ. Res. 74, 376–382. doi: 10.1161/01.res.74.3.376

CrossRef Full Text | Google Scholar

Macdonald, G. Z., Button, D. C., Drinkwater, E. J., and Behm, D. G. (2014). Foam rolling as a recovery tool after an intense bout of physical activity. Med. Sci. Sports Exerc. 46, 131–142. doi: 10.1249/MSS.0b013e3182a123db

PubMed Abstract | CrossRef Full Text | Google Scholar

MacDonald, G. Z., Penney, M. D., Mullaley, M. E., Cuconato, A. L., Drake, C. D., Behm, D. G., et al. (2013). An acute bout of self-myofascial release increases range of motion without a subsequent decrease in muscle activation or force. J. Strength Cond. Res. 27, 812–821. doi: 10.1519/JSC.0b013e31825c2bc1

PubMed Abstract | CrossRef Full Text | Google Scholar

Maillard, S. M., Jones, R., Owens, C., Pilkington, C., Woo, P., Wedderburn, L. R., et al. (2004). Quantitative assessment of MRI T2 relaxation time of thigh muscles in juvenile dermatomyositis. Rheumatology 43, 603–608. doi: 10.1093/rheumatology/keh130

PubMed Abstract | CrossRef Full Text | Google Scholar

Mastropietro, A., Porcelli, S., Cadioli, M., Rasica, L., Scalco, E., Gerevini, S., et al. (2018). Triggered intravoxel incoherent motion MRI for the assessment of calf muscle perfusion during isometric intermittent exercise. NMR Biomed. 31:e3922. doi: 10.1002/nbm.3922

PubMed Abstract | CrossRef Full Text | Google Scholar

McAnulty, S. R., McAnulty, L. S., Nieman, D. C., Morrow, J. D., Utter, A. C., Henson, D. A., et al. (2003). Influence of carbohydrate ingestion on oxidative stress and plasma antioxidant potential following a 3 h run. Free Radic. Res. 37, 835–840. doi: 10.1080/1071576031000136559

PubMed Abstract | CrossRef Full Text | Google Scholar

Messer, D. J., Bourne, M. N., Williams, M. D., Al Najjar, A., and Shield, A. J. (2018). Hamstring muscle use in women during hip extension and the nordic hamstring exercise: a functional magnetic resonance imaging study. J. Orthop. Sports Phys. Ther. 48, 607–612. doi: 10.2519/jospt.2018.7748

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogai, R., Yamane, M., Matsumoto, T., and Kosaka, M. (2008). Effects of petrissage massage on fatigue and exercise performance following intensive cycle pedalling. Br. J. Sports Med. 42, 834–838. doi: 10.1136/bjsm.2007.044396

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogura, A., Sotome, H., Asai, A., and Fuju, A. (2020). Evaluation of capillary blood volume in the lower limb muscles after exercise by intravoxel incoherent motion. Radiol. Med. 125, 474–480. doi: 10.1007/s11547-020-01163-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, T., Masuhara, M., and Ikuta, K. (2014). Acute effects of self-myofascial release using a foam roller on arterial function. J. Strength Cond. Res. 28, 69–73. doi: 10.1519/JSC.0b013e31829480f5

PubMed Abstract | CrossRef Full Text | Google Scholar

Overgaard, K., Fredsted, A., Hyldal, A., Ingemann-Hansen, T., Gissel, H., and Clausen, T. (2004). Effects of running distance and training on Ca2+ content and damage in human muscle. Med. Sci. Sports Exerc. 36, 821–829. doi: 10.1249/01.mss.0000126468.65714.60

PubMed Abstract | CrossRef Full Text | Google Scholar

Pablos, A., Ceca, D., Jorda, A., Rivera, P., Colmena, C., Elvira, L., et al. (2020). Protective effects of foam rolling against inflammation and notexin induced muscle damage in rats. Int. J. Med. Sci. 17, 71–81. doi: 10.7150/ijms.37981

PubMed Abstract | CrossRef Full Text | Google Scholar

Paniagua, O. A., Bryant, M. B., and Panza, J. A. (2001). Role of endothelial nitric oxide in shear stress-induced vasodilation of human microvasculature: diminished activity in hypertensive and hypercholesterolemic patients. Circulation 103, 1752–1758. doi: 10.1161/01.cir.103.13.1752

CrossRef Full Text | Google Scholar

Phillips, J., Diggin, D., King, D. L., and Sforzo, G. A. (2021). Effect of varying self-myofascial release duration on subsequent athletic performance. J. Strength Cond. Res. 35, 746–753. doi: 10.1519/JSC.0000000000002751

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero-Moraleda, B., González-García, J., Cuéllar-Rayo, Á, Balsalobre-Fernández, C., Muñoz-García, D., and Morencos, E. (2019). Effects of vibration and non-vibration foam rolling on recovery after exercise with induced muscle damage. J. Sports Sci. Med. 18, 172–180.

Google Scholar

Schleip, R., Duerselen, L., Vleeming, A., Naylor, I. L., Lehmann-Horn, F., Zorn, A., et al. (2011). Strain hardening of fascia: static stretching of dense fibrous connective tissues can induce a temporary stiffness increase accompanied by enhanced matrix hydration. J. Bodyw. Mov. Ther. 16, 94–100. doi: 10.1016/j.jbmt.2011.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharafi, A., Chang, G., and Regatte, R. R. (2017). Bi-component T1ρ and T2 relaxation mapping of skeletal muscle in-Vivo. Sci Rep. 7:14115. doi: 10.1038/s41598-017-14581-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegel, A. J., Stec, J. J., Lipinska, I., Van Cott, E. M., Lewandrowski, K. B., Ridker, P. M., et al. (2001). Effect of marathon running on inflammatory and hemostatic markers. Am. J. Cardiol. 88:A9. doi: 10.1016/s0002-9149(01)01909-9

CrossRef Full Text | Google Scholar

Škarabot, J., Beardsley, C., and Štirn, I. (2015). Comparing the effects of self-myofascial release with static stretching on ankle range-of-motion in adolescent athletes. Int. J. Sports Phys. Ther. 10, 203–212.

Google Scholar

Soares, R. N., Inglis, E. C., Khoshreza, R., Murias, J. M., and Aboodarda, S. J. (2020). Rolling massage acutely improves skeletal muscle oxygenation and parameters associated with microvascular reactivity: the first evidence-based study. Microvasc Res. 132:104063. doi: 10.1016/j.mvr.2020.104063

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, H., Chang, N. J., Wu, W. L., Guo, L. Y., and Chu, I. H. (2016). Acute effects of foam rolling. static stretching, and dynamic stretching during warm-ups on muscular flexibility and strength in young adults. J. Sport Rehabil. 26, 469–477. doi: 10.1123/jsr.2016-0102

PubMed Abstract | CrossRef Full Text | Google Scholar

Takayama, F., Aoyagi, A., Shimazu, W., and Nabekura, Y. (2017). Effects of marathon running on aerobic fitness and performance in recreational runners one week after a race. J. Sports Med. 2017:9402386. doi: 10.1155/2017/9402386

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, B., and Shan, G. (2016). Biomechanical modeling as a practical tool for predicting injury risk related to repetitive muscle lengthening during learning and training of human complex motor skills. Springerplus. 12:441. doi: 10.1186/s40064-016-2067-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Weerapong, P., Hume, P. A., and Kolt, G. S. (2005). The mechanisms of massage and effects on performance, muscle recovery and injury prevention. Sports Med. 35, 235–256. doi: 10.2165/00007256-200535030-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: recreational marathon runners, microvascular perfusion, intravoxel incoherent motion, inflammatory edema, T2 mapping

Citation: Shu D, Zhang C, Dai S, Wang S, Liu J and Ding J (2021) Acute Effects of Foam Rolling on Hamstrings After Half-Marathon: A Muscle Functional Magnetic Resonance Imaging Study. Front. Physiol. 12:723092. doi: 10.3389/fphys.2021.723092

Received: 10 June 2021; Accepted: 15 September 2021;
Published: 06 October 2021.

Edited by:

Hamdi Chtourou, University of Sfax, Tunisia

Reviewed by:

José Ramón Alvero Cruz, University of Malaga, Spain
Jochen Zange, Helmholtz Association of German Research Centers (HZ), Germany
David George Behm, Memorial University of Newfoundland, Canada

Copyright © 2021 Shu, Zhang, Dai, Wang, Liu and Ding. 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: Jianping Ding, jpding@vip.sina.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.