Gastrocnemius Medialis Contractile Behavior Is Preserved During 30% Body Weight Supported Gait Training

Rehabilitative body weight supported gait training aims at restoring walking function as a key element in activities of daily living. Studies demonstrated reductions in muscle and joint forces, while kinematic gait patterns appear to be preserved with up to 30% weight support. However, the influence of body weight support on muscle architecture, with respect to fascicle and series elastic element behavior is unknown, despite this having potential clinical implications for gait retraining. Eight males (31.9 ± 4.7 years) walked at 75% of the speed at which they typically transition to running, with 0% and 30% body weight support on a lower-body positive pressure treadmill. Gastrocnemius medialis fascicle lengths and pennation angles were measured via ultrasonography. Additionally, joint kinematics were analyzed to determine gastrocnemius medialis muscle–tendon unit lengths, consisting of the muscle's contractile and series elastic elements. Series elastic element length was assessed using a muscle–tendon unit model. Depending on whether data were normally distributed, a paired t-test or Wilcoxon signed rank test was performed to determine if body weight supported walking had any effects on joint kinematics and fascicle–series elastic element behavior. Walking with 30% body weight support had no statistically significant effect on joint kinematics and peak series elastic element length. Furthermore, at the time when peak series elastic element length was achieved, and on average across the entire stance phase, muscle–tendon unit length, fascicle length, pennation angle, and fascicle velocity were unchanged with respect to body weight support. In accordance with unchanged gait kinematics, preservation of fascicle–series elastic element behavior was observed during walking with 30% body weight support, which suggests transferability of gait patterns to subsequent unsupported walking.

Introduction

Orthopedic and neurological rehabilitation regimens often involve patients performing gait training with body weight support (BWS) in an attempt to retrain “natural” walking gait function. Whilst overhead suspension systems are largely employed to promote gait rehabilitation from neurologic disorders (Apte et al., 2018), lower-body positive pressure (LBPP) treadmills are frequently used following orthopedic injuries to re-expose patients to walking whilst bearing progressively greater proportions of their body weight (Quigley et al., 2000; Mishra, 2015). In order to restore gait function, movement patterns should be as similar, and thus transferable to daily activities, as possible albeit with a reduction of lower limb muscle and joint forces (Cutuk et al., 2006). Studies assessing LBPP have demonstrated that whilst ground reaction forces are reduced (Eastlack et al., 2005; Cutuk et al., 2006; Grabowski, 2010), gait kinematics are largely preserved (Apte et al., 2018).

During normal walking, mechanical energy is largely conserved due to the pendulum-like exchange between potential and kinetic energy (Cavagna et al., 2000). Despite this, additional mechanical work by the muscle–tendon unit (MTU) is required to sustain the movement of the body's center of mass. However, walking with BWS reduces the total mechanical energy of the center of mass, and thus presumably requires less force and work from the MTU to vertically support and accelerate the body (Cavagna et al., 2000; Pavei et al., 2015). In fact, significant reductions in the metabolic cost of locomotion have been observed (Farley and McMahon, 1992; Grabowski et al., 2005; Pavei et al., 2015). Furthermore, reductions of knee joint contact forces (Patil et al., 2013), ankle joint moments (Lewek, 2011; Goldberg and Stanhope, 2013; Fischer and Wolf, 2015) ankle joint angular momentum (McGowan et al., 2008) and ankle joint power (Lewek, 2011) have been reported during unloading. Despite reduced kinetic and metabolic requirements for vertical body support and forward acceleration, LBPP (unless BWS is >75%) has been reported to not induce significant differences in spatio-temporal gait parameters such as cadence, stride duration (Grabowski, 2010; Patil et al., 2013) and stride length (Quigley et al., 2000; Cutuk et al., 2006; Patil et al., 2013), nor range of ankle (Quigley et al., 2000; Cutuk et al., 2006) and knee (Eastlack et al., 2005; Cutuk et al., 2006) joint motion. In addition, whilst muscle activity patterns appear unchanged, lower limb muscle activity is reduced during LBPP-treadmill gait (Quigley et al., 2000; Liebenberg et al., 2011; Fischer et al., 2015) with the plantar flexor muscles being particularly susceptible to manipulations of body weight (McGowan et al., 2008). This demonstrates their critical role in human locomotion by providing the majority of the force necessary for vertical body weight support and horizontal propulsion (Neptune et al., 2001; Anderson and Pandy, 2003; McGowan et al., 2008). To gain a better understanding of the plantar flexors' response to different locomotor tasks, ultrasound imaging is a convenient technique to visualize architectural changes, which help to draw conclusions about muscle function.

