Next Article in Journal
Reasons and Determinants of BoNT-A Treatment Discontinuation in Patients Living with Spasticity: A 10-Year Retrospective Analysis
Next Article in Special Issue
Botulinum Toxin Intervention in Cerebral Palsy-Induced Spasticity Management: Projected and Contradictory Effects on Skeletal Muscles
Previous Article in Journal
Does Ultrasound Guidance Improve the Effectiveness of Neurotoxin Injections in Patients with Cervical Dystonia? (A Prospective, Partially-Blinded, Clinical Study)
Previous Article in Special Issue
Reduced Cross-Sectional Muscle Growth Six Months after Botulinum Toxin Type-A Injection in Children with Spastic Cerebral Palsy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Short-Term Impact of Botulinum Neurotoxin-A on Muscle Morphology and Gait in Children with Spastic Cerebral Palsy

1
Department of Rehabilitation Sciences, KU Leuven, 3001 Leuven, Belgium
2
Department of Rehabilitation Sciences, Ghent University, 9000 Ghent, Belgium
3
Inkendaal Rehabilitation Hospital, 1602 Vlezenbeek, Belgium
4
Rehabilitation Research Group, Vrije Universiteit Brussel, 1090 Brussels, Belgium
5
Department of Development and Regeneration, KU Leuven, 3000 Leuven, Belgium
6
Department of Orthopedic Surgery, University Hospitals Leuven, 3000 Leuven, Belgium
7
Clinical Motion Analysis Laboratory, University Hospitals Leuven, Pellenberg, 3212 Leuven, Belgium
*
Author to whom correspondence should be addressed.
Toxins 2022, 14(10), 676; https://doi.org/10.3390/toxins14100676
Submission received: 26 August 2022 / Revised: 19 September 2022 / Accepted: 21 September 2022 / Published: 29 September 2022
(This article belongs to the Special Issue Botulinum Toxin in the Management of Children with Cerebral Palsy)

Abstract

:
Children with spastic cerebral palsy (SCP) are often treated with intramuscular Botulinum Neurotoxin type-A (BoNT-A). Recent studies demonstrated BoNT-A-induced muscle atrophy and variable effects on gait pathology. This group-matched controlled study in children with SCP compared changes in muscle morphology 8–10 weeks post-BoNT-A treatment (n = 25, median age 6.4 years, GMFCS level I/II/III (14/9/2)) to morphological changes of an untreated control group (n = 20, median age 7.6 years, GMFCS level I/II/III (14/5/1)). Additionally, the effects on gait and spasticity were assessed in all treated children and a subgroup (n = 14), respectively. BoNT-A treatment was applied following an established integrated approach. Gastrocnemius and semitendinosus volume and echogenicity intensity were assessed by 3D-freehand ultrasound, spasticity was quantified through electromyography during passive muscle stretches at different velocities. Ankle and knee kinematics were evaluated by 3D-gait analysis. Medial gastrocnemius (p = 0.018, −5.2%) and semitendinosus muscle volume (p = 0.030, −16.2%) reduced post-BoNT-A, but not in the untreated control group, while echogenicity intensity did not change. Spasticity reduced and ankle gait kinematics significantly improved, combined with limited effects on knee kinematics. This study demonstrated that BoNT-A reduces spasticity and partly improves pathological gait but reduces muscle volume 8–10 weeks post-injections. Close post-BoNT-A follow-up and well-considered treatment selection is advised before BoNT-A application in SCP.
Key Contribution: Botulinum Neurotoxin type-A successfully reduced lower limb spasticity and improved ankle kinematics, but induced muscle volume reductions in the medial gastrocnemius and semitendinosus.

1. Introduction

Spasticity occurs in approximately 60–80% of the children with cerebral palsy [1,2]. Therefore, spastic cerebral palsy (SCP) is the most common subtype. Spasticity is thought to be a result of imbalanced signals from the central nervous system resulting in impaired gross motor function, such as pathological gait, as well as in discomfort and pain [3,4]. Graham et al. described four stages of musculoskeletal pathology in SCP [5].
Botulinum Neurotoxin type-A (BoNT-A) is widely used during the first stage of musculoskeletal SCP pathology, which is characterized by hypertonia, as an effective treatment in reducing the hyperactive stretch reflex by temporarily blocking the release of Acetylcholine at the neuromuscular junction [6,7,8]. By reducing the hyperactive stretch reflex, BoNT-A aims to increase range of motion (ROM), reduce contracture development and improve function and gait [8,9,10,11]. Additionally, patient comfort, physical activity and social participation may improve [12]. Yet, the etiology of contractures is very complicated and not fully understood, and spasticity is definitely not the only contributing factor [13,14,15,16,17]. Still, BoNT-A injections remain the first line of treatment to prevent the child from developing to the second stage of contractures, which can only be surgically corrected [5].
As BoNT-A induces a local, temporary paralysis in the muscle, recent concerns have been raised regarding the impact of repeated BoNT-A injections on muscle morphology in young children with SCP [18,19,20]. A review by Multani et al. summarized evidence of muscle atrophy, fatty infiltration, loss of contractile tissue and increased connective tissue after BoNT-A treatment [18]. However, it should be noted that the highest levels of muscle atrophy and quality reduction with only partial recovery of muscle morphology and function at 12 months after BoNT-A injection have been reported in animal studies or studies including healthy volunteers [21,22,23]. Yet, knowledge on BoNT-A-impact on muscle morphology in a pediatric SCP population is important, specifically since muscle growth is already impaired in SCP from an early age on [24]. Moreover, multiple studies showed that muscle volume and muscle quality can be negatively influenced by BoNT-A [25,26,27,28,29,30]. It has also been shown that SCP muscles have a reduced number of satellite cells, suggesting that the regenerative capacity of the muscle might not be sufficient to ensure efficient muscle recovery, for example to overcome the potential paralysis-induced atrophy following BoNT-A injections [17,31,32].
Furthermore, early retrospective studies in young children with SCP who received one to five BoNT-A treatment sessions following an integrated approach (i.e., BoNT-A injections combined with stretching casts if relevant, intensive physical therapy, and orthotic management) showed beneficial effects on gait pathology at five to ten years of age (defined by three-dimensional gait analysis in 60 children with SCP) [33]. Moreover, BoNT-A treatment resulted in the reduced frequency of surgical interventions and postponed the need for surgery (defined in 424 children with SCP) [34]. The positive effects of BoNT-A, following this integrated approach, on gross motor function, including gait, have been confirmed by more recent studies [35,36]. A cohort study of 28 children with SCP showed that children who received repeated integrated BoNT-A treatment at a young age, improved or maintained their gross motor function at an older age [35]. However, these benefits have been recently questioned, in other studies that showed limited reductions in pathological gait [37,38,39]. BoNT-A-induced muscle weakness has been suggested as the potential cause for the lack of functional improvement, whereas SCP is already characterized by reduced muscle strength [19].
Despite the early reported positive effects of BoNT-A on gait and the potential to postpone more invasive treatments, the reported negative effects are concerning. Due to its temporary effect, BoNT-A injections have to be repeated. Different research centers [40,41,42], as well as previous BoNT-A treatment recommendations, reported repeated injections every six months or even less, while other centers reported intervals between repeated BoNT-A sessions of more than one year. Interestingly, the most recent recommendations advise an interval of at least one year [18,26,43]. Previous animal [22] and human [21,26] studies indeed showed that muscles are not recovered after six months. Multiple experts have asserted the need to revise treatment protocols as irreversible muscle atrophy may occur and the detrimental effects might become more severe after multiple injections [44]. Clearly, there is an urgent need for a critical reflection on the benefits and harms of BoNT-A in young children with SCP, which should be based on a comprehensive assessment of a multidimensional set of outcomes [18,23]. These outcomes should not only include clinical measures but also measures of the muscle morphology and objective measures of spasticity, as well as functional outcomes, such as gait analysis. A quantitative objective approach is important, as previous studies suggested that the benefits of BoNT-A can be overestimated when using subjective outcomes [18,45,46]. Even though the number of studies investigating the impact of BoNT-A in SCP has been rising in the past decade, previous research had some limitations. Small sample sizes, without control groups, were often included, subjective outcomes or isolated gait parameters were used and the focus was limited to only one muscle.
Therefore, the primary aim of the current group-matched controlled study was to investigate the impact of BoNT-A injections administered following an integrated therapeutic approach (combined with stretching casts if relevant, intensive physical therapy, and orthotic management) on the muscle morphology of two lower-limb muscles in children with SCP, aged between 3–11 years. The secondary aim was to confirm the impact of BoNT-A on spasticity and gait. The objective measures of the muscle morphology of the medial gastrocnemius (MG) and the distal compartment of the semitendinosus (ST) by three-dimensional freehand ultrasound (3DfUS), instrumented spasticity assessments and three-dimensional gait analysis (3DGA) were used to quantify the effects of BoNT-A. It was hypothesized that BoNT-A would reduce the muscle volume and negatively affect muscle quality. Furthermore, it was expected that multilevel integrated BoNT-A treatment would reduce the hyperactive stretch reflex and significantly improve gait kinematics of the ankle and knee in the sagittal plane.

2. Results

2.1. Participants

For the current study, the participants of two sub-studies on the effect of integrated BoNT-A treatment on muscle morphology and gait were pooled. This resulted in a total sample of 45 children (29 boys, 16 girls), with a median age of 7.1 years (5.7–9.1), who all had routine physical therapy and the use of orthoses as usual care. Twenty-five children, with an indication for BoNT-A injections in the ST and/or MG were included in the intervention group, while 20 children without an indication for BoNT-A were included in the untreated control group. The indication for BoNT-A treatment was set a by a multidisciplinary team based on clinical examination including assessments of ROM, spasticity, selectivity and strength, and the longitudinal evaluations of the pathological gait. Participants’ baseline characteristics and BoNT-A treatment details are summarized in Table 1.
All children in the intervention group underwent an assessment of muscle morphology by 3DfUS and received 3DGA to assess the level of gait pathology before and 8–10 weeks after the BoNT-A intervention. For one of the two sub-studies, the participants (n = 14, median age of 7.0 years (5.3–9.6), 6 boys and 8 girls, bilateral/unilateral 7/7, GMFCS level I/II/III 9/4/1 received an additional instrumented spasticity assessment at the same time points. The median interval between the baseline assessment (3DfUS + instrumented spasticity assessment) and integrated BoNT-A treatment was 0.2 (0.0–0.6) months. Due to practical reasons and the clinical follow-up, the pre-3DGA could not always be combined with the baseline assessments (3DfUS + instrumented spasticity assessment) required for this study. Therefore, the median interval between the pre-3DGA and the BoNT-A treatment was 1.1 (0.5–2.6) months. The median interval between the BoNT-A treatment and post-treatment evaluation (3DfUS + instrumented spasticity assessment + 3DGA) was 2.1 (1.9–2.4) months.
For the primary aim, 20 children with SCP, receiving no intervention and continuing usual care, were included in the control group. These children underwent 3DfUS assessments at baseline and after a period of 8–10 weeks, with a median interval of 2.7 (1.9–2.8) months but no 3DGA nor instrumented spasticity assessment.
The duration of the interval between BoNT-A treatment and the post-assessment in the intervention group and the interval between the baseline and post-assessment of the control group was not significantly different (p = 0.740).
As part of the integrated approach, 20 children received lower-leg casts immediately following the BoNT-A injections. These casts were applied for at least 10 days. Of the 20 children, 15 received additional removable upper-leg casts (as part of the integrated approach). Two children had full leg casts for at least ten days. Three children did not receive casting.
At baseline, there was no difference in age (p = 0.144) between the intervention and control group. Body mass (p = 0.010) and height (p = 0.021) were significantly different. Therefore, the muscle morphology parameters were normalized to body dimension (body mass and height).

2.2. Muscle Morphology

At baseline, no significant differences in normalized muscle volume (normalized to the product of body mass × height) or echogenicity intensity of the MG and ST were found between the intervention and control group.
The normalized muscle volume and echogenicity intensity before and after BoNT-A treatment are presented in Figure 1, Figure 2 and in the supplementary data, Table S1. The normalized muscle volume of both the MG (p = 0.018) and the distal compartment of the ST (p = 0.030) showed a significant reduction of 5.2% and 16.2% (supplementary data, Table S1) after BoNT-A treatment, respectively (Figure 1). No differences in the echogenicity intensity of the MG (p = 0.407) and ST (p = 0.145) were found within the intervention group (Figure 2).
No changes in muscle volume or echogenicity intensity were reported in the control group (supplementary data, Table S1). The changes in MG (p = 0.118) and ST (p = 0.069) muscle volume post-BoNT-A were not significantly different from the evolution in muscle volume in the control group.
Absolute muscle growth rates are presented in Table 2. Muscle growth rates calculated over the period between the two assessments were not significantly different between the intervention and control group for both muscles (MG p = 0.346, ST p = 0.226).
All muscle morphology parameters, including differences’ scores and between-group comparisons are summarized in Table S1 of the supplementary data.

2.3. Gait

The kinematic and spatial-temporal parameters are summarized in Table 3.
The ankle angle at initial contact (p = 0.002), the maximum dorsiflexion angle during both stance (p = 0.004) and swing phase (p < 0.001) significantly improved (while accounting for multiple testing, Bonferroni correction). The maximum dorsiflexion angle during swing exceeded the previously reported standard error of measurement (SEM) of 4.7 degrees (supplementary data, Table S1) [47]. No significant improvements were found in knee kinematics in the sagittal plane following BoNT-A treatment. However, there was a tendency towards more knee flexion during stance (p = 0.065), whereas knee flexion at initial contact seemed slightly reduced (p = 0.048).
The gait profile score (GPS) showed a significant decrease (p = 0.005) after BoNT-A treatment, indicating less gait pathology. The median improvement of 1.76 degrees exceeded the minimally clinically important difference (MCID) of 1.6 degrees [48]. Cadence, walking velocity and stride length remained unchanged.
The statistical non-parametric mapping (SnPM) analyses of the full kinematic ankle and knee waveforms revealed significant improvements in ankle kinematics, located in two clusters. One large cluster (p < 0.001) reached from 19–100% of the gait cycle and exceeded the SEM (ankle dorsi-plantar flexion, 3.3 degrees) [49]. Another small cluster (p = 0.0165, 0–1.99%) at initial contact did not exceed 3% of the gait cycle and was therefore not considered clinically relevant [50]. The knee kinematics showed a small cluster at the end of swing from 92–96% of the gait cycle, exceeding 3% of the gait cycle and the SEM (knee flexion-extension, 2.4 degrees) [49].
An overlay of the pre- and post-treatment kinematic waveforms for the ankle and knee is presented in Figure 3.

2.4. Instrumented Spasticity Assessment

The instrumented spasticity assessment performed in 14 children of the intervention group showed a significant decrease in muscle activation of the MG (p = 0.019) and medial hamstrings (including the ST) (p = 0.041) during the high-velocity passive stretches. All results are summarized in Table 4.
The difference in EMG signal between the low- and high-velocity stretches, a secondary outcome, showed a borderline significant decrease after integrated BoNT-A treatment in the ST (p = 0.060), exceeding the SEM (Table S4, supplementary data) [51].

2.5. Clinical Examination and Range of Motion

The results of the clinical examination are summarized in the supplementary data, Table S3. Maximal dorsiflexion (p = 0.302) and popliteal angle (unilateral) (p = 0.564) remained unchanged after BoNT-A injections.

3. Discussion

3.1. Summary of Findings

The primary aim of this study was to investigate the effect of integrated BoNT-A treatment on muscle morphology of the MG and ST. The second aim was to provide additional evidence regarding the effects on spasticity and gait. To achieve the study aims, an integrated set of objective assessments was used. We hypothesized that BoNT-A injections would result in reduced spasticity and improved gait but would induce muscle atrophy. The current results revealed that integrated BoNT-A treatment successfully reduced the hyperactive stretch reflex but also reduced muscle volume of the MG and distal compartment of the ST 8–10 weeks after BoNT-A injections. No changes in muscle quality were observed. Ankle kinematics during gait improved significantly, whereas effects on knee kinematics were limited. The study hypotheses were thus partly confirmed.

3.2. Comparison with Literature

3.2.1. The Impact of BoNT-A on Muscle Volume

The current data showed a significant reduction (median change) in the normalized muscle volume of both the MG and ST of approximately 5.2% and 16.2%, 8–10 weeks after BoNT-A injections, which were applied as part of an integrated approach, respectively. Muscle growth rates, calculated over the period between the two assessments, were not significantly different from the control group. Several studies have reported similar rates of muscle atrophy after BoNT-A injections [21,22,23,25,26,27,28,29,30,52].
Alexander et al. reported reductions in MG volume, as absolute value and normalized to bone length at four (5.9%), 13 (9.4%) and 25 (6.8%) weeks, following the first single-level BoNT-a treatment to the MG. The greatest reduction, 13 weeks post-BoNT-A (9.4%), was slightly higher compared to the current data (5.2% median change, 8–10 weeks post-BoNT-A). Besides the MG atrophy, the (untreated) soleus muscle became hypertrophic in the sample of Alexander et al., resulting in an unchanged total volume of the plantar flexors. They assumed that the hypertrophy potentially developed to maintain plantar flexor strength and function. The reduction in volume of the MG combined with the hypertrophy of the soleus was also described by Williams et al. after five weeks post-BoNT-A [29]. They found an average reduction of 4.5% in the gastrocnemius (medial and lateral head). Five of those children received injections in the medial hamstrings as well, resulting in a volume reduction of 5.9%. All children in their study had a history of previous injections and received a BoNT-A injection in the MG as part of a multi-level treatment.
A recent study by De Beukelaer et al. was the first to quantify the BoNT-A effect on muscle morphology in an intervention group compared to a control group of children who did not receive BoNT-A treatment. They showed a reduction in anatomical cross-sectional area six months after the first BoNT-A injection, which was not observed in the control group of children, who remained BoNT-A naive during this period [26]. They did not find any significant between group differences in muscle volume following BoNT-A.
The decrease in ST volume reported in the current study is higher compared to the MG volume reduction and higher in comparison to previous studies of the distal lower limb muscles. The level of atrophy in the ST is more in line with a previous study by Van Campenhout et al. showing approximately 20.5% (79.5% of the initial volume) reduction in volume of the psoas two months after BoNT-A treatment [30], suggesting that proximal lower limb muscles might be more affected by BoNT-A. However, these findings should be confirmed in future studies including both distal and proximal lower limb muscles.
The current findings are in contrast with the results reported by Barber et al. [40], who showed an increase in absolute MG volume, reflecting an unchanged MG growth rate 12 months after BoNT-A. Williams et al. combined BoNT-A with strength training and also showed a significant increase in muscle volume of the hamstrings and plantar flexors (next to other lower limb muscles) six months post treatment [53]. However, it should be noted that improvements in absolute muscle volume observed in the latter studies were most likely a result of natural growth.
The differences in findings between the previously mentioned studies in SCP can be explained by the differences in study population, sample size, used methods, duration of follow-up, treatment protocol and content of usual care.
Besides the immediate impact of BoNT-A on muscle volume as described in the current prospective study, there are indications in the literature that previous BoNT-A injections in children who are not BoNT-A naïve may influence muscle morphology. A retrospective study by Schless et al. compared muscle volumes between children naïve to BoNT-A and children with a history of at least three previous BoNT-A treatments (median of five previous interventions) and showed that the absolute muscle volume was roughly 29% smaller in the children with a history of previous BoNT-A injections [25]. This may also further explain differences between the different studies. Longitudinal studies are required to further evaluate the impact of BoNT-A on muscle morphology and monitor the muscle-specific recovery profiles.
It should also be noted that the described alterations in muscle volume of children with SCP are less extreme compared to animal studies and studies with healthy adults [21,22,52]. This is potentially due to the aggressive treatment protocols that were applied in those studies.

3.2.2. The Impact of BoNT-A on Muscle Morphology—Echogenicity Intensity

In the current study, no changes in echogenicity intensity were found, suggesting that integrated BoNT-A treatment does not negatively impact muscle quality in the short-term period. This is in contrast with a previous retrospective study by our research group showing that the MG echogenicity intensity of children with a history of BoNT-A was significantly higher compared to children who remained BoNT-A naïve, pointing towards long-lasting alterations in the muscles as a result of BoNT-A [25]. A more recent study by De Beukelaer et al. showed a tendency towards higher echogenicity intensity values in children with a history of previous injections compared to BoNT-A naïve children at baseline [26]. The fact that no alterations in echogenicity intensity were found in the current study could be explained by the fact that only the short-term effects of a single BoNT-A treatment session or short follow-up period were investigated. The higher echogenicity intensity values, as described in previous studies, suggest a reduced muscle quality, which may be caused by increased connective tissue, loss of contractile elements, fibrosis and fatty infiltration [13,14,26,27]. Indeed, multiple animal studies reported fibrosis or collagen deposition, increased lipid accumulation and loss of contractile elements after BoNT-A injections [22,44,54,55] and echogenicity intensity has been linked to these factors, contributing to declined muscle quality and contracture development [56,57,58]. Moreover, a causal relation between muscle fibrosis and non-responsiveness to BoNT-A has been suggested [43,59]. Yet, since echogenicity intensity represents an indirect measure of muscle quality and cannot decompose the contribution of different components, such as increased fat infiltration and connective tissue, the earlier retrospective results in humans should be interpreted carefully. The validity of echogenicity intensity is still not fully understood. It has been related to muscle composition in several pathologies; however, validation studies in SCP are lacking [58,60,61].

3.2.3. The Impact of BoNT-A on Muscle/Tendon Length, Range of Motion and Spasticity

The secondary morphological outcomes (Table S1, supplementary data) showed an increase in muscle belly length (p = 0.009) and total muscle-tendon unit length (p = 0.011) of the MG after BoNT-A injections, while these length data did not change for the ST. However, the maximum ankle dorsiflexion angle defined by goniometry during the clinical examination did not change post-BoNT-A. A previous investigation showed a similar increase in MG muscle belly length combined with unchanged maximum ankle dorsiflexion already at two weeks post-BoNT-A injection [62]. Additionally, the resting ankle angle during the 3DfUS shifted towards increased dorsiflexion (median increase of 10 (12.5–0) degrees dorsiflexion, p = 0.003) after integrated BoNT-A treatment. This could explain the increase in muscle belly and total muscle-tendon unit length. Additionally, a previous observational study by Malaiya et al. did not find an association between MG muscle belly length and maximum dorsiflexion angle in a cohort of 16 children with unilateral SCP [63].
For the instrumented spasticity assessments, the average EMG signal while stretching the muscle at high-velocity was significantly reduced post-BoNT-A compared to the baseline condition. This resulted in a (tendency) towards a reduced difference in EMG between high and low stretch velocity (which exceeded the SEM). These results confirmed previous findings that BoNT-A injections in the MG and ST successfully reduces the hyperactive stretch reflex [62,64,65].

3.2.4. The Impact of BoNT-A on Gait

In contrast to the studies by Hastings-Ison, the current study found significant improvements in the overall gait pathology defined by the GPS [37,38]. The SnPM analyses of the full kinematic waveforms showed one large relevant ankle cluster, as well as one small cluster in the knee during swing. These improvements were confirmed by beneficial changes in gait features of the ankle (significant for ankle at initial contact and maximal dorsiflexion in stance and swing) and knee (tendencies for knee angle at initial contact and knee ROM). These findings are in line with a previous study [66]. Hastings-Ison et al. did not find any alteration in maximum ankle dorsiflexion or knee kinematics in the total group. However, they showed improved ankle dorsiflexion and increased knee flexion during mid stance in bilaterally involved children but not in unilaterally involved participants [37]. The current data were not separately analyzed for children with uni- and bilateral involvement due to lack of power. We previously performed a systematic review on the BoNT-A effect on gait and compared the BoNT-A effect on gait features with the effect defined by statistical parametric mapping in a group of 51 children with SCP who received multilevel BoNT-A treatment [66]. The latter study showed significant improvements at the ankle as well as at the knee joint, and these beneficial effects were also more pronounced at the ankle compared to the knee joint. Overall, we conclude that BoNT-A improved the gait of the children with SCP at 8–10 weeks post-BoNT-A.

3.3. Frequency of BoNT-A Administration

As highlighted by the study of Multani and colleagues, the BoNT-A injection frequency is an important matter of debate [18,67]. Studies by Kanovsky and Hastings-Ison showed that BoNT-A treatment once per year was effective with fewer adverse events compared to multiple injections per year [41,42]. These findings were supported by animal studies showing more muscle weakness and atrophy, increased loss of contractile tissue and exacerbated fibrosis after more frequent injections [23,44,54]. As the muscle might need more time to recover, even if the clinical effect of tone reduction is dissolved, the duration between injections should be based on assessments of the muscle, preferably at both the microscopic and macroscopic level [23,44,52].

3.4. Combined Outcome Assessments and Clinical Implication

This is the first study that combined the assessments of muscle morphology and gait, in order to critically review the benefits and potential harm of integrated BoNT-A treatment. The observed clinical benefits included reduced spasticity of the plantar flexors and knee flexors and improved ankle motion during gait but had rather limited effects on knee kinematics. On the other hand, the muscle volume was reduced 8–10 weeks post-BoNT-A treatment. These pros and cons should be weighted before the selection of a (repeated) BoNT-A treatment for young children with SCP. The recovery profile of the individual patient should be objectively monitored. The current results revealed muscle atrophy in the short-term period following BoNT-A injections and other studies highlighted that muscle recovery is still ongoing six months post-BoNT-A [26,27,40]. These findings suggest that muscle atrophy lasts longer than the clinical effects, such as tone reduction [23,52]. An animal study by Minamoto et al. hypothesized that the effect and potential negative effects of BoNT-A are much stronger if injected into a muscle which has not completely recovered [44]. The current evidence implies that BoNT-A injections should be provided with an interval of at least 12 months. As stated earlier, the first BoNT-A injections seem to be most distinctive in both the beneficial functional effects and the harmful effects on the muscle. Therefore, muscle recovery is especially crucial after the first injection and an interval of at least 12 months before the second injection is advised.

3.5. Future Perspectives

The observed muscle atrophy, after approximately two months post-BoNT-A in the current study and six months post-BoNT-A in a previous study [26], represents the rather short-term effect of BoNT-A on muscle morphology. Hence, no conclusion can be drawn on the long-term effects. More research is needed to define the time course of muscle recovery. Future studies should investigate the long-term effects and include the microscopic and molecular analyses of intrinsic muscle properties and functional evaluations in order to critically appraise the benefit–harm balance of BoNT-A treatment. Thereby, more research on dosage dependent effects, the frequency of BoNT-A injections and treatment history is relevant. In order to gain more insights into the patient- and or even muscle-specific response to BoNT-A and in order to properly identify responders and non-responders, the objective measures of muscle spasticity as well as morphology and function, such as gait [37,39,64,68], should be combined in future studies and implemented in clinical practice. These findings may improve efficient treatment planning.
Most studies, including the current study, focused on the effects of BoNT-A in ambulatory children. Therefore, future studies should also include participants with less functional capacity in order to evaluate the benefits and potential harms in non-ambulatory children.

3.6. Limitations

Some limitations have to be acknowledged. This study included a non-randomized control group. As BoNT-A has been proven to be effective in reducing spasticity and time is running out for growing children with spasticity, as this may lead to impaired musculoskeletal development, it was considered unethical to temporarily withdraw treatment from patients. Therefore, randomization procedures were not possible. We minimalized bias by group-matching age and GMFCS level. The control group presented with slightly lower spasticity levels compared to the intervention group at baseline, which was not surprising since these children had no indication for BoNT-A treatment. No control data were available for gait and spasticity assessments. However, gait and spasticity parameters were evaluated using previously reported reliability values [47,48,51]. In the current study, BoNT-A was applied as part of an integrated approach, whereby BoNT-A injections were combined with casting, orthotic wear and intensive physical therapy. Hence, no conclusions can be drawn on the contribution of the isolated treatments.
Due to the deep origin of the ST, only the distal compartment could be visualized by 3DfUS, even though both compartments were targeted by BoNT-A [69]. Muscles were only assessed in resting length, therefore no conclusion on alterations in extensibility can be drawn.
The assessments were performed by three different assessors. However, because all assessors were well-trained (i.e., minimum of two years of experience), regular comparisons of performance to ensure standardization were made and because the pre- and post-assessments were always performed by the same assessor, we do not believe that this could have influenced our results.
This study was ongoing during the onset of the COVID-19 pandemic. As the majority of participants was included after the first COVID-19 outbreak and the application of an additional inclusion criterion stating that the physical therapy had to be stable for at least three months, we believe that the pandemic did not influence the outcomes of the current investigation.

4. Conclusions

This study including a comprehensive set of objective measures showed that BoNT-A injections, when administered as part of an integrated approach (i.e., BoNT-A injections combined with stretching casts if relevant, intensive physical therapy and orthotic management), results in effective tone reduction, with positive effects on ankle kinematics but limited effects on knee kinematics during gait, combined with muscle atrophy, 8 to 10 weeks post-BoNT-A. There is a need for evidence-based guidelines rather than guidelines based on clinical experience in order to further evaluate the impact of BoNT-A. Close follow-ups with an integrated set of objective measures are required to gain more insight into the patient-specific response and recovery profile. Furthermore, patient-specific goals and expectations should be carefully evaluated and contemplated before BoNT-A treatment is planned in young children with SCP.

5. Materials and Methods

This study included participants of two different projects. Both projects were approved by the local ethics committee of the University Hospitals Leuven, Belgium (s59945 TAMTA project, s62187 and s62645 3D-MMAP project) (Ethische Commissie Onderzoek KU Leuven/UZ Leuven). Written informed consent was obtained from the parents or legal guardian. Participants were included between May 2019 and June 2022.
This trial is registered at ClinicalTrials.gov, identifier: NCT05126693.

5.1. Participants

Children aged between 3–11 years, GMFCS levels I-III with SCP, scheduled for BoNT-A in the MG and/or ST, as part of a multilevel treatment, were included. Participants were recruited via two parallel prospective studies. Both studies investigated the effect of integrated BoNT-A treatment on muscle morphology and gait, while one study also evaluated the effects on spasticity (s59945 TAMTA project). The selection criteria for the participants were similar in both studies. The exclusion criteria were: any signs of dystonia or ataxia in the lower limbs, severe muscle weakness or poor selectivity in the plantar flexors or knee flexors and cognitive problems that would impede the measurements, intrathecal baclofen pump, previous muscle surgery in the lower limbs or selective dorsal rhizotomy, previous bony corrections in the lower limb within two years before assessments, previous BoNT-A treatment within six months before the baseline assessment or within six months before the pre-BoNT-A 3DGA. Additionally, the physical therapy frequency had to be consistent for at least three months before inclusion. The participants were measured twice, before and 8–10 weeks after BoNT-A treatment to evaluate the short-term effects. All children (of both studies) received 3DfUS assessments and a 3DGA with the clinical examination of the lower limbs. A subset of children (i.e., the participants of one of the two studies) underwent an additional instrumented spasticity assessment in order to evaluate the BoNT-A-induced tone reduction. During the follow-up, all children received usual care, including intensive physical therapy, orthotic management and oral tone reduction if needed.
An age- and functional severity-matched (group-matching based on age and GMFCS level) control group, meeting the same inclusion criteria as the intervention group, but not receiving any intervention apart from their usual care, was included. The participants of the control group were recruited via Cerebral Palsy Reference Center of the University Hospitals Leuven, Belgium and the Inkendaal Rehabilitation Hospital Vlezenbeek, Belgium. All children eligible for inclusion were measured (n = 30, Figure 4); however, the final data selection was based on accurate matching with the participants in the intervention group. The participants in the control group underwent 3DfUS measurements only, with an interval of 8–10 weeks between the two assessments. Participant characteristics are summarized in Table 1, showing no structural differences, as assessed by a Mann–Whitney U test, in the maximum dorsiflexion angle, popliteal angle and spasticity levels assessed by the modified Ashworth scale.
An overview of the recruited and eventually included subjects and data availability is summarized in Figure 4.

5.2. Procedure

The most affected leg, based on the most recent clinical examination (highest modified Ashworth score for the ankle and knee joints), was assessed. The measures of body mass, height, total leg length and lower leg length were taken. A standardized anamnesis was performed to request treatment history and usual care details.
Muscle selection for BoNT-A injection was primarily based on the most recent clinical examination and multidisciplinary evaluation of the 3DGA data. BoNT-A, Onabotulinum toxin A (Botox®, Allergan, Diegem, Belgium), injections were applied by a pediatric orthopedic surgeon under general (mask) anesthesia, guided by ultrasonography or palpation [30,69]. A dose of 100 units of BoNT-A was diluted in 5 mL of saline. The BoNT-A dilution was injected at multiple sites, in areas with the largest concentration of motor endplates in the muscles selected for treatment [69]. A maximum of 50 units was injected per site.
BoNT-A is suggested to be more effective when part of an integrated approach, i.e., combined with casting, orthotic wear and physical therapy [68]. Therefore, lower leg casts and removable upper leg casts were applied as indicated. The majority of children used day and night orthoses (Table 1) and all children received intensive physical therapy in the 8–10 weeks following BoNT-A treatment.
All additional study assessments (3DfUS, clinical examination, instrumented spasticity assessment and 3DGA) were performed by three trained assessors. The pre- and post-assessments were always performed by the same assessor.

5.3. Assessments

5.3.1. Muscle Morphology

The muscle morphology of the MG and the distal compartment of the ST muscle were assessed by 3DfUS [70]. This technique, combining B-mode ultrasonography with three-dimensional motion tracking, is proven to be valid and reliable in muscle morphology assessments in children with SCP in previous studies [71,72,73].
Ultrasound assessments were performed while the participant was placed in a prone position. The lower leg was supported by a triangular cushion to allow slight knee flexion (approximately 25 degrees) and avoid bi-articular stretch on the MG and ST. The ankle was in resting position without fixation. The joint angles were measured by goniometry.
Large amounts of acoustic transmission gel were used to assure smooth movement of the transducer and optimal contact between the transducer and skin. An additional curved gel pad, the Portico, was used to reduce muscle deformation during the acquisition [74]. Muscle images were recorded by Stradwin software (version 6.0; Mechanical Engineering, Cambridge University, Cambridge, UK) while moving the transducer from muscle origin to tendinous insertion. The acquisitions were repeated in case of bad quality images, movement or muscle contraction or when the bony landmarks or muscle borders were missed.
After data collection, the muscle images were processed by a trained assessor using the Stradwin software (version 6.0; Mechanical Engineering, Cambridge University, Cambridge, UK). Muscle volume (in millimeters) was estimated by the cubic planimetry technique based on segmentations that were manually drawn along the inside of the muscle borders in 5–10% of the transverse 2D muscle images [75]. Additionally, echogenicity intensity, as a proxy of muscle quality, was based on the average of the interpolated reconstruction of the muscle and expressed as a gray value on an 8-bit scale (values between 0–255), whereby higher values (whiter images) indicated reduced muscle quality [56,58]. Muscle belly lengths were calculated as the Euclidean distance between the muscle origin (MG: posterior surface of medial femoral condyle, ST: ischial tuberosity) and the muscle-tendon junction. The tendon length of the MG was calculated between the muscle–tendon junction and tendinous insertion. The MTU length of the MG was calculated by adding the muscle belly to the tendon length.
In order to compare children with different body dimensions, muscle volume was normalized to the product of body mass × height and muscle and tendon lengths were expressed as a percentage of total leg length.
The absolute muscle-specific growth rates for muscle volume during the period between the two assessments were calculated by the following formulas [26,40]:
growth   rate ,   muscle   volume = muscle   volume   mL   POST   muscle   volume   mL   PRE   age   months   POST   age   months   PRE
The primary outcomes included the normalized muscle volume and echogenicity intensity of the MG and distal compartment of the ST and the absolute muscle volume growth rates. The secondary outcomes, muscle, tendon and muscle-tendon unit lengths and normalized muscle volume growth rates are included in the supplementary data (Table S1).

5.3.2. Gait Analysis

Gait, overground barefoot walking at a self-selected speed, was assessed by 3DGA in a 10 to 12 camera (Oxford Metrics, Oxford, UK) gait lab (Vicon Motion Systems, Oxford, UK). Twelve reflective markers were placed on anatomical landmarks according the Vicon Plug-In-Gait lower limb model [76]. Ground reaction forces were measured with force plates (Advanced Mechanical Technology Inc., Watertown, MA, USA) that were embedded in the walkway. During the 3DGA, kinematic data (joint angles for ankle, knee and hip and the segmental orientation of the pelvis and foot), kinetic data (joint moments and power) and the surface EMG data of the medial gastrocnemius, soleus, tibialis anterior, rectus femoris, vastus lateralis, the medial and lateral hamstrings and gluteus medius were collected for each leg. The kinematic data of the most affected leg (same as for the 3DfUS and instrumented spasticity assessments) and a selection of spatial-temporal parameters were included in this current investigation but the kinetic and EMG data were not included in any of the analyses.
Gait cycles were identified and manually marked using the Vicon Nexus software (Oxford Metrics, Oxford, UK). All trials were exported and imported to a custom-made Matlab script (The MathWorks, Natick, MA, USA, version 2020b). All gait cycles were judged on quality. Trials were excluded in case of artefacts, gaps in data or indications of inaccurate marker placement. Furthermore, the knee varus–valgus angle was evaluated and trials with absolute knee varus–valgus ROM exceeding 15 degrees, a knee valgus angle exceeding −10 degrees, or trials that were not representative of the participants’ gait pattern were excluded [50]. As many good-quality steps as possible were used to create average trials for the ankle and knee joint kinematics.
The GPS, and specific kinematic gait parameters that were found to be sensitive to detect BoNT-A-induced changes, were extracted per joint [66]. Additionally, the spatial-temporal parameters cadence, walking velocity and stride length were calculated. Finally, the full kinematic waveforms were exported for pre–post comparison.
The primary outcomes included the kinematic waveforms, BoNT-A-sensitive kinematic gait parameters [66], GPS and spatial-temporal parameters. The secondary outcomes of elaborate kinematic parameters are summarized in the supplementary data (Table S2).

5.3.3. Spasticity Assessment

A subset of participants (n = 14) underwent an instrumented spasticity assessment in order to objectively evaluate the change in the hyperactive stretch reflex during passive stretch [77]. This method has been previously described by Bar-On and colleagues [77] and is proven sensitive enough to detect treatment-induced changes [62,64,65].
Four low-velocity (over 5 seconds) passive stretches of the plantar flexors and medial hamstrings were performed over the full ROM. Next, four full ROM high-velocity passive stretches were performed, whereby the stretch was performed as fast as possible. The interval between the stretch repetitions was at least 7 seconds to avoid post-activation depression. During the stretches, muscle activation was recorded by surface EMG (Zerowire Cometa, Milan, Italy)
The raw EMG signals were filtered by a 6th order zero-phase Butterworth bandpass filter (20–500 Hz). The root mean square envelope of the EMG signal (RMS–EMG) was computed after application of a low-pass 30 Hz 6th order zero-phase Butterworth filter.
All data analyses were carried out with MATLAB Software (The MathWorks, Natick, MA, USA, version 2020b). The quality of the individual stretches were checked using custom-made software. Stretch repetitions were excluded in case of inconsistent stretch velocities, poor surface EMG data, out-of-plane movement or indications of insufficient participant relaxation (constant EMG activity of the antagonist).
In order to quantify the alteration in muscle tone after BoNT-A treatment, the average RMS–EMG expressed in millivolts was calculated as the square root of the area underneath the RMS–EMG time curve, divided by the duration of the time interval, during the high-velocity stretch. The interval started 200 milliseconds prior to the time point where maximum stretch velocity was reached and ended at the time corresponding to 90% of the full ROM. This value was referred to as the muscle activation during a high-velocity stretch. All BoNT-A-induced changes were compared to previous reported SEM or BoNT-A-induced values [51,62,64,65,77].
The average RMS–EMG expressed in millivolts during a high-velocity stretch was used as the primary outcome to assess alterations in muscle tone. Additionally, the increase during the high-velocity stretch was calculated by subtracting the muscle activation during a low-velocity stretch, which was calculated following the same method. This value was called the change in muscle activation and was included as a secondary outcome. This secondary outcome together with the RMS–EMG during a low-velocity stretch are enclosed in the supplementary data (Table S4).

5.3.4. Clinical Examination

A standardized clinical examination including assessments of ROM, spasticity, selectivity and strength was performed during the pre- and post-BoNT-A assessments of the intervention group and during the pre-assessment of the control group. These data were used to describe and compare both groups at baseline. The maximal dorsiflexion angle and popliteal angle (unilateral) were assessed by goniometry. The level of spasticity was assessed by the modified Ashworth scale and modified Tardieu scale [78,79]. Muscle strength was assessed by manual muscle testing and selectivity was assessed by the Selective Control Assessment of the Lower Extremity [80,81]. The evaluations of strength and selectivity were not included in the current investigation.

5.4. Statistical Analysis

The normality of data distribution was checked by a Kolmogorov–Smirnov test (p < 0.05). Since most parameters were not normally distributed, non-parametric statistical tests were performed.
Participant characteristics, maximal dorsiflexion angle (measured with the knee extended), popliteal angle (unilateral), and normalized muscle morphology parameters of the intervention and control group were compared at baseline by a Mann–Whitney U test (Table 1).
A Wilcoxon signed rank test was used to evaluate (within) treatment effects in muscle morphology, instrumented spasticity assessment outcomes and spatial-temporal parameters of the 3DGAs (alpha = 0.05). A Bonferroni correction was applied for the gait kinematic parameters per joint (ankle: alpha 0.05/4 = 0.0125, knee: alpha 0.05/3 = 0.0167). Gait parameters were considered meaningful if statistically significant changes were noted that were exceeding the SEM or MCID [47,48].
Additionally, median difference scores (xpost−xpre) for muscle morphology parameters were calculated in both the intervention and control group and between-group differences were tested by a Mann–Whitney U test. Additionally, the changes were expressed as a percentage change of the baseline values for both groups, which are summarized in Table S1 of the supplementary data.
Paired t-tests of the ankle and knee kinematics over the entire gait cycle using SnPM analysis were performed to identify treatment-induced changes in these sagittal plane motions. The averaged ankle and knee kinematic waveforms, consisting of 101 points were imported to the SPM1d software (version 0.4, available for download at http://www.spm1d.org/, accessed on 15 August 2022) in Matlab (The MathWorks, Natick, MA, USA, version 2020b). Clusters were formed if the critical threshold, based on random field theory, was crossed [50,82,83] (alpha = 0.05). In case of significant clusters, the p-value, details on the location and the duration (percentages of the entire gait cycle) were reported. Clusters were considered meaningful when their duration exceeded 3% of the gait cycle and when the differences between the means of the pre- and post-BoNT-A joint kinematics exceeded the respective SEMs [49,50].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins14100676/s1, Table S1: Muscle morphology parameters of both the intervention and control group; Table S2: Results of the gait kinematics, gait profile and gait variable score and the spatial-temporal parameters, intervention group; Table S3: Results of the clinical examination in the intervention group; Table S4: Results of the instrumented spasticity assessment of a subgroup (n = 14) in the intervention group.

Author Contributions

Conceptualization, N.P., M.D., A.V.C., C.V.d.B., P.C. and K.D.; methodology, N.P., E.P. and K.D.; investigation, N.P., E.P., B.H., N.D.B. and L.S.; data curation, N.P., E.P. and N.D.B.; writing—original draft preparation, N.P.; writing—review and editing, N.P., E.P., B.H., N.D.B., L.S., M.D., C.V.d.B., P.C., H.F., A.V.C. and K.D.; visualization, N.P.; supervision M.D., C.V.d.B., P.C., H.F., A.V.C. and K.D.; project administration, N.P., E.P., B.H., N.D.B. and L.S.; funding acquisition, K.D. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Flemish Organization for Scientific Research Foundation, Belgium (FWO-TBM/TAMTA-T005416N and FWO-G0B4619N) and Internal KU Leuven grant, Belgium, C24/18/103, Internal funding of KU Leuven Biomedical Sciences Group: Fund for Translational Biomedical Research, 2019. E.P. was supported by the Flemish Organization for Scientific Research Foundation, Belgium (grant number T003116N).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of University Hospitals Leuven/KU Leuven (S59945, approved 18 July 2017; S62187, approved 8 February 2019 and S62645, approved 7 May 2019).

Informed Consent Statement

All subjects, parents or legal guardians gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author. The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Acknowledgments

The authors wish to thank the clinical team of the Clinical Motion Analysis Laboratory, University Hospitals Leuven, Pellenberg, Belgium, for their assistance in gait analyses data collection. We wish to thank Ester Huyghe for her help in ultrasound data processing. The research team is very grateful to all participants in this study.

Conflicts of Interest

None of the authors has any conflict of interest to disclose.

Abbreviations

3DfUS = three-dimensional freehand ultrasound; 3DGA = three-dimensional gait analysis; AFO = ankle foot orthoses; BoNT-A = Botulinum Neurotoxin type-A; cm = centimeters; COVID-19 = Coronavirus Disease 2019; Deg = degrees; DF = dorsiflexion; EMG = electromyography; GC = gait cycle; GMFCS = Gross motor function classification system; GPS = gait profile score; IC = initial contact; ISA = instrumented spasticity assessment; Kg = kilograms; m = meters; MAS = modified Ashworth score; max = maximum; MCID = minimally clinically important difference; MEH = medial hamstrings; MG = medial gastrocnemius; ml = milliliters; n = number; n.a. = not applicable; n.s. = not significantly different; PF = plantar flexors; RMS = root mean square; ROM = range of motion; SCP = spastic cerebral palsy; SEM = standard error of measurement; SnPM = statistical non-parametric mapping; ST = semitendinosus; µV = microvolts.

References

  1. Stavsky, M.; Mor, O.; Mastrolia, S.A.; Greenbaum, S.; Than, N.G.; Erez, O. Cerebral Palsy—Trends in Epidemiology and Recent Development in Prenatal Mechanisms of Disease, Treatment, and Prevention. Front. Pediatr. 2017, 5, 7. [Google Scholar] [CrossRef] [PubMed]
  2. Sellier, E.; Platt, M.J.; Andersen, G.L.; Krägeloh-Mann, I.; De La Cruz, J.; Cans, C.; Surveillance of Cerebral Palsy Network. Decreasing Prevalence in Cerebral Palsy: A Multi-Site European Population-Based Study, 1980 to 2003. Dev. Med. Child Neurol. 2016, 58, 85–92. [Google Scholar] [CrossRef] [PubMed]
  3. Nielsen, J.B.; Petersen, N.T.; Crone, C.; Sinkjaer, T. Stretch Reflex Regulation in Healthy Subjects and Patients with Spasticity. Neuromodul. Technol. Neural Interface 2005, 8, 49–57. [Google Scholar] [CrossRef] [PubMed]
  4. Lance, J.W. Spasticity: Disordered Motor Control; Feldman, R.G., Young, R.R., Koella, W.P., Eds.; Year Book Medical Publishers: Chicago, IL, USA, 1980. [Google Scholar]
  5. Graham, H.K.; Thomason, P.; Willoughby, K.; Hastings-Ison, T.; Van Stralen, R.; Dala-Ali, B.; Wong, P.; Rutz, E. Musculoskeletal Pathology in Cerebral Palsy: A Classification System and Reliability Study. Children 2021, 8, 252. [Google Scholar] [CrossRef]
  6. Cocco, A.; Albanese, A. Recent Developments in Clinical Trials of Botulinum Neurotoxins. Toxicon 2018, 147, 77–83. [Google Scholar] [CrossRef]
  7. Jankovic, J. Botulinum Toxin: State of the Art. Mov. Disord. 2017, 32, 1131–1138. [Google Scholar] [CrossRef]
  8. Blumetti, F.C.; Belloti, J.C.; Tamaoki, M.J.; Pinto, J.A. Botulinum Toxin Type A in the Treatment of Lower Limb Spasticity in Children with Cerebral Palsy. Cochrane Database Syst. Rev. 2019, 10, CD001408. [Google Scholar] [CrossRef]
  9. Cosgrove, A.P.; Corry, I.S.; Graham, H.K. Botulinum Toxin in the Management of the Lower Limb in Cerebral Palsy. Dev. Med. Child Neurol. 1994, 36, 386–396. [Google Scholar] [CrossRef] [PubMed]
  10. Love, S.C.; Novak, I.; Kentish, M.; Desloovere, K.; Heinen, F.; Molenaers, G.; O’Flaherty, S.; Graham, H.K. Botulinum Toxin Assessment, Intervention and After-Care for Lower Limb Spasticity in Children with Cerebral Palsy: International Consensus Statement. Eur. J. Neurol. 2010, 17 (Suppl. S2), 9–37. [Google Scholar] [CrossRef] [PubMed]
  11. Molenaers, G.; Fagard, K. Botulinum Toxin A Treatment of the Lower Extremities in Children with Cerebral Palsy. J. Child. Orthop. 2013, 7, 383–387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bania, T.A.; Taylor, N.F.; Baker, R.J.; Graham, H.K.; Karimi, L.; Dodd, K.J. Gross Motor Function Is an Important Predictor of Daily Physical Activity in Young People with Bilateral Spastic Cerebral Palsy. Dev. Med. Child Neurol. 2014, 56, 1163–1171. [Google Scholar] [CrossRef]
  13. Gough, M.; Shortland, A.P. Could Muscle Deformity in Children with Spastic Cerebral Palsy Be Related to an Impairment of Muscle Growth and Altered Adaptation? Dev. Med. Child Neurol. 2012, 54, 495–499. [Google Scholar] [CrossRef] [PubMed]
  14. Barrett, R.S.; Lichtwark, G.A. Gross Muscle Morphology and Structure in Spastic Cerebral Palsy: A Systematic Review. Dev. Med. Child Neurol. 2010, 52, 794–804. [Google Scholar] [CrossRef] [PubMed]
  15. Smith, L.R.; Lee, K.S.; Ward, S.R.; Chambers, H.G.; Lieber, R.L. Hamstring Contractures in Children with Spastic Cerebral Palsy Result from a Stiffer Extracellular Matrix and Increased In Vivo Sarcomere Length. J. Physiol. 2011, 589 Pt 10, 2625–2639. [Google Scholar] [CrossRef] [PubMed]
  16. Hägglund, G.; Wagner, P. Spasticity of the Gastrosoleus Muscle Is Related to the Development of Reduced Passive Dorsiflexion of the Ankle in Children with Cerebral Palsy. Acta Orthop. 2011, 82, 744–748. [Google Scholar] [CrossRef] [PubMed]
  17. Mathewson, M.A.; Lieber, R.L. Pathophysiology of Muscle Contractures in Cerebral Palsy. Phys. Med. Rehabil. Clin. N. Am. 2015, 26, 57–67. [Google Scholar] [CrossRef]
  18. Multani, I.; Manji, J.; Ison, T.H.; Khot, A.; Graham, K. Botulinum Toxin in the Management of Children with Cerebral Palsy. Pediatr. Drugs 2019, 21, 261–281. [Google Scholar] [CrossRef]
  19. Gough, M.; Fairhurst, C.; Shortland, A. Botulinum Toxin and Cerebral Palsy: Time for Reflection? Dev. Med. Child Neurol. 2005, 47, 709. [Google Scholar] [CrossRef]
  20. Walhain, F.; Desloovere, K.; Declerck, M.; Van Campenhout, A.; Bar-On, L. Interventions and Lower-Limb Macroscopic Muscle Morphology in Children with Spastic Cerebral Palsy: A Scoping Review. Dev. Med. Child Neurol. 2021, 63, 274–286. [Google Scholar] [CrossRef]
  21. Schroeder, A.S.; Ertl-Wagner, B.; Britsch, S.; Schroder, J.M.; Nikolin, S.; Weis, J.; Muller-Felber, W.; Koerte, I.; Stehr, M.; Berweck, S.; et al. Muscle Biopsy Substantiates Long-Term MRI Alterations One Year after a Single Dose of Botulinum Toxin Injected into the Lateral Gastrocnemius Muscle of Healthy Volunteers. Mov. Disord. 2009, 24, 1494–1503. [Google Scholar] [CrossRef]
  22. Fortuna, R.; Horisberger, M.; Vaz, M.A.; Herzog, W. Do Skeletal Muscle Properties Recover Following Repeat Onabotulinum Toxin A Injections? J. Biomech. 2013, 46, 2426–2433. [Google Scholar] [CrossRef] [PubMed]
  23. Multani, I.; Manji, J.; Tang, M.J.; Herzog, W.; Howard, J.J.; Graham, H.K. Sarcopenia, Cerebral Palsy, and Botulinum Toxin Type A. JBJS Rev. 2019, 7, e4. [Google Scholar] [CrossRef] [PubMed]
  24. Herskind, A.; Ritterband-Rosenbaum, A.; Willerslev-Olsen, M.; Lorentzen, J.; Hanson, L.; Lichtwark, G.; Nielsen, J.B. Muscle Growth Is Reduced in 15-Month-Old Children with Cerebral Palsy. Dev. Med. Child Neurol. 2016, 58, 485–491. [Google Scholar] [CrossRef]
  25. Schless, S.; Cenni, F.; Bar-On, L.; Hanssen, B.; Kalkman, B.; O’brien, T.; Aertbeliën, E.; Van Campenhout, A.; Molenaers, G.; Desloovere, K. Medial Gastrocnemius Volume and Echo-intensity after Botulinum Neurotoxin A Interventions in Children with Spastic Cerebral Palsy. Dev. Med. Child Neurol. 2019, 61, 783–790. [Google Scholar] [CrossRef]
  26. De Beukelaer, N.; Weide, G.; Huyghe, E.; Vandekerckhove, I.; Hanssen, B.; Peeters, N.; Uytterhoeven, J.; Deschrevel, J.; Maes, K.; Corvelyn, M.; et al. Reduced Cross-Sectional Muscle Growth Six Months after Botulinum Toxin Type-A Injection in Children with Spastic Cerebral Palsy. Toxins 2022, 14, 139. [Google Scholar] [CrossRef]
  27. Alexander, C.; Elliott, C.; Valentine, J.; Stannage, K.; Bear, N.; Donnelly, C.J.; Shipman, P.; Reid, S. Muscle Volume Alterations after First Botulinum Neurotoxin A Treatment in Children with Cerebral Palsy: A 6-Month Prospective Cohort Study. Dev. Med. Child Neurol. 2018, 60, 1165–1171. [Google Scholar] [CrossRef]
  28. Park, E.S.; Sim, E.; Rha, D.; Jung, S. Architectural Changes of the Gastrocnemius Muscle after Botulinum Toxin Type A Injection in Children with Cerebral Palsy. Yonsei Med. J. 2014, 55, 1406–1412. [Google Scholar] [CrossRef]
  29. Williams, S.A.; Reid, S.; Elliott, C.; Shipman, P.; Valentine, J. Muscle Volume Alterations in Spastic Muscles Immediately Following Botulinum Toxin Type-A Treatment in Children with Cerebral Palsy. Dev. Med. Child Neurol. 2013, 55, 813–820. [Google Scholar] [CrossRef] [PubMed]
  30. Van Campenhout, A.; Verhaegen, A.; Pans, S.; Molenaers, G. Botulinum Toxin Type A Injections in the Psoas Muscle of Children with Cerebral Palsy: Muscle Atrophy after Motor End Plate-Targeted Injections. Res. Dev. Disabil. 2013, 34, 1052–1058. [Google Scholar] [CrossRef] [PubMed]
  31. Kinney, M.C.; Dayanidhi, S.; Dykstra, P.B.; McCarthy, J.J.; Peterson, C.A.; Lieber, R.L. Reduced Skeletal Muscle Satellite Cell Number Alters Muscle Morphology after Chronic Stretch but Allows Limited Serial Sarcomere Addition. Muscle Nerve 2017, 55, 384–392. [Google Scholar] [CrossRef] [Green Version]
  32. Dayanidhi, S.; Dykstra, P.B.; Lyubasyuk, V.; McKay, B.R.; Chambers, H.G.; Lieber, R.L. Reduced Satellite Cell Number in Situ in Muscular Contractures from Children with Cerebral Palsy. J. Orthop. Res. 2015, 33, 1039–1045. [Google Scholar] [CrossRef] [PubMed]
  33. Desloovere, K.; Molenaers, G.; De Cat, J.; Pauwels, P.; Van Campenhout, A.; Ortibus, E.; Fabry, G.; De Cock, P. Motor Function Following Multilevel Botulinum Toxin Type A Treatment in Children with Cerebral Palsy. Dev. Med. Child Neurol. 2006, 49, 56–61. [Google Scholar] [CrossRef] [PubMed]
  34. Molenaers, G.; Desloovere, K.; Fabry, G.; De Cock, P. The Effects of Quantitative Gait Assessment and Botulinum Toxin a on Musculoskeletal Surgery in Children with Cerebral Palsy. J. Bone Jt. Surg. Am. 2006, 88, 161–170. [Google Scholar] [CrossRef]
  35. Valentine, J.; Davidson, S.A.; Bear, N.; Blair, E.; Paterson, L.; Ward, R.; Forbes, D.; Elliott, C. A Prospective Study Investigating Gross Motor Function of Children with Cerebral Palsy and GMFCS Level II after Long-Term Botulinum Toxin Type A Use. BMC Pediatr. 2020, 20, 7. [Google Scholar] [CrossRef] [PubMed]
  36. Fonseca, P.R.; Calhes Franco de Moura, R.; Galli, M.; Santos Oliveira, C. Effect of Physiotherapeutic Intervention on the Gait after the Application of Botulinum Toxin in Children with Cerebral Palsy: Systematic Review. Eur. J. Phys. Rehabil. Med. 2018, 54, 757–765. [Google Scholar] [CrossRef]
  37. Hastings-Ison, T.; Sangeux, M.; Thomason, P.; Rawicki, B.; Fahey, M.; Graham, H.K. Onabotulinum Toxin-A (Botox) for Spastic Equinus in Cerebral Palsy: A Prospective Kinematic Study. J. Child. Orthop. 2018, 12, 390–397. [Google Scholar] [CrossRef]
  38. Schwartz, M.H.; Ries, A.J.; Georgiadis, A.G. Short-Term Causal Effects of Common Treatments in Ambulatory Children and Young Adults with Cerebral Palsy: Three Machine Learning Estimates. Sci. Rep. 2022, 12, 7818. [Google Scholar] [CrossRef]
  39. Read, F.A.; Boyd, R.N.; Barber, L.A. Longitudinal Assessment of Gait Quality in Children with Bilateral Cerebral Palsy Following Repeated Lower Limb Intramuscular Botulinum Toxin-A Injections. Res. Dev. Disabil. 2017, 68, 35–41. [Google Scholar] [CrossRef]
  40. Barber, L.; Hastings-Ison, T.; Baker, R.; Kerr Graham, H.; Barrett, R.; Lichtwark, G. The Effects of Botulinum Toxin Injection Frequency on Calf Muscle Growth in Young Children with Spastic Cerebral Palsy: A 12-Month Prospective Study. J. Child. Orthop. 2013, 7, 425–433. [Google Scholar] [CrossRef]
  41. Kaňovský, P.; Baresš, M.; Severa, S.; Richardson, A. Long-Term Efficacy and Tolerability of 4-Monthly versus Yearly Botulinum Toxin Type A Treatment for Lower-Limb Spasticity in Children with Cerebral Palsy. Dev. Med. Child Neurol. 2009, 51, 436–445. [Google Scholar] [CrossRef]
  42. Hastings-Ison, T.; Blackburn, C.; Rawicki, B.; Fahey, M.; Simpson, P.; Baker, R.; Graham, K. Injection Frequency of Botulinum Toxin A for Spastic Equinus: A Randomized Clinical Trial. Dev. Med. Child Neurol. 2016, 58, 750–757. [Google Scholar] [CrossRef] [PubMed]
  43. Heinen, F.; Desloovere, K.; Schroeder, A.S.; Berweck, S.; Borggraefe, I.; van Campenhout, A.; Andersen, G.L.; Aydin, R.; Becher, J.G.; Bernert, G.; et al. The Updated European Consensus 2009 on the Use of Botulinum Toxin for Children with Cerebral Palsy. Eur. J. Paediatr. Neurol. 2010, 14, 45–66. [Google Scholar] [CrossRef]
  44. Minamoto, V.B.; Suzuki, K.P.; Bremner, S.N.; Lieber, R.L.; Ward, S.R. Dramatic Changes in Muscle Contractile and Structural Properties after 2 Botulinum Toxin Injections. Muscle Nerve 2015, 52, 649–657. [Google Scholar] [CrossRef]
  45. Hastings-Ison, T.; Graham, K. Lanbotulinumtoxin-A: Safe for Children with Cerebral Palsy but Is It Effective? Dev. Med. Child Neurol. 2022, 1–2. [Google Scholar] [CrossRef]
  46. Delgado, M.R.; Tilton, A.; Russman, B.; Benavides, O.; Bonikowski, M.; Carranza, J.; Dabrowski, E.; Dursun, N.; Gormley, M.; Jozwiak, M.; et al. AbobotulinumtoxinA for Equinus Foot Deformity in Cerebral Palsy: A Randomized Controlled Trial. Pediatrics 2016, 137, e20152830. [Google Scholar] [CrossRef]
  47. Klejman, S.; Andrysek, J.; Dupuis, A.; Wright, V. Test-Retest Reliability of Discrete Gait Parameters in Children with Cerebral Palsy. Arch. Phys. Med. Rehabil. 2010, 91, 781–787. [Google Scholar] [CrossRef] [PubMed]
  48. Baker, R.; McGinley, J.L.; Schwartz, M.; Thomason, P.; Rodda, J.; Graham, H.K. The Minimal Clinically Important Difference for the Gait Profile Score. Gait Posture 2012, 35, 612–615. [Google Scholar] [CrossRef]
  49. Kainz, H.; Graham, D.; Edwards, J.; Walsh, H.P.J.; Maine, S.; Boyd, R.N.; Lloyd, D.G.; Modenese, L.; Carty, C.P. Reliability of Four Models for Clinical Gait Analysis. Gait Posture 2017, 54, 325–331. [Google Scholar] [CrossRef] [PubMed]
  50. Papageorgiou, E.; Simon-Martinez, C.; Molenaers, G.; Ortibus, E.; Van Campenhout, A.; Desloovere, K. Are Spasticity, Weakness, Selectivity, and Passive Range of Motion Related to Gait Deviations in Children with Spastic Cerebral Palsy? A Statistical Parametric Mapping Study. PLoS ONE 2019, 14, e0223363. [Google Scholar] [CrossRef]
  51. Schless, S.; Desloovere, K.; Aertbeliën, E.; Molenaers, G.; Huenaerts, C.; Bar-On, L. The Intra- and Inter-Rater Reliability of an Instrumented Spasticity Assessment in Children with Cerebral Palsy. PLoS ONE 2015, 10, e0131011. [Google Scholar] [CrossRef] [Green Version]
  52. Koerte, I.K.; Schroeder, A.S.; Fietzek, U.M.; Borggraefe, I.; Kerscher, M.; Berweck, S.; Reiser, M.; Ertl-Wagner, B.; Heinen, F. Muscle Atrophy beyond the Clinical Effect after a Single Dose of Onabotulinumtoxina Injected in the Procerus Muscle: A Study with Magnetic Resonance Imaging. Dermatol. Surg. 2013, 39, 761–765. [Google Scholar] [CrossRef]
  53. Williams, S.A.; Elliott, C.; Valentine, J.; Gubbay, A.; Shipman, P.; Reid, S. Combining Strength Training and Botulinum Neurotoxin Intervention in Children with Cerebral Palsy: The Impact on Muscle Morphology and Strength. Disabil. Rehabil. 2013, 35, 596–605. [Google Scholar] [CrossRef]
  54. Fortuna, R.; Vaz, M.A.; Youssef, A.R.; Longino, D.; Herzog, W. Changes in Contractile Properties of Muscles Receiving Repeat Injections of Botulinum Toxin (Botox). J. Biomech. 2011, 44, 39–44. [Google Scholar] [CrossRef]
  55. Pingel, J.; Nielsen, M.S.; Lauridsen, T.; Rix, K.; Bech, M.; Alkjaer, T.; Andersen, I.T.; Nielsen, J.B.; Feidenhansl, R. Injection of High Dose Botulinum-Toxin A Leads to Impaired Skeletal Muscle Function and Damage of the Fibrilar and Non-Fibrilar Structures. Sci. Rep. 2017, 7, 14746. [Google Scholar] [CrossRef]
  56. Pitcher, C.A.; Elliott, C.M.; Panizzolo, F.A.; Valentine, J.P.; Stannage, K.; Reid, S.L. Ultrasound Characterization of Medial Gastrocnemius Tissue Composition in Children with Spastic Cerebral Palsy. Muscle Nerve 2015, 52, 397–403. [Google Scholar] [CrossRef]
  57. Svane, C.; Forman, C.R.; Rasul, A.; Nielsen, J.B.; Lorentzen, J. Muscle Contractures in Adults with Cerebral Palsy Characterized by Combined Ultrasound-Derived Echo Intensity and Handheld Dynamometry Measures. Ultrasound Med. Biol. 2022, 48, 694–701. [Google Scholar] [CrossRef]
  58. Paris, M.T.; Mourtzakis, M. Muscle Composition Analysis of Ultrasound Images: A Narrative Review of Texture Analysis. Ultrasound Med. Biol. 2021, 47, 880–895. [Google Scholar] [CrossRef] [PubMed]
  59. Picelli, A.; Bonetti, P.; Fontana, C.; Barausse, M.; Dambruoso, F.; Gajofatto, F.; Girardi, P.; Manca, M.; Gimigliano, R.; Smania, N. Is Spastic Muscle Echo Intensity Related to the Response to Botulinum Toxin Type a in Patients with Stroke? A Cohort Study. Arch. Phys. Med. Rehabil. 2012, 93, 1253–1258. [Google Scholar] [CrossRef] [PubMed]
  60. Pillen, S.; Tak, R.O.; Zwarts, M.J.; Lammens, M.M.Y.; Verrijp, K.N.; Arts, I.M.P.; van der Laak, J.A.; Hoogerbrugge, P.M.; van Engelen, B.G.M.; Verrips, A. Skeletal Muscle Ultrasound: Correlation Between Fibrous Tissue and Echo Intensity. Ultrasound Med. Biol. 2009, 35, 443–446. [Google Scholar] [CrossRef] [PubMed]
  61. Taniguchi, M.; Fukumoto, Y.; Yagi, M.; Yamagata, M.; Kobayashi, M.; Yamada, Y.; Kimura, M.; Ichihashi, N. Enhanced Echo Intensity and a Higher Extracellular Water-to-Intracellular Water Ratio Are Helpful Clinical Signs for Detecting Muscle Degeneration in Patients with Knee Osteoarthritis. Clin. Rheumatol. 2021, 40, 4207–4215. [Google Scholar] [CrossRef]
  62. Peeters, N.; Van Campenhout, A.; Hanssen, B.; Cenni, F.; Schless, S.-H.; Van den Broeck, C.; Desloovere, K.; Bar-On, L. Joint and Muscle Assessments of the Separate Effects of Botulinum NeuroToxin-A and Lower-Leg Casting in Children with Cerebral Palsy. Front. Neurol. 2020, 11, 210. [Google Scholar] [CrossRef]
  63. Malaiya, R.; McNee, A.E.; Fry, N.R.; Eve, L.C.; Gough, M.; Shortland, A.P. The Morphology of the Medial Gastrocnemius in Typically Developing Children and Children with Spastic Hemiplegic Cerebral Palsy. J. Electromyogr. Kinesiol. 2007, 17, 657–663. [Google Scholar] [CrossRef] [PubMed]
  64. Bar-On, L.; Aertbeliën, E.; Van Campenhout, A.; Molenaers, G.; Desloovere, K. Treatment Response to Botulinum Neurotoxin-A in Children with Cerebral Palsy Categorized by the Type of Stretch Reflex Muscle Activation. Front. Neurol. 2020, 11, 378. [Google Scholar] [CrossRef] [PubMed]
  65. Bar-On, L.; Aertbeliën, E.; Molenaers, G.; Van Campenhout, A.; Vandendoorent, B.; Nieuwenhuys, A.; Jaspers, E.; Hunaerts, C.; Desloovere, K. Instrumented Assessment of the Effect of Botulinum Toxin-A in the Medial Hamstrings in Children with Cerebral Palsy. Gait Posture 2014, 39, 17–22. [Google Scholar] [CrossRef] [PubMed]
  66. Nieuwenhuys, A.; Papageorgiou, E.; Pataky, T.; De Laet, T.; Molenaers, G.; Desloovere, K. Literature Review and Comparison of Two Statistical Methods to Evaluate the Effect of Botulinum Toxin Treatment on Gait in Children with Cerebral Palsy. PLoS ONE 2016, 11, e0152697. [Google Scholar] [CrossRef] [PubMed]
  67. Kahraman, A.; Seyhan, K.; Değer, Ü.; Kutlutürk, S.; Mutlu, A. Should Botulinum Toxin A Injections Be Repeated in Children with Cerebral Palsy? A Systematic Review. Dev. Med. Child Neurol. 2016, 58, 910–917. [Google Scholar] [CrossRef]
  68. Franki, I.; Bar-On, L.; Molenaers, G.; Van Campenhout, A.; Craenen, K.; Desloovere, K.; Feys, H.; Pauwels, P.; De Cat, J.; Ortibus, E. Tone Reduction and Physical Therapy: Strengthening Partners in Treatment of Children with Spastic Cerebral Palsy. Neuropediatrics 2020, 51, 89–104. [Google Scholar] [CrossRef]
  69. Van Campenhout, A.; Molenaers, G. Localization of the Motor Endplate Zone in Human Skeletal Muscles of the Lower Limb: Anatomical Guidelines for Injection with Botulinum Toxin. Dev. Med. Child Neurol. 2011, 53, 108–119. [Google Scholar] [CrossRef]
  70. Cenni, F.; Monari, D.; Schless, S.-H.; Bar-On, L.; Bruyninckx, H.; Desloovere, K. Combining Motion Analysis and Ultrasound to Analyse Muscles in Children with Cerebral Palsy. Gait Posture 2018, 42, S2–S3. [Google Scholar] [CrossRef]
  71. Cenni, F.; Schless, S.H.; Bar-On, L.; Aertbeliën, E.; Bruyninckx, H.; Hanssen, B.; Desloovere, K. Reliability of a Clinical 3D Freehand Ultrasound Technique: Analyses on Healthy and Pathological Muscles. Comput. Methods Programs Biomed. 2018, 156, 97–103. [Google Scholar] [CrossRef]
  72. Cenni, F.; Schless, S.-H.; Bar-On, L.; Molenaers, G.; Van Campenhout, A.; Aertbeliën, E.; Bruyninckx, H.; Hanssen, B.; Desloovere, K. Can in Vivo Medial Gastrocnemius Muscle–Tendon Unit Lengths Be Reliably Estimated by Two Ultrasonography Methods? A Within-Session Analysis. Ultrasound Med. Biol. 2018, 44, 110–118. [Google Scholar] [CrossRef] [PubMed]
  73. Hanssen, B.; De Beukelaer, N.; Schless, S.-H.; Cenni, F.; Bar-On, L.; Peeters, N.; Molenaers, G.; Van Campenhout, A.; Van den Broeck, C.; Desloovere, K. Reliability of Processing 3-D Freehand Ultrasound Data to Define Muscle Volume and Echo-Intensity in Pediatric Lower Limb Muscles with Typical Development or with Spasticity. Ultrasound Med. Biol. 2021, 47, 2702–2712. [Google Scholar] [CrossRef] [PubMed]
  74. Cenni, F.; Schless, S.; Monari, D.; Bar-on, L.; Aertbeliën, E.; Bruyninckx, H.; Hanssen, B.; Desloovere, K. An Innovative Solution to Reduce Muscle Deformation during Ultrasonography Data Collection. J. Biomech. 2018, 77, 194–200. [Google Scholar] [CrossRef] [PubMed]
  75. Treece, G.M.; Prager, R.W.; Gee, A.H.; Berman, L. Fast Surface and Volume Estimation from Non-Parallel Cross-Sections, for Freehand Three-Dimensional Ultrasound. Med. Image Anal. 1999, 3, 141–173. [Google Scholar] [CrossRef]
  76. Davis, R.B.; Õunpuu, S.; Tyburski, D.; Gage, J.R. A Gait Analysis Data Collection and Reduction Technique. Hum. Mov. Sci. 1991, 10, 575–587. [Google Scholar] [CrossRef]
  77. Bar-On, L.; Aertbeliën, E.; Wambacq, H.; Severijns, D.; Lambrecht, K.; Dan, B.; Huenaerts, C.; Bruyninckx, H.; Janssens, L.; Van Gestel, L.; et al. A Clinical Measurement to Quantify Spasticity in Children with Cerebral Palsy by Integration of Multidimensional Signals. Gait Posture 2013, 38, 141–147. [Google Scholar] [CrossRef] [PubMed]
  78. Bohannon, R.W.; Smith, M.B. Interrater Reliability of a Modified Ashworth Scale of Muscle Spasticity. Phys. Ther. 1987, 67, 206–207. [Google Scholar] [CrossRef]
  79. Li, F.; Wu, Y.; Li, X. Test-Retest Reliability and Inter-Rater Reliability of the Modified Tardieu Scale and the Modified Ashworth Scale in Hemiplegic Patients with Stroke. Eur. J. Phys. Rehabil. Med. 2014, 50, 9–15. [Google Scholar]
  80. Hislop, H.J.; Montgomery, J. Daniels and Worthingham’s Muscle Testing: Techniques of Manual Examination, 8th ed.; Saunders, Elsevier: St. Louis, MO, USA, 2007. [Google Scholar]
  81. Fowler, E.G.; Staudt, L.A.; Greenberg, M.B.; Oppenheim, W.L. Selective Control Assessment of the Lower Extremity (SCALE): Development, Validation, and Interrater Reliability of a Clinical Tool for Patients with Cerebral Palsy. Dev. Med. Child Neurol. 2009, 51, 607–614. [Google Scholar] [CrossRef]
  82. Pataky, T.C.; Vanrenterghem, J.; Robinson, M.A. The Probability of False Positives in Zero-Dimensional Analyses of One-Dimensional Kinematic, Force and EMG Trajectories. J. Biomech. 2016, 49, 1468–1476. [Google Scholar] [CrossRef]
  83. Gómez-Pérez, C.; Font-Llagunes, J.M.; Martori, J.C.; Vidal Samsó, J. Gait Parameters in Children with Bilateral Spastic Cerebral Palsy: A Systematic Review of Randomized Controlled Trials. Dev. Med. Child Neurol. 2019, 61, 770–782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Box plots with individual data points of normalized muscle volume of the medial gastrocnemius and distal compartment of the semitendinosus before and after BoNT-A treatment in the BoNT-A intervention group and at the two assessments in the control untreated group. * alpha < 0.05, BoNT-A = Botulinum Neurotoxin-A, kg = kilograms, m = meters, ml = milliliters, n.s. = not significantly different.
Figure 1. Box plots with individual data points of normalized muscle volume of the medial gastrocnemius and distal compartment of the semitendinosus before and after BoNT-A treatment in the BoNT-A intervention group and at the two assessments in the control untreated group. * alpha < 0.05, BoNT-A = Botulinum Neurotoxin-A, kg = kilograms, m = meters, ml = milliliters, n.s. = not significantly different.
Toxins 14 00676 g001
Figure 2. Box plots with individual data points of echogenicity intensity of the medial gastrocnemius and distal compartment of the semitendinosus before and after BoNT-A treatment in the BoNT-A intervention group and at the two assessments in the untreated control group. BoNT-A = Botulinum Neurotoxin-A, n.s. = not significantly different.
Figure 2. Box plots with individual data points of echogenicity intensity of the medial gastrocnemius and distal compartment of the semitendinosus before and after BoNT-A treatment in the BoNT-A intervention group and at the two assessments in the untreated control group. BoNT-A = Botulinum Neurotoxin-A, n.s. = not significantly different.
Toxins 14 00676 g002
Figure 3. Overlay of the kinematic waveforms of the ankle and knee joint during gait, pre- and post-intervention. The green waveform represents the pre-treatment kinematics, and the blue waveform represents the post-treatment joint kinematics. The black bar at the bottom of the image indicates the location of the significant clusters, exceeding 3% of the gait cycle and exceeding the standard error of measurement. Deg = degrees. t = critical threshold, calculated by the random field theory.
Figure 3. Overlay of the kinematic waveforms of the ankle and knee joint during gait, pre- and post-intervention. The green waveform represents the pre-treatment kinematics, and the blue waveform represents the post-treatment joint kinematics. The black bar at the bottom of the image indicates the location of the significant clusters, exceeding 3% of the gait cycle and exceeding the standard error of measurement. Deg = degrees. t = critical threshold, calculated by the random field theory.
Toxins 14 00676 g003
Figure 4. Flow chart of all recruited and finally included subjects in both the intervention and untreated control group, together with a summary of data availability. 3DfUS = three-dimensional freehand ultrasound, 3DGA = three-dimensional gait analysis, ISA = instrumented spasticity assessment, MEH = medial hamstrings, MG = medial gastrocnemius, n = number, PF = plantar flexors, ST = semitendinosus.
Figure 4. Flow chart of all recruited and finally included subjects in both the intervention and untreated control group, together with a summary of data availability. 3DfUS = three-dimensional freehand ultrasound, 3DGA = three-dimensional gait analysis, ISA = instrumented spasticity assessment, MEH = medial hamstrings, MG = medial gastrocnemius, n = number, PF = plantar flexors, ST = semitendinosus.
Toxins 14 00676 g004
Table 1. Participant characteristics of the included participants of the intervention and control group.
Table 1. Participant characteristics of the included participants of the intervention and control group.
Intervention Group (n = 25)Control Group (n = 20)Between Group
p-Value
Age (years) 16.4 (4.6–8.5)7.6 (6.2–9.4)0.144
Gender (n)boys = 13, girls = 12boys = 16, girls = 4
GMFCS (n)I = 14, II = 9, III = 2I = 14, II = 5, III = 1
Paralysis (n)bilateral = 15
unilateral = 10
bilateral = 10
unilateral = 10
Body height (cm) 1114.5 (106.5–130.5)127.3 (115.7–133.8)0.021
Body mass (kg) 120.0 (17.0–23.9)24.4 (21.1–30.7)0.010
Maximal dorsiflexion angle, knee extended 15.0 (0.0–10.0)10.0 (0.0–10.0)0.267
Popliteal angle, unilateral 1−55.0 (−42.5–−62.5)−55.0 (−50.0–−63.8)0.737
MAS PF, knee extended 12.0 (1.6–3.0)1.5 (1.4–1.6)
MAS PF, knee 90 degrees 12.8 (1.5–2.0)1.5 (1.4–1.5)
MAS hamstrings 11.5 (1.0–1.9)1.0 (0.8–1.5)
BoNT-A MG (units/kg), n = 24 12.7 (2.2–2.9)n.a.
BoNT-A ST (units/kg), n = 21 11.7 (1.4–2.0)n.a.
BoNT-A treatment history (n)BoNT-A naive = 10
1–3 previous BoNT-A = 11
>4 previous BoNT-A = 4
BoNT-A naive = 9
1–3 previous BoNT-A = 10
>4 previous BoNT-A = 1
Orthotic management, AFO’s during the day (n)Frequently used (>50% of the day) = 20
Not frequently used (≤50% of the day) = 3
Not used = 2
Frequently used (>50% of the day) = 13
Not frequently used (≤50% of the day) = 2
Not used = 2
Unknown = 3
Oral tone reduction (n)Yes = 3
No = 22
Yes = 0
No = 20
1 Values are presented as medians with 25th and 75th percentiles, AFO = ankle foot orthoses, BoNT-A = Botulinum Neurotoxin type-A, cm = centimeters, GMFCS = Gross motor function classification system, kg = kilograms, MAS = modified Ashworth score, MG = medial gastrocnemius, n = number, n.a. = not applicable, PF = plantar flexors, ST = semitendinosus.
Table 2. Muscle-specific growth rates calculated over the period between the two assessments for both the intervention and control group.
Table 2. Muscle-specific growth rates calculated over the period between the two assessments for both the intervention and control group.
Intervention GroupControl GroupBetween-Group
p-Value
Medial gastrocnemius growth rate (mL/month) 10.00
(−1.34–0.97)
0.17
(−0.23–0.87)
0.346
Semitendinosus, distal compartment growth rate (mL/month) 1−0.49
(−1.45–0.58)
0.02
(−0.39–0.48)
0.226
1 Growth rates are presented as medians with corresponding 25th and 75th percentiles. mL = milliliters.
Table 3. Sagittal plane gait kinematic parameters, gait profile score and spatial-temporal parameters of the intervention group during the pre- and post-treatment assessments and difference scores.
Table 3. Sagittal plane gait kinematic parameters, gait profile score and spatial-temporal parameters of the intervention group during the pre- and post-treatment assessments and difference scores.
Gait Intervention Group
Kinematic Parameters—Ankle 1Pre BoNT-APost BoNT-ADifference Scorep-value, treatment-induced/within group
ankle ROM in sagittal plane (total GC)
(degrees)
25.85
(19.83–34.55)
24.18
(21.51–29.13)
−1.74
(−9.04–3.07)
0.143
ankle angle in sagittal plane at IC (degrees)−3.83
(−11.13–2.14)
−2.25
(−5.25–4.94)
3.56
(0.47–6.36)
0.002 **
maximum ankle angle (max DF) in sagittal plane during stance phase (degrees)9.47
(0.70–14.52)
12.89
(7.37–17.76)
3.99
(−0.32–9.99)
0.004 **
maximum ankle angle (max DF) in sagittal plane during swing phase (degrees)−1.71
(−9.63–3.12)
2.53
(−3.66–7.03)
5.34
(1.16–9.70)
<0.001 **
Kinematic parameters—Knee 1pre BoNT-Apost BoNT-Adifference scorep-value, treatment-induced/within group
knee ROM in sagittal plane (degrees)53.50
(45.13–62.03)
54.84
(46.85–60.12)
1.94
(−3.79–8.84)
0.158
knee angle in sagittal plane at IC (degrees)26.45
(20.93–35.43)
25.53
(17.26–30.53)
−1.69
(−9.19–1.84)
0.048
minimal knee angle in sagittal plane during stance (degrees)5.06
(0.79–14.05)
6.74
(−1.22–10.09)
−2.44
(−7.61–2.23)
0.065
Gait profile score 1pre BoNT-Apost BoNT-Adifference scorep-value, treatment-induced/within group
Gait profile score for angles (degrees)9.80
(7.41–11.33)
8.30
(6.71–9.41)
−1.76
(−2.51 – −0.15)
0.005 *
Spatial-temporal parameters 1 pre BoNT-Apost BoNT-Adifference scorep-value, treatment-induced/within group
Cadence (number of steps/second) 2.25
(1.95–2.43)
2.19
(1.91–2.35)
−0.05
(−0.28–0.10)
0.156
Walking velocity (meter/second)0.89
(0.67–1.10)
0.85
(0.73–1.05)
0.01
(−0.10–0.10)
0.911
Stride length (meters)0.81
(0.66–0.96)
0.79
(0.69–1.04)
0.02
(−0.03–0.08)
0.108
1 Pre, post and difference scores are presented as median and 25th and 75th percentiles. A Bonferroni correction was applied for the kinematic parameters per joint (ankle: alpha 0.05/4 = 0.0125, knee: alpha 0.05/3 = 0.0167), * alpha < 0.05, ** significant different after Bonferroni correction: alpha < 0.0125 for the ankle and alpha < 0.0167 for the knee. DF = dorsiflexion, GC = gait cycle, IC = initial contact, max = maximum, ROM = range of motion.
Table 4. Outcomes of the instrumented spasticity assessment for a subgroup (n = 14) of the BoNT-A intervention group.
Table 4. Outcomes of the instrumented spasticity assessment for a subgroup (n = 14) of the BoNT-A intervention group.
Instrumented Spasticity Assessment—Treated MusclesIntervention Group (n = 14)
Parameterpre BoNT-Apost BoNT-ADifference Post-Prep-Value, Treatment-Induced/within Group
Medial gastrocnemius, EMG—absolute high RMS–EMG (µV)11.52
(4.87–22.07)
4.91
(3.67–9.45)
−4.73
(−15.74–1.31)
0.019 *
Medial hamstrings, EMG—absolute high RMS–EMG (µV)12.07
(5.77–16.04)
7.98
(3.68–10.58)
−2.82
(7.64–0.21)
0.041 *
* alpha < 0.05, EMG = electromyography, RMS = root mean square, µV = microvolts.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peeters, N.; Papageorgiou, E.; Hanssen, B.; De Beukelaer, N.; Staut, L.; Degelaen, M.; Van den Broeck, C.; Calders, P.; Feys, H.; Van Campenhout, A.; et al. The Short-Term Impact of Botulinum Neurotoxin-A on Muscle Morphology and Gait in Children with Spastic Cerebral Palsy. Toxins 2022, 14, 676. https://doi.org/10.3390/toxins14100676

AMA Style

Peeters N, Papageorgiou E, Hanssen B, De Beukelaer N, Staut L, Degelaen M, Van den Broeck C, Calders P, Feys H, Van Campenhout A, et al. The Short-Term Impact of Botulinum Neurotoxin-A on Muscle Morphology and Gait in Children with Spastic Cerebral Palsy. Toxins. 2022; 14(10):676. https://doi.org/10.3390/toxins14100676

Chicago/Turabian Style

Peeters, Nicky, Eirini Papageorgiou, Britta Hanssen, Nathalie De Beukelaer, Lauraine Staut, Marc Degelaen, Christine Van den Broeck, Patrick Calders, Hilde Feys, Anja Van Campenhout, and et al. 2022. "The Short-Term Impact of Botulinum Neurotoxin-A on Muscle Morphology and Gait in Children with Spastic Cerebral Palsy" Toxins 14, no. 10: 676. https://doi.org/10.3390/toxins14100676

APA Style

Peeters, N., Papageorgiou, E., Hanssen, B., De Beukelaer, N., Staut, L., Degelaen, M., Van den Broeck, C., Calders, P., Feys, H., Van Campenhout, A., & Desloovere, K. (2022). The Short-Term Impact of Botulinum Neurotoxin-A on Muscle Morphology and Gait in Children with Spastic Cerebral Palsy. Toxins, 14(10), 676. https://doi.org/10.3390/toxins14100676

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop