Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gait Patterns in Patients with Hereditary Spastic Paraparesis

  • Mariano Serrao ,

    Contributed equally to this work with: Mariano Serrao, Martina Rinaldi

    mariano.serrao@uniroma1.it

    Affiliations Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza, Latina, Italy, Rehabilitation Centre, Policlinico Italia, Rome, Italy

  • Martina Rinaldi ,

    Contributed equally to this work with: Mariano Serrao, Martina Rinaldi

    Affiliation Department of Engineering, Roma TRE University, Rome, Italy

  • Alberto Ranavolo,

    Affiliation Department of Occupational and Environmental Medicine, Epidemiology and Hygiene, INAIL, Monte Porzio Catone, Rome, Italy

  • Francesco Lacquaniti,

    Affiliations Centre of Space Bio-Medicine, University of Rome Tor Vergata, Rome, Italy, Laboratory of Neuromotor Physiology, Istituto Di Ricovero e Cura a Carattere Scientifico Santa Lucia Foundation, Rome, Italy, Department of Systems Medicine, University of Rome Tor Vergata, Rome, Italy

  • Giovanni Martino,

    Affiliations Centre of Space Bio-Medicine, University of Rome Tor Vergata, Rome, Italy, Laboratory of Neuromotor Physiology, Istituto Di Ricovero e Cura a Carattere Scientifico Santa Lucia Foundation, Rome, Italy

  • Luca Leonardi,

    Affiliation Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza, Latina, Italy

  • Carmela Conte,

    Affiliation Fondazione Don Gnocchi, Milan, Italy

  • Tiwana Varrecchia,

    Affiliation Department of Engineering, Roma TRE University, Rome, Italy

  • Francesco Draicchio,

    Affiliation Department of Occupational and Environmental Medicine, Epidemiology and Hygiene, INAIL, Monte Porzio Catone, Rome, Italy

  • Gianluca Coppola,

    Affiliation G.B. Bietti Foundation-IRCCS, Department of Neurophysiology of Vision and Neurophthalmology, Rome, Italy

  • Carlo Casali,

    Affiliation Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza, Latina, Italy

  • Francesco Pierelli

    Affiliations Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza, Latina, Italy, IRCCS, Neuromed, Pozzilli, Isernia, Italy

Abstract

Background

Spastic gait is a key feature in patients with hereditary spastic paraparesis, but the gait characterization and the relationship between the gait impairment and clinical characteristics have not been investigated.

Objectives

To describe the gait patterns in hereditary spastic paraparesis and to identify subgroups of patients according to specific kinematic features of walking.

Methods

We evaluated fifty patients by computerized gait analysis and compared them to healthy participants. We computed time-distance parameters of walking and the range of angular motion at hip, knee, and ankle joints, and at the trunk and pelvis. Lower limb joint moments and muscle co-activation values were also evaluated.

Results

We identified three distinct subgroups of patients based on the range of motion values. Subgroup one was characterized by reduced hip, knee, and ankle joint range of motion. These patients were the most severely affected from a clinical standpoint, had the highest spasticity, and walked at the slowest speed. Subgroup three was characterized by an increased hip joint range of motion, but knee and ankle joint range of motion values close to control values. These patients were the most mildly affected and had the highest walking speed. Finally, subgroup two showed reduced knee and ankle joint range of motion, and hip range of motion values close to control values. Disease severity and gait speed in subgroup two were between those of subgroups one and three.

Conclusions

We identified three distinctive gait patterns in patients with hereditary spastic paraparesis that correlated robustly with clinical data. Distinguishing specific features in the gait patterns of these patients may help tailor pharmacological and rehabilitative treatments and may help evaluate therapeutic effects over time.

Introduction

Hereditary spastic paraparesis is a heterogeneous group of inherited neurodegenerative disorders characterized by retrograde degeneration of the corticospinal axonal fibers [1]. Lower limb spasticity, usually more prominent than muscle weakness, is the key clinical feature in patients with hereditary spastic paraparesis [2] and impairs walking ability, autonomy, and quality of life [3,4]. No treatment is known to reduce disease progression, but antispastic drugs and physiotherapy [58] may help reduce the functional impairment of gait. Quantifying and typifying the specific gait disorder in hereditary spastic paraparesis is crucial to designing individual pharmacological and rehabilitative treatments. Most descriptions of paraparetic gait are based on qualitative clinical observations [1,2,5], [917]. Some studies have quantitatively evaluated gait impairment in hereditary spastic paraparesis patients, revealing several gait abnormalities of reduced step length, increased step width, reduced range of motion (RoM) at the knee joint [1820], impaired knee torque and stiffness [19,20], and decreased activity of the rectus femoris muscle [19]. Despite the great relevance of such quantitative assessments, they remain generic without reflecting the wide clinical heterogeneity of gait disorders in hereditary spastic paraparesis patients. Spasticity of the lower limb muscles represents the most important clinical sign of hereditary spastic paraparesis, but it affects different patients to different extents [5,18,20]. Individual differences in spasticity should translate into corresponding biomechanical features of gait; specifically, more spastic patients should have more reduced RoMs during walking [18,20]. We hypothesized that the individual kinematic behavior of patients with hereditary spastic paraparesis could be used to identify distinct subgroups of patients, and that these subgroups would exhibit different levels of limb spasticity. Our aims were as follows: i) to perform a comprehensive analysis of kinematics, kinetics and sEMG (surface electromyography) in adult patients with hereditary spastic paraparesis, and ii) to identify specific gait patterns in subgroups of patients categorized according to their kinematic behavior.

