Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 26 February 2021
Sec. Pediatric Neurology
This article is part of the Research Topic Cerebral Palsy: New Developments View all 17 articles

Muscle Characteristics in Pediatric Hereditary Spastic Paraplegia vs. Bilateral Spastic Cerebral Palsy: An Exploratory Study

  • 1KU Leuven Department of Rehabilitation Sciences, Leuven, Belgium
  • 2Clinical Motion Analysis Laboratory, University Hospitals Leuven, Leuven, Belgium
  • 3Department of Rehabilitation Medicine, Amsterdam University Medical Center, Amsterdam, Netherlands
  • 4Department of Orthopedics, University Hospitals Leuven, Leuven, Belgium
  • 5KU Leuven Department of Development and Regeneration, Leuven, Belgium

Hereditary spastic paraplegia (HSP) is a neurological, genetic disorder that predominantly presents with lower limb spasticity and muscle weakness. Pediatric pure HSP types with infancy or childhood symptom onset resemble in clinical presentation to children with bilateral spastic cerebral palsy (SCP). Hence, treatment approaches in these patient groups are analogous. Altered muscle characteristics, including reduced medial gastrocnemius (MG) muscle growth and hyperreflexia have been quantified in children with SCP, using 3D-freehand ultrasound (3DfUS) and instrumented assessments of hyperreflexia, respectively. However, these muscle data have not yet been studied in children with HSP. Therefore, we aimed to explore these MG muscle characteristics in HSP and to test the hypothesis that these data differ from those of children with SCP and typically developing (TD) children. A total of 41 children were retrospectively enrolled including (1) nine children with HSP (ages of 9–17 years with gross motor function levels I and II), (2) 17 age-and severity-matched SCP children, and (3) 15 age-matched typically developing children (TD). Clinically, children with HSP showed significantly increased presence and severity of ankle clonus compared with SCP (p = 0.009). Compared with TD, both HSP and SCP had significantly smaller MG muscle volume normalized to body mass (p ≤ 0.001). Hyperreflexia did not significantly differ between the HSP and SCP group. In addition to the observed pathological muscle activity for both the low-velocity and the change in high-velocity and low-velocity stretches in the two groups, children with HSP tended to present higher muscle activity in response to increased stretch velocity compared with those with SCP. This exploratory study is the first to reveal MG muscle volume deficits in children with HSP. Moreover, high-velocity-dependent hyperreflexia and ankle clonus is observed in children with HSP. Instrumented impairment assessments suggested similar altered MG muscle characteristics in pure HSP type with pediatric onset compared to bilateral SCP. This finding needs to be confirmed in larger sample sizes. Hence, the study results might indicate analogous treatment approaches in these two patient groups.

Introduction

Hereditary spastic paraplegia (HSP) is described as a genetic, heterogeneous disorder leading to axonal degeneration of the spinal pathways and mostly classified in either pure or complex (i.e., additional symptoms rather than pyramidal signs) types (1, 2). In children with the pure HSP type presenting symptoms already from infancy or childhood, bilateral spasticity and muscle weakness in the lower limbs contribute to gait deviations (3, 4). This predominant clinical picture is similarly observed in children with bilateral spastic cerebral palsy (SCP) (2, 4). Hence, the symptomatic treatment management in HSP is often analogous to the common treatment in SCP. Conservative treatment includes orthotic management and regular physiotherapy, which incorporates strengthening, stretching, and gait rehabilitation. This approach is commonly combined with spasticity management [i.e., oral medication, intrathecal baclofen therapy, and/or intramuscular botulinum neurotoxin (BTX) injections; (57)].

Caused by upper motor neuron lesions in the developing brain, children with SCP present with neuromotor symptoms already within the first 2 years of life. Therefore, these children have an atypical development resulting in altered muscle function and pathological gait (5). In the last decade, there has been growing research interest in the muscle morphology of children with SCP, with studies showing evidence for decreased muscle volume (MV) and muscle belly length of the medial gastrocnemius (MG) muscle in comparison with typically developing (TD) children (8). The calf muscle is most often selected as a target muscle because of its pathological involvement in the clinical presentation and consequently functional importance (9). Muscle morphology is considered related to muscle function in the sense that decreased MV has been indicated as one of the major contributors to muscle weakness (10, 11). Moreover, altered MG muscle morphology combined with the manifestation of lower limb spasticity has been observed already from early ages (12, 13).

Muscle imaging using 3D-freehand ultrasound (3DfUS) is a reliable and valid method to determine macroscopic alterations of muscle morphology in children with SCP (1416). In addition, instrumented assessment of hyperreflexia provides reliable, objective parameters to evaluate the stretch of hyperreflexia by combining surface electromyography (sEMG) signals with joint motion analysis during manually applied muscle stretches at low and high velocity (17, 18). This quantitative approach enables the assessment of the neural component of muscle hyper-resistance, which is further indicated as hyperreflexia (19). Moreover, using such instrumentation, different muscle activation patterns in reaction to increased stretch velocity have been reported (20). Treatment response may differ between these phenotypes, highlighting the importance to use instrumented impairment assessment for treatment guidance (21, 22).

In more than half of the children who are clinically suspected of having HSP, the genetic diagnosis is not confirmed by means of whole exome sequencing or HSP gene panels (23). The latter considers the pathogenic mutations in the spastic paraplegia genes (SPG) known to date (24). Moreover, since the clinical picture is often suggestive of SCP, other markers for a good differential diagnosis are a priority (2). Previous research has already pointed out possible indications for the HSP condition based on the gait deviations, whereby these results could contribute to the diagnostic process of children with HSP (25, 26). In the same sense, an objective and quantitative investigation of both hyperreflexia and potential alterations in muscle morphology in children with HSP might help to differentiate them from SCP.

Whereas, the clinical presentation indicates an overlap in children with HSP and SCP, it is unclear if similarities are also presented in terms of muscle characteristics. These insights might help to distinguish HSP and SCP and, thus, potentially provide markers for more accurate diagnosis, which may support the treatment management in the HSP population and optimize treatment plans toward delineated patient-specific approaches. However, to date, muscle characteristics in children with HSP have not yet been studied.

Therefore, the first aim of this retrospective study is to explore the MG muscle characteristics in children with HSP using instrumented impairment assessments. Second, we tested the hypothesis that the muscle morphology and hyperreflexia in children with HSP differ from those of children with SCP.

Materials and Methods

Participants

This study retrospectively selected 41 participants who were recruited via the Clinical Motion Analysis Laboratory in collaboration with the pediatric–orthopedic department of the University Hospitals Leuven. Children underwent one or more assessments including clinical examination, 3DfUS and instrumented assessment of hyperreflexia. Written informed consent was signed by all parents of the children. The study assessments were performed as part of different ongoing research projects and were approved by the Ethical Committee of the University Hospitals Leuven and KU Leuven (S56977 and S59945).

In the previous described research projects, children with a genetically confirmed diagnosis of HSP were eligible for inclusion. In the absence of a genetic diagnosis after performing extensive genetic workup, a clinical HSP diagnosis was made by qualified professionals based on (a) neurological findings, such as spasticity and weakness in the plantar flexor muscles, (b) the age of onset, and if available, (c) a positive family history of gait disturbances, and (d) normal neuro-imaging, electromyography measures, and metabolic investigations [Supplementary Table 1; (1, 2, 27)]. The latter is performed to ensure the exclusion of other disorders (2, 28). As a result, from the population of children with HSP who received regular clinical follow-up between May 2017 and July 2020 at the University Hospitals Leuven, 17 patients were assessed as part of ongoing research projects (S56977 and S59945).

For this retrospective study, the following inclusion criteria for children with either clinically or genetically confirmed HSP were applied: (a) aged between 8 and 17 years; (b) level I or II on the gross motor function classification system (GMFCS); (c) uncomplicated or pure type of HSP; (d) no BTX injections in the previous 6 months; (e) no orthopedic surgery in the previous 2 years; (f) no history of soft-tissue musculoskeletal- or neurosurgery (e.g., selective dorsal rhizotomy), and (g) no intrathecal baclofen pump. This resulted in a group of nine children with HSP, including six boys and three girls. Mutations in SP-designated genes were confirmed in four patients (two children with SPG3A, one with SPG4, and one with SPG56). For the other five participants, genetic analysis could not indicate mutations in the SP-designated genes. The flowchart of the HSP patient selection and inclusion process is presented in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Flowchart of the patient selection process for the children with hereditary spastic paraplegia (HSP). N, number.

Based on the age, GMFCS levels and gender of this HSP group, a cohort of children with bilateral SCP were group matched. These children were selected from a retrospective database (S59945) that was established during an ongoing research project on muscle characteristics in children with SCP. To aim for a maximum number of children that were eligible from the retrospective database, group matching was performed per available assessment. As a result, 15 children with SCP were included with 3DfUS data, and nine children with SCP were included with hyperreflexia data. Seven children with SCP were measured for both assessments. Supplementary Figure 1 gives a schematic overview of the study samples per assessment. To confirm the diagnosis of bilateral SCP, these children showed abnormal neuro-imaging and the clinical presentation of plantar flexor muscle spasticity.

The muscle morphology of these two patient groups was compared with a reference group of 15 TD children, which was retrospectively included based on the available 3DfUS data in the Clinical Motion Analysis Laboratory database (S59945). This cohort was group matched following the age and gender of the HSP group. No reference of instrumently assessed hyperreflexia data were used for the TD group, since the muscle activation during passive muscle stretches in TD children is negligible (17).

Data Collection and Processing

Assessments were carried out by the different experienced examiners who were involved in the research projects, whereas all data processing was performed by the same experienced researcher (ND). As defined by the clinical examination, the most affected leg was measured in the children with HSP and SCP, whereas a random leg was selected in the TD children. Anthropometric data including body mass, body length, and leg lengths were measured. The latter was measured from the lower border of the anterior superior iliac spine to the lower border of the medial malleolus.

Clinical Assessments

Clinical examination was performed including the passive range of motion to maximal ankle dorsiflexion (ROM), plantarflexor spasticity assessed with the Modified Ashworth Scale (MAS) and Modified Tardieu Angle (MTA), and plantarflexor strength using the Medical Research Council grade scale (2931). ROM, MAS, and MTA were measured with the knee extended. The presence and severity of the clonus was clinically graded based on how quickly a rhythmic, oscillating reflex appears, and disappears while stretching the ankle towards the dorsiflexion position. The clonus grading is as follows: 0, no clonus; 1, clonus after multiple stretches; 2, quick clonus, quick stop; 3, quick clonus, slow stop, and 4, quick clonus, no stop. Clinical scores including the Functional Mobility Scale (FMS), walking ability, self-selected gait speed, presence of pain during gait, and bladder and bowel function were retained from the clinical assessments, which were performed at the time of or within 1 year of the study assessment, as part of the routine clinical follow-up. The FMS describes the functional mobility by scoring the walking ability over three distances (i.e., 5–50–500 m) on a six-point scale for which lower scores indicate more need for assistance (32).

3D-Freehand Ultrasound

To aim for a comfortable and resting joint position, the children were positioned in prone with a triangular cushion placed under the lower leg whereby the ankle was placed over the edge of this cushion. This resting position provided ~25° of knee flexion and 30° of plantar flexion. Muscle morphology was assessed using 3DfUS, combining a 2D ultrasound system and motion-tracking system of three fixed optical cameras, with the same instrumentation settings as described by Schless et al. (33). US images of the MG muscle were acquired in a transverse orientation, starting from the medial femoral condyle to the distal end of the calcaneus. Both data collection and processing were completed using STRADWIN software (version 6.0; Mechanical Engineering, Cambridge University, Cambridge, UK). The following processing guidelines to assign landmarks in the transverse plane were applied: (a) MG muscle origin was defined as the most superficial part of the medial femoral condyle, (b) the muscle–tendon junction (MTJ) was defined as the first image after the last image with visible muscle belly mass, and (c) the first image of the calcaneus bone was defined as the tendon insertion. Data processing started with segmenting the created 3D-reconstruction by manually outlining the cross-sectional areas along the inside of the muscle borders starting from the defined origin to the last image before the MTJ. A linear interpolation between the outlined borders was applied to compute the MV (in ml). In addition, MV was normalized to body mass (ml/kg). Muscle belly lengths (ML) were defined from the origin to the MTJ, and the tendon lengths (TL) were defined from the MTJ to the tendon insertion. The distance between the muscle origin and insertion was specified as the muscle–tendon unit complex length (MTUL). The same assessor assigned the landmarks twice, and the average of the length scores was calculated. ML, TL, and MTUL data were reported in millimeter (mm). Both ML and MTU L were normalized for leg length (i.e., nML and nMTUL) and therefore expressed as percentages.

Instrumented Assessment of Hyperreflexia

A previously described measurement protocol was applied in which muscle activation and joint motion analysis data were simultaneously collected during manually applied ankle rotations to stretch the MG muscle at low (LV) and high (HV) velocity (17). Muscle activity was registered using sEMG data (Zerowire, Cometa, Milan, IT) whereby circular Ag/Ag Cl electrodes (diameter 2 cm) were attached on the MG and tibialis anterior muscle bellies following the SENIAM guidelines (34). Inertial measurement units were placed on the shank and foot to analyze the ankle joint motions. Participants were assessed in supine position on an examination table, while the lower leg rested on a support to avoid contact of the MG muscle sEMG sensors with the table. The foot was fixated in a custom-made orthotic with an attached force sensor. Data of the latter were not reported in the current study. While holding on the force sensor with the shank in a fixed position, the examiner manually rotated the ankle from maximal plantarflexion to maximal dorsiflexion position. Four passive rotations over the full ankle ROM were performed at LV (i.e., in 5 s) and then four rotations were repeated at HV (i.e., performed as fast as possible). During each velocity trial, the participants were asked to remain relaxed, and a rest period of 7 s between the four repetitions was taken to minimize post-activation depression.

Offline data processing was performed using a custom-made Matlab software (Mathworks, R2017b). Data quality checks for velocity consistency and quality of the EMG signals were performed as described by Bar-On et al. (17). ROM and the maximum angular velocity (Vmax) were calculated and averaged over the repetitions per velocity trial. Hyperreflexia parameters were based on the registered MG muscle activity whereby first, the raw sEMG data were filtered, and second, root-mean square (RMS)-EMG (μV) was calculated, as previously described for both steps (17). Hyperreflexia was defined as the average RMS-EMG over the duration of the maximal velocity zone. This time zone was restricted to 200 ms prior to Vmax up to 90% of the full ROM. In addition, this hyperreflexia parameter was normalized to the peak RMS-EMG value of the obtained maximum voluntary isometric contractions (MVIC) and expressed as percentage. Both non-normalized and normalized hyperreflexia parameters were averaged over the stretch repetitions for the LV trial as well as for the change between the high and low velocity trials (change HV–LV). The latter is important to investigate the velocity-dependent nature of muscle hyperreflexia. In addition, visualization of normalized sEMG data across the LV and HV stretches over the ROM were used to classify responses according to previously described muscle activation patterns (i.e., velocity- or length-dependent pattern). More details on the definition of hyperreflexia parameters and muscle activation patterns can be found in earlier publications (17, 20).

Statistical Analysis

Statistical analyses were performed using SPSS (IBM SPPS Statistics version 27), and graphs were designed in GraphPad Prism 9 and Matlab software. Normality of the data was both visually checked and tested with the Shapiro–Wilk test. Due to the small sample size and non-normally distributed data, descriptive statistics included median and interquartile ranges and non-parametric test were applied. Independent-sample Kruskal–Wallis test was carried out to explore differences in mean ranks of age, anthropometric data, and muscle morphology between the three groups. Group differences were tested with post-hoc Mann–Whitney U (MWU) test. Bonferroni adjustment for multiple testing was applied after these post-hoc MWU analyses (p < 0.017). Comparisons of the clinical examination data, performance, and hyperreflexia parameters between the HSP and SCP groups were performed with MWU analyses, whereby the level of significance was set at p < 0.05.

Results

Participant characteristics and anthropometric data of the nine children in the HSP group were compared to (a) data of 15 participants in the SCP group and 15 participants of the TD group defined by 3DfUS, and (b) data of nine participants in the SCP group, defined by instrumented assessment of hyperreflexia (Table 1). No differences were found between the groups for age, body mass, body length, and leg length. Results of the clinical examination for ankle ROM, spasticity, and muscle strength did not differ significantly between the HSP group (n = 8, missing data for one child) and all children with SCP (n = 17) (Table 2). For the HSP group (n = 8), only one child had low MAS-scores (MAS 1+) and seven children had high MAS-scores (four MAS of 2; three MAS of 3), whereas in the SCP group (n = 17), 10 children had low MAS-scores (two MAS 1 and eight MAS 1+) and seven children had high MAS-scores (two MAS 2; five MAS 3). The ankle clonus scores differed significantly between the HSP and the total SCP group (p < 0.009). Specifically, only one child with HSP did not present a clonus, whereas the majority of children with SCP (i.e., 11 out of 17) had no clonus. All children in the HSP and SCP groups had regular physiotherapy with a median of 120 min per week and 100 min per week, respectively. One third of the HSP children frequently used day orthoses, whereas 10 out of the 17 children with SCP used orthoses during the day. For both groups, more than 50% of the children were never administered with BTX injections (Supplementary Table 2). The HSP and SCP cohort presented with median FMS scores of 6-6-5 and 6-5-5, respectively. More than half of the children in both groups were able to walk for more than 3,000 m, whereas the self-selected walking speed was slightly or even moderately reduced. Especially in the SCP group, 12 children out of 17 presented with a walking speed of more than 10 s over a 10-m distance. Pain during gait was reported by only two and three children in the HSP and SCP groups, respectively (Supplementary Table 3).

