3.1 Study population
In total, 34 children aged 6.0 to 12.6 years were enrolled in the study, with 16 and 18 participants randomised into 20Hz and 25Hz groups, respectively (Figure 1). The demographic characteristics of the study population are presented in Table 1. Four participants withdrew from the study (Figure 1): one soon after the control assessment for having a semi-elective surgery scheduled; one after 8 weeks and two after 12 weeks of VT due to lack of time to perform the sessions. Please note that the study was partially conducted while New Zealand was under COVID-19 lockdown restrictions (30). This scenario markedly impacted our ability to recruit participants and perform the assessments, also affecting the ability of some participants to undergo VT in school settings. Three participants completed 20 weeks of VT but were unable to attend the final 20VT assessment (T3).
Table 1. Characteristics of study participants
Parameters
|
|
Group 20 Hz
|
Group 25 Hz
|
n
|
|
16
|
18
|
Age (years)
|
|
9.5 [4.5, 11.7]
|
9.2 [6.9, 10.4]
|
Sex
|
Females
|
6 (37%)
|
7 (39%)
|
Ethnicity
|
NZ European
|
13 (81%)
|
11 (61%)
|
|
Māori
|
2 (13%)
|
4 (22%)
|
|
Other
|
1 (6%)
|
3 (17%)
|
GMFCS
|
Level I
|
6 (37%)
|
7 (39%)
|
|
Level II
|
7 (44%)
|
8 (44%)
|
|
Level III
|
3 (19%)
|
3 (17%)
|
CP type
|
Spastic
|
13 (81%)
|
13 (72%)
|
|
Dystonic
|
3 (19%)
|
2 (11%)
|
|
Ataxic
|
nil
|
1 (6%)
|
|
Unknown
|
nil
|
2 (11%)
|
CP, Cerebral palsy; GMFCS, Gross Motor Function Classification System.
Age data are median [quartile 1, quartile 3] and categorical data are n (%).
3.2 Effects of VT frequency and duration
There were no observed differences between the two VT frequencies (i.e., 20 Hz vs 25 Hz) on study outcomes, but there were some differences associated with VT duration (i.e., 12 vs 20 weeks). As a result, study outcomes are reported for the overall pairwise differences (Control period vs 12VT and Control vs 20VT), except for outcomes with an observed effect of VT duration, for which differences between 12VT and 20VT are also reported.
3.3 Mobility (primary outcome)
The results of mobility outcomes are presented in Table 2. For the primary outcome, 20VT (but not 12VT) led to improvements in the 6MWT, with participants covering additional 23 meters according to both ITT (p=0.007) and PPA (p=0.011) analyses (Table 2), with distance milestones reached progressively faster (Figure 2). Participants also showed improvements in the 10MWT (Table 2), with an increase of 0.18 m/s in gait speed after 20 VT (p=0.047).
Table 2. Mobility parameters outcomes
Parameters
|
n
|
Control period
|
12VT
|
20VT
|
12VT vs Control
|
20VT vs Control
|
20VT vs 12VT
|
6MWT ITT (m)
|
34
|
406 (391, 422)
|
407 (391, 423)
|
429 (413, 446)
|
0 (-12, 13)
|
23 (6, 39)**
|
22 (9, 36)††
|
6MWT PPA (m)
|
29
|
421 (404, 438)
|
422 (405, 439)
|
444 (427, 462)
|
1 (-12, 15)
|
23 (6, 41)*
|
22 (8, 37)††
|
10MWT (m/s)
|
25
|
2.27 (2.12, 2.42)
|
2.36 (2.21, 2.51)
|
2.46 (2.29, 2.62)
|
0.09 (-0.05, 0.26)
|
0.18 (0.00, 0.37)*
|
0.10 (-0.06, 0.25)
|
6MWT, 6-minute walk test; 10MWT, 10-meter walk test; 12VT, assessment after 12 weeks of vibration therapy; 20VT, assessment after 20 weeks of vibration therapy; ITT, intention-to-treat analyses; PPA, per-protocol analysis.
Data at each assessment are the adjusted means and 95% confidence intervals (CI), while differences between assessments are the adjusted mean differences and 95% CI; all values were derived from linear mixed models based on repeated measures including the participant's GMFCS level, randomisation group (20 Hz / 25Hz), and the baseline value of the outcome as a covariate.
n is the number of participants at baseline; the number of participants who completed a given assessment is provided in Additional file 1.
p-values for statistically significant differences (at p<0.05) between two given assessments are shown in bold; *p < 0.05 and **p < 0.01 for pairwise differences compared to the Control period; ††p<0.01 for a difference between 12VT and 20VT.
3.4 Gross motor function
Both 12VT and 20VT led to improvements in gross motor function, measured by both GMFM-D and GMFM-E (Table 3). GMFM-E scores improved by 2.5 points after 12VT (+3.5%; 95% CI 2.3, 4.7%; p<0.0001) and by 3.7 points after 20VT (+5.1%; 95% CI 3.6, 6.6%; p<0.0001) (Table 3). For GMFM-D, there was a significant interaction between GMFCS level and assessment (p=0.0009), indicating a differential response to VT. Among participants with GMFCS level I and II, GMFM-D scores increased by 0.8 points after 12VT (+2.1%; 95% CI 0.6, 3.6%; p=0.008) and by 1.3 points after 20VT (+3.4%; 95% CI 1.4, 5.4%; p=0.001) (Table 3). There was a greater increase in GMFM-D scores among children with GMFCS level III of 3.0 points after 12VT (+7.7%; 95% CI 4.1, 11.3%; p<0.0001) and 5.0 points after 20VT (+12.8%; 95% CI 8.5, 17.1%; p<0.0001) (Table 3).
After 20VT, the extra 8 weeks of VT lead to an additional 1.2-point increase in GMFM-E scores (+1.6%; 95% CI 0.3, 2.9%; p=0.011), as well as a 2.0-point increase in GMFM-D scores for participants with GMFCS level III (+5.1%; 95% CI 1.5, 8.7%; p=0.006) (Table 3).
Table 3. Gross motor function outcomes
Parameters, units
|
n
|
Control period
|
12VT
|
20VT
|
12VT vs
Control
|
20VT vs
Control
|
20 VT vs
12VT
|
GMFM-D (level I-II), score
|
25
|
32.8 (31.9, 33.7)
|
33.6 (32.7, 34.5)
|
34.1 (33.2, 35.0)
|
0.8 (0.2, 1.4)**
|
1.3 (0.5, 2.1)**
|
0.5 (-0.1, 1.2)
|
GMFM-D (level III), score
|
5
|
26.7 (23.4, 30.1)
|
29.7 (26.4, 33.1)
|
31.7 (28.4, 35.1)
|
3.0 (1.6, 4.4)****
|
5.0 (3.3, 6.7)****
|
2.0 (0.6, 3.4)††
|
GMFM-E, score
|
30
|
49.2 (47.6, 50.8)
|
51.7 (50.1, 53.3)
|
52.9 (51.3, 54.5)
|
2.5 (1.6, 3.3)****
|
3.7 (2.6, 4.8)****
|
1.2 (0.3, 2.2)†
|
12VT, assessment after 12 weeks of vibration therapy; 20VT, assessment after 20 weeks of vibration therapy; GMFM-D, gross motor function measure dimension D; GMFM-E, gross motor function measure dimension E.
Data at each assessment are the adjusted means and 95% confidence intervals (CI), while differences between assessments are the adjusted mean differences and 95% CI; GMFM-E values were derived from linear mixed models based on repeated measures including the participant's GMFCS level, randomisation group (20 Hz / 25Hz), and the baseline value of the outcome as a covariate; GMFM-D values were similarly derived but with the inclusion of an interaction term between GMFCS level and assessment.
n is the number of participants at baseline; the number of participants who completed a given assessment is provided in Additional file 1.
Statistically significant differences (at p<0.05) between assessments are shown in bold; *p < 0.05, **p < 0.01 and ****p < 0.0001 for pairwise differences compared to the Control period; †p<0.05 and ††p<0.01 for differences between 12VT and 20VT.
Body composition
VT led to no observed changes in anthropometry (i.e., height, weight, and BMI z-scores), lean mass, or fat mass (Table 4). In contrast, spine aBMD z-scores increased by 0.14 after 20VT (p=0.015), with a 1.5-g improvement also seen in spine BMC (L1-L4) after 12VT (p=0.046) that was not detected after 20VT (p=0.09) (Table 4). There were no observed changes in TBLH aBMD, TBLH BMC, or spine aBMD (Table 4).
Table 4. Anthropometry and body composition outcomes
Parameters (units)
|
n
|
Control period
|
12VT
|
20VT
|
12VT vs Control
|
20VT vs Control
|
Anthropometry
|
Height
(z-score)
|
30
|
-0.63 (-0.72, -0.55)
|
-0.60 (-0.69, -0.52)
|
-0.58 (-0.67, -0.49)
|
0.03 (-0.05, 0.12)
|
0.06 (-0.05, 0.16)
|
Weight
(z-score)
|
30
|
-0.50 (-0.60, -0.39)
|
-0.51 (-0.62, -0.40)
|
-0.40 (-0.52, -0.29)
|
-0.02 (-0.14, 0.10)
|
0.10 (-0.05, 0.24)
|
BMI
(z-score)
|
30
|
-0.19 (-0.35, -0.03)
|
-0.23 (-0.40, -0.06)
|
-0.11 (-0.28, 0.07)
|
-0.04 (-0.25, 0.17)
|
0.08 (-0.15, 0.32)
|
aBMD
|
TBLH
(g/cm2)
|
30
|
0.631 (0.612, 0.650)
|
0.633 (0.624, 0.641)
|
0.633 (0.614, 0.651)
|
0.002 (-0.018, 0.022)
|
0.002 (-0.033, 0.036)
|
TBLH
(z-score)
|
30
|
-0.97 (-1.06, -0.89)
|
-0.94 (-1.02, -0.85)
|
-0.93 (-1.01, -0.84)
|
0.03 (-0.03, 0.11)
|
0.05 (-0.05, 0.14)
|
Spine L1-L4 (g/cm2)
|
29
|
0.685 (0.662, 0.707)
|
0.691 (0.682, 0.701)
|
0.699 (0.676, 0.723)
|
0.007 (-0.018, 0.031)
|
0.015 (-0.029, 0.058)
|
Spine L1-L4 (z-score)
|
27
|
-0.71 (-0.79, -0.63)
|
-0.63 (-0.72, -0.54)
|
-0.57 (-0.66, -0.48)
|
0.08 (-0.02, 0.18)
|
0.14 (0.03, 0.26)*
|
BMC
|
TBLH (g)
|
30
|
741 (714, 768)
|
763 (749, 776)
|
772 (745, 800)
|
21 (-7, 50)
|
31 (-18, 81)
|
Spine L1-L4 (g)
|
29
|
21.7 (20.4, 23.0)
|
23.2 (22.6, 23.7)
|
23.9 (22.5, 25.2)
|
1.5 (0.0, 2.9)*
|
2.2 (-0.3, 4.7)
|
Lean mass
|
Legs (kg)
|
30
|
5.9 (5.6, 6.3)
|
6.0 (5.9, 6.2)
|
6.1 (5.8, 6.5)
|
0.1 (-0.3, 0.4)
|
0.2 (-0.5, 0.8)
|
Fat mass
|
Total (%)
|
30
|
28.6 (28.3, 28.9)
|
28.8 (28.5, 29.2)
|
28.8 (28.4, 29.2)
|
0.2 (-0.3, 0.7)
|
0.2 (-0.3, 0.7)
|
12VT, assessment after 12 weeks of vibration therapy; 20VT, assessment after 20 weeks of vibration therapy; aBMD, areal bone mineral density; BMC, bone mineral content; BMI, body mass index; TBLH, total body less head.
Data at each assessment are the adjusted means and 95% confidence intervals (CI), while differences between assessments are the adjusted mean differences and 95% CI; all values were derived from linear mixed models based on repeated measures including the participant's GMFCS level, randomisation group (20 Hz / 25Hz), the baseline value of the outcome, as well as the number of days elapsed from baseline (except for anthropometric outcomes).
n is the number of participants at baseline; the number of participants who completed a given assessment is provided in Additional file 1.
Statistically significant differences (at p<0.05) between assessments are shown in bold, where *p < 0.05 indicates a pairwise difference compared to the Control period.
3.6 Muscle function
Tests on the Leonardo Mechanography plate showed that VT improved some parameters of muscle function (Table 5). The maximum velocity rise in the chair-rise test increased by 0.14 m/s (≈17%) after 20VT (p=0.021; Table 5). The maximum force in the single two-leg jump test increased by 0.30 N/kg after 12VT (p=0.0005), with eight extra weeks of VT leading to an additional 0.15 N/kg improvement (95% CI 0.02, 0.28 N/kg; p=0.022), so that force increased by 0.46 N/kg after 20VT (p=0.022) (Table 5). However, there were no observed changes in double-leg balance (Table 5), or in muscle strength measured by HHD (Additional file 2).
Table 5. Muscle function and physical activity outcomes
Parameters (units)
|
n
|
Control
|
12VT
|
20VT
|
12VT vs
Control
|
20VT vs
Control
|
CRT
|
Forcemax (N/kg)
|
8
|
0.52 (-0.21, 1.25)
|
0.71 (0.49, 0.93)
|
0.82 (0.04, 1.61)
|
0.19 (-0.71, 1.09)
|
0.30 (-1.21,1.82)
|
Velocity rise max (m/s)
|
8
|
0.83 (0.76, 0.90)
|
0.89 (0.80, 0.98)
|
0.97 (0.89, 1.06)
|
0.06 (-0.06, 0.18)
|
0.14 (0.03, 0.26)*
|
STLJT
|
Forcemax , (N/kg)
|
19
|
0.43 (0.29, 0.57)
|
0.74 (0.67, 0.80)
|
0.89 (0.75, 1.03)
|
0.30 (0.14, 0.47)***
|
0.46 (0.19, 0.73)**
|
Velocitymax (m/s)
|
19
|
1.46 (1.30, 1.63)
|
1.42 (1.26, 1.59)
|
1.46 (1.29, 1.62)
|
-0.04 (-0.13, 0.05)
|
-0.01 (-0.13, 0.12)
|
Double leg balance
|
Elliptical area (cm2)
|
26
|
1.38 (1.07, 1.79)
|
1.34 (1.02, 1.75)
|
1.18 (0.90, 1.56)
|
0.97 (0.74, 1.27)
|
0.86 (0.61, 1.20)
|
12VT, assessment after 12 weeks of vibration therapy; 20VT, assessment after 20 weeks of vibration therapy; CRT, chair rising test; STLJT, Single two-leg jump test.
Data at each assessment are the adjusted means and 95% confidence intervals (CI), while differences between assessments are the adjusted mean differences (aMD) and 95% CI; all values were derived from linear mixed models based on repeated measures including the participant's GMFCS level, randomisation group (20 Hz / 25Hz), and the baseline value of the outcome, as well as the number of days elapsed from baseline for CRT Forcemax. Note the values for double-leg balance were log-transformed for analyses and back-transformed for reporting in this table, so the aMD represents the proportional difference compared to the Control period.
n is the number of participants at baseline; the number of participants who completed a given assessment is provided in Additional file 1.
Statistically significant differences (at p<0.05) between assessments are shown in bold; *p < 0.05, **p < 0.01, and *** p < 0.001 for pairwise differences compared to the Control period.
3Health-related quality of life outcomes
CPQOL was filled out by parents or caregivers of 25 participants at all assessments (Additional file 1). After 20VT, there was a 7-point improvement in general health scores (95% CI 2, 12; p=0.0095), but there were no other observed effects of VT on health-related quality of life (Additional file 3).
-
Compliance and side effects
Overall, participants had a high level of compliance with the prescribed VT protocol, after both 12 weeks [median=99% (Q1=92%, Q3=100%)] and 20 weeks [median=99% (Q1=89%, Q3=100%)]. Only 5 out of 32 (16%) and 6 out of 30 (20%) participants who attended the 12-week and 20-week assessments, respectively, completed less than 80% of prescribed VT sessions. The main reported reasons for missing sessions were lack of time and being unwell. Parents/caregivers of participants who performed sessions at home highlighted a higher need for child encouragement after 12 weeks of VT, as participants appeared to be gradually losing interest in performing VT at home. Notably, 50% of the participants were undergoing VT during COVID-19 lockdowns, which might have also changed their level of routine activities, and impacted their mental health and behaviour. In turn, these might have negatively affected their motivation to continue VT (31).
VT was well-tolerated with no severe adverse events reported. Nine participants (30%) reported occasional mild itchiness in the calf and ankle areas during VT sessions, which resolved within approximately 30 seconds to 2 minutes after cessation of the VT session. Two of these participants also complained about occasional warmth and redness of the skin on ankle area quickly resolved after the VT session.