Materials and Methods

Subjects

We recruited fifty patients with hereditary spastic paraparesis (twenty women and thirty men, mean age 47.70 ± 16.06 years, height 1.64 ± 0.11 m, weight: 75.97 ± 18.51 kg, disease duration 17.65 ± 12.50 years). All patients included in the study were able to walk without assistance or walking aids on a level surface. A defined molecular diagnosis of hereditary spastic paraparesis was applied to thirty patients. Of these, twenty-two patients had spastic paraplegia (SPG) type four (mutations in SPAST), two patients had SPG3A (mutations in ATL1), one patient had SPG5 (mutations in CYP7B1), two patients had SPG7 (mutations in the PGN), and three patients had SPG31 (mutations in REEP1). Twenty patients did not have a molecular diagnosis at the time of examination, but all patients unequivocally showed either a recessive (eight patients) or dominant (twelve patients) inheritance pattern. None of the patients showed any involvement of neurological systems other than the pyramidal one (e.g. cerebellar or sensory deficits). All patients were evaluated independently by two experienced neurologists (C.C. and F.P.) who assessed cognitive functions, cranial nerves, muscle tone, muscle strength, joint coordination, tendon reflexes, and sensory function.

The severity of the disease was rated using the Spastic Paraplegia Rating Scale (SPRS). The spasticity of hip and knee joint muscles was scored by the Modified Ashworth scale included in SPRS as a spasticity-related subscale [21]. Table 1 summarizes the clinical features and genotypes of all patients. Twelve out of fifty patients were assuming oral antispastic drugs (baclofen or tizanidine) since 4–6 years, All patients were clinically stable at the time of the study evaluation. Indeed, their clinical assessment (SPRS) did not change over the last six months prior to the study. At the time of the evaluation, all patients were undergoing physical therapy, which included lower limb and stretching exercises, balance, and gait training.

The control group was fifty healthy subjects (twenty-three women and twenty-seven men, mean age 49.12 ± 11.76 years, height 1.68 ± 0.07 m, weight 70.83 ± 13.22 kg).

All participants provided written informed consent before taking part in the study, which complied with the Helsinki Declaration and had local ethics committee approval (ICOT- Sapienza, Polo Pontino).

Gait analysis

Kinematic data were recorded at 300 Hz using an optoelectronic motion analysis system (SMART-D System, BTS, Milan, Italy) consisting of eight infrared cameras spaced around the walkway. In accordance with a validated biomechanical model, twenty-two reflective spherical markers (15 mm in diameter) were attached on the anatomical landmarks in accordance with a validated biomechanical model [22], using double-adhesive tape in such a way as to prevent them from falling out of place during the test. In detail, the markers were placed over the cutaneous projections of the spinous processes of the seventh cervical vertebra and sacrum and bilaterally over acromion, anterior superior iliac spine, great trochanter, lateral femoral condyle, fibula head, lateral malleoli and metatarsal head. In addition to markers directly applied to the skin, sticks, or wand, varying in length from 7 to 10 cm, placed at 1/3 of the length of the body segment (femur and leg) were used. Anthropometric data were collected for each subject [23].

Ground reaction forces were acquired by two dynamometric platforms (Kistler 9286B, Winterthur, Switzerland), attached to each other in the longitudinal direction but displaced by 0.2 m in the lateral direction (sampling rate 1200 Hz).

Surface myoelectric signals were recorded at 1000 Hz using a 16-channel wireless system (FreeEMG300 System, BTS, Milan, Italy). After skin preparation, bipolar Ag/AgCl surface electrodes (H124SG Kendall ARBO, Donau, Germany) were placed over the muscle belly in the direction of the muscle fibers according to the European Recommendations for Surface Electromyography [24] and the atlas of muscle innervation zones [25]. Bipolar electrodes, eight in total, were placed on the right side of the body of each subject on the tibialis anterior (TA), gastrocnemius lateralis (LG), gastrocnemius medialis (MG), vastus lateralis (VL), vastus medialis (VM), rectus femoris (RF), biceps femoris (BF), and semitendinosus (ST). Acquisition of kinematic, kinetic, and electromyographic data was integrated and synchronized.

Experimental Procedure

Patients and controls were asked to walk barefoot at a comfortable, self-selected speed along a walkway approximately 10 m in length while looking forward. Because we were interested in natural locomotion, only general, qualitative instructions were provided. Before the recording session, subjects practiced for a few minutes to familiarize themselves with the procedure. The starting position was adjusted to ensure that the right foot always landed at least on one of the two force platforms embedded in the middle of the pathway. Given that typical walking speeds were slow in these patients, we instructed the control subjects to also walk at a low but comfortable speed. In this way, the general characteristics of gait could be compared between the groups without any potential bias due to speed differences (see below speed matching procedure).