TABLE 1
www.frontiersin.org

Table 1. Participant characteristics and anthropometric data in the HSP, SCP, and TD group.

TABLE 2
www.frontiersin.org

Table 2. Median (minimum and maximum values) of all clinical examination results in the HSP and SCP group.

Compared with the TD reference group, both children with HSP and SCP differed significantly in the MG muscle volume normalized to the body mass (p ≤ 0.001) with the greatest deficits in the SCP group, namely, 32%. Post-hoc MWU test revealed only a significantly smaller absolute MV in the children with SCP compared with TD children (p = 0.004; Supplementary Figure 2). Muscle lengths were not significantly different between the three groups, although normalized muscle belly lengths (i.e., ML/MTUL and nML) tended to be lower in both patient groups compared with the TD group, with more alterations in the SCP group (Table 3).

TABLE 3
www.frontiersin.org

Table 3. Median (and IQR) of medial gastrocnemius muscle morphology outcomes in the HSP, SCP, and TD group.

Ankle joint rotations during HV stretch were performed over a median ROM of 50.7° (15.1) and 52.8° (18.4) with median maximal velocity of 81.1°/s (88.4) and 122.9°/s (45.6) in the HSP and SCP groups, respectively. The between-group differences in ankle joint rotations were not significant. Both non-normalized and normalized hyperreflexia parameters for LV and for the change between HV and LV were not significantly different between the HSP and SCP groups (Table 4). Higher hyperreflexia values, especially for change in HV–LV, tended to be present in the HSP group, compared with the SCP group. Eight children with HSP were categorized in a pure high-velocity-dependent pattern, whereas this pattern was observed in only six children with SCP. Additionally, muscle activity in one child of the HSP and three children of the SCP group was categorized as mixed high-velocity-dependent activation pattern. Figure 2 shows an example joint motion and muscle activation data during the instrumented assessment of hyperreflexia in a child with HSP and age-matched child with bilateral SCP.

TABLE 4
www.frontiersin.org

Table 4. Median (and IQR) of performance- and hyperreflexia-parameters in the HSP and SCP group.

FIGURE 2
www.frontiersin.org

Figure 2. Example illustration of (A) ankle joint angle vs. time graph, (B) ankle angular velocity vs. time graph, (C) root-mean square-electromyography (RMS-EMG) vs. time graph, and (D) 3D bar graph of average normalized RMS-EMG across three position zones and the two velocities of the medial gastrocnemius (MG) muscle, in a child with hereditary spastic paraplegia (HSP) and an age-matched child with spastic cerebral palsy (SCP) (20). Graphs (A–C) presented both a slow (dashed line) and high (continuous line) velocity stretch, whereby 0 s represent the time at the maximal velocity. In addition, the RMS-EMG vs. time graph (C) for the child with HSP indicates a clear ankle clonus, with 24 oscillations over a time period of 4 s. For the SCP case, this graph (C) indicates a catch during the high velocity stretch. In graph (D), the position zones indicate 10–90% of range of motion (ROM) divided in three equal parts. Velocity 1 and velocity 2 represent the low-velocity and high-velocity muscle stretches, respectively. High-velocity-dependent activation patterns are presented in both the HSP and SCP case.

Discussion

The purpose of this exploratory study was to investigate whether the MG exhibits muscle-related differences in two patient groups who are known to present with similar clinical impairments, despite their different etiology. To date, HSP is unexplored in terms of muscle morphology and instrumented assessment of hyperreflexia. In this study, only children presenting the pure or uncomplicated HSP type (i.e., solely pyramidal signs) with a pediatric onset were retrospectively included. This restriction minimized the possibility of a progressive disease course as seen in complicated and late onset HSP types (35). In the HSP group of this study, the mean age of symptom onset was 2.7 years, ranging from 2 months to 9 years old (Supplementary Table 1). The confirmed gene mutations in three of the four patients are described as the most common cause of autosomal dominant HSP. The mutations in these SPG3A and SPG4 cases and their HSP-associated proteins atlastin-1 and spastin, respectively, result in dysfunctional axonal transport with consequently corticospinal motor neuron degeneration (2, 3638). Moreover, children with SPG3A and SPG4 are known to present similarities in upper motor neuron signs that are commonly observed in children with bilateral SCP (39). In line with the current literature, the included child with the SPG56 subtype, which is described as recessive–autosomal HSP, is presenting as pure form with an early onset [i.e., at the age of 2 years; (40, 41)]. The involved gene mutation results in altered mitochondrial architecture and function (42). To further ensure differences in etiology between the two groups, only children with SCP showing clear, abnormal neuro-imaging findings (e.g., periventricular leukomalacia) were included (43). In the current study, the results of the clinical examination showed comparable levels of ankle mobility, spasticity, and muscle strength between both patient groups. However, the ankle clonus score was significantly higher in children with HSP compared with those with SCP, indicating increased presence and severity with mostly a quick clonus appearance with either slow or no stop (Supplementary Table 4).

Muscle Morphology

Altered MG muscle morphology in children with HSP was demonstrated by significantly smaller MV normalized to body mass of 30%, compared to TD. Furthermore, no significant differences were observed between the HSP and SCP groups, which indicate similarities in altered MG muscle morphology. Although the children in the SCP group tended to be lower in body mass and shorter in length than the children in the HSP group, normalization to this anthropometric data did not reveal significant differences in muscle volumes, nor in muscle lengths between both groups. However, a trend towards longer muscle belly lengths normalized to leg length was observed in the HSP group compared with the SCP group.

