Participant characteristics
Table 1 summarizes the participant characteristics. Median age was 68.5 (interquartile range 63-74.25) and median disease duration was 11 years (interquartile range 7.75–13.25). One of the sixteen participants was excluded from analysis because of data loss due to a detached ECG sensor. Median heart rate during the resting periods was 83.9 bpm (interquartile range 79.6–88.1 bpm); median heart rate variability during rest, defined as the coefficient of variation, was 1.6% (interquartile range 1.0-2.2%).
Table 1
participant characteristics
| Median | Interquartile range |
Age | 68.5 | 63-74.25 |
Gender (% male) | 69 | |
Disease duration | 11 | 7.75–13.25 |
Daily levodopa dosage (mg) | 1120 | 500–1300 |
MDS-UPDRS part III | 23.5 | 19-26.75 |
H&Y | 2 | 2–3 |
MMSE | 27 | 26–29 |
FAB | 16 | 14.75–17.25 |
NFOGQ | 19.5 | 16.5-23.25 |
HR rest (bpm) | 83.5 | 78.4–87.2 |
HRV rest (%) | 1.5 | 1.0-2.1 |
MDS-UPDRS part III: Movement Disorders Society Unified Parkinson Disease Rating scale part III, H&Y: Hoehn and Yahr Staging Scale, MMSE: mini-mental state examination (range 0–30), FAB: Frontal Assessment Battery (range 0–18), NFOGQ: New Freezing of Gait Questionnaire (range 0–28), HR rest: heart rate during the first 10 seconds of standing preceding each round, HRV rest: heart rate variability defined as the coefficient of variation during the first 10 seconds of standing preceding each round.
FOG annotations
In total, 4683 gait events and 685 freezing episodes were observed with a high interrater agreement (kappa correlation coefficient 0.90; Spearman correlation for number and total duration of FOG of 0.84 and 0.89 respectively). The median duration of the FOG events was 5.8 seconds (interquartile range 3.4–11.5 sec.).
Figure 2.B-D shows the distribution of FOG episodes over the participants for each FOG type, triggering event, and DT condition. 81% of all freezing was annotated as trembling, but was often seen in combination with shuffling (6%), as reflected by a low interrater agreement for both types (in 86% of the shuffling events, the raters disagreed on the FOG type). Because both likely share a similar pathophysiological substrate that is different from akinesia (31, 32), we decided to compare the combination of trembling-shuffling with akinesia. Notably, akinesia -accounting for 13% of the events- was mainly present in two of the participants (see Fig. 2.B, orange). Turning, as expected from the literature (33–35), was most provocative for FOG (73% of the events), followed by narrow passages (19%) and starting hesitation (6%). The first two, turning and narrow passages, were considered for further analysis. The number of FOG events was comparable over the three DT conditions: 36% of the episodes occurred during noDT, 35% during cDT, and 29% during mDT. After exclusion of trials that were close (< 6 seconds) to other FOG episodes, 406 FOG events, 377 stops, and 1446 normal gait events were included in the analysis.
Freezing index (FI)
Figure 3.A shows the overall time course of the FI during FOG, normal gait events and stopping. Overall, the FI increased over time during FOG (orange), while it remained at the same level during normal gait events (green). However, the FI also increased during stopping (blue).
The first model (model 1.a), evaluating the increase of the FI during preFOG and FOG compared to baseline levels, showed that the FI increased significantly over time (main effect of time: FI(time): p < 0.001), but that its performance was dependent on the type of FOG (significant interaction effect between time and type: FI(time*type): p < 0.001) (Table 2). Figure 4.A-C shows the results of the post-hoc analyses of this model. The FI increased significantly for trembling during preFOG (post-hoc FI(preFOG−baseline): B = 0.39, p = 0.018) and during FOG compared to baseline (post-hoc FI(FOG−baseline): B = 1.87, p < 0.001), but not for akinesia (post-hoc FI(FOG−baseline): B = 0.48, p = 0.270) (Fig. 4.A). No significant interaction effects were found for trigger or for DT. The significant main effects for type (p = 0.024), for trigger (p = 0.033), and for DT (p < 0.001) indicate that the FI was different at baseline for these conditions. Overall, these results show that the FI performed as expected (i.e., significant increase compared to baseline) to detect the trembling type of FOG in all tested circumstances, but not for the akinetic type.
Figure 5.A-C shows the results the second model (model 2.a) that compared the FI during FOG, normal gait events, and stops under all conditions. For the ideal FOG detector, the green (FOG – normal gait event), and blue (FOG - stopping) point ranges would appear above the zero threshold line. Significant interaction effects were found between condition (i.e., FOG-normal gait event, or FOG-stop) and all three factors (type, trigger and DT) (FI(condition*type): p < 0.001; FI(condition*trigger): p < 0.001; FI(condition*DT): p = 0.012), indicating that the performance of the FI differed significantly between the conditions (Table 2).
Firstly, the distinctiveness of the FI was dependent on the type of FOG (Fig. 5.A). More specifically, for trembling FOG, the difference between FOG and a normal gait event was significant (post-hoc FI(FOG−normal gait event): B = 1.92, p < 0.001), but not between FOG and stopping (post-hoc FI(FOG−stop): B = 0.37, p = 0.559). In contrast, for akinetic FOG, the difference between FOG and a normal gait event was smaller (post-hoc FI(FOG−normal gait event): B = 0.58, p = 0.019), and between FOG and stopping even reversed, i.e. the FI during akinetic FOG was lower than during stopping (post-hoc FI(FOG−stop): B=-1.21, p < 0.001). Even though this reversed relationship was significant, this has no clinical advantage because the FI should rise above a certain threshold to identify a FOG episode. A significant main effect of type of the second model indicated again that FI during akinesia was lower than during trembling (p < 0.001) (Table 2).
Secondly, the relation of the FI between FOG and a normal gait event depended on the trigger (Fig. 5.B). In a narrow passage, the FI was higher during an event with FOG than without FOG (post-hoc FI(FOG−normal gait event): B = 2.08, p < 0.001); but for turning, the relationship became smaller (post-hoc FI(FOG−normal gait event): B = 0.42, p = 0.002).
Thirdly, the difference in FI between FOG and a normal gait event varied between DT conditions, with slightly smaller differences between FOG and a normal gait event during cDT and mDT than during noDT (Fig. 5.C).
Taken together, these results show that the FI performs well to differentiate FOG from normal gait events (although only marginally significant for akinetic FOG, or when comparing FOG turns with normal turns), but not to differentiate FOG from stopping.
Heart rate
Figure 3.B shows that overall, heart rate increased during FOG compared to baseline levels, while a clear drop in heart rate was seen during stopping, which is in line with previous studies (14, 36, 37). However, in contrast to previous findings (14), in our results heart rate increased during normal gait events.
Indeed, the first model (model 1.b) showed a significant main effect of time (HR(time): p = 0.006), with post-hoc analysis revealing a significant increase in heart rate during FOG compared to baseline (HR(FOG−baseline): B = 0.29, p = 0.001), but not during preFOG (HR(preFOG−baseline): B = 0.69, p = 0.253) (Table 2). The increase in heart rate was dependent on type (HR(time*type): p < 0.001), and on DT condition (HR(time*DT): p = 0.011), but not on the triggering event (HR(time*trigger): p = 0.312).
Namely, Fig. 4.D shows a clear and significant increase in heart rate for trembling during preFOG (post-hoc HR(preFOG−baseline): B = 0.36, p = 0.029) and FOG periods (post-hoc HR(FOG−baseline): B = 1.26 p < 0.001), while the increase during akinesia was not significant in these periods (post-hoc HR(preFOG−baseline): B = 0.27, p = 0.729 and post-hoc HR(FOG−baseline): B = 0.12, p = 0.896 respectively). Considering the DT conditions, heart rate increased significantly during FOG for the noDT and the cDT condition (post-hoc HR(FOG−baseline): B = 0.69, p = 0.009 and HR(FOG−baseline): B = 1.13, p < 0.001), but not for the mDT condition (HR(FOG−baseline): B = 0.26, p = 0.525) (Fig. 4.I).
Of less interest, a significant main effect of the first model was found for trigger (HR(trigger): p = 0.006) and DT (HR(DT): p < 0.001), meaning that the heart rate levels for these conditions differed at baseline (Table 2). Namely, heart rate was higher during baseline levels of narrow passages than during turning (post-hoc HR(narrow passage − turn): B = 1.47, p = 0.007), and heart rate during baseline levels of mDT was higher than during noDT (post-hoc HR(mDT−noDT): B = 1.35, p = 0.020) or cDT (HR(mDT−cDT): B = 2.20, p < 0.001).
Overall, these results confirm that heart rate increases before and during trembling FOG, but stays at a similar level during akinetic FOG. Furthermore, the heart rate during the mDT condition had significantly higher baseline levels and did not increase significantly during FOG.
The second model (model 2.b) investigated the differences in heart rate changes during FOG with stopping and normal gait events (Fig. 5.D-F). Of major interest, a significant interaction effect was found between condition and type (HR(condition*type): p = 0.005) (Table 2, Fig. 5.D). For both trembling and akinesia, heart rate change was significantly different from stopping (post-hoc HR(FOG−stop): B = 2.41, p < 0.001 and HR(FOG−stop): B = 1.14, p = 0.004 respectively). However, heart rate change during a normal gait event was not different from trembling (post-hoc HR(FOG−normal gait event): B = 0.18, p = 0.473), or from akinesia (post-hoc HR(FOG−normal gait event): B=-0.86, p = 0.058).
Moreover, a significant interaction effect was observed between condition and DT with a slightly smaller difference in heart rate change between FOG and normal gait event for the noDT condition than the other two conditions (Fig. 5.F). In addition, a significant main effect of the second model was seen for type and DT (HR(type): p = 0.005; HR(DT): p = 0.0003) indicating differences in heart rate change for these factors during the FOG condition (Table 2). More specifically, heart rate change was greater during trembling compared to akinesia (post-hoc HR(trembling−akinesia): B = 0.16, p = 0.0014); and greater during cDT than noDT (post-hoc HR(cDT−noDT): B = 0.46, p < 0.0001).
In conclusion, heart rate change was significantly different between FOG and stopping of any type and triggering situation, but not between FOG and a normal gait event.
Table 2
results of the type III F test (with Kenward-Roger adjustment for degrees of freedom) for the fixed effects of the mixed models
| Model 1.a: FI ~ type*time + trigger*time + DT*time + (1|participant) | Model 1.b: heart rate ~ type*time + trigger*time + DT*time + (1│participant) + (1|trial) |
| F | Df | Df.res | p value | F | Df | Df.res | p value |
(Intercept) | 170.70 | 1 | 24 | < 0.001* | 1355.21 | 1 | 14.907 | < 0.001* |
time | 15.00 | 2 | 1189 | < 0.001* | 5.13 | 2 | 802 | 0.006* |
trigger | 4.55 | 1 | 1195 | 0.033* | 7.64 | 1 | 468 | 0.006* |
type | 5.13 | 1 | 1180 | 0.024* | 2.48 | 1 | 426 | 0.116 |
DT | 9.27 | 2 | 1191 | < 0.001* | 10.76 | 2 | 475 | < 0.001* |
time*type | 10.35 | 2 | 1189 | < 0.001* | 7.27 | 2 | 802 | 0.001* |
time*trigger | 0.52 | 2 | 1189 | 0.595 | 1.17 | 2 | 802 | 0.312 |
time*DT | 2.33 | 4 | 1189 | 0.054 | 3.28 | 4 | 802 | 0.011* |
| Model 2.a: FI ~ type*condition + trigger*condition + DT*condition + (1│participant) | Model 2.b: heart rate change ~ type*condition + trigger*condition + DT*condition + (1|participant) |
| F | Df | Df.res | p value | F | Df | Df.res | p value |
(Intercept) | 325.54 | 1 | 24 | < 0.001* | 2.91 | 1 | 292 | 0.089 |
condition | 71.69 | 2 | 2204 | < 0.001* | 60.15 | 2 | 2167 | < 0.001* |
trigger | 0.91 | 1 | 2205 | 0.340 | 2.75 | 1 | 2199 | 0.097 |
DT | 1.99 | 2 | 2202 | 0.137 | 5.92 | 2 | 2212 | 0.003* |
type | 32.67 | 1 | 2212 | < 0.001* | 7.82 | 1 | 1687 | 0.005* |
condition*trigger | 43.10 | 2 | 2202 | < 0.001* | 2.43 | 2 | 2209 | 0.088 |
condition*DT | 3.24 | 4 | 2201 | 0.012* | 2.86 | 4 | 2210 | 0.022* |
condition*type | 17.19 | 2 | 2207 | < 0.001* | 5.26 | 2 | 1985 | 0.005* |
*p < 0.05 ; F = F values; Df = degrees of freedom; Df.res = residual degrees of freedom; FI = freezing index, DT = dual-task |