Ten trials per patient were recorded, instead healthy subjects were evaluated for a total of 15 trials (10 trials self selected speed and 5 trials slow walking). To avoid muscle fatigue, blocks of three trials were separated by a one-minute rest period.

Speed matching procedure

Walking speed was matched between groups as follows: we considered only those control group subjects whose mean walking speed fell within the range identified by patients’ mean walking speed ± SD [26]. Unpaired two-sample t-test was used to investigate differences in walking velocity between patients and controls. In this way, the mean speed values were not statistically different between groups (patients 2.40 ± 1.29 km/h; controls 2.63±0.71 km/h, p = 0.283).

Data Analysis

Kinematic, kinetic and electromyographic data were normalized to the duration of the gait cycle and interpolated to 201 samples using a polynomial procedure. Gait cycle was defined as the time between two successive foot contacts of the same leg. In this study, heel strike and toe-off events were determined by maximum and minimum of limb angle excursions. Limb angle was calculated as the angle between a vertical axis from the greater trochanter and a vector drawn from the greater trochanter to lateral malleolus projected on the sagittal plane: a 0° limb angle means that the leg was positioned vertically under the body; positive angles denote flexion (i.e. limb positioned in front of the vertical axis) and negative angles denote extension (i.e. limb positioned behind the vertical axis) [2729]. When subjects stepped on the force platforms, these kinematic criteria were verified by comparison with foot strike and lift-off measured from a threshold crossing event in the vertical force: stance phase was defined as the interval during which the vertical reaction force exceeded 7% of body weight. In general, the difference between the time events measured from kinematics and kinetics was no more than 3% [29] and kinematic criterion proved to be very robust in both healthy subejcts [27] and neurological patients [28]. The raw sEMG signals were band-pass filtered using a zero-lag third-order Butterworth filter (20–450 Hz), rectified, and low-pass filtered with a zero-lag fourth-order Butterworth filter (10 Hz). For each individual, the sEMG signal from each muscle was normalized to its peak value across all trials [29].

Time-distance, kinematic, kinetic and sEMG parameters were evaluated after preprocessing procedures.

Time-distance parameters.

The following time-distance gait parameters were calculated for each subject: walking speed (km/h), stance duration (% gait cycle), swing duration (% gait cycle), first and second double support duration (% gait cycle), cadence (step/s), step length (% limb length), and step width (% limb length).

Kinematic parameters.

We computed the anatomical angles for hip, knee, and ankle joints (in the sagittal plane), and trunk and pelvis (frontal, sagittal, and transverse plane). From these variables, we derived the RoM at each joint or segment, defined as the difference between the maximum and minimum value during the gait cycles.

Kinetic parameters.

Net internal joint moments (MomentAnkle, MomentKnee, MomentHip) were calculated with an inverse dynamics approach [30] and were normalized to the subject’s body weight. Joint moment curves were used to calculate the angular impulse (AI), i.e., the area under the joint moment curve within a specific time interval [31,32].

Angular impulse quantifies the total contribution of a joint moment to the production of movement and accounts for different gait adaptations (e.g., changes in walking speed) more accurately than peak moment values and it is defined as: where M is the flexor-extensor moment of the joint of interest, and Δt is the time interval used to calculate the integral. These angular impulses were hip extensor angular impulse during the first double support subphase (AI1stDS_Hip); hip flexor angular impulse during the second double support subphase (AI2ndDS_Hip); knee first and second extensor angular impulse (AI1st_Knee and AI2nd_Knee respectively) during the stance phase; ankle dorsiflexor angular impulse during the first double support subphase (AI1stDS_Ankle); ankle plantar flexor angular impulse during the mid-stance subphase (AIMidStance_Ankle); and ankle plantar flexor angular impulse during the second double support subphase (AI2ndDS_Ankle) (S1 Fig). We also evaluated the moment of support (MS) as follow: calculated as the sum of the moments MH, MK, and MA, which refer to the total curves of MomentHip MomentKnee and MomentAnkle, respectively.

Specifically, we considered the area (MSArea) within the gait cycles and the values of the two peaks of the curve (MSPeak1 and MSPeak2).

sEMG parameters. From the processed EMG signals, we calculated the simultaneous activation by considering the time-varying multi-muscle co-activation function (TMCf) proposed by Ranavolo and colleagues [33]: where d(t) is the mean of the differences, N is the number of muscles considered in the analysis, and EMGi is the sEMG signal of ith muscle. For each subject, data over individual strides were calculated and then averaged across cycles.

As co-activation indices, we considered the area of the TMCf (TMCfArea) within the gait cycles. We calculated the TMCf and TMCfArea by considering knee (RF-VL-VM vs BF-ST) and ankle (MG-LG vs TA) antagonistic muscles (TMCfKnee, TMCfArea_Knee, TMCfAnkle, TMCfArea_Ankle, respectively).

Patients’ subgroups classification.