In agreement with previous research, muscle volumes in children with SCP were significantly smaller than the muscle volumes of the TD reference group (8). In the current study, average muscle volume deficits of 21% in the HSP group and 35% in the SCP group were observed when compared with the TD data. Interestingly, these differences are larger than the previously reported minimal detectable changes of 9.9% for MG muscle volumes (15). In the same sense, Noble et al. previously reported deficits of 38% in MG muscle volume for children with bilateral CP, with the same age ranges (i.e., older than 9 years) and similar levels of ambulatory ability [i.e., GMFCS levels I and II; (44)]. No significant differences were found with respect to the muscle lengths. Nonetheless, in comparison with data of the TD peers, absolute muscle belly lengths were 4 and 12% shorter in the HSP and SCP groups, respectively. These results were in line with previous research reporting 10% shorter muscle lengths in the paretic limb of children with unilateral SCP compared with data of TD children (45). However, only limited 3D ultrasound data of MG muscle lengths in children older than age 12 years have been reported so far (9). In conclusion, altered muscle volumes were observed in children with HSP and showed similarities to the muscle alterations in children with SCP.

Hyperreflexia

Hyperreflexia parameters were not significantly different between the HSP and SCP groups. Pathological muscle activity with increasing stretch velocity was observed in both patient groups, with a tendency of more hyperreflexia in the HSP group. This velocity-dependent hyperreflexia was confirmed by the results of the classification in hyperreflexia muscle activation patterns. Most of the children with HSP were classified with pure high-velocity-dependent muscle activation pattern. The amount of normalized EMG over increased velocity (i.e., EMG-change HV–LV) during MG muscle stretches was comparable to previous reported MG data. Specifically, the SCP group showed a median EMG-change of 3.96% (7.26) when comparing muscle activity at HV to LV stretches, which is in line with previously reported median scores of 2.57% (3.25) of children with the same spasticity scores (MAS 1.5) (17). The children in the HSP group presented MAS scores of 2 and normalized EMG-change of 7.88%. Previous research reported similar EMG change scores (median of 7.25%) in a cohort of children with SCP presenting MAS scores of 2.3 (17).

Clinical Implications

The objective instrumented investigation of MG muscle characteristics in the HSP group may help to gain more insight in the clinical impairments of children with HSP and to reveal potential markers for accurate differential diagnosis to SCP. In line with the findings of Harding et al., our results suggest the potential predominance of hyperreflexia rather than muscle weakness with the latter indirectly indicated based on the muscle volume results (1). Smaller muscle volumes are likely to generate less active force and are consequently indicative for decreased muscle strengths (46, 47). In the HSP group, muscle volumes tended to be less altered in comparison with the SCP group. This observation was combined with a tendency of higher hyperreflexia values, suggesting the predominance of spasticity in HSP. Yet, this trend needs to be confirmed in larger study samples.

Noticeably, it is important to investigate hyperreflexia using an objective instrumented assessment while performing quantitative analysis of the neurophysiological muscle response during passive muscle stretches in order to describe the neural contributor of hyper-resistance. Manual clinical examinations, such as MAS and MTA, are unable to discriminate among neural or non-neural contributors (19, 22, 48). The observed tendency of more hyperreflexia in the HSP group in comparison with rather similar MAS and MTA scores between the HSP and SCP groups confirms this statement.

The occurrence of clonus (i.e., repetitive oscillatory muscle contractions following stretch) is described as an upper motor neuron lesion sign and is generally presented in combination with hyperreflexia (49, 50). In an additional exploration, normalized RMS-EMG graphs of the high velocity stretches were visually inspected for the number of oscillations over the time period for which a minimal RMS-EMG activity was observed. Figure 2 represents an HSP case with a clear clonus, indicating 24 oscillations over a time period of 4 s. Both the results of the clinical examination scores for clonus and number of oscillations over time are presented per participant in Supplementary Table 4. The presence of clonus in particular in the HSP group is equally observed according to these two approaches. Moreover, in the children with HSP showing clinically a quick clonus with no stop, an increasing number of oscillations was suggested in comparison with children with a quick stop. In contrast, in the SCP group, only three children showed a clonus in the clinical examination, whereas the phenomenon was confirmed by EMG data for only two children with SCP. However, more research should be conducted to provide valid and reliable methods to quantify clonus with the described assessment of hyperreflexia in both children with HSP and SCP. In the context of the current study, hyperreflexia and clonus can be suggested as clinical markers to distinguish the two disorders.

Treatment approaches in children with HSP are symptomatic and generally aimed to reduce muscle spasticity. Treatment management mainly involves physiotherapy and tone-reducing medication. However, a recent review highlighted the lack of high-level studies, like randomized control trials, that provide sufficient evidence to promote these treatments in children with HSP (6). Yet, the hyperreflexia data of the current study might indicate tone reduction as an important treatment goal in children with HSP. Despite limited data on the effectiveness of BTX injections on focal spasticity and function in children with HSP (6, 7), the presence of high-velocity-dependent phenotypes indicates the potential positive outcomes of BTX treatment. Indeed, in children with SCP, previous research on the effect of BTX demonstrated more pronounced decreases in muscle tone in children with high-velocity-dependent muscle activation patterns compared with length-dependent muscle activation patterns (21). Apart from daily orthotic management and regular physiotherapy as the general treatment approach, the use of BTX is widely applied as tone-reducing treatment in children with bilateral SCP (51). Since the muscle characteristics were found to be similar, it is reasonable to recommend analogous treatment approaches in children with HSP and SCP.

Limitations and Future Perspectives

