Lower complexity of motor primitives ensures robust control of high-speed human locomotion (Q11977)

Walking and running are mechanically and energetically different locomotion modes. For selecting one or another, speed is a parameter of paramount importance. Yet, both are likely controlled by similar low-dimensional neuronal networks that reflect in patterned muscle activations called muscle synergies. Here, we investigated how humans synergistically activate muscles during locomotion at different submaximal and maximal speeds. We analysed the duration and complexity (or irregularity) over time of motor primitives, the temporal components of muscle synergies. We found that the challenge imposed by controlling high-speed locomotion forces the central nervous system to produce muscle activation patterns that are wider and less complex relative to the duration of the gait cycle. The motor modules, or time-independent coefficients, were redistributed as locomotion speed changed. These outcomes show that robust locomotion control at challenging speeds is achieved by modulating the relative contribution of muscle activations and producing less complex and wider control signals, whereas slow speeds allow for more irregular control. In this supplementary data set we made available: a) the metadata with anonymized participant information, b) the raw EMG, c) the touchdown and lift-off timings of the recorded limb, d) the filtered and time-normalized EMG, e) the muscle synergies extracted via NMF and f) the code to process the data, including the scripts to calculate the Higuchis fractal dimension (HFD) of motor primitives. In total, 180 trials from 30 participants are included in the supplementary data set. The file metadata.dat is available in ASCII and RData format and contains: Code: the participants code Group: the experimental group in which the participant was involved (G1 = walking and submaximal running; G2 = submaximal and maximal running) Sex: the participants sex (M or F) Speeds: the type of locomotion (W for walking or R for running) and speed at which the recordings were conducted in 10*[m/s] Age: the participants age in years Height: the participants height in [cm] Mass: the participants body mass in [kg] PB: 100 m-personal best time (for G2). The RAW_DATA.RDataR list consists of elements of S3 class EMG, each of which is a human locomotion trial containing cycle segmentation timings and raw electromyographic (EMG) data from 13 muscles of the right-side leg. Cycle times are structured as data frames containing two columns thatcorrespond to touchdown (first column) and lift-off (second column).Raw EMG data sets are also structured as data frames with one row for each recorded data pointand 14 columns. The first column contains the incremental time in seconds. The remaining 13 columns contain the raw EMG data, named with the following muscle abbreviations:ME = gluteus medius, MA = gluteus maximus, FL = tensor fasci lat, RF = rectus femoris, VM = vastus medialis, VL = vastus lateralis, ST = semitendinosus, BF = biceps femoris, TA = tibialis anterior, PL = peroneus longus, GM = gastrocnemius medialis, GL = gastrocnemius lateralis, SO = soleus. Please note that the following trials include less than 30 gait cycles (the actual number shown between parentheses): P16_R_83 (20), P16_R_95 (25), P17_R_28 (28), P17_R_83 (24), P17_R_95 (13), P18_R_95 (23), P19_R_95 (18), P20_R_28 (25), P20_R_42 (27), P20_R_95 (25), P22_R_28 (23), P23_R_28(29), P24_R_28 (28), P24_R_42 (29), P25_R_28 (29), P25_R_95 (28), P26_R_28 (29), P26_R_95 (28), P27_R_28 (28), P27_R_42 (29), P27_R_95 (24), P28_R_28 (29), P29_R_95 (17).All the other trialsconsist of 30 gait cycles.Trials are named like P20_R_20, where the characters P20 indicate the participant number (in this example the 20th), the character R indicate the locomotion type (W=walking, R=running), and the numbers 20 indicate the locomotion speed in 10*m/s (in this case the speed is 2.0 m/s). The filtered and time-normalized emg data is named, following the same rules, like FILT_EMG_P03_R_30. Old versions not compatible with the R package musclesyneRgies The files containing the gait cycle breakdown are available in RData format, in the file named CYCLE_TIMES.RData. The files are structured as data frames with as many rows as the available number of gait cycles and two columns. The first column named touchdown contains the touchdown incremental times in seconds. The second column named stance contains the duration of each stance phase of the right foot in seconds. Each trial is saved as an element of a single R list. Trials are named like CYCLE_TIMES_P20_R_20, where the characters CYCLE_TIMES indicate that the trial contains the gait cycle breakdown times, the characters P20 indicate the participant number (in this example the 20th), the character R indicate the locomotion type (W=walking, R=running), and the numbers 20 indicate the locomotion speed in 10*m/s (in this case the speed is 2.0 m/s). Please note that the following trials include less than 30 gait cycles (the actual number shown between parentheses): P16_R_83 (20), P16_R_95 (25), P17_R_28 (28), P17_R_83 (24), P17_R_95 (13), P18_R_95 (23), P19_R_95 (18), P20_R_28 (25), P20_R_42 (27), P20_R_95 (25), P22_R_28 (23), P23_R_28(29), P24_R_28 (28), P24_R_42 (29), P25_R_28 (29), P25_R_95 (28), P26_R_28 (29), P26_R_95 (28), P27_R_28 (28), P27_R_42 (29), P27_R_95 (24), P28_R_28 (29), P29_R_95 (17). The files containing the raw, filtered and the normalized EMG data are available in RData format, in the files named RAW_EMG.RData and FILT_EMG.RData. The raw EMG files are structured as data frames with as many rows as the amount of recorded data points and 13 columns. The first column named time contains the incremental time in seconds. The remaining 12 columns contain the raw EMG data, named with muscle abbreviations that follow those reported above. Each trial is saved as an element of a single R list. Trials are named like RAW_EMG_P03_R_30, where the characters RAW_EMG indicate that the trial contains raw emg data, the characters P03 indicate the participant number (in this example the 3rd), the character R indicate the locomotion type (see above), and the numbers 30 indicate the locomotion speed (see above). The filtered and time-normalized emg data is named, following the same rules, like FILT_EMG_P03_R_30. The files containing the muscle synergies extracted from the filtered and normalized EMG data are available in RData format, in the files named SYNS_H.RData and SYNS_W.RData. The muscle synergies files are divided in motor primitives and motor modules and are presented as direct output of the factorisation and not in any functional order. Motor primitives are data frames with 6000 rows and a number of columns equal to the number of synergies (which might differ from trial to trial) plus one. The rows contain the time-dependent coefficients (motor primitives), one column for each synergy plus the time points (columns are named e.g. time, Syn1, Syn2, Syn3, where Syn is the abbreviation for synergy). Each gait cycle contains 200 data points, 100 for the stance and 100 for the swing phase which, multiplied by the 30 recorded cycles, result in 6000 data points distributed in as many rows. This output is transposed as compared to the one discussed in the methods section to improve user readability. Each set of motor primitives is saved as an element of a single R list. Trials are named like SYNS_H_P12_W_07, where the characters SYNS_H indicate that the trial contains motor primitive data, the characters P12 indicate the participant number (in this example the 12th), the character W indicate the locomotion type (see above), and the numbers 07 indicate the speed (see above). Motor modules are data frames with 12 rows (number of recorded muscles) and a number of columns equal to the number of synergies (which might differ from trial to trial). The rows, named with muscle abbreviations that follow those reported above, contain the time-independent coefficients (motor modules), one for each synergy and for each muscle. Each set of motor modules relative to one synergy is saved as an element of a single R list. Trials are named like SYNS_W_P22_R_20, where the characters SYNS_W indicate that the trial contains motor module data, the characters P22 indicate the participant number (in this example the 22nd), the character W indicates the locomotion type (see above), and the numbers 20 indicate the speed (see above). Given the nature of the NMF algorithm for the extraction of muscle synergies, the supplementary data set might show non-significant differences as compared to the one used for obtaining the results of this paper. The files containing the HFD calculated from motor primitives are available in RData format, in the file named HFD.RData. HFD results are presented in a list of lists containing, for each trial, 1) the HFD, and 2) the interval time k used for the calculations. HFDs are presented as one number (mean HFD of the primitives for that trial), as are the interval times k. Trials are named like HFD_P01_R_95, where the characters HFD indicate that the trial contains HFD data, the characters P01 indicate the participant number (in this example the 1st), the character R indicates the locomotion type (see above), and the numbers 95 indicate the speed (see above). All the code used for the pre-processing of EMG data, the extraction of muscle synergies and the calculation of HFD is available in R format. Explanatory comments are profusely present throughout the script muscle_synergies.R.