In order to classify patients according to their kinematic behavior, we used a z-score with a one tailed z-test for statistical significance [34]. Thus, we chose a z-score of mean±1.5*SD (93% percentile) of the joint RoMs of the control group as the threshold for subgrouping patients with hereditary spastic paraparesis. This z-score threshold is considered as a fairly selective score used in several research fields [3538]. According to this criterion, each patient joint RoM could be either reduced (below threshold), increased (above threshold), or not significantly different from the values of the healthy controls. Thus, three subgroups of patients were identified. Subgroup one was patients with a statistically significant reduction of RoM at hip, knee, and ankle joints; subgroup two was patients with knee and ankle joint RoMs significantly reduced, but hip joint RoM not significantly different from the control value; and subgroup three was patients with hip joint RoM significantly increased, but ankle and knee joint RoMs not significantly different from the control values (Fig 1).

thumbnail
Fig 1. Patients’ subgroups classification according to lower limb joint kinematic behavior.

The threshold of mean±1.5*SD of the joint RoMs of the control group is reported. Each patient joint RoM could be either reduced (below the threshold), increased (above the threshold) or close to the values of healthy controls. It is possible to note that the subgroups’ division corresponds also to the severity of the disease scored by SPRS scale (higher values of SPRS correspond to higher disease’s severity). Each triangle represents a patient with different colors according to the SPRS scores. Each circle represents a qualitative characterization of subgroup of patients with different color shades corresponding, in the three spatial dimensions, to the SPRS scores.

https://doi.org/10.1371/journal.pone.0164623.g001

Statistical analysis

The Kolmogorov–Smirnov and Shapiro-Wilk tests were used to analyze the normal distribution of the data. Unpaired two-sample t-test or the Mann-Whitney test (two-tailed) were used for between-group differences in the demographic characteristics, time-distance parameters, joint kinematics, joint kinetics and sEMG values. Cohen's d values were also evaluated to estimate the effect size for the comparison between the two means. A multivariate ANOVA was used to compare demographic and clinical parameters (age, gender, disease onset and duration, Ashworth and SPRS scores) between subgroups of patients. One-way ANOVA was used to evaluate the differences in gait variables between the subgroups. Post hoc analyses (with Bonferroni’s corrections) were performed when significant differences were found with the ANOVA. Descriptive statistics included means ± SD, and significance level was set at p<0.05. Lower (LB) and upper (UB) bound of 95% confidence interval are reported for SPRS.

Results

Demographic and clinical characteristics

There were no significant differences between the group of patients with hereditary spastic paraparesis and the control group regarding gender, age, weight and height (all p>0.05). Multivariate analysis showed a significant main effect of the subgroup (F(14,84) = 4.468, p<0.001). Specifically, we observed significantly lower values of hip-spasticity Ashworth-score in both subgroups two (mean ± SD: 1.25 ± 0.86) and three (0.82 ± 0.88) as compared with subgroup one (2.23 ± 1.15) (p = 0.017 and p<0.001, respectively); significantly lower values of knee-spasticity Ashworth-score in subgroup three (0.88 ± 0.93) than both subgroup one (2.47 ± 1.12) and subgroup two (2.06 ± 0.68) (p<0.001 and p = 0.002, respectively); significantly lower values of total SPRS score in both subgroup two (16.07 ± 6.48, LB: 12.62, UB: 19.52) and subgroup three (5.64 ± 5.36, LB: 2.89, UB: 8.40) than subgroup one (21.33 ± 7.17, LB: 17.65, UB: 25.02), as well as in subgroup three than subgroup two (p = 0.049, p<0.001, p<0.001, respectively). No significant differences between subgroups were found for any other variable.

Time-distance parameters

When comparing the whole sample of patients with the healthy participants, no significant differences were found in any time-distance parameters, except for step width and step length, whose values were significantly increased and reduced, respectively, in patients compared with controls (Table 2). A significant effect of patients’ subgroup was found, using one-way ANOVA, on most of the time-distance parameters. These were the main effect of walking speed (F(2,47) = 6.703, p = 0.003); stance duration (F(2,47) = 4.923, p = 0.011); swing duration (F(2,47) = 4.900, p = 0.012); second double support duration (F(2,47) = 4.551, p = 0.016); and step length (F(2,47) = 11.173, p<0.001). Post-hoc analysis revealed significantly higher values of walking speed in subgroup three than in subgroup one, lower stance duration in subgroup three than in subgroup one, higher swing duration in both subgroups two and three than in subgroup one, lower second double support duration in subgroup three than in subgroup one and higher step length in subgroup three than in both subgroups one and two and in subgroup two than in subgroup one (Fig 2A).

thumbnail
Fig 2. Time-distance and joint and trunk kinematic parameters in HSP subgroups.

(A) Mean values (±SD) of time distance parameters. (B) (B) Mean (with SDs in light colors) kinematic plot of joint angular displacements during the gait cycle. (C) Mean values (±SD) of range of angular motion (RoM). Mean values of healthy controls for both time-distance and kinematic parameters, are reported in each bar graph (dotted line) and plot (black line). Asterisks indicate significant differences among the three subgroups at post hoc analysis (* p<0.05, ** p<0.001).

https://doi.org/10.1371/journal.pone.0164623.g002

thumbnail
Table 2. Time-distance, kinematic, kinetic and sEMG mean±SD data in patients and controls.

https://doi.org/10.1371/journal.pone.0164623.t002

Kinematic parameters

Significant lower values in knee and ankle RoMs and significant higher values in trunk lateral bending, flexion-extension, and rotation RoMs and pelvis rotation RoM were found in patients than in controls (Table 2). A significant main effect of the subgroup was found, using one-way ANOVA, on hip (F(2,47) = 33.747, p<0.001), knee (F(2,47) = 38.555, p<0.001), ankle (F(2,47) = 14.043, p<0.001), and pelvis tilt (F(2,47) = 4.328, p = 0.019) RoMs. Post-hoc analysis revealed significant higher values in hip RoM in both subgroups two and three than in subgroup one, higher values of knee RoM in both subgroups two and three than in subgroup one and in subgroup three than in subgroup two, higher values of ankle RoM in subgroup three than in both subgroups one and two, and lower values of pelvis tilt RoM in both subgroups two and three than in subgroup one (Fig 2B and 2C).

Kinetic parameters

Significant differences were found only for AI1st_Kneewhose value was higher in patients than controls (Table 2). A significant effect of the subgroup was found, using one-way ANOVA, on AI1stDS_Hip (main effect, F(2,47) = 3.517, p = 0.043). Post-hoc analysis showed lower values of this parameter in subgroup three than subgroup one (Fig 3A and 3B).

thumbnail
Fig 3. Joint kinetic and muscles parameters in HSP subgroups.

(A) Mean (with SDs in light colors) kinetic plot of joint moments (hip, knee and ankle) and support moment. The patterns are normalized to body weight and plotted vs. normalized stance. (B) Mean values (±SD) of kinetic parameters. (C) Mean TMCfArea (±SD) calculated for ankle and knee muscles (TMCfArea_Ankle and TMCfArea_Knee). (D) Mean (with SDs in light colors) plot of co-activation of ankle and knee joint antagonist muscles during the gait cycle. Mean values of healthy controls for kinetic and muscles parameters are reported in each bar graph (dotted line) and plot (black line). Asterisks indicate significant differences among the three subgroups at post hoc analysis (* p<0.05).

https://doi.org/10.1371/journal.pone.0164623.g003

sEMG parameters

Significant higher values in TMCfArea_Ankle were found in patients than controls (Table 2). No significant effect of the subgroup was found, using one-way ANOVA, on sEMG parameters (Fig 3C and 3D).

Discussion

We investigated the gait patterns in patients with hereditary spastic paraparesis by performing a comprehensive analysis of all time-distance, kinematic, kinetic, and sEMG parameters. In particular, our study was aimed at identifying specific subgroups of patients according to their kinematic behavior. Our assumption herein was that the decrease in the joint RoMs reflected the presence and extent of spasticity, and thus the primary deficit characterizing the gait of patients with hereditary spastic paraparesis. Few studies have previously investigated the gait in adults or children with hereditary spastic paraparesis [1820]. In line with these previous studies, we found an abnormal gait pattern characterized by reduced step length, increased step width, and reduced RoM at the knee joint in the whole sample of patients as compared with the control group. Furthermore, we found increased trunk RoM in all three spatial planes, increased pelvic tilt, increased hip joint torques (AI1st_Knee), reduced ankle joint RoM, and increased co-activation of muscles acting at the ankle joint. In addition to these general biomechanical characteristics of gait, one would expect some differential characteristics in distinct subgroups of patients according to clinical involvement of the pyramidal tract, given that patients with hereditary spastic paraparesis exhibit different degrees of severity both within and between families [1]. Thus, some specific biomechanical features may not emerge because they are hidden within their global walking strategy. Compared to previous studies, we enrolled a greater sample of patients with HSP (fifty in our study compared with twenty-two [18], nine [19] and twenty [20]) and performed an overall analysis of time-distance, kinematic (upper and lower body), kinetic and sEMG parameters. This allowed us to identify subgroups of patients and to define a global picture of walking strategies adopted by them. When subgrouping patients according to the hip, knee and ankle joint kinematic behavior, three clear gait patterns emerged. The gait pattern of subgroup one was characterized by reduced RoMs at hip, knee and ankle joints. Patients of this subgroup were the most severely affected (highest SPRS score) (Fig 1), and walked at the slowest speed. The gait pattern was characterized by the highest stance and second double support durations and the shortest swing duration and step length (Fig 2A). Such gait pattern reflects on one hand the reduced gait speed; on the other hand, the attempt to increase the most stable configuration duration (bipedal support), aimed at maintaining the dynamic balance. Furthermore, in these patients, we observed increased values of RoM for pelvis tilt and hip extensor angular impulse during the first double support subphase (Figs 2C and 3A and 3B). The former result might be due to spasticity and contracture of hip muscles as reported in neurological disorders with lower limb spasticity [20], [3941]. The last result indicates that HSP patients, although they have reduced hip RoMs, need to greatly involve the hip joint for weight acceptance increasing the internal torques. Interestingly, this finding further reinforces the notion that spasticity predominates on muscle weakness in the most severely involved patients [42,43]. The gait pattern of subgroup three was characterized by increased hip joint RoM and knee and ankle joint RoMs close to control values. These patients were the most mildly affected (lowest SPRS score) (Fig 1) and showed the highest walking speed. Their gait pattern included the highest swing duration and step length, the shortest stance and second double support duration values (Fig 2A) and the lowest pelvis tilt RoM and hip extensor angular impulse (during the first double support sub phase) than the other subgroups (Figs 2C and 3A and 3B). From Fig 2, it is possible to note that patients of this subgroup showed a gait pattern which was close to that of healthy controls in terms of time-distance parameters. It is also important to note that, with respect to healthy subjects, this subgroup of patients showed increased trunk RoM in all three spatial planes. Considering the very low SPRS score in this subgroup, this result suggests that the compensatory mechanisms represented by the increased trunk movements and hip RoM are the most important biomechanical features characterizing the gait disorders from the early phase of the disease. Patients of subgroup two had characteristics between those of subgroups one and three, in terms of disease severity (Fig 1) and gait speed, and showed hip joint RoM close to controls but decreased knee and ankle joint RoMs. In particular, their gait pattern was characterized by intermediate values with respect to the other two subgroups in terms of step length, swing duration and pelvis tilt RoM (Fig 2).

As regards the sEMG, we observed significantly increased co-activation of antagonist muscles acting at the ankle for the whole group of HSP patients compared with healthy controls (Table 2). When analyzing patients’ subgroups, a trend, although not statistically significant, for higher co-activation values of knee and ankle antagonist muscles was also observed (Fig 3C and 3D). Such a finding reflects the inability of the CNS to selectively activate lower limb joint muscles and may be explained by inefficient mechanisms of reciprocal inhibition [44] and the supraspinal and spinal plastic neuronal changes associated with the development of spasticity [45,46]. In general the different patterns of gait disturbance seem to correlate fairly well with the different degree of disease’s severity among patients, as measured by the SRPS score, with mild (SPRS < ten), moderate (SPRS < twenty) and more severe presentation (SPRS > twenty). In addition, since we chose to study only patients with pure pyramidal signs, regardless of the genetic form (namely SPG3A, 4, 5, 7 and 31), we can safely assume that the identified three different gait patterns based on lower limb kinematic behavior, also reflect the different degree of pyramidal tract involvement in individual patients. We think that identifying specific gait patterns [41] in patients with hereditary spastic paraparesis may be useful in: i) improving our understanding on gait disorder in hereditary spastic paraparesis by sorting out the most meaningful gait features from the complexity of locomotion; ii) recognizing specific abnormalities and their impact on clinical decision-making; and iii) individualizing rehabilitative treatment and better evaluating its effects over the time.

Supporting Information

S1 Data. Main data underlying the findings described in the manuscript.

https://doi.org/10.1371/journal.pone.0164623.s001

(XLSX)

S1 Fig. Method of angular impulse evaluation on hip, knee, and ankle joint moments curves.

https://doi.org/10.1371/journal.pone.0164623.s002

(TIF)

Author Contributions

  1. Conceptualization: MS MR AR.
  2. Data curation: MR MS AR.
  3. Formal analysis: MR MS GM.
  4. Funding acquisition: FP.
  5. Investigation: MR MS.
  6. Methodology: MS MR AR C. Conte GC.
  7. Project administration: FP FL.
  8. Resources: MS MR LL C. Casali FP.
  9. Software: MR GM.
  10. Supervision: MS C. Casali FP.
  11. Validation: MS AR FL.
  12. Visualization: MS C. Casali FP.
  13. Writing – original draft: MS MR AR.
  14. Writing – review & editing: MS MR AR GM TV C. Casali FP FL FD.

References

  1. 1. Lo Giudice T, Lombardi F, Santorelli FM, Kawarai T, Orlacchio A. Hereditary spastic paraplegia: clinical-genetic characteristics and evolving molecular mechanisms. Exp Neurol. 2014; 261: 518–539. pmid:24954637
  2. 2. Faber I, Servelhere KR, Martinez AR, D'Abreu A, Lopes-Cendes I, França MC Jr. Clinical features and management of hereditary spastic paraplegia. Arq Neuropsiquiatr. 2014; 72(3): 219–226. pmid:24676440
  3. 3. Klimpe S, Schüle R, Kassubek J, Otto S, Kohl Z, Klebe S, et al. Disease severity affects quality of life of hereditary spastic paraplegia patients. Eur J Neurol. 2012; 19(1): 168–171. pmid:21631647
  4. 4. Orsucci D, Petrucci L, Ienco EC, Chico L, Simi P, Fogli A, et al. Hereditary spastic paraparesis in adults. A clinical and genetic perspective from Tuscany. Clin Neurol Neurosurg. 2014; 120: 14–19. pmid:24731568
  5. 5. Fink JK. Hereditary spastic paraplegia. Curr Neurol Neurosci Rep. 2006; 6(1): 65–76. pmid:16469273
  6. 6. Zhang Y, Roxburgh R, Huang L, Parsons J, Davies TC. The effect of hydrotherapy treatment on gait characteristics of hereditary spastic paraparesis patients. Gait Posture. 2014; 39(4): 1074–1079. pmid:24556467
  7. 7. Heetla HW, Halbertsma JP, Dekker R, Staal MJ, van Laar T. Improved gait performance in a patient with hereditary spastic paraplegia after a continuous intrathecal baclofen test infusion and subsequent pump implantation: a case report. Arch Phys Med Rehabil. 2015; 96(6): 1166–1169. pmid:25626112
  8. 8. Bertolucci F, Di Martino S, Orsucci D, Ienco EC, Siciliano G, Rossi B, et al. Robotic gait training improves motor skills and quality of life in hereditary spastic paraplegia. NeuroRehabilitation. 2015; 36(1): 93–99. pmid:25547770
  9. 9. Adams RD, Victor M, Ropper AH, Samuels MA, Klein J, editors. Principles of Neurology, 10th ed, New York, McGraw-Hill Education Medical, 2015.
  10. 10. Rowland LP, Pedley TA, editors. Neurologia de Merritt, 12nd ed, Nederland, Wolters Kluwer Health, 2011.
  11. 11. Martinuzzi A, Montanaro D, Vavla M, Paparella G, Bonanni P, Musumeci O, et al. Clinical and Paraclinical Indicators of Motor System Impairment in Hereditary Spastic Paraplegia: A Pilot Study. PLoS One. 2016;11(4). pmid:27077743
  12. 12. Henson BJ, Zhu W, Hardaway K, Wetzel JL, Stefan M, Albers KM, et al. Transcriptional and post-transcriptional regulation of SPAST, the gene most frequently mutated in hereditary spastic paraplegia. PLoS One. 2012;7(5). pmid:22574173
  13. 13. Wedding IM, Koht J, Tran GT, Misceo D, Selmer KK, Holmgren A, et al. Spastic paraplegia type 7 is associated with multiple mitochondrial DNA deletions. PLoS One. 2014;9(1). pmid:24466038
  14. 14. Beetz C, Schüle R, Deconinck T, Tran-Viet KN, Zhu H, Kremer BP. REEP1 mutation spectrum and genotype/phenotype correlation in hereditary spastic paraplegia type 31. Brain. 2008; 131(4): 1078–86. pmid:18321925
  15. 15. Goizet C, Boukhris A, Durr A, Beetz C, Truchetto J, Tesson C. CYP7B1 mutations in pure and complex forms of hereditary spastic paraplegia type 5. Brain. 2009; 132(6): 1589–600. pmid:19439420
  16. 16. Solowska JM, Baas PW. Hereditary spastic paraplegia SPG4: what is known and not known about the disease. Brain. 2015; 138(9): 2471–84.
  17. 17. Wilkinson PA, Crosby AH, Turner C, Bradley LJ, Ginsberg L, Wood NW. A clinical, genetic and biochemical study of SPG7 mutations in hereditary spastic paraplegia. Brain. 2004; 127(5): 973–80. pmid:14985266
  18. 18. Klebe S, Stolze H, Kopper F, Lorenz D, Wenzelburger R, Volkmann J, et al. Gait analysis of sporadic and hereditary spastic paraplegia. J Neurol. 2004; 251(5): 571–578. pmid:15164190
  19. 19. Piccinini L, Cimolin V, D'Angelo MG, Turconi AC, Crivellini M, Galli M. 3D gait analysis in patients with hereditary spastic paraparesis and spastic diplegia: a kinematic, kinetic and EMG comparison. Eur J Paediatr Neurol. 2011; 15(2): 138–145. pmid:20829081
  20. 20. Marsden J, Ramdharry G, Stevenson V, Thompson A. Muscle paresis and passive stiffness: key determinants in limiting function in Hereditary and Sporadic Spastic Paraparesis. Gait Posture. 2012; 35(2): 266–271. pmid:22050971
  21. 21. Schüle R, Holland-Letz T, Klimpe S, Kassubek J, Klopstock T, Mall V, et al. The Spastic Paraplegia Rating Scale (SPRS): a reliable and valid measure of disease severity. Neurology. 2006; 67(3): 430–434. pmid:16894103
  22. 22. Davis RB III, Õunpuu S, Tyburski D, Gage JR. A gait analysis data collection and reduction technique. Hum Mov Sci. 1991; 10: 575–587.
  23. 23. Winter DA, editor. Biomechanics of Human Movement, 2nd ed, New York, Wiley and Sons, 1979.
  24. 24. Hermens HJ, Freriks B, Disselhorst-Klug C, Rau G. Development of recommendations for SEMG sensors and sensor placement procedures. J Electromyogr Kinesiol. 2000; 10(5): 361–374. pmid:11018445
  25. 25. Barbero M, Merletti R, Rainoldi A, editors. Atlas of Muscle Innervation Zones: Understanding Surface Electromyography and Its Applications, New York, Springer, 2012. https://doi.org/10.1007/978-88-470-2463-2
  26. 26. Mari S, Serrao M, Casali C, Conte C, Martino G, Ranavolo A, et al. Lower limb antagonist muscle co-activation and its relationship with gait parameters in cerebellar ataxia. Cerebellum. 2014; 13: 226–236. pmid:24170572
  27. 27. Borghese NA, Bianchi L, Lacquaniti F. Kinematic determinants of human locomotion. J Physiol. 1996; 494(3): 863–879.
  28. 28. Vasudevan EV, Torres-Oviedo G, Morton SM, Yang JF, Bastian AJ. Younger is not always better: development of locomotor adaptation from childhood to adulthood. J Neurosci. 2011; 31(8): 3055–3065. pmid:21414926
  29. 29. Martino G, Ivanenko YP, Serrao M, Ranavolo A, d'Avella A, Draicchio F, et al. Locomotor patterns in cerebellar ataxia. J Neurophysiol. 2014; 112: 2810–2821. pmid:25185815
  30. 30. Vaughan C, Brian D, O’Connor J. Dynamics of Human Gait, 2nd ed. Champaign, IL: Human Kinetics Publishers; 1992.
  31. 31. De Vita P, Lassiter T Jr, Hortobagyi T, Torry M. Functional knee brace effects during walking in patients with anterior cruciate ligament reconstruction. Am J Sports Med. 1998; 26: 778–784. pmid:9850778
  32. 32. Don R, Serrao M, Vinci P, Ranavolo A, Cacchio A, Ioppolo F, et al. Foot drop and plantar flexion failure determine different gait strategies in Charcot-Marie-Tooth patients. Clin Biomech. 2007; 22(8): 905–916. pmid:17686557
  33. 33. Ranavolo A, Mari S, Conte C, Serrao M, Silvetti A, Iavicoli S, et al. New muscle co-activation index for biomechanical load evaluation in work activities. Ergonomics. 2015; 58(6): 966–979. pmid:25555042
  34. 34. van der Waal JM, Terwee CB, van der Windt DA, Bouter LM, Dekker J. The impact of non-traumatic hip and knee disorders on health-related quality of life as measured with the SF-36 or SF-12. A systematic review. Qual Life Res. 2005; 14(4): 1141–1155. pmid:16041909
  35. 35. Meyer AC, Boscardin WJ, Kwasa JK, Price RW. Is it time to rethink how neuropsychological tests are used to diagnose mild forms of HIV-associated neurocognitive disorders? Impact of false-positive rates on prevalence and power. Neuroepidemiology. 2013; 41(3–4): 208–216. pmid:24157541
  36. 36. Wasser K, Pilgram-Pastor SM, Schnaudigel S, Stojanovic T, Schmidt H, Knauf J, et al. New brain lesions after carotid revascularization are not associated with cognitive performance. J Vasc Surg. 2011; 53(1): 61. pmid:20875716
  37. 37. France CR, Rhudy JL, McGlone S. Using normalized EMG to define the nociceptive flexion reflex (NFR) threshold: further evaluation of standardized NFR scoring criteria. Pain. 2009; 145(1–2): 211–218. pmid:19595510
  38. 38. Weinstein G, Beiser AS, Decarli C, Au R, Wolf PA, Seshadri S. Brain imaging and cognitive predictors of stroke and Alzheimer disease in the Framingham Heart Study. Stroke. 2013; 44(10): 2787–2794. pmid:23920020
  39. 39. Mao Y, Chen P, Li L, Li L, Huang D. Changes of pelvis control with subacute stroke: A comparison of body-weight- support treadmill training coupled virtual reality system and over-ground training. Technol Health Care. 2015; 23 Suppl 2: S355–S364. pmid:26410502
  40. 40. De Quervain IA, Simon SR, Leurgans S, Pease WS, McAllister D. Gait pattern in the early recovery period after stroke. J Bone Joint Surg Am. 1996; 78(10): 1506–1514. pmid:8876578
  41. 41. Roche N, Pradon D, Cosson J, Robertson J, Marchiori C, Zory R. Categorization of gait patterns in adults with cerebral palsy: a clustering approach. Gait Posture. 2014; 39(1): 235–240. pmid:23948331
  42. 42. McDermott C, White K, Bushby K, Shaw P. Hereditary spastic paraparesis: a review of new developments. J Neurol Neurosurg Psychiatry. 2000; 69(2):150–160. Review. pmid:10896685
  43. 43. Schüle R, Wiethoff S, Martus P, Karle KN, Otto S, Klebe S, et al. Hereditary Spastic Paraplegia -clinico-genetic lessons from 608 patients. Ann Neurol. 2016. pmid:26856398
  44. 44. Meunier S, Pol S, Houeto JL, Vidailhet M. Abnormal reciprocal inhibition between antagonist muscles in Parkinson’s disease. Brain J Neurol. 2000; 123(5): 1017–1026. pmid:10775546
  45. 45. Simmons RW, Richardson C. Peripheral regulation of stiffness during arm movements by coactivation of the antagonist muscles. Brain Res. 1988; 473: 134–140. pmid:3208115
  46. 46. Falconer K, Winter DA. Quantitative assessment of co-contraction at the ankle joint in walking. Electromyogr Clin Neurophysiol. 1985; 25: 135–149. pmid:3987606