It should be noted that some limitations are present in the current study. The main drawback of this retrospective study is the small sample size, which compromises the power of the statistical tests and may prevent generalization of the observations. Due to the nature of the single-center, retrospective study design, we were limited to the available datasets for each subgroup. In addition, specific inclusion criteria for both the HSP and SCP groups were applied, aiming to reduce the clinical heterogeneous presentation. Despite these efforts, the muscle data presented high variability, although normalization to anthropometric data successfully decreased this heterogeneity. This finding suggests that children of the different cohorts should be matched for anthropometric data rather than for age. Notably, HSP is a rare disease with an estimated prevalence of 1:10,000, and not all case were genetically confirmed (52). The restriction to only include pure HSP types enlarges the rareness of the disease. Future studies should aim for a prospective case-control design and could be carried out as a multi-center study to increase the size of the dataset.

This study only addressed the MG muscle. The restriction to the MG as muscle of interest was fixed due to (a) the available muscle-specific datasets and (b) the ability to compare our results with previous publications. Hamstring muscle spasticity and weakness have been described as clinical impairments in HSP (4, 35). Since muscle activation patterns were found to be muscle specific, future studies should investigate these muscle-specific characteristics. In the current study, the observed high-velocity-dependent pattern of the MG muscle confirmed this previously reported MG-specific pattern (20).

Furthermore, the current study focused on the neurophysiological response to passive muscle stretch. Future studies should further explore the specific contribution of both neural and non-neural components of the hyper-resistance against passive muscle stretch using instrumented assessments in children with HSP. In addition, these assessments should be used to consider the effectiveness of frequently applied tone-reducing treatment approaches in children with HSP.

Conclusion

This is the first study that explored the MG muscle characteristics in children with HSP using instrumented impairment assessments. Muscle volume deficits, high-velocity hyperreflexia, and ankle clonus were presented in children with HSP. In general, instrumented assessments of the altered muscle morphology and hyperreflexia suggested similar MG muscle-related impairments in pediatric pure HSP types with pediatric onset and in children with bilateral SCP. Hence, the use of analogous treatment approaches might be supported by the study results. Due to the limited statistical power caused by the limitations in sample size, the study results should be confirmed in larger cohorts.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethical Committee of the University Hospitals Leuven and KU Leuven (S56977 and S59945). Written informed consent was provided by the participants' parents/legal guardians.

Author Contributions

This study was designed by ND, LB-O, EO, KD, and AV. ND, BH, and NP contributed to the data collection of the retrospective database. EO and AV evaluated the eligibility of the participants. ND was responsible for the data processing and conducted all presented analyses. ND, LB-O, KD, and AV contributed to the interpretation of the results and were involved in the critical revision and editing of the manuscript that was written by ND. All authors approved the final version of the manuscript and agreed to be accountable for the content of the work. All authors have had complete access to the study data throughout the study.

Funding

ND was supported by Internal Funds of KU Leuven (C24/18/103). This work was further supported by the TBM grant (TAMTA-T005416N) and post-doctoral grant (12R4215N) from the Research Foundation Flanders (FWO), Belgium, from La Fondation Motrice (contract 2016/8) and from the Netherlands Organization for Scientific Research (NWO-016.186.144).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank Francesco Cenni and Simon-Henri Schless for their contribution to this work and the development of the 3DfUS method. We would like to thank Liesbeth Keldermans and Valerie Race for providing all information related to the genetic work-up. This research team was very grateful to all children and parents who participated in this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.635032/full#supplementary-material

References

1. Harding AE. Classification of the hereditary ataxias and paraplegias. Lancet. (1983) 321:1151–5. doi: 10.1016/S0140-6736(83)92879-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Salinas S, Proukakis C, Crosby A, Warner TT. Hereditary spastic paraplegia: clinical features and pathogenetic mechanisms. Lancet Neurol. (2008) 7:1127–38. doi: 10.1016/S1474-4422(08)70258-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Klebe S, Stolze H, Kopper F, Lorenz D, Wenzelburger R, Volkmann J, et al. Gait analysis of sporadic and hereditary spastic paraplegia. J Neurol. (2004) 251:571–8. doi: 10.1007/s00415-004-0366-7

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fink JK. Hereditary spastic paraplegia: clinical principles and genetic advances. Semin Neurol. (2014) 34:293–305. doi: 10.1055/s-0034-1386767

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Graham HK, Rosenbaum P, Paneth N, Dan B, Lin J. Cerebral palsy. Nat Rev. (2016) 2:15082. doi: 10.1038/nrdp.2015.82

CrossRef Full Text | Google Scholar

6. Bellofatto M, De Michele G, Iovino A, Filla A, Santorelli FM. Management of hereditary spastic paraplegia: a systematic review of the literature. Front Neurol. (2019) 10:3. doi: 10.3389/fneur.2019.00003

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rousseaux M, Launay MJ, Kozlowski O, Daveluy W. Botulinum toxin injection in patients with hereditary spastic paraparesis. Eur J Neurol. (2007) 14:206–12. doi: 10.1111/j.1468-1331.2006.01617.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Barrett RS, Lichtwark GA. Gross muscle morphology and structure in spastic cerebral palsy: a systematic review. Dev Med Child Neurol. (2010) 52:794–804. doi: 10.1111/j.1469-8749.2010.03686.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Williams I, Reid L, Stott NS, 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:263–73. doi: 10.1111/dmcn.14714

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pitcher CA, Elliott CM, Williams SA, Licari MK, Kuenzel A, Shipman PJ, et al. Childhood muscle morphology and strength: alterations over six months of growth. Muscle Nerve. (2012) 46:360–6. doi: 10.1002/mus.23326

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Reid SL, Pitcher CA, Williams SA, Licari MK, Valentine JP, Shipman PJ, et al. Does muscle size matter? The relationship between muscle size and strength in children with cerebral palsy. Disabil Rehabil. (2015) 37:579–84. doi: 10.3109/09638288.2014.935492

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Herskind A, Ritterband-Rosenbaum A, Willerslev-Olsen M, Lorentzen J, Hanson L, Lichtwark G, et al. Muscle growth is reduced in 15-month-old children with cerebral palsy. Dev Med Child Neurol. (2016) 58:485–91. doi: 10.1111/dmcn.12950

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Willerslev-Olsen M, Choe Lund M, Lorentzen J, Barber L, Kofoed-Hansen M, Nielsen JB. Impaired muscle growth precedes development of increased stiffness of the triceps surae musculotendinous unit in children with cerebral palsy. Dev Med Child Neurol. (2018) 60:672–9. doi: 10.1111/dmcn.13729

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Cenni F, Monari D, Desloovere K, Aertbeliën E, Schless S-H, Bruyninckx H. The reliability and validity of a clinical 3D freehand ultrasound system. Comput Methods Programs Biomed. (2016) 136:179–87. doi: 10.1016/J.CMPB.2016.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Cenni F, Schless SH, Bar-On L, Aertbeliën E, Bruyninckx H, Hanssen B, et al. Reliability of a clinical 3D freehand ultrasound technique: analyses on healthy and pathological muscles. Comput Methods Programs Biomed. (2018) 156:97–103. doi: 10.1016/j.cmpb.2017.12.023