Ultrasonic visualization of muscle fascicle behavior during locomotion without BWS has not only demonstrated the importance of the storage and release of elastic energy in the Achilles tendon for running and walking (Fukunaga, 2001; Lichtwark et al., 2007), but also that the plantar flexor muscles modulate their behavior depending on gait type, and speed (Farris and Sawicki, 2012). In fact, increased walking speeds have been shown to increase gastrocnemius medialis (GM) shortening velocities (Farris and Sawicki, 2012), and to shorten soleus fascicles (Lai et al., 2015), thereby impairing plantar flexors force generation due to shifting the force–length–velocity relationship toward less favorable contractile conditions (Neptune and Sasaki, 2005; Arnold et al., 2013). However, it is unknown whether walking with BWS modulates fascicle and series elastic element (SEE) behavior to meet the reduced locomotor demands (Richter et al., 2017). Knowledge of any changes in GM's muscle architecture (primarily fascicle length and pennation angle) in addition to fascicle shortening velocity, which affect the force–length–velocity relationships, would facilitate inference of the mechanisms determining mechanical power generation when BWS is applied. Whereas, preservation of fascicle contraction behavior concurrent with preservation of gait kinematics would support the validity of rehabilitative gait training with BWS.

30% BWS is typically recommended for rehabilitative re-introduction to walking and running, due to the preservation of kinematic and spatio-temporal gait parameters (Fischer and Wolf, 2015; Apte et al., 2018) in addition to muscle activation patterns (Neal et al., 2016; Hansen et al., 2017). As during early postoperative rehabilitation patients usually start with recovering their walking function, the present study focuses on walking with BWS. Increasing BWS is known to result in walk-to-run transitions occurring at slower absolute walking speeds (Kram et al., 1997), but similar Froude number, a dimensionless number embedding gait speed, leg lengths and gravitational acceleration (in the present paper expressed as BWS) (Kram et al., 1997; Labini et al., 2011). Thus, to obtain mechanically equivalent speeds (i.e., a similar walking speed relative to the preferred walk-to-run transition speed) at different BWS levels, walking speeds should be adjusted to the same Froude number (Kram et al., 1997; Minetti, 2001; Vaughan and O'Malley, 2005), which requires a reduction in absolute walking speed.

Therefore, the aim of the present study was to determine via ultrasonography GM's fascicle–SEE behavior during walking at mechanically equivalent speeds, namely 75% of the preferred walk-to-run transition speed (PTS), on a LBPP treadmill, with, and without 30% BWS.

It was hypothesized that during walking with BWS (i.e., where forces acting on the SEE are reduced) peak SEE length decreases and is compensated for by longer fascicles and/or smaller pennation angles, rather than by a shorter MTU as ankle and knee joint kinematics are reported to be preserved.

Materials and Methods

Participants

Eight healthy male volunteers (mean ± standard deviation: 31.9 ± 4.7 years, 178.4 ± 5.7 cm heights, 94.2 ± 5.6 cm leg lengths, and 73.5 ± 7.3 kg body masses) with treadmill running experience provided informed written consent to participate in this observational study, which received approval from the “Ärztekammer Nordrhein” Ethical Committee of Düsseldorf, Germany. The study was conducted in the Physiology Laboratory of the Institute of Aerospace Medicine at the German Aerospace Center in Cologne, where all participants underwent a standard medical examination. Exclusion criteria included any cardiovascular, musculoskeletal or neurological disorders within the previous 2 years in addition to any lower limb surgery that may affect MTU behavior.

Study Design and Experimental Protocol

Participants attended the laboratory on a single occasion and walked on an Anti-Gravity Treadmill (AlterG; AlterG^®, M320, Fremont, USA; Figure 1), an LBPP treadmill, with 0% BWS and thereafter with ~30% BWS (recommended load for rehabilitative gait training; Hesse, 2008; Fischer and Wolf, 2015). Before each trial, participants familiarized themselves until they have acclimatized to the BWS level and the predefined walking speed (~4 min). After another 2 min accommodation time given to produce reproducible gait kinematics (Karamanidis et al., 2003) and hence a total warm-up time of ~6 min, which is further required for the Achilles tendon to achieve a relatively stable steady-state behavior (Hawkins et al., 2009), data were collected for 30 s. Blinding of participants was not applicable due to the nature of the experimental set-up.

FIGURE 1

Figure 1. Experimental set-up. Participant walking on the lower-body positive pressure treadmill (the AlterG) with an ultrasound transducer attached to the midbelly of the gastrocnemius medialis muscle and electrogoniometers (added in green to accentuate placement) to measure ankle and knee joint angles.

TABLE 1

Table 1. Means and standard deviations of kinematic outcome measures while participants walked at 75% of their preferred walk-to-run transition speed with 0% and 30% body weight support.

Walking speeds were defined as 75% of the preferred walk-to-run transition speed (PTS) expressed as a Froude number (PTS_FR). PTS_FR was estimated by fitting an exponential regression equation $(PT{S}_{FR}(a)=1.183e^{-5.952a}+0.4745)$ with a least-squares method (r² = 0.99) to the experimental data of Kram et al. (1997) using the resulting acceleration (a) as the independent variable. Hence, for a = 0.7 g (g = 9.81 m·s⁻²), a PTS_FR value of 0.49 was obtained. By accounting for the participants' leg lengths (l), measured from the greater trochanter to the ground, the individual $PTS(a)=\sqrt{PT{S}_{FR}(a)·a·l}$ expressed in meters per second was determined resulting in walking speeds of 1.58 ± 0.05 m·s⁻¹ at 0% BWS, and 1.34 ± 0.04 m·s⁻¹ at 30% BWS.

The AlterG was enclosed within a sealed height-adjustable chamber, which allowed air pressure to increase inside the chamber and generated an additional vertical buoyant force to produce controlled and stable BWS levels. A seal between the participant and the chamber was created through a neoprene kayak-type skirt that could be zipped into the aperture of the chamber.

Joint Kinematics

Knee and ankle joint angles were recorded using a twin-axis (Penny and Giles Biometrics Ltd., Blackwood Gwent, UK) and a custom-made potentiometer based electrogoniometer, respectively. The end blocks of the knee electrogoniometer were placed along the leg from the greater trochanter to the lateral femur epicondyle and along the leg from the lateral epicondyle of the femur to the lateral malleolus. The end blocks of the ankle electrogoniometer were placed along the leg from the lateral femur epicondyle to the lateral malleolus and from the lateral malleolus to the most distal end of the fifth metatarsal. Before each walking trial, a reference measurement was taken in the anatomical neutral position to define the 0° joint angles. Data were sampled at a frequency of 1,500 Hz via the TeleMyo 2400 G2 Telemetry System (Noraxon USA., Inc., Scottsdale, USA) and MyoResearch XP software (Master Edition 1.08.16).

Spatio-Temporal Parameters

To determine gait cycle events and thereby define stance phases of the left leg, participant plantar pressure was measured (83 Hz) via insoles (Novel GmbH, loadsol^® version 1.4.60, Munich, Germany). Touchdown and toe-off were automatically detected using a 20 N threshold for 0.1 s via a custom-made script (MATLAB R2018a, MathWorks, Inc., Natick, United States). Insole and electrogoniometer signals were time-synchronized via recording of a rectangular pulse generated by pressing on a custom-made pedal.

GM Muscle Fascicle Length and Pennation Angle

Real-time B-mode ultrasound (Prosound α7, ALOKA, Tokyo, Japan) captured at 73 Hz using a T-shaped 6-cm linear array transducer (13 MHz) was performed over the midbelly of the left GM muscle. Transducer position was standardized by determining the intersection of the mediolateral and proximodistal midline of the GM and aligning the transducer longitudinally to the fascicles, while transducer movement was minimized by using a custom-made cast, which was secured with elastic Velcro. Ultrasound recordings and electrogoniometer signals were time-synchronized via a rectangular pulse generated by a hand switch, which was recorded synchronously through the electrocardiography channel of the ultrasound and the MyoResearch software. A semi-automatic tracking algorithm (UltraTrack Software, version 4.2; Farris and Lichtwark, 2016) was used to quantify muscle fascicle length (distance between the insertion of the fascicles into the superficial and the deep aponeuroses) and pennation angles (angle between the fascicle and the deep aponeurosis) during the stance phase. Manual correction of the digitized fascicle and the deep aponeurosis, defined as a second fascicle, were performed where appropriate. If the field of view of the transducer was not sufficiently wide to capture the entire fascicle, the missing portion was estimated via manual extrapolation based on the assumption that the fascicle and the aponeurosis extend linearly. Ultrasonography has been frequently used in dynamic conditions (Cronin and Lichtwark, 2013) and has been demonstrated to provide reliable measures of GM fascicle lengths and pennation angles (Aggeloussis et al., 2010; Van Hooren et al., 2020).

SEE and MTU Lengths

Series elastic element length was estimated using an MTU model by subtracting muscle fascicle lengths multiplied by the cosine of their pennation angles from the MTU lengths (Fukunaga, 2001). Muscle–tendon unit length was calculated via a linear regression equation (Hawkins and Hull, 1990), using participant's shank length data (the distance from the lateral malleolus to the lateral femur epicondyle) in addition to recorded knee and ankle joint angles.

Data Processing

For each participant, and each outcome, eight consecutive stance phases (touchdown to toe-off) of the left foot per condition were analyzed and averaged using custom-made scripts (MATLAB R2018a, MathWorks, Inc., Natick, United States). Prior to being resampled to 101 data points per stance phase (to represent data as a percentage), fascicle lengths and pennation angles were smoothed with a five-point moving average. Electrogoniometer signals were smoothed with a fifth-order Butterworth low-pass filter, and a 10-Hz cut-off frequency. Fascicle velocity was calculated as the time derivative of its length using the central difference method (Robertson et al., 2013).

Based on the ultrasound and joint-angle recordings SEE length, MTU length, fascicle length, pennation angle, and fascicle velocity were determined at the time when peak SEE length was achieved, and thus force acting on the SEE is presumably at its greatest. Furthermore, average values across the stance phase were determined. Overall fascicle shortening was calculated by subtracting the minimum from maximum fascicle length. Knee and ankle joint range of motion were defined as the delta between their respective minimum and maximum joint angles. Additionally, the difference in knee and ankle joint angles between touchdown to the time of first local maximum and maximum dorsiflexion, were defined as knee flexion and ankle dorsiflexion, respectively. Knee and Ankle joint angles at touchdown and toe-off as well as ground-contact times were determined. To estimate the level of BWS achieved by applying LBPP, average plantar forces over the stance phase were determined and expressed as percentage of the average plantar forces when walking without BWS.

Statistical Analysis

Distribution normality was assessed using the Shapiro–Wilk normality test. If normally distributed, a two-tailed paired t-test was performed, whereas if not, a non-parametric Wilcoxon (matched-pairs) signed rank test was used to compare conditions (30 vs. 0% BWS). All tests were performed in GraphPad Prism (v 7.04) with a significance level of α = 0.05. Effect Sizes (d_z) were calculated using the G^*Power software version 3.1.9.4 (Faul et al., 2007). Thresholds of 0.2, 0.5 and 0.8 were defined as small, moderate and large effects between the two comparison groups (Cohen, 1988).

Results

Participants walking with 30% BWS generated significantly lower average plantar forces (−194 ± 32 N, P < 0.001, d_z = −6.07) corresponding to 68 ± 4% of the average plantar forces when walking without BWS which did not differ significantly from the target of 70% (P = 0.223). Ground-contact times were 0.03 ± 0.03 s longer when walking with 30% BWS, however, the effect was statistically not significant (P = 0.078, d_z = 0.80) (Table 1).

Figure 2 presents the averages and standard errors of joint angles and muscle–SEE outcomes time normalized to a single stance phase for participants walking with 0 and 30% BWS.

FIGURE 2

Figure 2. Sample average and standard error of plantar forces (A), knee joint angle (B), ankle joint angle (C), muscle-tendon unit length (D), series elastic element length (E), fascicle length (F), pennation angle (G), and fascicle velocity (H) for participants walking at 75% of their preferred walk-to-run transition speed with 0% body weight support (black line) and 30% body weight support (red line) during the entire stance phase. The solid lines represent the sample average, and the corresponding shaded areas represent the standard error of measurement.

No statistically significant differences in knee and ankle joint angles at touchdown (P = 0.164, d_z = −0.55; P = 0.635, d_z = −0.18), at toe-off (P = 0.848, d_z = −0.07; P = 0.641, d_z < 0.01) and at the time of the peak SEE length (P = 0.461, d_z = 0.42; P = 0.742, d_z = 0.04) were observed (Table 1). Furthermore, knee and ankle joint range of motion (P = 0.860, d_z = 0.06; P = 0.844, d_z = −0.04), knee flexion (P = 0.347, d_z = −0.36) and ankle dorsiflexion (P = 0.204, d_z = 0.50) were unaffected by walking with 30% BWS (Table 1).

Walking with 30% BWS had no effect on peak SEE length (P = 0.976, d_z = −0.01) (Figure 3A). Furthermore, at the time when peak SEE length was reached, no statistically significant differences from 0% BWS were observed for MTU length (P = 0.641, d_z= −0.04), fascicle length (P = 0.890, d_z = −0.05), pennation angle (P = 0.945, d_z = −0.03) and fascicle velocity (P = 0.576, d_z = −0.21) (Figures 3A–C).

FIGURE 3

Figure 3. Series elastic element length (A), muscle-tendon unit length (A), fascicle length (A), pennation angle (B), and fascicle velocity (C) at the time of the peak series elastic element length as presented as boxplots for participants walking without body weight support (black box) and 30% body weight support (red box). The lower and upper parts of the box represent the first and third quartile, respectively. The length of the whisker represent the minimum and maximum values. The horizontal line in the box represents the statistical median of the sample; +, the mean of the sample; °, individual data points.

No statistically significant differences of the average values across the entire stance phase were also observed for SEE length (P = 0.945, d_z = 0.05), MTU length (P = 0.641, d_z = 0.01), fascicle length (P = 0.790, d_z = −0.10), pennation angle (P = 0.641, d_z = 0.16) and fascicle velocity (P = 0.148, d_z = 0.51) between 0 and 30% BWS walking (Table 2). Furthermore, overall fascicle shortening did not differ between conditions (P = 0.313, d_z = −0.43) (Table 2).

TABLE 2

Table 2. Means and standard deviations of gastrocnemius medialis muscle and SEE outcome measures while participants walked at 75% of their preferred walk-to-run transition speed with 0 and 30% body weight support.

Discussion

Effects of Walking With 30% BWS on Contractile and Series Elastic Element Behavior

During the walking trials, participants were successfully unloaded by 30% of their body weight, as the average plantar forces actually achieved by inducing LBPP did not differ significantly from the target average plantar forces. The main findings were that walking with 30% BWS did not significantly affect joint kinematics. Furthermore, in contrast to the hypotheses, walking with 30% BWS induced no statistically significant differences from 0% BWS in peak SEE length as well as MTU length, fascicle length, and pennation angle neither at the time of the peak SEE length, nor on average across the stance phase. Also, in contrast with the hypotheses, no statistically significant effect of 30% BWS was found on fascicle shortening velocity at the time of the peak SEE length, nor on average across stance, despite a reduction in absolute walking speed (albeit same Froude number). These findings are further supported by the overall small effect sizes.

Previous studies and simulation models have shown that the GM force-length-velocity behavior shifts with gait type and speed to meet the varying locomotor demands (Farris and Sawicki, 2012; Arnold et al., 2013). However, in the present study, fascicle length and pennation angle were unchanged when walking with 30% BWS, which implies that the GM remains operating on a similar part of the force–length relationship, thereby preserving its force generation ability (Arnold et al., 2013). Moreover, GM fascicle velocity has been reported to decrease with decreasing walking speed, thereby increasing GM's force generation ability (Farris and Sawicki, 2012; Arnold et al., 2013). In fact, in the present study average fascicle velocity was 5.0 ± 10 mm·s⁻¹ slower when walking at a slightly slower speed at 30 vs. 0% BWS (−0.24 m·s⁻¹) reaching a moderate effect (d_z = 0.51), however, high variability may have contributed to it failing to reach statistical significance.

It has been reported that the ankle plantarflexion moment decreases with increasing BWS (Lewek, 2011; Goldberg and Stanhope, 2013; Fischer and Wolf, 2015). In fact, the present results suggests a reduction in average plantar force by almost 200 N whilst ankle joint kinematics were largely preserved when walking with 30% BWS suggesting that ankle joint moment was reduced. Interestingly, this did not affect peak SEE length, which incorporates the length of the free tendon and aponeurosis. As aponeurosis stiffness varies upon contractile conditions (e.g., reduced muscle activity results in lower orthogonal muscle expansion linked to lower transverse strain) (Azizi and Roberts, 2009), SEE length can remain similar despite a reduction in ankle joint moment.

However, as in the present observational study MTU interaction was only modeled for the GM, which accounts for a modest fraction (~17%) of the physiologic cross-sectional area of the plantar flexor muscles (Ward et al., 2009), changes that influence the ankle joint moment might not be fully reflected. Furthermore, joint moments were not determined and SEE length was not measured directly but estimated using an MTU model. Thus, if tuning of the mechanical properties of the SEE actually causes preserved SEE and fascicle kinematics warrants further study.

Implications for Rehabilitative Body Weight Supported Gait Training

Maintenance of joint kinematics and GM behavior may facilitate rehabilitative gait training by preserving “natural” movement patterns, despite joint loads and related pain being reduced (Eastlack et al., 2005; Cutuk et al., 2006). Preserved fascicle's operating range suggests that the stimuli exerted on the muscle remain the same and thus help to maintain optimum fascicle length for force production, which is key for locomotor recovery. Furthermore, the maintenance of SEE strain, as possibly achieved by an increased aponeurosis strain, might help to prevent degeneration and maintain function of the aponeurosis despite external unloading. Patients who may benefit from LBPP gait training during their early postoperative rehabilitation include not only those with tendon, ligament and meniscus repairs but also joint replacements or fractures (Eastlack et al., 2005). However, the increased aponeurosis strain, which is required to compensate for the decreased free tendon strain (and thus to maintain SEE strain), could pose a potential risk to patients after Achilles tendon rupture if the rupture does not exclusively affect the free tendon. Therefore, BWS rehabilitation should be individualized to the specific pathological characteristics of patients, depending on the impaired biological tissues that require unloading, e.g., rehabilitation after total knee arthroplasty vs. ankle or Achilles tendon injury. Based on the current findings, further studies including different patient groups are required.

The present data, are not only in agreement with a recent systematic review, which concluded that spatio-temporal and kinematic gait parameters can be preserved with up to 30% BWS (Apte et al., 2018), but extends this view to preserved muscle–SEE mechanics. In fact, healthy individuals appear able to retain normal walking kinematics even when unloaded by up 50% BWS (Van Hedel et al., 2006; Awai et al., 2017). The absence of any effects when BWS was increased from 0 to 30% suggests that the modulation of fascicle–SEE behavior does not develop linearly with increasing BWS but is determined by a certain threshold, however if this threshold is below or above 50% BWS remains to be determined. Additionally, if non-LBPP BWS systems, such as overhead suspension harnesses, therapist-assisted waist belts or robotic-assisted gait-training devices, are also able to preserve GM behavior warrants further study. Nevertheless, the present observational study supports the recommendation (Fischer and Wolf, 2015) for LBPP-induced 30% BWS in rehabilitative gait training. Finally, it should be noted that walking speed was intentionally reduced with increasing BWS via the adjustment to the same Froude number to obtain mechanically equivalent walking speeds (Vaughan and O'Malley, 2005). Thus, the observation that the neural system appears to largely preserve GM overall contraction behavior in addition to joint kinematics suggests that the approach of producing comparable gait patterns across the different walking conditions was successful and should be considered for future gait rehabilitation.

Conclusions

This is the first study to examine in vivo GM fascicle–SEE behavior during walking at 30% BWS, frequently employed in gait rehabilitation, at 75% PTS on an LBPP treadmill. The present findings reveal that during walking with 30% BWS fascicle–SEE behavior was largely preserved, in contrast to the hypothesis. Thus, the present study not only supports the contention made in previous studies that walking with the recommended therapeutic dose of 30% BWS largely retains spatio-temporal and joint kinematic characteristics but extends this to GM fascicle and SEE mechanics. This may be advantageous during rehabilitative gait training with BWS as it indicates transferability of gait patterns to subsequent unsupported walking.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The study involving human participants were reviewed and approved by Ärztekammer Nordrhein Ethical Committee of Düsseldorf, Germany. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

BB, TW, and DG conceptualized research. BB, JA, TW, KM, JR, DG, and KA designed research. CR, BB, BS, JA, and AS acquired data. CR, BB, BS, AS, and KA analyzed data. CR, BB, JR, DG, and KA interpreted data. CR and DG drafted manuscript. BB, JA, TW, KM, JR, DG, and KA revised manuscript. All authors approved manuscript.

Funding

This study was funded by the European Space Agency's (ESA) Space Medicine Team (HRE-OM) of the European Astronaut Center (EAC) in Cologne, Germany and the University of Applied Sciences Aachen obtained funding by the Federal Ministry for Economic Affairs and Energy (50WB1728). Spaceship EAC provided the lower body positive pressure treadmill (AlterG) used in this study.

Conflict of Interest

DG and TW are employed by KBR on behalf of the European Space Agency.

The remaining 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.

Abbreviations

AlterG, AlterG® Anti-Gravity TreadmillTM; BWS, body weight support; g, earth's gravitational acceleration; GM, gastrocnemius medialis; LBPP, lower-body positive pressure; MTU, muscle–tendon unit; PTS, preferred walk-to-run transition speed; SEE, series elastic element.

References

Aggeloussis, N., Giannakou, E., Albracht, K., and Arampatzis, A. (2010). Reproducibility of fascicle length and pennation angle of gastrocnemius medialis in human gait in vivo. Gait Posture 31, 73–77. doi: 10.1016/j.gaitpost.2009.08.249

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, F. C., and Pandy, M. G. (2003). Individual muscle contributions to support in normal walking. Gait Posture 17, 159–169. doi: 10.1016/S0966-6362(02)00073-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Apte, S., Plooij, M., and Vallery, H. (2018). Influence of body weight unloading on human gait characteristics: a systematic review. J. Neuroeng. Rehabil. 15:53. doi: 10.1186/s12984-018-0414-7

CrossRef Full Text | Google Scholar

Arnold, E. M., Hamner, S. R., Seth, A., Millard, M., and Delp, S. L. (2013). How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds. J. Exp. Biol. 216, 2150–2160. doi: 10.1242/jeb.075697

PubMed Abstract | CrossRef Full Text | Google Scholar

Awai, L., Franz, M., Easthope, C., Vallery, H., Curt, A., and Bolliger, M. (2017). Preserved gait kinematics during controlled body unloading. J. Neuroeng. Rehabil. 14:25. doi: 10.1186/s12984-017-0239-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Azizi, E., and Roberts, T. J. (2009). Biaxial strain and variable stiffness in aponeuroses. J. Physiol. 587, 4309–4318. doi: 10.1113/jphysiol.2009.173690

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavagna, G., Willems, P., and Heglund, N. (2000). The role of gravity in human walking: pendular energy exchange, external work and optimal speed. J. Physiol. 528, 657–668. doi: 10.1111/j.1469-7793.2000.00657.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences. New York, NY: Academic.

Google Scholar

Cronin, N. J., and Lichtwark, G. (2013). The use of ultrasound to study muscle-tendon function in human posture and locomotion. Gait Posture 37, 305–312. doi: 10.1016/j.gaitpost.2012.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Cutuk, A., Groppo, E. R., Quigley, E. J., White, K. W., Pedowitz, R. A., and Hargens, A. R. (2006). Ambulation in simulated fractional gravity using lower body positive pressure: cardiovascular safety and gait analyses. J. Appl. Physiol. 101, 771–777. doi: 10.1152/japplphysiol.00644.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Eastlack, R. K., Hargens, A. R., Groppo, E. R., Steinbach, G. C., White, K. K., and Pedowitz, R. A. (2005). Lower body positive-pressure exercise after knee surgery. Clin. Orthop. Relat. Res. 431, 213–219. doi: 10.1097/01.blo.0000150459.92012.f7

PubMed Abstract | CrossRef Full Text | Google Scholar

Farley, C. T., and McMahon, T. A. (1992). Energetics of walking and running: insights from simulated reduced-gravity experiments. J. Appl. Physiol. 73, 2709–2712. doi: 10.1152/jappl.1992.73.6.2709

PubMed Abstract | CrossRef Full Text | Google Scholar

Farris, D. J., and Lichtwark, G. A. (2016). UltraTrack: software for semi-automated tracking of muscle fascicles in sequences of B-mode ultrasound images. Comput. Methods Programs Biomed. 128, 111–118. doi: 10.1016/j.cmpb.2016.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Farris, D. J., and Sawicki, G. S. (2012). Human medial gastrocnemius force-velocity behavior shifts with locomotion speed and gait. Proc. Natl. Acad. Sci. U.S.A. 109, 977–982. doi: 10.1073/pnas.1107972109

PubMed Abstract | CrossRef Full Text | Google Scholar

Faul, F., Erdfelder, E., Lang, A.-G., and Buchner, A. (2007). G^* Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods 39, 175–191. doi: 10.3758/BF03193146

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischer, A. G., Debbi, E. M., and Wolf, A. (2015). Effects of body weight unloading on electromyographic activity during overground walking. J. Electromyogr. Kinesiol. 25, 709–714. doi: 10.1016/j.jelekin.2015.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischer, A. G., and Wolf, A. (2015). Assessment of the effects of body weight unloading on overground gait biomechanical parameters. Clin. Biomech. 30, 454–461. doi: 10.1016/j.clinbiomech.2015.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukunaga, T. (2001). In vivo behaviour of human muscle tendon during walking. Proc. R. Soc. Land. B 268, 229–233. doi: 10.1098/rspb.2000.1361

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, S. R., and Stanhope, S. J. (2013). Sensitivity of joint moments to changes in walking speed and body-weight-support are interdependent and vary across joints. J. Biomech. 46, 1176–1183. doi: 10.1016/j.jbiomech.2013.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowski, A., Farley, C. T., and Kram, R. (2005). Independent metabolic costs of supporting body weight and accelerating body mass during walking. J. Appl. Physiol. 98, 579–583. doi: 10.1152/japplphysiol.00734.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowski, A. M. (2010). Metabolic and biomechanical effects of velocity and weight support using a lower-body positive pressure device during walking. Arch. Phys. Med. Rehabil. 91, 951–957. doi: 10.1016/j.apmr.2010.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, C., Einarson, E., Thomson, A., Whiteley, R., and Witvrouw, E. (2017). Hamstring and calf muscle activation as a function of bodyweight support during treadmill running in ACL reconstructed athletes. Gait Posture 58, 154–158. doi: 10.1016/j.gaitpost.2017.07.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkins, D., and Hull, M. (1990). A method for determining lower extremity muscle-tendon lengths during flexion/extension movements. J. Biomech. 23, 487–494. doi: 10.1016/0021-9290(90)90304-L

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawkins, D., Lum, C., Gaydos, D., and Dunning, R. (2009). Dynamic creep and pre-conditioning of the Achilles tendon in-vivo. J. Biomech. 42, 2813–2817. doi: 10.1016/j.jbiomech.2009.08.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Hesse, S. (2008). Treadmill training with partial body weight support after stroke: a review. NeuroRehabilitation 23, 55–65. doi: 10.3233/NRE-2008-23106

PubMed Abstract | CrossRef Full Text | Google Scholar

Karamanidis, K., Arampatzis, A., and Brüggemann, G.-P. (2003). Symmetry and reproducibility of kinematic parameters during various running techniques. Med. Sci. Sports Exerc. 35, 1009–1016. doi: 10.1249/01.MSS.0000069337.49567.F0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kram, R., Domingo, A., and Ferris, D. P. (1997). Effect of reduced gravity on the preferred walk-run transition speed. J. Exp. Biol. 200, 821–826.

PubMed Abstract | Google Scholar

Labini, F. S., Ivanenko, Y. P., Cappellini, G., Gravano, S., and Lacquaniti, F. (2011). Smooth changes in the EMG patterns during gait transitions under body weight unloading. J. Neurophysiol. 106, 1525–1536. doi: 10.1152/jn.00160.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, A., Lichtwark, G. A., Schache, A. G., Lin, Y. C., Brown, N. A., and Pandy, M. G. (2015). In vivo behavior of the human soleus muscle with increasing walking and running speeds. J. Appl. Physiol. 118, 1266–1275. doi: 10.1152/japplphysiol.00128.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewek, M. D. (2011). The influence of body weight support on ankle mechanics during treadmill walking. J. Biomech. 44, 128–133. doi: 10.1016/j.jbiomech.2010.08.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Lichtwark, G., Bougoulias, K., and Wilson, A. (2007). Muscle fascicle and series elastic element length changes along the length of the human gastrocnemius during walking and running. J. Biomech. 40, 157–164. doi: 10.1016/j.jbiomech.2005.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Liebenberg, J., Scharf, J., Forrest, D., Dufek, J. S., Masumoto, K., and Mercer, J. A. (2011). Determination of muscle activity during running at reduced body weight. J. Sports Sci. 29, 207–214. doi: 10.1080/02640414.2010.534806

PubMed Abstract | CrossRef Full Text | Google Scholar

McGowan, C. P., Neptune, R. R., and Kram, R. (2008). Independent effects of weight and mass on plantar flexor activity during walking: implications for their contributions to body support and forward propulsion. J. Appl. Physiol. 105, 486–494. doi: 10.1152/japplphysiol.90448.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Minetti, A. E. (2001). Invariant aspects of human locomotion in different gravitational environments. Acta Astronaut 49, 191–198. doi: 10.1016/S0094-5765(01)00098-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mishra, D. K. (2015). Summary of Clinical Research Involving The AlterG Anti-Gravity Treadmill^®. Available online at: https://static1.squarespace.com/static/5a77c1f76f4ca3157898f698/t/5a78edadf9619a2745217400/1517874616880/Summary+of+Clinical+Research+Involving+the+AlterG.pdf (accessed November 20, 2020).

Google Scholar

Neal, M., Fleming, N., Eberman, L., Games, K., and Vaughan, J. (2016). Effect of body-weight-support running on lower-limb biomechanics. J. Orthop. Sports Phys. Ther. 46, 784–793. doi: 10.2519/jospt.2016.6503

PubMed Abstract | CrossRef Full Text | Google Scholar

Neptune, R., Kautz, S., and Zajac, F. (2001). Contributions of the individual ankle plantar flexors to support, forward progression and swing initiation during walking. J. Biomech. 34, 1387–1398. doi: 10.1016/S0021-9290(01)00105-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Neptune, R. R., and Sasaki, K. (2005). Ankle plantar flexor force production is an important determinant of the preferred walk-to-run transition speed. J. Exp. Biol. 208, 799–808. doi: 10.1242/jeb.01435

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, S., Steklov, N., Bugbee, W. D., Goldberg, T., Colwell, C. W., and D'lima, D. D. (2013). Anti-gravity treadmills are effective in reducing knee forces. J. Orthop. Res. 31, 672–679. doi: 10.1002/jor.22272

PubMed Abstract | CrossRef Full Text | Google Scholar

Pavei, G., Biancardi, C. M., and Minetti, A. E. (2015). Skipping vs. running as the bipedal gait of choice in hypogravity. J. Appl. Physiol. 119, 93–100. doi: 10.1152/japplphysiol.01021.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Quigley, E., Noh, H., Groppo, E., Cutuk, A., Pedowitz, R., and Hargens, A. (2000). Gait mechanics using a lower body positive pressure chamber for orthopaedic rehabilitation. Trans. Orthop. Res. Soc. 25, 828. Available online at: http://www.ors.org/Transactions/46/0828.pdf

Google Scholar

Richter, C., Braunstein, B., Winnard, A., Nasser, M., and Weber, T. (2017). Human biomechanical and cardiopulmonary responses to partial gravity–a systematic review. Front. Physiol. 8:583. doi: 10.3389/fphys.2017.00583

PubMed Abstract | CrossRef Full Text | Google Scholar

Robertson, G. E., Caldwell, G. E., Hamill, J., Kamen, G., and Whittlesey, S. (2013). Research Methods in Biomechanics. Champaign, IL: Human Kinetics. doi: 10.5040/9781492595809

CrossRef Full Text | Google Scholar

Van Hedel, H., Tomatis, L., and Müller, R. (2006). Modulation of leg muscle activity and gait kinematics by walking speed and bodyweight unloading. Gait Posture 24, 35–45. doi: 10.1016/j.gaitpost.2005.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Hooren, B., Teratsias, P., and Hodson-Tole, E. F. (2020). Ultrasound imaging to assess skeletal muscle architecture during movements: a systematic review of methods, reliability, and challenges. J. Appl. Physiol. 128, 978–999. doi: 10.1152/japplphysiol.00835.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaughan, C. L., and O'Malley, M. J. (2005). Froude and the contribution of naval architecture to our understanding of bipedal locomotion. Gait Posture 21, 350–362. doi: 10.1016/j.gaitpost.2004.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, S. R., Eng, C. M., Smallwood, L. H., and Lieber, R. L. (2009). Are current measurements of lower extremity muscle architecture accurate? Clin. Orthop. Relat. Res. 467, 1074–1082. doi: 10.1007/s11999-008-0594-8

PubMed Abstract | CrossRef Full Text | Google Scholar