1. Introduction
Skeletal muscle changes are frequently manifested in children with spastic cerebral palsy (SCP) which interact with the neuromotor functioning throughout development[
1,
2,
3]. However, the underlying mechanisms of this potential bilateral causations or associated presence remain unknown. Additionally, neuromotor impairments, such as spasticity, manifest from young ages, progress with age and result in heterogenous clinical presentations of children with SCP[
1,
4,
5,
6,
7]. To classify the children based on the functional abilities and limitations in daily life, the Gross Motor Function Classification System Expanded and Revised (GMFCS-E&R) is frequently used [
8]. The 5 GMFCS levels are defined using standardized descriptions for which the lower levels indicate higher motor abilities. To reduce and manage the impact of neuromotor impairments on the daily motor functioning, children with SCP are frequently treated at the muscular level (e.g., botulinum neurotoxin injection or casting) from young ages [
9].
In typical development, morphological muscle growth is triggered by body size increases and is responding to the demands through protein synthesis[
10,
11,
12]. However, for children with SCP, previous research highlighted lower increases in medial gastrocnemius (MG) muscle volume (MV) with increasing age, which was unrelated to the rate of skeletal growth and was already observed at the age of 15 months [
13]. Taking body length and mass into account, the MVs were found to reduce by 3% over a 6-month follow-up period in 2-5 years (y) old children with SCP [
14]. Over the age ranges of 2 to 12 years, cross-sectional investigations identified normalized medial gastrocnemius MV deficits of 22%-57% and shorter normalized muscle belly lengths (ML) with deficits of 11-19% compared to age-matched typically developing (TD) children[
15,
16,
17,
18,
19]. Since the muscle is a very plastic tissue, the observed muscular changes in SCP could reflect adaptations to the functional demands including altered muscle tone and limited motor functioning [
12]. Indeed, previous cross-sectional studies demonstrated increasing muscle pathology, such as deficits in MV, anatomical cross-sectional area (CSA) and ML, in relation to higher GMFCS levels [
8,
13,
18,
20,
21]. Yet, research to date has not determined the impact of the severity in motor impairments on the trajectory of altered muscle growth. Overall, these previous cross-sectional observations highlighted significant alterations in the SCP muscle from a young age, which were suggested to be multifactorial in origin, for example due to impaired neural activation, inflammation and altered or non-use of the muscle [
3,
12,
22,
23,
24].
Up to now, cross-sectional investigations have already suggested altered associations between muscle morphology and age in SCP versus TD children [
13,
16]. These previous studies assumed a linear trajectory of morphological muscle growth with age. However, much uncertainty remains in the variability of the trajectories for altered muscle growth over the time-span from infancy to school-ages and between different levels of severity in children with SCP [
13,
16,
25]. Further, MV is frequently reported apart from the estimates of muscle growth in the cross-sectional dimension (e.g., CSA) and longitudinal dimension (e.g., ML). To delineate the muscle growth trajectories, repeated assessments over time per individual are required. These longitudinal datasets generate child-specific profiles per morphological measure. A better understanding of this time course may shed light on opportunities for the planning, development and promotion of muscle growth strategies, whether or not combined with the conventional therapeutic strategies.
The current prospective longitudinal study investigated morphological MG muscle growth over a wide age-span, including infants (6mo-2y), preschool (2-5y) and school-aged SCP-children (6-11y). The first objective was to model the muscle growth trajectory with respect to age by using repeated muscle assessments. Growth in the MG muscle morphology was expressed for different absolute and normalized outcomes. Second, we aimed to compare the trajectories in muscle growth between the GMFCS levels I and II-III. We hypothesized gradually increased morphological muscle outcomes with increasing age, with lower muscle growth rates in children of the higher GMFCS levels.
3. Results
Patient characteristics at the time of baseline assessment are summarized in
Table 1. At the end of the follow-up, muscle data was collected until 10.3 years for GMFCS I and 11.1 years for GMFCS II-III (
Figure 1). Only 9 children changed from GMFCS level during the follow-up for which the re-classification was performed around 2-3 years of age. The study sample received standardized clinical care such as regular physical therapy, orthotic devices and state-of-the-art medical and orthopedic services (e.g., serial casting and/or BoNT-A injections when indicated) within a multidisciplinary clinical setting (
Table S2). The results of the piecewise models for the subgroups in GMFCS levels are summarized in
Table 2,
Figure 2 and
Table S3. The comparisons of slopes and breakpoints between the groups are presented in
Table 3.
For the GMFCS Level I group, the MV increased with 12.8 ml/year (β
1, p<0.0001) up to 2.1 years of age (c
1). After these infant-ages, an MV increase of 5.7 ml/year was found until the age of 7.8y (β
2,
p<0.0001 and c
2, respectively). From this second breakpoint at these older ages, a growth rate of 3.1 ml/year was found (β
3,
p=0.0024) which was significantly lower compared to growth rates for the younger ages (β
1 vs. β
3 and β
2 vs. β
3,
p<0.05,
Table S3B). On the other hand, after an increase in nMV of 0.83 ml/kg⋅m per year until 2.1y, the nMV showed a non-significant rate of -0.03 ml/kg⋅m per year until the age of 8.0y (β
2, p=0.0933 and c
2). After this second breakpoint, a rate of -0.12 ml/kg⋅m per year was found (β
3,
p<0.0001) for which this linear trend showed significantly more decline compared to the muscle growth rates before 2 years and between 2-8 years old (β
1 vs. β
3 and β
2 vs. β3,
p<0.05,
Table S3B). Similar to the growth rate in MV, the trajectory of CSA growth showed two breakpoints (age 2.2y and 6.7y, respectively) with only in the first and second linear trend, a significant yearly increase in CSA (β
1=158.4mm
2/year, p<0.0001; β
2=37.0 mm
2/year,
p<0.0001, respectively). The nCSA trajectory showed after an increase of 8.88 mm
2/kg per year only 1 breakpoint at the age of 2.1y followed by decrease in nCSA per year (c
1 and β
2= -0.96 mm
2/kg per year,
p<0.0001, respectively).
After a significant MV increase with a rate of 4.5 ml/y (β1, p<0.0001), children in the GMFCS level II-III group showed only a breakpoint at the age of 9.1y (c1) which was followed by a rate of -0.8 ml/y (β2, p=0.7154). The nMV showed increases of 0.28 ml/kg⋅m per year until the breakpoint at the age of 2.7y (c1) and followed by decreases of 0.07 ml/kg⋅m per year (β2, p=0.0003). The trajectories for absolute and normalized CSA were modelled with two breakpoints, around approximately 2 and 9 years of age. During the (pre)school ages, nCSA decreased with significantly higher decline after the age of 9.1y (β2 vs β3, p=0.0424).
Before the age of 2 years, children with GMFCS level I showed significantly higher yearly increases for both absolute and normalized MV and CSA compared to level II-III (β
1, GMFCS-I vs. β
1, GMFCSII-III,
Table 3,
p<0.05). After the age of 2 years, absolute MV and CSA were significantly increasing in both GMFCS subgroups until the age of 6-9 years old. However, GMFCS level I tended to have a higher absolute MV growth rate compared to GMFCS level II-III (β
2 , GMFCS-I vs. β
1, GMFCSII-III =1.16ml/y, p=0.0712). Despite the significant earlier second breakpoint for GMFCS levels I (Δ1.31y for MV,
p=0.0326 and Δ2.66y for CSA,
p=0.0033), no significant differences were found in the MV and CSA growth rate at these oldest ages between the GMFCS subgroups (β
3, GMFCS-I vs β
2, GMFCS-II-IIII,
p=0.1069 and β
3, GMFCS-I vs. β
3, GMFCSII-III,
p=0.0915, respectively). The trajectory for nMV and nCSA showed decreased rates from the age of 2 years in both GMFCS groups, with a significant higher decline in nCSA in children with GMFCS level II-III after the age of 6-9y compared to level I (β
2, GMFCS-I vs.β
3, GMFCS-II-III,
p=0.0334).
Absolute ML increased with age, with significantly lower ML growth rate after the age of 5.11y in the GMFCS level I and after the age of 2.1y in GMFCS levels II-III (β1 vs. β2, p<0.05). Normalized ML showed significant different rates with increasing age between the GMFCS-groups, for which the GMFCS levels I increased with 0.33%/y until 4.88y (β1, p=0.0125) followed by a rate of -0.10%/year (β2, p=0.1977), whereas the levels II-III increased in nML with 0.12% per year from infancy to school-ages, without a breakpoint (β1, p=0.0096).
4. Discussion
In this prospective longitudinal follow-up between 6 months and 11 years of age, the trajectory of morphological muscle growth for children with SCP were observed as piecewise profiles with increasing age, indicating linear trends interrupted with breakpoints. Therefore, our hypothesis of gradually increases in muscle morphology with increasing age could only be partially accepted as expressed by the changes in observed yearly muscle growth rates before and after the breakpoints. Indeed, MV and CSA significantly increased before the age of 2 years but was followed by slower increases until the age of 6-7 years for GMFCS level I and 9 years for GMFCS level II-III. Normalized MV and CSA also increased during infancy but showed already reduced growth rates after the age of 2 years. Furthermore, from the teenage years, both the absolute and normalized MV and CSA outcomes reduced per year which is in contrast with the hypothesis of a linear trajectory of muscle growth with increasing age. Absolute ML increased over the entire age-range with a breakpoint at the age of 5 years for GMFCS level I and 2 years for GMFCS level II-III whereas nML showed less increase after this age for GMFCS level I. Lower muscle growth rates with increasing age in children of the higher GMFCS levels was hypothesized which was only confirmed for the absolute morphological growth before the age of 2 years and in teenagers. Next, the normalized MV and CSA outcomes reduced per year with a more pronounced decrease in CSA normalized to changes in body dimensions for GMFCS level II-III. Further, ML trajectories were different between the GMFCS levels and more specific, after normalization for skeletal growth.
These results revealed higher rate of muscle growth in early years of life, with lower growth rates for children with GMFCS level II-III compared to GMFCS level I. During these youngest ages, the current MV increased with 12.8 ml per year for GMFCS level I, which was slightly higher compared to the median growth rate of 10.3 ml per year, computed for the 6 months to 2 years old TD children of our retrospective cross-sectional database (
Figure S2 & Table S1). Only one previous study reported a cross-section) al growth rate for 8 months to 5 years old children, with 8.16 ml/y for TD and 3.84 ml/y for CP children [
13]. The latter is comparable to the growth rate of 4.5ml/y in the current GMFCS level II-III children suggesting an early onset of more pronounced growth deficits in children with more severe impairments. Furthermore, the nMV rates of 0.83 ml/(kg⋅m) and 0.28 ml/(kg⋅m) per year were lower than our cross-sectional TD growth rate of 1.17 ml/(kg⋅m) per year. No growth rates for nMV during early ages have been previously reported. Nevertheless, growth rate calculated on cross-sectional datasets pretended linear growth whereas the current results revealed non-linear trajectories of muscle growth with increasing age. However, caution is warranted with the interpretation of observed increases in normalized muscle data. In case of optimal normalization for changes in body sizes, growth rates close to zero are expected, representing harmonization between muscle and anthropometric growth [
39]. The validity of ratio-scaling in growing TD children, over a wide age range, has not yet been defined. The body composition in children is assumed to change with less amount of fat mass after the age of 2 years. Hence, normalization by taking body mass into account is considered invalid during infancy [
40]. Further, significantly lower intercepts for nMV, nCSA and nML indicated that children with less motor abilities develop more muscle pathology already early in life (
Table 3). These observations suggest early mechanisms hampering the muscle growth, such as the contribution of neural alterations as result of the brain lesion as well as altered patterns of muscle use. However, limited muscle data on TD infants is available to compare the current suggestion of an early onset of muscle alterations, especially in GMFCS levels II-III. Further research involving also longitudinal assessments of TD infants is required to explore the impact of normalization techniques and provide valid references to investigate the onset and development of muscle alterations in SCP infants.
Interestingly, from 2 years of age, similar muscle growth trajectories in (normalized) MV and CSA for the GMFCS subgroups were found. The average MV rates of 5.1 ml/y and 4.5 ml/y in GMFCS level I and level II-III respectively, were comparable to earlier reported MV growth rates based on 6 months and 12 months follow-up studies in 2-5 years old CP children (6.0 ml/y and 6.6 ml/y, respectively) [
14,
41]. In line with the cross-sectional ML deficits for SCP compared to TD children, the current average ML growth rates of 14.2 mm/y and 6.8 mm/y in the GMFCS level I and 18.8 mm/y and 8.6 mm/y in the GMFCS level II-III children were lower than the current cross-sectional median ML rate of 28.0 mm/y for 6mo-9y old TD children (
Table S1). Furthermore, hampered morphological muscle growth is more clearly highlighted by investigating the changes in normalized MV and CSA compared to the investigation of the absolute parameters. It should be noted that altered anthropometric growth indicated by shorter body length and lower body mass, have already been shown in children with SCP compared to TD peers [
42]. These growth deficits were found to be associated with increasing age, as well as with bilateral involvement and more severe gross motor impairments [
43]. Although anthropometric growth in children with minor motor impairments is close to age-matched TD peers [
42], the significant negative normalized muscle growth rates observed in the current study for the GMFCS level I group showed that the muscle growth was not in accordance with the skeletal growth. Since muscle morphological growth is assumed to be triggered by anthropometric growth and patterns of muscle use which might be defined by the gross motor functional abilities, the observed early presence of muscle pathology for GMFCS level II-III and increasing deficits after accounting for normalization to changes in body dimensions suggesting that pathology-related biomechanical triggers contribute to hampered muscle growth.
From the age of 6-9 years, hampered muscle growth was presented in both GMFCS-groups. Despite the lower number of datapoints at these older ages for children with GMFCS levels I, the level of gross motor mobility in daily life is suggested to further contribute to hampered growth as shown by the tendency of lower CSA growth and significantly lower nCSA growth for GMFCS level II-III compared to level I. These findings are in line with cross-sectional observations of significantly smaller MV and CSA in GMFCS levels II compared to level I for 5-12 years old children [
18]. However, the current data set did not further distinguish the growth trajectories between GMFCS levels II and III. While the underlying mechanisms of these changes in trajectory between the GMFCS levels are currently unclear, they might be attributed to reduced levels of physical activity, higher incidence of nutritional problems and increased (secondary) musculoskeletal alterations, especially from 6-9 years [
11].
The current findings also suggested a different trajectory for the growth in cross-sectional and longitudinal dimension of the MG muscle, i.e., changes in CSA and ML, respectively. For CSA, a growth trajectory with an early breakpoint at 2 years of age was observed in both GMFCS groups, whereas the trajectory of ML showed a breakpoint at the age of 4-5y in GMFCS level I and at the age of 2y in GMFCS level II-III. The trajectories in CSA are corresponding nicely to the changes in MV, suggesting that alterations in cross-sectional MG muscle growth represent the dominant contribution to the change in overall MG muscle size. This is in line with previous findings indicating that reduction in physiological CSA, a determinant for force generation, is determined by reduced MV rather than by changes in fascicle length [
44,
45]. Hence, to minimize the efforts for clinicians, the assessment of only the mid-belly anatomical CSA can be suggested as primary outcome for monitoring the muscle status and treatment follow-up.
With increasing age, evolution in neuromotor symptoms can be expected due to the natural history in children with SCP. Interestingly, breakpoint models describing changes in development at specific ages were already found in longitudinal follow-up studies for spasticity and lower limb ROM [
5,
6,
46]. These previous studies showed that spasticity is increasing to the age of 4 years, followed by a decrease each year until the age of 12 years [
6]. Decreasing ankle joint ROM was reported up to 5 years of age, followed by further decrease with age in GMFCS level I and II, while levels III showed increased ROM [
5]. The current study showed similar non-linear models for the ML outcomes in GMFCS level I with a breakpoint at the age of 5.11 years. The nML changed from 0.33% per year to -0.10% per year for GMFCS level I whereas a constant increase of 0.12%/y was presented for GMFCS level II-III. Therefore, the restricted ROM for GMFCS level I might result from reduced growth of muscle length. Furthermore, these similarities in longitudinal observations of neuromotor impairment and muscle development including the breakpoint at 5 years, support the hypothesis that formation of contractures is triggered by different factors, of which spasticity is probably not the dominant one [
46,
47].
Over the last years, 3DfUS has already been extensively used mainly to describe muscle deficits in children with SCP compared to TD children [
48]. These cross-sectional investigations included diverse age-ranges and GMFCS levels resulting in variable results of deficits (for example 22% for age-range of 2-5 years and 41% 5-12years old children) [
16,
19]. However, age-specific observations of hampered muscle growth for SCP children were missing in these previous studies. Our findings, derived from longitudinal models, emphasized the need of repeated assessments to accurately delineate muscle growth. Furthermore, the current muscle growth trajectories can be used to monitor the status of muscle pathology in individual children with SCP. The current results may support the clinical decision making of therapy selection, goals and planning at specific ages and for each GMFCS level with the aim to stimulate the muscle growth. For example, the application of strength training, considering the appropriate age and cognitive functioning, might be beneficial to maintain the muscle size relative to the skeletal growth and might be further combined with age-and child-specific prescription of protein intake to increase the muscle size[
49,
50,
51]. Beside interventions, the current observation of early muscle alterations for higher GMFCS levels and the breakpoint in growth trajectory at a very young age highlight the opportunity for prevention strategies aiming to maintain muscle size and lengths comparable to TD children and persevere muscle growth during childhood, e.g., intensive physical therapy with stimulation of lower leg movements and mechanical loading and a nutritional plan. Previous studies focusing already on the impact of BoNT-A injections, a frequently applied tone-reducing treatment, demonstrated hampered muscle growth in response to the first BoNT-A treatment [
14,
52]. This post-treatment effect on muscle growth was only assessed after 6 months for which the interference of BoNT-A injection with the trajectory of muscle changes is yet unclear. Future intervention studies could use the currently modelled muscle growth trajectories as a reference to assess the effect of treatment on muscle morphology and could finetune the timing of treatment in relation to the potential hampered growth.
This study has some limitations to consider. First, the age-matched TD data could not be included in the models due to the lack of longitudinal assessments. Yet, visual inspection of cross-sectionally measured typical muscle morphology allowed to judge the overall level of alterations in the data of children with SCP. We included children aged between 6 months and 9 years at baseline resulting in limited available data points at the beginning and end of the age-range. More longitudinal data before and after this range is relevant to further enrich the muscle growth trajectories during the entire childhood. Combined with longitudinal data of TD children, a focused analysis of the muscle growth before 2 years of age and after 6 years of age would be interesting for future research. Second, this unique longitudinal database is limited to provide trajectories of muscle growth for the specific diversity in SCP phenotypes. In this study, only the GMFCS levels representing ambulant children with SCP were included. The GMFCS level II and III were merged to ensure sufficient power for the data-analysis and to provide as much as homogeneity as possible. The current dataset was, however, limited in datapoints after the age of 6 years to describe the trajectories of muscle growth per GMFCS level. Hereto, we only descriptively explored the muscle growth per GMFCS level with the individual observed profiles (
Figure S1), indicating similar trajectories in morphological muscle growth. Nevertheless, it is important to further consider the potential heterogeneity in muscle growth based on the specificity in gross motor abilities and limitations through daily life per GMFCS level. Future studies should aim for more participants after the age of 6 years and with equally distributed number of children over each of the 3 ambulatory GMFCS levels in order to distinct the models for muscle growth per GMFCS level. Next, the muscle growth associated to the SCP topography i.e., unilateral versus bilateral SCP was also not specifically investigated. Anthropometric growth, lower limb strength and gait were found to be less involved in unilateral compared to bilateral SCP, whereas only one investigation showed more muscle growth deficits for this SCP motor type [
53,
54]. These findings suggest the need to describe, in future research, separate muscle growth trajectories for the children with uni- and bilateral SCP in interaction with GMFCS level. As a first exploration, the observed individual trajectories per motor subtype for the current GMFCS subgroups are provided in
Figure S3. Visual inspection suggests less muscle growth in children with bilateral SCP, especially for GMFCS level II-III of which 65% of the children were bilaterally involved. Further, each muscle assessment during the follow-up was not accompanied with evaluation of motor impairments and functioning such as assessment of strength, spasticity and selective motor control as well as the gross motor function measure (GMFM). These clinical results combined with muscle growth data could create a comprehensive monitoring of SCP muscle pathology, which may help to understand the underlying mechanisms of altered muscle growth and eventually help to improve patient-specific treatment management. Of note, in the current study, children who received BoNT-treatment were allowed for further follow-up in case the BoNT-A treatment occurred minimal 10 months prior to the assessment. The current database included 24 children who received BoNT-A injection during the follow-up time. Our previous research highlighted that the muscle recovery is still ongoing at 6 months post BoNT-A injections [
14]. Despite the 10-month interval between the BoNT-A treatment and muscle morphology assessments, which is already longer than the frequently used criteria of an interval of 6 months, there is no guarantee for full recovery of muscle growth at the follow-up assessment. Future studies should investigate the prolonged impact of treatments. Long-term follow-up assessments are highly needed to further understand if BoNT-A treatment is a confounding factor on the currently established trajectories of muscle growth.
Author Contributions
Conceptualization, N.D.B., E.O., A.V.C. and K.D.; methodology, N.D.B, I.V, G.M. and K.D.; formal analysis, N.D.B., I.V and G.M.; investigation, N.D.B., E.H., N.P and B.H.; resources, E.O., A.V.C., and K.D.; data curation, N.D.B., E.H., N.P. and B.H.; writing—original draft preparation, N.D.B.; writing—review and editing, N.D.B., I.V., E.H., G.M., N.P., B.H., E.O., A.V.C. and K.D.; visualization, N.D.B., I.V. and G.M.; supervision, E.O., A.V.C., and K.D.; project administration, N.D.B. and K.D. All authors have read and agreed to the published version of the manuscript.