CrossRef Full Text | Google Scholar

16. 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:384–91. doi: 10.1111/joa.12927

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bar-On L, Aertbeliën E, Wambacq H, Severijns D, Lambrecht K, Dan B, et al. A clinical measurement to quantify spasticity in children with cerebral palsy by integration of multidimensional signals. Gait Posture. (2013) 38:141–7. doi: 10.1016/j.gaitpost.2012.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Schless S-H, 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. doi: 10.1371/journal.pone.0131011

PubMed Abstract | CrossRef Full Text | Google Scholar

19. van den Noort JC, Bar-On L, Aertbeliën E, Bonikowski M, Braendvik SM, Broström EW, et al. European consensus on the concepts and measurement of the pathophysiological neuromuscular responses to passive muscle stretch. Eur J Neurol. (2017) 24:981–e38. doi: 10.1111/ene.13322

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bar-On L, Aertbeliën E, Molenaers G, Desloovere K. Muscle activation patterns when passively stretching spastic lower limb muscles of children with cerebral palsy. PLoS ONE. (2014) 9:e91759. doi: 10.1371/journal.pone.0091759

PubMed Abstract | CrossRef Full Text | Google Scholar

21. 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. doi: 10.3389/fneur.2020.00378

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bar-On L, Harlaar J, Desloovere K. Spasticity assessment in cerebral palsy. In: Miller F, Bachrach S, Lennon N, O'Neil ME, editors. Cerebral Palsy. Cham: Springer (2020). p. 585–600. doi: 10.1007/978-3-319-74558-9_40

CrossRef Full Text | Google Scholar

23. Schüle R, Wiethoff S, Martus P, Karle KN, Otto S, Klebe S, et al. Hereditary spastic paraplegia: clinicogenetic lessons from 608 patients. Ann Neurol. (2016) 79:646–58. doi: 10.1002/ana.24611

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Shribman S, Reid E, Crosby AH, Houlden H, Warner TT. Hereditary spastic paraplegia: from diagnosis to emerging therapeutic approaches. Lancet Neurol. (2019) 18:1136–46. doi: 10.1016/S1474-4422(19)30235-2

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wolf SI, Braatz F, Metaxiotis D, Armbrust P, Dreher T, Dö Derlein L, et al. Gait analysis may help to distinguish hereditary spastic paraplegia from cerebral palsy. Gait Posture. (2011)33:556–61. doi: 10.1016/j.gaitpost.2011.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Bickley C, Mitchell K, Scott A, Bury M, Oyelami M. Familiarity with Hereditary Spastic Paraplegia (HSP) and differentiation of upper body gait characteristics between children with HSP and spastic diplegic cerebral palsy. Phys Occup Ther Pediatr. (2020) 41:99–113. doi: 10.1080/01942638.2020.1819934

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Fink JK. Advances in the hereditary spastic paraplegias. Exp Neurol. (2003) 184:106–10. doi: 10.1016/j.expneurol.2003.08.005

CrossRef Full Text | Google Scholar

28. McDermott CJ, White K, Bushby K, Shaw PJ. Hereditary spastic paraparesis: a review of new developments. J Neurol Neurosurg Psychiatry. (2000) 69:150–60. doi: 10.1136/jnnp.69.2.150

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Bohannon RW, Smith MB. Interrater reliability of a modified Ashworth scale of muscle spasticity. Phys Ther. (1987) 67:206–7. doi: 10.1093/ptj/67.2.206

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Boyd RN, Graham HK. Objective measurement of clinical findings in the use of botulinum toxin type A for the management of children with cerebral palsy. Eur J Neurol. (1999) 6:23–s35. doi: 10.1111/j.1468-1331.1999.tb00031.x

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Matthews WB. Aids to the examination of the peripheral nervous system. J Neurol Sci. (1977) 33:299. doi: 10.1016/0022-510x(77)90205-2

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Graham HK, Harvey A, Rodda J, Nattrass GR, Pirpiris M. The functional mobility scale. J Pediatr Orthop. (2004) 24:514–520. doi: 10.1097/00004694-200409000-00011

CrossRef Full Text | Google Scholar

33. Schless SH, Cenni F, Bar-On L, Hanssen B, Kalkman B, O'brien T, et al. Medial gastrocnemius volume and echo-intensity after botulinum neurotoxin A interventions in children with spastic cerebral palsy. Dev Med Child Neurol. (2018) 61:783–90. doi: 10.1111/dmcn.14056

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Hermens HJ, Freriks B, Disselhorst-Klug C, Rau G. Development of recommendations for SEMG sensors and sensor placement procedures. J Electromyogr Kinesiol. (2000) 10:361–74. doi: 10.1016/s1050-6411(00)00027-4

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fink J. Hereditary spastic paraplegia: clinico-pathologic features and emerging molecular mechanisms. Acta Neuropathol. (2013) 126:307–28. doi: 10.1007/s00401-013-1115-8.Hereditary

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Botzolakis EJ, Zhao J, Gurba KN, Macdonald RL, Hedera P. The effect of HSP-causing mutations in SPG3A and NIPA1 on the assembly, trafficking, and interaction between atlastin-1 and NIPA1. Mol Cell Neurosci. (2011) 46:122–35. doi: 10.1016/j.mcn.2010.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Hazan J, Fonknechten N, Mavel D, Paternotte C, Samson D, Artiguenave F, et al. Spastin, a new AAA protein, is altered in the most frequent form of autosomal dominant spastic paraplegia. Nat Genet. (1999) 23:296–303. doi: 10.1038/15472

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Park SH, Zhu PP, Parker RL, Blackstone C. Hereditary spastic paraplegia proteins REEP1, spastin, and atlastin-1 coordinate microtubule interactions with the tubular ER network. J Clin Invest. (2010) 120:1097–110. doi: 10.1172/JCI40979

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Rainier S, Sher C, Reish O, Thomas D, Fink JK. De novo occurrence of novel SPG3A/atlastin mutation presenting as cerebral palsy. Arch Neurol. (2006) 63:445–7. doi: 10.1001/archneur.63.3.445

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Klebe S, Stevanin G, Depienne C. Clinical and genetic heterogeneity in hereditary spastic paraplegias: from SPG1 to SPG72 and still counting. Rev Neurol. (2015) 171:505–30. doi: 10.1016/j.neurol.2015.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Durand CM, Dhers L, Tesson C, Tessa A, Fouillen L, Jacqueré S, et al. CYP2U1 activity is altered by missense mutations in hereditary spastic paraplegia 56. Hum Mutat. (2018) 39:140–51. doi: 10.1002/humu.23359

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Tesson C, Nawara M, Salih MAM, Rossignol R, Zaki MS, Al Balwi M, et al. Alteration of fatty-acid-metabolizing enzymes affects mitochondrial form and function in hereditary spastic paraplegia. Am J Hum Genet. (2012) 91:1051–64. doi: 10.1016/j.ajhg.2012.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Krägeloh-Mann I, Horber V. The role of magnetic resonance imaging in furthering understanding of the pathogenesis of cerebral palsy. Dev Med Child Neurol. (2007) 49:948. doi: 10.1111/j.1469-8749.2007.00948.x

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Noble JJ, Fry NR, Lewis AP, Keevil SF, Gough M, Shortland AP. Lower limb muscle volumes in bilateral spastic cerebral palsy. Brain Dev. (2014) 36:294–300. doi: 10.1016/j.braindev.2013.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Malaiya R, McNee AE, Fry NR, Eve LC, Gough M, Shortland AP. The morphology of the medial gastrocnemius in typically developing children and children with spastic hemiplegic cerebral palsy. J Electromyogr Kinesiol. (2007) 17:657–63. doi: 10.1016/J.JELEKIN.2007.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fukunaga T, Miyatani M, Tachi M, Kouzaki M, Kawakami Y, Kanehisa H. Muscle volume is a major determinant of joint torque in humans. Acta Physiol Scand. (2001) 172:249–55. doi: 10.1046/j.1365-201X.2001.00867.x

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Mockford M, Caulton JM. The pathophysiological basis of weakness in children with cerebral palsy. Pediatr Phys Ther. (2010) 22:222–33. doi: 10.1097/PEP.0b013e3181dbaf96

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Bar-On L, Van Campenhout A, Desloovere K, Aertbeliën E, Huenaerts C, Vandendoorent B, et al. Is an instrumented spasticity assessment an improvement over clinical spasticity scales in assessing and predicting the response to integrated botulinum toxin type A treatment in children with cerebral palsy? Arch Phys Med Rehabil. (2014) 95:515–23. doi: 10.1016/J.APMR.2013.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Hidler JM, Rymer WZ. A simulation study of reflex instability in spasticity: origins of clonus. IEEE Trans Rehabil Eng. (1999) 7:327–40. doi: 10.1109/86.788469

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Boyraz I, Uysal H, Sarman H, Ethica P, Hospital I. Clonus: definition, mechanism, treatment. Med Glas. (2015) 12:19–26.

Google Scholar

51. Multani I, Manji J, Ison TH, Khot A, Graham K. Botulinum toxin in the management of children with cerebral palsy. Pediatr Drugs. (2019) 21:261–81. doi: 10.1007/s40272-019-00344-8

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Ruano L, Melo C, Silva MC, Coutinho P. The global epidemiology of hereditary ataxia and spastic paraplegia: a systematic review of prevalence studies. Neuroepidemiology. (2014) 42:174–83. doi: 10.1159/000358801

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hereditary spastic paraplegia, cerebral palsy, instrumented impairment assessments, ultrasound, muscle morphology, muscle volume, hyperreflexia, spasticity

Citation: De Beukelaer N, Bar-On L, Hanssen B, Peeters N, Prinsen S, Ortibus E, Desloovere K and Van Campenhout A (2021) Muscle Characteristics in Pediatric Hereditary Spastic Paraplegia vs. Bilateral Spastic Cerebral Palsy: An Exploratory Study. Front. Neurol. 12:635032. doi: 10.3389/fneur.2021.635032

Received: 29 November 2020; Accepted: 26 January 2021;
Published: 26 February 2021.

Edited by:

Pasquale Striano, University of Genoa, Italy

Reviewed by:

Filippo M. Santorelli, Fondazione Stella Maris (IRCCS), Italy
Alessandra Tessa, Fondazione Stella Maris (IRCCS), Italy
Sergio Lerma, Autonomous University of Madrid, Spain

Copyright © 2021 De Beukelaer, Bar-On, Hanssen, Peeters, Prinsen, Ortibus, Desloovere and Van Campenhout. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Nathalie De Beukelaer, nathalie.debeukelaer@kuleuven.be

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.