Preprint
Article

Morphological Medial Gastrocnemius Muscle Growth in Ambulant Children with Spastic Cerebral Palsy: A Prospective Longitudinal Study

Altmetrics

Downloads

140

Views

70

Comments

0

Submitted:

09 January 2023

Posted:

12 January 2023

You are already at the latest version

Alerts
Abstract
Only cross-sectional studies have demonstrated muscle deficits in children with spastic cerebral palsy (SCP). The impact of functional impairments on altered muscle growth remains unclear. This prospective longitudinal study modelled the morphological muscle growth in 87 children with SCP (age range 6 months to 11 years, Gross Motor Function Classification System [GMFCS] level I/II/III=47/22/18). Repeated ultrasound assessments were performed during 2-year (y) follow-up with an interval of minimal 6 months. Three-dimensional freehand ultrasound was applied to assess medial gastrocnemius muscle volume (MV), mid-belly cross-sectional area (CSA) and muscle belly length (ML). Non-linear mixed models compared trajectories of (normalized) muscle growth between GMFCS-I versus GMFCS-II&III. MV and CSA growth trajectories showed piecewise model with 2 breakpoints, with the highest growth before 2y and negative growth rates after 6-9y. Before 2y, children with GMFCS-II&III already showed lower growth rates compared to GMFCS-I. From 2-9y, the growth rates did not differ between GMFCS-levels. After 9y, a more pronounced reduction in normalized CSA was observed in GMFCS-II&III. Different trajectories in ML growth were shown between the GMFCS-level subgroups. These longitudinal trajectories could be used to monitor the SCP muscle pathology during childhood and to optimize treatment planning and goals aiming to stimulate muscle growth
Keywords: 
Subject: Medicine and Pharmacology  -   Pediatrics, Perinatology and Child Health

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.

2. Materials and Methods

2.1. Study design & participants

Via the Clinical Motion Analysis Laboratory (CMAL-Pellenberg) and Cerebral Palsy Reference Center of the University Hospitals Leuven, children with a diagnosis of SCP confirmed by a pediatric neurologist, were recruited. A prospective longitudinal research protocol of 2 years with intervals of 6 months up to 2 years resulting in 2 to 5 assessments per child, was conducted between April 2018 and March 2022. SCP-children aged between 6 months and 9 years at baseline, with bi- or unilateral involvement and with GMFCS-E&R levels I to III, defined by the criteria per age-band [26], were included. Considering potential change in GMFCS level over time, which could be expected by consequence of motor development, especially at the youngest ages, the GMFCS level at the last assessment was used to define subgroups. These subgroups, i.e., GMFCS level I vs level II & III, were selected to maximize the number of observations per group and to reduce heterogeneity [7,27]. For the latter, children of GMFCS level I might be more homogenous in their possibilities and limitations in gross motor mobility during daily life compared to both children with GMFCS level II and children with level III. Furthermore, it has been described that children with GMFCS level I are more stable in motor function in comparison to children with levels II and III, and this stability is found to be already established by the age of four years. The following exclusion criteria were applied at baseline and during follow-up: (a) Botulinum neurotoxin type A (BoNT-A) treatment 10 months prior to the assessment (b) serial casting 3 months prior to assessment and (c) history of orthopedic and/or neurosurgery. To illustrate muscle deficits in the children with SCP, a retrospective dataset of cross-sectional assessments in age-matched TD children were selected from the established reference database at the CMAL. This database included TD children with no history of neurological, neuromuscular or orthopedic disorders. Criteria for selecting the TD children for the current study were: (a) aged between 6 months and 9 years old and (b) assessments collected with similar measurement equipment and protocol as applied for the SCP participants. The study protocol was approved by the Ethical Commission UZ/KU Leuven, Leuven, Belgium (S59945, S62187, S62645) and registered (NCT05197764). Written informed consent was obtained for all participants.

2.2. Measurements

The three-dimensional freehand ultrasound (3DfUS), combining a conventional 2D US device and a motion tracking system, was used. This technique has proven to be valid and reliable for measuring healthy and pathological muscles, already from early ages[13,28,29,30]. The US device settings were kept constant throughout the study period. Details on the 3DfUS equipment and measurement protocol are described in previous studies [19,20]. While positioned in prone, the most affected leg was measured in all SCP participants whereas a random leg was selected for TD children. The most affected leg was defined according to the MG spasticity and ankle joint range of motion (ROM) measured during a standardized clinical examination (Modified Ashworth Scale and goniometry, respectively) at the time of the baseline assessment [31]. The assessed leg remained the same for all follow-up assessments. Only 3 experienced researchers were involved in the assessments with only 2 researchers for the same participant. Of these 3 researchers, only 2 were also involved in the processing for which 1 only researcher performed the further data-analysis. The researchers were blinded for the GMFCS levels and timing of the assessment during the processing of the ultrasound images. All involved researchers were extensively trained in the CMAL research group following predefined guidelines and regular meetings to maximize similarities in the prescribed acquisition and processing workflow [32]. The inter-rater and intra-rater inter-session reliability was previously demonstrated with an intra-class correlation of minimal 0.95 [33] Furthermore, the Portico (i.e., concave-shaped holder combined with concave gel pad) was used to minimize inter-acquire differences in probe pressure [34]. STRADWIN software (Version 6.0, Mechanical Engineering, Cambridge, UK) was used for both data acquisition and processing. MV (in ml) was defined by manually drawing muscle segmentations alongside the inner muscle border, starting from the medial femoral condyle until the last image of the muscle belly before approaching the muscle tendon junction. Subsequently, a cubic planimetry technique was applied for interpolation of the defined segmentations [35]. The mid-belly segmented anatomical CSA was extracted and is further referred to as the CSA (in mm2). The ML of the MG muscle (in mm) was defined as the Euclidean distance between the muscle origin at the most superficial aspect of the medial femoral condyle and the muscle tendon junction. Accounting for inter-individual anthropometric variability and in line with previous studies, ratio scaling was applied for which the MV was normalized to the product of body mass and length (nMV, ml/(kg⋅m)), CSA to body mass (nCSA, mm2/kg) and ML to body length (nML, %) [14,36].

2.3. Statistical analysis

Participant characteristics and muscle outcomes at the baseline assessment were summarized by descriptive statistics reporting frequencies and median (interquartile ranges, IQ1-IQ3). Non-normality was confirmed with Shapiro-Wilk’s test. Two subgroups considering the motor impairments were used: subgroup 1, GMFCS level I and subgroup 2, GMFCS level II and III. To model the change in muscle morphology with respect to age, non-linear mixed models were used [37,38]. These models accommodate (a) unbalanced longitudinal datasets, resulting from variable spacing of follow-up intervals and different number of follow-up assessment among the children with SCP, and (b) the correlation between the repeated assessments taking the variance between and within the children with SCP into account. Fixed effects were (a) the age (expressed in years) at the time of assessments within each child, representing the time-effect and (b) covariate, representing GMFCS with its specific subgroups to describe within-group and compare between-group changes in the muscle growth. Random intercepts were added to model the variability in starting point between the children with SCP. To capture the changes in the trajectory of morphological muscle growth with age, piecewise models were allowed. The following workflow was used to construct these models: first, the average trend (i.e., mean structure) was explored by performing Loess regressions and by plotting the observed individual longitudinal trajectories. These explorations suggested piecewise trends, i.e., constatation of linear trends interrupted by breaking points. Second, residual trajectories calculated from the Loess regressions and the observed variance function (i.e., the change of the squared residuals over time) were explored, to define the random-effect structure. A random intercept was selected based on the exploration closely approaching a constant variance. Estimated starting values from the graphical explorations were used in the logistic piecewise models to obtain the predicted intercept, regression coefficients, breakpoints at specific ages, variance of the random intercept and measurement error.
The observed morphological outcomes included the MV, nMV, CSA, nCSA, ML and nML. The model with the jth observation in child i for GMFCS group g (GMFCS I=I and GMFCS II-III=II-III to estimate the response was described as follow:
Preprints 67218 i001
With α0 = intercept; a1i = random intercept; c1 = first breakpoint; c2 = second breakpoint; β1= regression coefficient of the slope before first breakpoint; β2= regression coefficient of the slope after first breakpoint and before second breakpoint; β3= regression coefficient of the slope after second breakpoint; ε(1)ijg = measurement error.
An F-test was used to test (a) if the slopes and the breakpoints differ from zero, and (b) if the slopes before and after the breakpoint differ from each other, confirming the non-linear trajectory. After formulating the piecewise regression model, so-called empirical Bayes estimates were calculated and used to assess the potential presence of outliers, i.e., patients with an “exceptional” starting point and evolution over time. Further, the F-test was used to compare the regression coefficients and breakpoints within and between the GMFCS groups. P-values were set at <0.05.
To enhance interpretation of the intercept, the age at all assessments was subtracted with 2 years, and therefore, the intercept represented the response’s value at the age of 2 years. The regression coefficients of the slope represent the magnitude of response in morphological muscle growth per year. The observed individual outcomes, observed individual trajectories, predicted individual trajectories and predicted average trajectory per morphological muscle outcome were plotted. The observed individual trajectories were also presented per GMFCS level, providing descriptive results per GMFCS level (Figure S1). The individual observed outcomes of 102 TD children were visualized with grey dots combined with boxplots for median (interquartile ranges) per 1-year age-groups, illustrating a cross-sectional reference dataset compared to the longitudinal trajectories of the SCP children. For the TD children, the median cross-sectional growth rate for a specific age range was calculated as the median of the muscle outcome divided by the median age. Illustrations of morphological growth trajectories in each GMFCS group of the children with SCP compared to a retrospective dataset of cross-sectional assessments in TD children were added in Figure S2 and Table S1. All analyses were performed in SAS® (Statistical Analysis Software version 9.4, SAS Institute Inc., Cary, NC, United States).

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 (c1). 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 c2, 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 c2). 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.4mm2/year, p<0.0001; β2=37.0 mm2/year, p<0.0001, respectively). The nCSA trajectory showed after an increase of 8.88 mm2/kg per year only 1 breakpoint at the age of 2.1y followed by decrease in nCSA per year (c1 and β2= -0.96 mm2/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.

5. Conclusions

To our knowledge, this was the first study that performed repeated muscle assessments over 2 years in an extended group of growing children with SCP and demonstrated piecewise model for MG muscle growth. The study sample is believed to be representative for the population of children receiving standardized clinical care in Belgium (Table S2). After high rates of muscle growth during first 2 years of life, the trajectory of absolute MV and CSA changed, resulting in a slower gradual growth until 6-9 years, and subsequent reduction in muscle size growth after 9 years. With increasing age, normalized MV and mid-belly CSA represented a decrease in muscle morphology after approximately 2 years of age. Children with GMFCS level II-III showed more muscle growth reduction compared to children with level I. This lower growth rate started already at a very young age. After infancy until approximately 9 years of age, the trajectory of muscle growth did not much differ between GMFCS levels. These longitudinal trajectories could be used to monitor the SCP muscle pathology during childhood and to optimize treatment planning and goals aiming to stimulate muscle growth.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Figure S1: Observed individual profiles for the (A) absolute and (B) normalized muscle morphology for (i) the total group of children with SCP (blue=GMFCS level I, dark red=level II and pink=level III), (ii) children with GMFCS level I, (iii) children with GMFCS level II and (iv) children with GMFCS level III. Figure S2: A Reference dataset for absolute morphological muscle outcomes and B: Reference dataset for normalized morphological muscle outcomes, Table S1: Muscle morphology of retrospective cross-sectional dataset with TD children. Table S2: Standard clinical care of all SCP participants. Table S3: A Fixed and random effects of the piecewise regressions for muscle morphology (GMFCS I n=47 and GMFCS II-III n=40).& B Results of the differences between slopes and breakpoints (GMFCS, I n=47, and GMFCS II-III n=40). Figure S3: Observed individual trajectories for the (i) absolute and (ii) normalized muscle morphology for children with unilateral SCP (n=51) (grey) and children with bilateral SCP (n=36) (green).

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.

Funding

ND was funded by an internal KU Leuven grant, Belgium, 3D-MMAP-C24/18/103. This work was further supported by Internal funding of KU Leuven Biomedical Sciences Group: Fund for Translational Biomedical Research, Belgium, 2019 and the Research Foundation Flanders (FWO), Belgium, G0B4619N. BH and NP were funded by a TBM grant (TAMTA-T005416N). IV was funded through a research fellowship from FWO (1188921N).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by Ethical 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

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors via an online repository, without undue reservation.

Acknowledgments

The authors wish to thank Heidi Devolder and Julie Uytterhoeven for assistance in patient recruitment and planning of the follow-up appointments. We would like to thank Jorieke Deschrevel for the support during the data acquisition and Tijl Dewit for the support of 3DfUS measures and database management.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Graham, H. K.; Rosenbaum, P.; Paneth, N.; Dan, B.; Lin, J.-P.; Damiano, Di. L.; Becher, J. G.; Gaebler-Spira, D.; Colver, A.; Reddihough, Di. S.; Crompton, K. E.; Lieber, R. L. Cerebral Palsy. Nat. Rev. Dis. Prim. 2016, 2(1), 15082. [Google Scholar] [CrossRef]
  2. Barrett, R. S.; Lichtwark, G. A. Gross Muscle Morphology and Structure in Spastic Cerebral Palsy: A Systematic Review. Dev. Med. Child Neurol. 2010, 52(9), 794–804. [Google Scholar] [CrossRef]
  3. Schless, S. H.; Cenni, F.; Bar-On, L.; Hanssen, B.; Goudriaan, M.; Papageorgiou, E.; Aertbeliën, E.; Molenaers, G.; Desloovere, K. Combining Muscle Morphology and Neuromotor Symptoms to Explain Abnormal Gait at the Ankle Joint Level in Cerebral Palsy. Gait Posture 2019, 68(December 2018), 531–537. [Google Scholar] [CrossRef]
  4. Rosenbaum, P. L.; Walter, S. D.; Hanna, S. E.; Palisano, R. J.; Russell, D. J.; Raina, P.; Wood, E.; Bartlett, D. J.; Galuppi, B. E. Prognosis for Gross Motor Function in Cerebral Palsy Creation of Motor Development Curves. JAMA 2002, 208(11), 1357–1363. [Google Scholar] [CrossRef]
  5. Nordmark, E.; Hägglund, G.; Lauge-Pedersen, H.; Wagner, P.; Westbom, L. Development of Lower Limb Range of Motion from Early Childhood to Adolescence in Cerebral Palsy: A Population-Based Study. BMC Med. 2009, 7. [Google Scholar] [CrossRef]
  6. Lindén, O.; Hägglund, G.; Rodby-Bousquet, E.; Wagner, P. The Development of Spasticity with Age in 4,162 Children with Cerebral Palsy: A Register-Based Prospective Cohort Study. Acta Orthop. 2019, 90(3), 286–291. [Google Scholar] [CrossRef]
  7. Burgess, A.; Reedman, S.; Chatfield, M. D.; Ware, R. S.; Sakzewski, L.; Boyd, R. N. Development of Gross Motor Capacity and Mobility Performance in Children with Cerebral Palsy: A Longitudinal Study. Dev. Med. Child Neurol. 2022, 64(5), 578–585. [Google Scholar] [CrossRef]
  8. Palisano, R. J.; Rosenbaum, P.; Bartlett, D.; Livingston, M. H. Content Validity of the Expanded and Revised Gross Motor Function Classification System. Dev. Med. Child Neurol. 2008, 50(10), 744–750. [Google Scholar] [CrossRef]
  9. Molenaers, G.; Fagard, K.; Van Campenhout, A.; Desloovere, K. Botulinum Toxin A Treatment of the Lower Extremities in Children with Cerebral Palsy. J.Child Orthop 2013, 7, 383–387. [Google Scholar] [CrossRef]
  10. Kruse, A.; Rivares, C.; Weide, G.; Tilp, M.; Jaspers, R. T. Stimuli for Adaptations in Muscle Length and the Length Range of Active Force Exertion—A Narrative Review. Front. Physiol. 2021, 12(October). [Google Scholar] [CrossRef] [PubMed]
  11. Verschuren, O.; Smorenburg, A. R. P.; Luiking, Y.; Bell, K.; Barber, L.; Peterson, M. D. Determinants of Muscle Preservation in Individuals with Cerebral Palsy across the Lifespan: A Narrative Review of the Literature. J. Cachexia. Sarcopenia Muscle 2018, 9(3), 453–464. [Google Scholar] [CrossRef]
  12. Gough, M.; Shortland, A. P. Early Muscle Development in Children with Cerebral Palsy. the Consequences for Further Muscle Growth, Muscle Function, and Long-Term Mobility; Roberta Shepherd, 2013. [CrossRef]
  13. 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(5), 485–491. [Google Scholar] [CrossRef]
  14. De Beukelaer, N. De; Weide, G.; Huyghe, E.; Vandekerckhove, I.; Hanssen, B.; Peeters, N.; Uytterhoeven, J.; Deschrevel, J.; Maes, K.; Corvelyn, M.; Costamagna, D.; Gayan-Ramirez, G.; Van Campenhout, A.; Desloovere, K. Reduced Cross-Sectional Muscle Growth Six Months after Botulinum Toxin Type-A Injection in Children with Spastic Cerebral Palsy. Toxins (Basel). 2022, 14(139). [Google Scholar] [CrossRef]
  15. 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(6), 657–663. [Google Scholar] [CrossRef]
  16. Barber, L.; Hastings-Ison, T.; Baker, R.; Barrett, R.; Lichtwark, G. Medial Gastrocnemius Muscle Volume and Fascicle Length in Children Aged 2 to 5 Years with Cerebral Palsy. Dev. Med. Child Neurol. 2011, 53(6), 543–548. [Google Scholar] [CrossRef]
  17. Oberhofer, K.; Stott, N. S.; Mithraratne, K.; Anderson, I. A. Subject-Specific Modelling of Lower Limb Muscles in Children with Cerebral Palsy. Clin. Biomech. 2010, 25(1), 88–94. [Google Scholar] [CrossRef]
  18. Pitcher, C. A.; Elliott, C. M.; Valentine, J. P.; Stannage, K.; Williams, S. A.; Shipman, P. J.; Reid, S. L. Muscle Morphology of the Lower Leg in Ambulant Children with Spastic Cerebral Palsy. Muscle and Nerve 2018, 58(6), 818–823. [Google Scholar] [CrossRef]
  19. Hanssen, B.; Peeters, N.; Vandekerckhove, I.; De Beukelaer, N.; Bar-On, L.; Molenaers, G.; Van Campenhout, A.; Degelaen, M.; Van den Broeck, C.; Calders, P.; Desloovere, K. The Contribution of Decreased Muscle Size to Muscle Weakness in Children With Spastic Cerebral Palsy. Front. Neurol. 2021, 12(July), 1–14. [Google Scholar] [CrossRef]
  20. Schless, S. H.; 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. 2018, 61(7), 783–790. [Google Scholar] [CrossRef]
  21. Noble, J. J.; Chruscikowski, E.; Fry, N. R. D.; Lewis, A. P.; Gough, M.; Shortland, A. P. The Relationship between Lower Limb Muscle Volume and Body Mass in Ambulant Individuals with Bilateral Cerebral Palsy. BMC Neurol. 2017, 17(1), 1–9. [Google Scholar] [CrossRef]
  22. Lieber, R. L.; Steinman, S.; Barash, I. A.; Chambers, H. Structural and Functional Changes in Spastic Skeletal Muscle. Muscle Nerve 2004, 29(5), 615–627. [Google Scholar] [CrossRef]
  23. Gundersen, K. Excitation-Transcription Coupling in Skeletal Muscle: The Molecular Pathways of Exercise. Biol. Rev. Camb. Philos. Soc. 2011, 86(3), 564. [Google Scholar] [CrossRef]
  24. Clowry, G. J. The Dependence of Spinal Cord Development on Corticospinal Input and Its Significance in Understanding and Treating Spastic Cerebral Palsy. Neuroscience and Biobehavioral Reviews. Pergamon January 1, 2007, pp 1114–1124. [CrossRef]
  25. Noble, J. J.; Fry, N. R.; Lewis, A. P.; Keevil, S. F.; Gough, M.; Shortland, A. P. Lower Limb Muscle Volumes in Bilateral Spastic Cerebral Palsy. 2014. [CrossRef]
  26. Gorter, J. W.; Tol, E. Van; Schie, P. Van; Ketelaar, M.; Palisano, R.; Rosenbaum, P.; Bartlett, D.; Livingston, M.; Palisano, R.; Rosenbaum, P.; Walter, S.; Russell, D.; Wood, E.; Galuppi, B. GMFCS – E & R Gross Motor Function Classification System Expanded and Revised. 2009, 2007–2010.
  27. Palisano, R. J.; Cameron, D.; Rosenbaum, P. L.; Walter, S. D.; Russell, D. Stability of the Gross Motor Function Classification System. Dev. Med. Child Neurol. 2006, 48(6), 424–428. [Google Scholar] [CrossRef]
  28. 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]
  29. Barber, L.; Alexander, C.; Shipman, P.; Boyd, R.; Reid, S.; Elliott, C. Validity and Reliability of a Freehand 3D Ultrasound System for the Determination of Triceps Surae Muscle Volume in Children with Cerebral Palsy. J. Anat. 2019, 234(3), 384–391. [Google Scholar] [CrossRef]
  30. Williams, S. A.; Bell, M.; Kim, H. K.; Salim Al Masruri, G.; Stott, N. S.; Fernandez, J.; Mirjalili, S. A. The Reliability and Validity of Triceps Surae Muscle Volume Assessment Using Freehand Three-dimensional Ultrasound in Typically Developing Infants. J. Anat. 2021, (No. July), 1–12. [Google Scholar] [CrossRef]
  31. Thomason, P.; Willoughby, K.; Graham, H. K. Orthopaedic Assessment. In Cerebral palsy: science and clinical practice; Mac Keith Press: London, UK, 2014; pp. 287–312. [Google Scholar]
  32. 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(9), 2702–2712. [Google Scholar] [CrossRef]
  33. Hanssen, B. Muscle Structure, Weakness, and Strengthening: A Complex Triangle for Children with Spastic Cerebral Palsy, KU Leuven, 2022.
  34. Cenni, F.; Schless, S. H.; 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]
  35. Treece, G.; Prager, R.; Gee, A.; Berman, L. Fast Surface and Volume Estimation from Non-Parallel Cross-Sections, for Freehand Three-Dimensional Ultrasound. Med. Image Anal. 1999, 3(2), 141–173. [Google Scholar] [CrossRef]
  36. Handsfield, G. G.; Meyer, C. H.; Hart, J. M.; Abel, M. F.; Blemker, S. S. Relationships of 35 Lower Limb Muscles to Height and Body Mass Quantified Using MRI. J. Biomech. 2014, 47(3), 631–638. [Google Scholar] [CrossRef]
  37. Verbeke, G.; Molenberghs, G. Linear Mixed Models in Practice: A SAS-Oriented Approach; Springer: New York, 1997. [Google Scholar]
  38. Molenberghs, G.; Verbeke, G. An Introduction to (Generalized) (Non-)Linear Mixed Models. In Explantory im respons models: A generalized linear and nonlinear approach; De Boeck, P., Wilson, M., Eds.; Springer-Verlag: New York, 2004; pp 111–153.
  39. Jaric, S.; Mirkov, D.; Markovic, G. Normalizing Physical Performance Tests for Body Size: A Proposal for Standardization. J. Strength Cond. Res. 2005, 19(2), 467–474. [Google Scholar] [CrossRef]
  40. Wells, J. C. K.; Davies, P. S. W.; Fewtrell, M. S.; Cole, T. J. Body Composition Reference Charts for UK Infants and Children Aged 6 Weeks to 5 Years Based on Measurement of Total Body Water by Isotope Dilution. Eur. J. Clin. Nutr. 2020, 74(1), 141–148. [Google Scholar] [CrossRef] [PubMed]
  41. 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(5), 425–433. [Google Scholar] [CrossRef]
  42. Day, S. M.; Strauss, D. J.; Vachon, P. J.; Rosenbloom, L.; Shavelle, R. M.; Wu, Y. W. Growth Patterns in a Population of Children and Adolescents with Cerebral Palsy. Dev. Med. Child Neurol. 2007, 49(3), 167–171. [Google Scholar] [CrossRef] [PubMed]
  43. de las Mercedes Ruiz Brunner, M.; Cuestas, E.; Heinen, F.; Schroeder, A. S. Growth in Infants, Children and Adolescents with Unilateral and Bilateral Cerebral Palsy. Sci. Rep. 2022, 12(1). [Google Scholar] [CrossRef]
  44. Bénard, M. R.; Harlaar, J.; Becher, J. G.; Huijing, P. A.; Jaspers, R. T. Effects of Growth on Geometry of Gastrocnemius Muscle in Children: A Three-Dimensional Ultrasound Analysis. J. Anat. 2011, 219(3), 388–402. [Google Scholar] [CrossRef] [PubMed]
  45. Lieber, R. L.; Fridén, J. Functional and Clinical Significance of Skeletal Muscle Architecture. Muscle Nerve 2000, 23(November), 1647–1666. [Google Scholar]
  46. 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: A Registry Analysis of 2,796 Examinations in 355 Children. Acta Orthop. 2011, 82(6), 744–748. [Google Scholar] [CrossRef]
  47. Mathewson, M. A.; Lieber, R. L. Pathophysiology of Muscle Contractures in Cerebral Palsy. Phys. Med. Rehabil. Clin. N. Am. 2015, 26(1), 57–67. [Google Scholar] [CrossRef]
  48. Williams, I.; Reid, L.; Stott, N. S.; Valentine, J.; Elliott, C. Measuring Skeletal Muscle Morphology and Architecture with Imaging Modalities in Children with Cerebral Palsy: A Scoping Review. Dev. Med. Child Neurol. 2020, 63(3), 263–273. [Google Scholar] [CrossRef]
  49. Park, E.-Y.; Kim, W.-H. Meta-Analysis of the Effect of Strengthening Interventions in Individuals with Cerebral Palsy. Res. Dev. Disabil. 2014, 35, 239–249. [Google Scholar] [CrossRef] [PubMed]
  50. Anker–van der Wel, I.; Smorenburg, A. R. P.; de Roos, N. M.; Verschuren, O. Dose, Timing, and Source of Protein Intake of Young People with Spastic Cerebral Palsy. Disabil. Rehabil. 2020, 42(15), 2192–2197. [Google Scholar] [CrossRef] [PubMed]
  51. Hanssen, B.; Peeters, N.; De Beukelaer, N.; Vannerom, A.; Peeters, L.; Molenaers, G.; Van Campenhout, A.; Deschepper, E.; Van den Broeck, C.; Desloovere, K. Progressive Resistance Training for Children with Cerebral Palsy: A Randomized Controlled Trial Evaluating the Effects on Muscle Strength and Morphology. Frontiers in Physiology. 2022. [CrossRef] [PubMed]
  52. 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(11), 1165–1171. [Google Scholar] [CrossRef] [PubMed]
  53. Damiano, D.; Abel, M.; Romness, M.; Oeffinger, D.; Tylkowski, C.; Gorton, G.; Bagley, A.; Nicholson, D.; Barnes, D.; Calmes, J.; Kryscio, R.; Rogers, S. Comparing Functional Profiles of Children with Hemiplegic and Diplegic Cerebral Palsy in GMFCS Levels I and II: Are Separate Classifications Needed? Dev. Med. Child Neurol. 2006, 48(10), 797–803. [Google Scholar] [CrossRef]
  54. Barber, L. A.; Read, F.; Lovatt Stern, J.; Lichtwark, G.; Boyd, R. N. Medial Gastrocnemius Muscle Volume in Ambulant Children with Unilateral and Bilateral Cerebral Palsy Aged 2 to 9 Years. Dev. Med. Child Neurol. 2016, 58(11). [Google Scholar] [CrossRef]
Figure 1. Representation of the number of the performed longitudinal assessment at the different ages for the group of children with spastic cerebral palsy (SCP) (n=87) of gross motor function classification system level I (blue) versus levels II-III (red). The x-as represents the ages at the time of the repeated assessment. On the y-as, each child is represented on one line, with each dot representing the performed assessments. The children were further ordered with increasing age. The group of GMFCS level I included 47 children with total of 130 assessments (average (range) of 2.8 (2-5) assessments per child and average time of 10 months (4-25 months) between the assessments. The group of GMFCS level II-III included 40 children with total of 104 assessments (average (range) of 2.6 (2-5) assessments per child and average time of 9 months (4-25 months) between the assessments.
Figure 1. Representation of the number of the performed longitudinal assessment at the different ages for the group of children with spastic cerebral palsy (SCP) (n=87) of gross motor function classification system level I (blue) versus levels II-III (red). The x-as represents the ages at the time of the repeated assessment. On the y-as, each child is represented on one line, with each dot representing the performed assessments. The children were further ordered with increasing age. The group of GMFCS level I included 47 children with total of 130 assessments (average (range) of 2.8 (2-5) assessments per child and average time of 10 months (4-25 months) between the assessments. The group of GMFCS level II-III included 40 children with total of 104 assessments (average (range) of 2.6 (2-5) assessments per child and average time of 9 months (4-25 months) between the assessments.
Preprints 67218 g001
Figure 2. Predicted average trajectories (bold solid line) for the (i) absolute and (ii) normalized muscle morphology in SCP children with GMFCS level I (blue) and GMFCS level II-III (red). Observed outcomes (dots) and individual predicted trajectories (lines) are presented. SCP, spastic cerebral palsy; GMFCS, gross motor function classification system; MV, muscle volume; ml, milliliter; CSA, anatomical cross-sectional area; mm, millimeter; ML, muscle belly length; n, normalized; kg, kilogram; m, meters.
Figure 2. Predicted average trajectories (bold solid line) for the (i) absolute and (ii) normalized muscle morphology in SCP children with GMFCS level I (blue) and GMFCS level II-III (red). Observed outcomes (dots) and individual predicted trajectories (lines) are presented. SCP, spastic cerebral palsy; GMFCS, gross motor function classification system; MV, muscle volume; ml, milliliter; CSA, anatomical cross-sectional area; mm, millimeter; ML, muscle belly length; n, normalized; kg, kilogram; m, meters.
Preprints 67218 g002
Table 1. Participant characteristics of the children with SCP at baseline assessment.
Table 1. Participant characteristics of the children with SCP at baseline assessment.
Group GMFCS I GMFCS II-III
Participants 47 40
Observations 130 104
Sex (boy/girls) 25 / 22 21 / 19
GMFCS level (I /II/ III) 47 / 0 / 0 0 / 22 / 18
Involvement (unilateral/bilateral) 37 / 10 14 (II, =10 & III=4) / 26 (II=12 & III=14)
History of BoNT-A injection 11 7
Antropometric & muscle outcomes
Age (y) 4.66 (2.25-6.83) 2.88 (1.47-4.81)
Body mass (kg) 18.0 (13.0-22.0) 12.5 (10.0-16.4)
Body length (cm) 104.7 (88.5-121.1) 90.0 (77.8-106.0)
MV (ml) 28.2 (16.7-38.4) 13.6 (7.4-25.5)
CSA (mm2) 335.6 (260.1-392.1) 205.5 (162.3-285.8)
ML (mm) 127.1 (99.2-144.9) 96.0 (82.8-127.1)
nMV (ml/kg⋅m) 1.45 (1.26-1.58) 1.12 (0.92-1.39)
nCSA (mm2/kg) 19.4 (16.1-22.1) 16.1 (13.8-19.0)
nML (%) 11.7 (11.3-12.6) 10.9 (9.8-11.9)
The frequencies are presented for the general characteristics. Anthropometric and muscle data are presented as median (interquartile 1 – interquartile 3). SCP, spastic cerebral palsy; GMFCS, gross motor function classification system; BoNT-A, botulinum neurotoxin type A; y, years; kg, kilogram; cm, centimeter; MV, muscle volume; ml, milliliter; CSA, anatomical cross-sectional area; mm, millimeter; ML, muscle belly length; n, normalized; kg, kilogram; m, meter.
Table 2. Fixed effects of the piecewise model for muscle morphology in GMFCS subgroups (level I, n=47; levels II-III, n=40).
Table 2. Fixed effects of the piecewise model for muscle morphology in GMFCS subgroups (level I, n=47; levels II-III, n=40).
Preprints 67218 i002Preprints 67218 i003
The asterisks (*) indicate significance level at p<0.05 for α and c. P-values in bold indicate significance level at p<0.05. The following symbols represent: α = outcome at the age of 2 years, β = rate of change in muscle outcome per year, and c = age (years) of the breakpoint. The numbers (in subscript) associated with symbols β and c refer to the order of the observed growth ratio and breakpoint, respectively, e.g. β1 as the first ratio before the c1, the first breakpoint. SCP, spastic cerebral palsy; GMFCS, gross motor function classification system; n, number; CI 95% confidence interval; MV, muscle volume; n, normalized; CSA, anatomical cross-sectional area; ML, muscle belly length.
Table 3. Results of the comparison in slopes and breakpoints between the GMFCS subgroups.
Table 3. Results of the comparison in slopes and breakpoints between the GMFCS subgroups.
Outcome GMFCS I vs II- III Δ p-value
MV β1, I vs. β1, II-III 8.24 0.0003
β2, I vs. β1, II-III 1.16 0.0712
β 3, I vs. β2, II-III 3.91 0.1069
c2, I vs c1, II-III 1.31 0.0326
nMV β1, I vs. β1, II-II 0.55 0.0007
β2, I vs. β2 II-III 0.04 0.1710
β 3, I vs. β2, II-III 0.05 0.1272
c1, I vs c1, II-III 0.58 0.0040
α1, I vs α1, II-III 0.32 0.0056
CSA β1, I vs. β1, II-III 69.8 0.0112
β2, I vs. β2, II-III 7.17 0.4328
β3, I vs. β3, II-III 35.2 0.0915
c1, I vs c1, II-III 0.09 0.3783
c2, I vs c2, II-III 2.66 0.0033
nCSA β1, I vs. β1, II-II 6.18 0.0039
β2, I vs. β2, II-III -0.09 0.8168
β2, I vs β3, II-III 1.56 0.0334
c1, I vs c1, II-III 0.61 0.1152
α0, I vs α0, II-III 3.97 0.0063
ML β1, I vs. β1, II-II -4.61 0.2482
β2, I vs. β2, II-III -1.79 0.1384
c1 vs c1 -3.01 <0.001
nML β1, I vs. β1, II-II 0.21 0.1311
β2, I vs. β1, III-II -0.22 0.0146
α0, I vs α0., II-III 0.91 0.0024
P-values in bold indicate significance level at p<0.05. The following symbols represent: α = outcome at the age of 2 years, β = rate of change in muscle outcome per year, c = age (years) of the breakpoint, Δ = difference score, I = GMFCS level I and II-III = GMFCS level II & III. The numbers (in subscript) associated with symbols β and c refer to the order of the observed growth ratio and breakpoint, respectively, e.g. β1 as the first ratio before the c1, the first breakpoint.; GMFCS, gross motor function classification system; ; MV, muscle volume; n, normalized; CSA, anatomical cross-sectional area; ML, muscle belly length.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated