- 1Department of Physiotherapy, HESAV School of Health Sciences, HES-SO University of Applied Sciences and Arts Western Switzerland, Lausanne, Switzerland
- 2Swiss BioMotion Lab, Department of Musculoskeletal Medicine, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
- 3Institute of Microengineering, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
Despite a large body of evidence demonstrating spinal movement alterations in individuals with chronic low back pain (CLBP), there is still a lack of understanding of the role of spinal movement behavior on LBP symptoms development or recovery. One reason for this may be that spinal movement has been studied during various functional tasks without knowing if the tasks are interchangeable, limiting data consolidation steps. The first objective of this cross-sectional study was to analyze the influence of the functional tasks on the information carried by spinal movement measures. To this end, we first analyzed the relationships in spinal movement between various functional tasks in patients with CLBP using Pearson correlations. Second, we compared the performance of spinal movement measures to differentiate patients with CLBP from asymptomatic controls among tasks. The second objective of the study was to develop task-independent measures of spinal movement and determine the construct validity of the approach. Five functional tasks primarily involving sagittal-plane movement were recorded for 52 patients with CLBP and 20 asymptomatic controls. Twelve measures were used to describe the sagittal-plane angular amplitude and velocity at the lower and upper lumbar spine as well as the activity of the erector spinae. Correlations between tasks were statistically significant in 91 out of 99 cases (0.31 ≤ r ≤ 0.96, all p < 0.05). The area under the curve (AUC) to differentiate groups did not differ substantially between tasks in most of the comparisons (82% had a difference in AUC of ≤0.1). The task-independent measures of spinal movement demonstrated equivalent or higher performance to differentiate groups than functional tasks alone. In conclusion, these findings support the existence of an individual spinal movement signature in patients with CLBP, and a limited influence of the tasks on the information carried by the movement measures, at least for the twelve common sagittal-plane measures analysed in this study. Therefore, this work brought critical insight for the interpretation of data in literature reporting differing tasks and for the design of future studies. The results also supported the construct validity of task-independent measures of spinal movement and encouraged its consideration in the future.
Introduction
Alterations in spinal movement have been suggested as one of the key physical factors in the persistence of chronic low back pain (CLBP) (Marras et al., 1995; O’Sullivan, 2005; Dubois et al., 2014), however, the understanding of spinal movement behavior in CLBP remains limited. The abundance of measures used to describe spinal movement, as outlined by two recent systematic reviews (Papi et al., 2018; Moissenet et al., 2021), is undoubtedly one of the reasons limiting a better understanding. Consequently, to advance the field, there is a need to determine if that many measures are needed or if it would be possible to focus on a selection of measures. While characterizing the movement in terms of angular amplitude, angular velocity and muscle activity appears appropriate for a comprehensive description (Laird et al., 2014; Papi et al., 2018; Moissenet et al., 2021), the necessity to test multiple functional tasks remains to be determined.
So far, CLBP spinal movement has been assessed during a range of different functional tasks (Shum et al., 2007; Christe et al., 2016b, 2020; Lima et al., 2018; Matheve et al., 2019). However, these prior works mainly assessed one task at a time, or when multiple tasks were assessed, the influence of the task on spinal movement was not analyzed (Papi et al., 2018; Moissenet et al., 2021). Yet, individual consistency in spinal angular amplitudes across different functional tasks has been demonstrated in pain-free participants (Alqhtani et al., 2015; Seerden et al., 2019), suggesting that it may also be the case for individuals with CLBP and for angular velocity and muscle activity measures. This possibility is particularly supported in the clinical setting, where consistent movement patterns are frequently observed in CLBP patients. If measures from different tasks were to carry similar information (for example, patient X moves with relatively lower flexion than the other individuals independently of the tasks analyzed), this would suggest the existence of “individual spinal movement signatures.” Clarifying this point appears essential, as it would help the interpretation of data in literature and the design of new experiments. In fact, on one side, knowing the extend of movement signatures would provide a basis for the comparison of studies testing different tasks and, on the other side, it would provide a rationale for the number and specificity of the tasks to include in future studies.
If individual spinal movement signatures were to exist in CLBP patients, this would question the possibility to develop task–independent measures of spinal movement. Such task–independent measures could produce more robust assessments and reduce the number of variables to deal with in statistical analyses, which would be beneficial for both the design and the interpretation of future studies. If possible, this simplified description of spinal movement could prove particularly useful to detangle the role of spinal movement alterations in CLBP development or recovery and inform rehabilitation principles (Wernli et al., 2020b; Schmid et al., 2021).
Therefore, the first objective of the study was to determine the influence of the functional tasks on the information carried by spinal movement measures. To this end, this study aimed: 1) to analyze the correlations among various functional tasks in patients with CLBP; 2) to compare the performance to differentiate patients with CLBP from asymptomatic controls among functional tasks. Based on prior research, measures corresponding to peaks and ranges of sagittal-plane lumbar angular amplitude and angular velocity, as well as maximal erector spinae activity were analyzed (Shum et al., 2005; Dankaerts et al., 2009; Hidalgo et al., 2012; Christe et al., 2016b, 2020). Following previous work in asymptomatic people and a pilot study in patients with CLBP (Alqhtani et al., 2015; Christe et al., 2016a; Seerden et al., 2019), it was hypothesized that the spinal movement measures would be positively correlated among the functional tasks. It was also hypothesized that there would be no relevant performance difference among functional tasks (Dankaerts et al., 2009; Laird et al., 2014; Christe et al., 2016b, 2020; Moissenet et al., 2021);
The second objective of the study was to determine the construct validity of task–independent measures of spinal movement obtained by grouping measures across multiple tasks. Specifically, we aimed: 3) to compare the performance of task–independent measures to differentiate patients with CLBP from asymptomatic controls to the performance of task-specific measures. It was hypothesized that the performance of task–independent measures would not be inferior compared to the performance of task-specific measures.
Materials and Methods
Design
This cross-sectional case-controlled study is reported according to the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) criteria (von Elm et al., 2007).
Participants and Setting
Recruitment took place in an interdisciplinary rehabilitation program (IRP) at the local university hospital. This IRP is a full time 3-week program that includes patients with difficulties to maintain their leisure and professional activity because of CLBP. Participants to the IRP were invited to take part in the study if they met the inclusion/exclusion criteria. Both males and females could participate if they had a diagnosis of non-specific LBP with or without leg pain for more than 3 months, a sufficient French level and an age from 18 to 65 years old. Exclusion criteria for the CLBP group in this study were the presence of a diagnosis of specific LBP and/or previous back surgery that limited spinal mobility (i.e., spinal fusion). Asymptomatic controls were a convenience sample recruited via emails and flyers. To be included, they had to have no history of LBP requiring third-party attention during the last 2 years. They were also excluded in the presence of any recent or current episode of LBP. Exclusion criteria for both groups included pregnancy, a body mass index (BMI) above 32 kg/m2 and other concomitant pain or condition that could compromise the evaluation of lumbar kinematics. The BMI cutoff was selected to limit the influence of body shape on lumbar kinematics and experimental complications, without compromising external validity and patients’ recruitment. The research was approved by the local Research Ethics Committee (CER-VD 2018-00188) and all participants signed an informed consent form before enrolment in the study.
Experimental Procedures
Participants were invited to the movement analysis laboratory at the university hospital for a measurement session before the IRP. First, participants completed three reliable and valid questionnaires to document mean pain intensity during the last 24 h, kinesiophobia and catastrophizing using the numeric pain rating scale (24h-NPRS), the Tampa Scale of Kinesiophobia (TSK) and the Pain Catastrophizing Scale (PCS), respectively (Sullivan, 1995; Vlaeyen et al., 1995; Chapman et al., 2011). Then, after having cleaned the skin with alcohol and shaved it, if necessary, two pairs of electrodes were placed bilaterally parallel to the erector spinae fibers 3 cm lateral to the L3 spinous process (Dupeyron et al., 2013; Wong et al., 2016). Participants then performed one submaximal voluntary contraction in crook lying as described by Dankaerts et al. (2004). Reflective markers were then attached to the participants lumbar region and pelvis following a previously described protocol (Seay et al., 2008; Christe et al., 2016b, 2017, 2020). Lumbar markers were placed on the spinous processes of L1, L3 and L5 with four additional markers attached between these markers on each side of the spine, at a distance of 5 cm. Pelvis markers were placed on the posterior superior iliac spines, anterior superior iliac spines and iliac crest tips. Marker trajectories and lumbar muscle activity were measured using an optoelectronic motion capture system with 14 cameras (Vicon, Oxford Metrics, Oxford, United Kingdom) and an electromyography system (Myon, Schwarzenberg, CH) recording synchronously at 120 and 1200 Hz, respectively.
Data collection started with the recording of a reference standing posture, where participants were standing upright and looking forward with arms at 60° of shoulder abduction. Then, five functional tasks were recorded in the same order for every participant to avoid varying remnant effects among participants as some were more exacerbating for participants: standing flexion, sit-to-stand, stepping-up on a 36 cm high step, picking-up a sponge from the floor and lifting a 4.5 kg box from the floor. All functional tasks were first demonstrated in a video with standardized instructions (Supplementary Appendix SA). Each functional activity was practiced between one and three times and then recorded three times, except for picking-up that was recorded ten times (for the purpose of another study, Christe et al., 2019). Following the video instruction, and before performing each task, participants rated on a zero to ten scale how much do they think the task to-be-performed is harmful for the back (perceived harm; 0: not harmful at all; 10: extremely harmful). After each task, they also rated their pain during the task with a numeric pain rating scale. At the end of the session, participants completed the French version of the Oswestry Disability Index (ODI) (Fairbank and Pynsent, 2000; Vogler et al., 2008).
Data Processing
Spinal kinematics were calculated based on a three-segment biomechanical model that includes the pelvis and the lower lumbar and upper lumbar spine (Christe et al., 2016b, 2017, 2020). Briefly, markers’ trajectories were used to calculate the orientation of anatomical frames embedded in each segment. The joint coordinate system (Grood and Suntay, 1983) was then used to calculate sagittal-plane joint angles at the lower lumbar (LLSa) and the upper lumbar (ULSa) joints. LLSa was defined as the angle between the lower lumbar segment (L3-L5 central and lateral markers) and the pelvis segment, while ULSa was defined as the angle between the upper lumbar segment (L1-L3 central and lateral markers) and the lower lumbar segment. Angles were low-pass filtered using a 15 Hz Butterworth filter. The amplitude of the angles during the reference standing posture were subtracted from the angle curves to limit the inter-individual variations in morphology. Angular velocity curves (LLSv and ULSv) were obtained by numerical differentiation of the angle curves.
Electromyography recordings were band-pass filtered using a Butterworth filter with cutoff frequencies at 20 and 450 Hz. Then, for both muscles, the minimal amplitude of the electromyography signals recorded during the entire session was identified and this minimal amplitude was subtracted from the signals. This operation defined a zero-value (0%) for the electromyography data. Next, for both muscles, the signals were scaled in order to have the amplitude recorded during the submaximal voluntary contraction in crook lying equal 100%. Submaximal contraction was chosen for the normalization because its reliability was shown to be superior to maximal contraction in CLBP patients (Dankaerts et al., 2004).
In order to extract the movement measures from the angular amplitude, angular velocity and muscle activity curves, first, the curves were time-normalized to 0–100% for each repetition of each task. The beginning and the end of the task were determined visually using strict criteria based on markers displacements (Christe et al., 2016b, 2020). Then, following the methodology in prior studies (Christe et al., 2016b, 2017, 2020), the curves were tested using the coefficient of multiple correlation (CMC) (Kadaba et al., 1989) and the curves presenting a characteristic pattern were described by discrete measures. In total, 12 measures were identified. They included the peak flexion angle and sagittal-plane range of motion (ROM) at the LLS (LLSaflexion and LLSarange) and at the ULS (ULSaflexion and ULSarange); the peak angular velocity in flexion, the peak angular velocity in extension and the range between velocity peaks at the LLS (LLSvflexion, LLSvextension and LLSvrange) and at the ULS (ULSvflexion, ULSvextension and ULSvrange) and; the peaks of erector spinae muscle activity during the first (EMGpeak1) and second (EMGpeak2) halves of the tasks. Not all tasks presented the characteristic features necessary to the extraction of the 12 measures. Indeed, ULSvflexion, ULSvextension, ULSvrange and EMGpeak2 were only present in flexion, picking-up and lifting. The measures were averaged over the repetitions in order to have only one value per participant and task. Finally, for EMGpeak1 and EMGpeak2, the maximal value observed between the left and right erector spinae muscles was kept for analysis.
Task–independent measures were calculated by averaging the measures obtained with the diverse tasks. The averaging was done independently for each participant and measure. To give similar weight to all the tasks, a Z-score transformation was applied to the measures before averaging over the tasks. The transformations were based on the means and standard deviations of the asymptomatic controls. Consequently, as illustrated in Figure 1, the task-independent measures were dimensionless: their values indicated how they were situated compared to the reference (asymptomatic) population. All calculations were performed with Matlab (R2019b, MathWorks, Inc, Natick, MA).
FIGURE 1. Illustration of the task-independent measures concept. Task-independent measures are expressed according to the mean and SD of the reference population. These dimensionless measures therefore indicate how they situate compared to the reference population.
Statistical Analysis
The normality of the measures was assessed visually using QQ plots and tested with the Shapiro-Wilk test (Ghasemi and Zahediasl, 2012). Extreme outliers were discarded from the analyses using a standard procedure (Portney and Watkins, 2000).
For the first aim, relationships among functional tasks in patients with CLBP were tested with Pearson correlations. These analyses were performed separately for each of the 12 measures. Correlation coefficients (r) were interpreted as small (0.1 ≤ r < 0.3), medium (0.3 ≤ r < 0.5) and large (r ≥ 0.5) (Cohen, 1988).
For the second aim, we conducted binary logistic regression models with the movement measures as the independent variable and the group as the dependent variable. Regression models were performed separately for each measure of each task. The performance to differentiate patients with CLBP from asymptomatic controls was primarily tested with the area under the curve (AUC) value. These values were categorized as poor (AUC<0.7), acceptable (0.7 ≤ AUC <0.8), excellent (0.8 ≤ AUC <0.9) or outstanding (AUC≥ 0.9) discriminations (Hosmer et al., 2013). Following these categories, a difference in AUC between tasks above 0.1 was considered indicative of a relevant difference in groups’ differentiation performance. For completeness, other usual statistics of logistic regression models were calculated: coefficient of determination (r2), sensitivity (Sn), specificity (Sp), positive likelihood ratios (LR+) and negative likelihood ratios (LR−). In addition, independent t-tests were conducted to determine if group differences were statistically significant and Cohen’s d effect sizes (ES) were computed to quantify the size of the differences between groups (Cohen, 1988). ES were interpreted as very small (0.01 ≤ ES < 0.2), small (0.2 ≤ ES < 0.5), moderate (0.5 ≤ ES < 0.8), large (0.8 ≤ ES < 1.2), very large (1.2 ≤ ES < 2.0) and huge (ES ≥ 2.0) (Cohen, 1988; Sawilowsky, 2009).
For the third aim regarding the construct validity of task–independent measures of spinal movement, binary logistic regression models were conducted for each of the task–independent measure. AUC values were interpreted as detailed above for the tasks’ comparison. Independent t-tests and ES were also calculated for task–independent measures. For completeness with the first aim, Pearson correlations were performed between task-independent and task-specific data for each measure. Statistical analyses were performed with SPSS (Version 25, IBM, NY, United States), using a significance level at α < 0.05.
Sample Size
To detect a correlation coefficient between functional tasks of r ≥ 0.4, as reported in prior asymptomatic and pilot studies (Alqhtani et al., 2015; Christe et al., 2016a; Bujang and Baharum, 2016), with a power of 0.8 and α error of 0.05, the minimum sample size in the patients’ group was 46. For the logistic regression, the usual recommendations were followed, indicating a minimum of 15 participants per group in the case of models with a single independent variable (Portney and Watkins, 2000). Five asymptomatic participants and six participants with CLBP were added to prevent insufficient power due to potential drop-out or corrupted movement data.
Results
Fifty-two patients with CLBP (sex: 63.5% male; age (mean ± SD): 40.0 ± 10.4 years old; BMI: 25.3 ± 3.3 kg/m2) and 20 asymptomatic controls (55% male; 38.2 ± 10.9 years old; 22.7 ± 2.8 kg/m2) were included in the study. The mean 24 h-NPRS, TSK, PCS and ODI scores of the patients were 5.6 ± 2.1, 44.3 ± 7.5, 25.2 ± 11.7 and 35.3 ± 11.2, respectively. Mean pain during movement was 4.6 ± 2.5 during flexion, 4.8 ± 2.7 during lifting, 4.3 ± 2.6 during picking-up, 2.7 ± 2.2 during stepping-up and 2.7 ± 2.5 during sit-to-stand for the patients. Mean perceived harm by the patients for each movement was 4.5 ± 3.5 for flexion, 6.3 ± 3.1 for lifting, 5.4 ± 3.2 for picking-up, 1.9 ± 2.4 for stepping-up and 2.6 ± 3.0 for sit-to-stand. Mean movement measures are reported in Supplementary Appendix SB for both groups. Movement data were available for at least 48 CLBP patients and 17 asymptomatic controls.
Correlations between tasks were statistically significant in 91 out of 99 (92%) cases (0.31 ≤ r ≤ 0.96, all p < 0.05) and 59 (60%) had a coefficient above 0.5 (Table 1; Figure 2). For angular amplitude measures, all but one correlation coefficients (39/40) were significant (p ≤ 0.01). All significant coefficients were above 0.3 and 28 (70%) exceeded 0.5. Correlation coefficients were significant and above 0.3 (p < 0.05) in 31/39 (79%) cases for angular velocity measures. The coefficients for angular velocity measures were larger than 0.5 in 16 (41%) cases. All correlations coefficients were significant and above 0.3 (p < 0.05) for muscle activity measures and 14/19 were above 0.5 (74%).
TABLE 1. Correlations in angular amplitude, angular velocity and muscle activity between different functional tasks. The darkness of the blue represents the strength of the correlation: white: r < 0.3; pale blue: 0.3 ≤ r < 0.5; blue: 0.5 ≤ r < 0.7; dark blue: r ≥ 0.7. Correlations in blue (r ≥ 0.3) all have a p-value <0.05. NA: variables not available as there were no characteristic pattern (see Data Processing). STS: sit-to-stand; LLS: Lower lumbar spine; ULS: Upper lumbar spine. EMGpeak1: first peak of maximal paraspinal muscles activity; EMGpeak2: second peak of maximal paraspinal muscles activity. Correlations between EMGpeak2 of flexion, lifting and picking-up and EMGpeak1 of stepping-up and sit-to-stand are reported in the last line (during stepping-up and sit-to-stand, there is only one peak of paraspinal muscle activity, see Data Processing).
FIGURE 2. Scatterplots of LLSaflexion among the five specific functional tasks as well as with respect to the task-independent LLSaflexion measure in patients with CLBP. Task-specific measures are in degree and task-independent measures are dimensionless.
The AUC with their 95% confidence interval (CI), as well as the other performance values, of each movement measure in each functional task are presented in Tables 2–4. For the angular amplitude measures, five comparisons between tasks out of 40 (13%) reported a difference of AUC above 0.1 (range 0.11–0.15). Performance to differentiate groups during sit-to-stand compared to flexion and lifting was higher for ULSarange (AUC of 0.74 compared to AUC of 0.63 and 0.59, respectively), but was lower for LLSarange (AUC of 0.6 compared to AUC of 0.75 for flexion and AUC of 0.72 for lifting). AUC was also lower during lifting (AUC of 0.66) compared to flexion (AUC of 0.77) for LLSaflexion. Regarding angular velocity measures, 7 comparisons out of 42 (17%) reported a difference of AUC above 0.1 (range 0.11–0.14). For LLSvflexion, differentiation performance was smaller during sit-to-stand (AUC of 0.76) compared to picking-up, lifting and flexion (0.87 ≤ AUC ≤0.90), and smaller during stepping-up (AUC of 0.79) compared to picking-up (AUC of 0.90). AUC during picking-up for LLSvextension (AUC of 0.88) was higher than during stepping-up (AUC of 0.75), and was also higher for LLSvrange (AUC of 0.92) compared to stepping-up (AUC of 0.79) and sit-to-stand (AUC of 0.81). Five out of 14 comparisons (36%) for muscle activity measures reported a difference in AUC above 0.1 (range 0.13–0.17). AUC during stepping-up (AUC of 0.51) was lower than during picking-up, lifting and flexion (0.64 ≤ AUC ≤0.68) for EMGpeak1. For EMGpeak2, flexion and lifting had higher AUC (AUC of 0.60 and 0.59, respectively) compared to picking-up (AUC of 0.47).
TABLE 2. Logistic regression models for angular amplitude measures. AUC: area under the curve [with its confidence interval (95%CI)]; r2: coefficient of determination; Sn: sensitivity; Sp: specificity; LR+: positive likelihood ratio; LR−: negative likelihood ratio. Symbols (°, §, #, *) indicate a difference of AUC >0.1 between two tasks or between a specific-task and the task–independent measure.
TABLE 3. Logistic regression models for angular velocity measures. AUC: area under the curve [with its confidence interval (95%CI)]; r2: coefficient of determination; Sn: sensitivity; Sp: specificity; LR+: positive likelihood ratio; LR−: negative likelihood ratio. Symbols (°, §, #, *) indicate a difference of AUC >0.1 between two tasks or between a specific-task and the task–independent measure.
TABLE 4. Logistic regression models for lumbar muscle activity measures. AUC: area under the curve [with its confidence interval (95%CI)]; r2: coefficient of determination; Sn: sensitivity; Sp: specificity; LR+: positive likelihood ratio; LR−: negative likelihood ratio. Symbols (°, §, #, *) indicate a difference of AUC >0.1 between two tasks or between a specific-task and the task–independent measure.
Results from independent t-tests and ES regarding the differences in movement measures between groups are reported in Figure 3 and Supplementary Appendix SB. Patients with CLBP moved with statistically significantly reduced LLSaflexion and LLSarange in all functional tasks compared to asymptomatic controls and all ES were moderate to large. At the ULS, ULSarange during picking-up, stepping-up and sit-to-stand were significantly reduced in patients with CLBP, with moderate to large ES for the three tasks. Angular velocity measures were all significantly reduced in patients with CLBP, and ES were at least large. For muscle activity measures, EMGpeak1 during flexion and lifting was significantly higher in patients with CLBP. ES ranged from 0.18 to 0.6 for EMGpeak1 and from 0.02 to 0.42 for EMGpeak2.
FIGURE 3. Effect sizes of the differences between patients with CLBP and asymptomatic controls. *: p-value<0.05 at t-tests. No ES means that the variable was not available because there was no characteristic pattern (see Data processing).
Task–independent measures reported AUC between 0.59 and 0.71, between 0.88 and 0.94 and between 0.55 and 0.64 for angular amplitude, angular velocity and muscle activity measures, respectively. Compared to each task, performance to differentiate groups for task–independent measures was superior (AUC difference >0.1) in 7/52 cases (13%) and was not different (AUC difference ≤0.1) in 45/52 cases (87%) (Tables 2–4). Correlation coefficients between task-independent and task-specific measures were large for angular amplitudes (range 0.63–0.95, all p < 0.001), angular velocities (range 0.74–0.92, all p < 0.001) and muscle activities (range 0.70–0.93, all p < 0.001) (Supplementary Appendix SC). Groups differences were statistically significant for LLSaflexion, LLSarange, ULSarange and for all angular velocity measures. ES for angular amplitude measures ranged from 0.33 to 0.97, for angular velocity measures from 1.54 to 2.36 and for muscle activity measures from 0.07 to 0.51 (Figure 3 and Supplementary Appendix SB).
Discussion
The correlation and capacity to differentiate patients from controls results indicated a limited influence of the tasks on the information carried by spinal movement measures and highlighted an individual spinal movement signature. This study also showed the construct validity of task–independent measures of spinal movement, and encouraged its consideration in future research. These important findings are discussed in the following sections.
Consistency of Spinal Movement in Different Tasks
Ninety-two percent of the measures were significantly correlated between different functional tasks in patients with CLBP, with correlation coefficients demonstrating at least a medium effect (r ≥ 0.3). The coefficients were even large in 70, 41 and 76% of angular amplitude, angular velocity and muscle activity measures, respectively. Correlation coefficients tended to be larger for peak flexion angle at the lower (LLSaflexion) and upper lumbar spine (ULSaflexion), which is consistent with what was found in asymptomatic individuals (Alqhtani et al., 2015; Seerden et al., 2019). Interestingly, the present study also demonstrated individual consistency in lumbar angular velocity and in the level of muscle activity, suggesting that the consistency across tasks is not limited to angular amplitudes as previously shown (Alqhtani et al., 2015; Seerden et al., 2019). Correlation coefficients were very large between analogous functional tasks, such as lifting and picking-up, but consistency was also observed between differing tasks, such as flexion and sit-to-stand. While the large standard deviation in most measures of spinal movement in patients with CLBP showed heterogeneity between participants, the correlations between tasks demonstrated consistency within patients.
These findings support an individual spinal movement signature in the sagittal plane in patients with CLBP, suggesting that each individual has a consistent spinal movement across tasks. These results questioned the need to analyze multiple tasks involving primarily sagittal-plane movement independently and the need to investigate new primarily sagittal-plane functional tasks in future studies, as those may lead to redundant data, therefore complexifying the procedure without gaining information to improve our understanding of spinal movement in CLBP (Papi et al., 2018; Moissenet et al., 2021). These findings also suggest that spinal movement is probably more influenced by individual factors than by the tasks. However, which individual factors are associated with the spinal movement signature is still unclear. Previous research in patients with LBP reported association between spinal kinematics and pain intensity, psychological characteristics, sex, age and BMI, among others (Mitchell et al., 2008; Arshad et al., 2019; Christe et al., 2021). However, these factors demonstrated only small associations with spinal movement, suggesting that other unknown factors are in play. Consequently, further research is needed to better understand which factors influence the spinal movement signature in patients with CLBP.
Difference in Performance to Differentiate Groups Between Functional Tasks
The functional tasks had little influence on the performance to differentiate patients with CLBP from asymptomatic controls. Less than 20% of the comparisons demonstrated a difference in AUC between tasks larger than 0.1. For angular amplitude measures, performance to differentiate patients from controls was consistent across the tasks in 79% of the comparisons. Independent t-tests, ES and other outcomes of models’ performance also supported these findings. When there was a difference, it was in the ranges of motion and not the peak amplitudes. However, differences were inconsistent among tasks, as sit-to-stand demonstrated higher performance at the ULSarange and poorer performance at the LLSarange compared to flexion and lifting. Regarding angular velocity measures, AUC did not differ substantially in 83% of the comparisons. When it differed, it consistently showed poorer performance during sit-to-stand and stepping-up. Although stepping-up and sit-to-stand also reported mostly very large effect sizes, their capacity to differentiate groups was smaller compared to flexion, lifting and picking-up. Regarding muscle activity measures, the capacity to differentiate groups was poor in all tasks and ES were moderate at most. The performance to differentiate groups was poorer during stepping-up for EMGpeak1 and during picking-up for EMGpeak2.
Globally, these observations indicate that assessing a range of functional tasks in the sagittal-plane may provide similar findings in terms of differentiating patients with CLBP from asymptomatic controls. The differences found between sit-to-stand and stepping-up compared to flexion, lifting and picking-up for angular velocity measures may be explained by the fact that participants with CLBP rated flexion, lifting and picking-up as more painful and more harmful for the back than stepping-up and sit-to-stand. The perceived harm might have led to increased pain-related fear and, together with the higher pain intensity, influenced angular velocity. Yet, these differences did not seem to consistently influence the angular amplitudes, questioning the effect of pain intensity and pain-related fear on lumbar angular amplitude. These findings are in line with a recent meta-analysis that showed a very small association between pain-related fear or pain intensity and spinal angular amplitudes, which was consistent across a wide range of tasks and measures of spinal angular amplitude (Christe et al., 2021). Therefore, although this study showed a redundancy among the tasks, it is possible that some functional tasks may be more sensitive, for example with respect to pain intensity or pain-related fear. In this regard, some authors suggested that selecting a specific task for each individual based on their identified limitations could be helpful when analyzing the relationships between spinal movement and patient-related outcomes (Wernli et al., 2020a). It is currently unknown if selecting one sagittal-plane task based on the individual limitation would be more appropriate to analyze such relationships and future studies should address this gap. If it would be the case, this would support the assessment of one specific primarily sagittal-plane task; the others tasks complexifying the procedure in vain (without brining supplementary information). If not, this might support the use of task–independent measures as discussed below.
Task–Independent Measures of Spinal Movement
In this study, we showed the possibility and potential of averaging spinal movement measures across different functional tasks. Averaging spinal movement measures across multiple tasks was particularly supported by the individual spinal movement signature found in this study. The performance to differentiate groups was higher or did not differ compared to individual tasks in all the measures. Furthermore, angular velocity task–independent measures showed a high performance to differentiate groups, with statistically significant differences between the groups and very large to huge ES. ES were also large for the angular amplitude measures at the lower lumbar spine. Therefore, these results support the construct validity of task–independent measures of spinal movement and its consideration in future research.
Using an “average” measure across different tasks may have some interest in future studies. The method we used is simple as it consisted in averaging the Z-scores of each of the task, which could be easily replicated with any other biomechanical model or measure. These task–independent measures could notably be more robust because they are not reliant on a single task.
Nevertheless, future research is strongly recommended to determine the value of task–independent measures. Based on this study, it is not known if task–independent measures can provide more information than task-specific measures, nor how many movements should be averaged. Furthermore, reliability of averaging spinal movement measures from different functional tasks remains to be tested. While the method of averaging tasks using Z-scores has the advantage of its simplicity, other more advanced methods to group spinal movement measures (i.e., machine learning methods) will also need to be investigated.
Capacity to Differentiate Patients With CLBP From Asymptomatic Controls
While it was not the objective of this study, the capacity of spinal movement measures to differentiate patients with CLBP from asymptomatic controls is worth discussing. First, patients with CLBP moved with reduced sagittal-plane lumbar amplitude and range of motion at the lower lumbar spine, in all functional tasks. ES were moderate to large. However, the low specificity and LR + suggested that small amplitudes are also frequent in asymptomatic controls. Second, angular velocity measures demonstrated very large ES in the majority of the tasks, which were always larger than the ES from angular amplitude or lumbar muscle activity measures. The capacity of angular velocity measures to differentiate patients from controls was even rated as outstanding for flexion, lifting and picking-up. The high sensitivity and very low LR− in all tasks showed that moving with high angular velocity was very rare in our sample of patients with CLBP, suggesting that moving at high angular velocity is very difficult with CLBP. Third, selected peaks of erector spinae activity demonstrated poor performance to differentiate groups. There were only two muscle activity measures that showed a statistically significant difference between the groups, and most ES were small. These findings are in agreement with previous studies analyzing erector spinae activity during dynamic tasks and reporting inconsistent results (Geisser et al., 2005). Yet, when a difference in muscle activity was observed between groups, it corresponded to higher levels of activity in patients with CLBP.
Overall, based on the present findings and previous reports (Shum et al., 2005; Laird et al., 2014; Christe et al., 2016b, 2020; Papi et al., 2018; Moissenet et al., 2021), reduced lumbar amplitude and angular velocity seems to be key characteristics of patients with CLBP. Furthermore, the consistent reduced sagittal-plane lumbar angular amplitude and velocity across all the functional tasks suggest that these spinal movement alterations generalize across a wide range of daily-life activities. Our results thus support the measurement of lumbar angular amplitude and angular velocity in any functional task in future studies. Nevertheless, there is an urgent need for well-conducted longitudinal studies to detangle if and how spinal kinematic changes are associated with patients’ changes in pain and disability (Wernli et al., 2020b; Schmid et al., 2021).
Limitations
This study has some limitations that are important to discuss. First, the asymptomatic population was small, despite a number of participants above the minimum indicated by the sample size calculation (Portney and Watkins, 2000). While the number of participants had certainly little influence on the task comparisons, the performance values from the logistic regression models should be confirmed with larger groups. Second, the findings may not be transferable to all patients with CLBP. Although our results are consistent with current knowledge in the field, patients with CLBP included in this study had high levels of disability, pain-related fear and catastrophizing that are common in patients participating to interdisciplinary rehabilitation programs. These individual factors may have influenced the large differences found between the groups. Third, the tasks were not assessed in random order. Therefore, the higher level of pain found during picking-up and lifting may be related to the fact that these two tasks were collected last (order effect). Fourth, the video recordings used to present the daily-life tasks may have influenced how participants performed the tasks, requesting caution when interpreting the findings in the context of movement behavior. Video recordings were used to give standardized instructions and avoid differences in the ways of completing the tasks. This was particularly important for picking-up and lifting, as these tasks can be performed in different ways (i.e., stoop or squat). While the instructions could have limited inter-individual variability, large variations were observed among individuals, suggesting that participants were not too constrained and could express their individual movement signature. On the other hand, one cannot exclude that more pronounced signatures could have been observed if the tasks would have been less standardized. Finally, this study focused on sagittal-plane lumbar biomechanics during functional tasks primarily involving sagittal-plane movements, because this corresponds to the alterations the most frequently reported in literature and because patients often complain of movement-related pain in primarily sagittal-plane activities. Therefore, it has yet to be determined if other functional tasks with larger solicitations in the frontal and/or transverse–planes, such as gait (Christe et al., 2017; Schmid et al., 2017), would display individual movement signatures in these other planes or even three-dimensionally.
Conclusion
This study showed that individuals with CLBP have consistent spinal movement in the sagittal plane across different functional tasks, supporting the existence of an individual biomechanical signature. Furthermore, the capacity to differentiate patients with CLBP from asymptomatic controls did not differ between functional tasks in most of the cases. Therefore, this study highlighted a redundancy among tasks, questioning the most appropriate measures to describe spinal movement behavior in the framework of CLBP. While further research will be necessary in this regard, this study showed the feasibility of task-independent measures, a promising approach towards an effective quantification of spinal movement.
Data Availability Statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Research Ethics Committee of VAUD (CER-VD 2018-00188). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
GC contributed to the conception of the study, collected the data, analyzed the data, interpreted the results and wrote the initial manuscript. CA processed all the data and edited the manuscript. BJ participated to the conception of the study, edited the manuscript and supervised the project. JF contributed to the conception of the study, interpreted the results, edited the manuscript and supervised the project. All authors approved the final version of the manuscript. Both BJ and JF supervised this research and should be considered as last authors.
Funding
This project was partly supported by a grant from the University of Applied Sciences and Arts Western Switzerland//HES-SO, Faculty of Health Science.
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.
Publisher’s Note
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.
Acknowledgments
We would like to thanks Valentina Pizzolato and Benjamin Coty for their important contribution in markers’ tracking and data management.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2021.767974/full#supplementary-material
References
Alqhtani, R. S., Jones, M. D., Theobald, P. S., and Williams, J. M. (2015). Correlation of Lumbar-Hip Kinematics Between Trunk Flexion and Other Functional Tasks. J. Manipulative Physiol. Ther. 38, 442–447. doi:10.1016/j.jmpt.2015.05.001
Arshad, R., Pan, F., Reitmaier, S., and Schmidt, H. (2019). Effect of Age and Sex on Lumbar Lordosis and the Range of Motion. A Systematic Review and Meta-Analysis. J. Biomech. 82, 1–19. doi:10.1016/j.jbiomech.2018.11.022
Bujang, M. A., and Baharum, N. (2016). Sample Size Guideline for Correlation Analysis. World J. Soc. Sci. 3, 37. doi:10.22158/wjssr.v3n1p37
Chapman, J. R., Norvell, D. C., Hermsmeyer, J. T., Bransford, R. J., DeVine, J., McGirt, M. J., et al. (2011). Evaluating Common Outcomes for Measuring Treatment Success for Chronic Low Back Pain. Spine. 36, S54–S68. doi:10.1097/BRS.0b013e31822ef74d
Christe, G., Crombez, G., Edd, S., Opsommer, E., Jolles, B. M., and Favre, J. (2021). Relationship Between Psychological Factors and Spinal Motor Behaviour in Low Back Pain: a Systematic Review and Meta-Analysis. Pain. 162, 672–686. doi:10.1097/j.pain.0000000000002065
Christe, G., Jolles, B. M., and Favre, J. (2019). “Variability of Lumbar Flexion Angle During a Pick up Task in Patients With Chronic Low Back Pain,” in World Conference on Low Back Pain (Antwerp).
Christe, G., Kade, F., Jolles, B. M., and Favre, J. (2017). Chronic Low Back Pain Patients Walk With Locally Altered Spinal Kinematics. J. Biomech. 60, 211–218. doi:10.1016/j.jbiomech.2017.06.042
Christe, G., Redhead, L., Jolles-Haeberli, B., and Favre, J. (2016a). Correlation of Regional Lumbar Kinematics Between Forward Flexion and Functional Activities in Patients With Chronic Low Back Pain. Man. Ther. 25, e78. doi:10.1016/j.math.2016.05.127
Christe, G., Redhead, L., Legrand, T., Jolles, B. M., and Favre, J. (2016b). Multi-segment Analysis of Spinal Kinematics During Sit-To-Stand in Patients With Chronic Low Back Pain. J. Biomech. 49, 2060–2067. doi:10.1016/j.jbiomech.2016.05.015
Christe, G., Rochat, V., Jolles, B. M., and Favre, J. (2020). Lumbar and Thoracic Kinematics during Step‐up: Comparison of Three‐Dimensional Angles between Patients With Chronic Low Back Pain and Asymptomatic Individuals. J. Orthop. Res. 38, 1248–1256. doi:10.1002/jor.24575
Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences. Hillsdale, N.J.: L. Erlbaum Associates. Available at: http://www.amazon.com/books/dp/0805802835.
Dankaerts, W., O’Sullivan, P. B., Burnett, A. F., Straker, L. M., and Danneels, L. A. (2004). Reliability of EMG Measurements for Trunk Muscles During Maximal and Sub-Maximal Voluntary Isometric Contractions in Healthy Controls and CLBP Patients. J. Electromyogr. Kinesiol. 14, 333–342. doi:10.1016/j.jelekin.2003.07.001
Dankaerts, W., O’Sullivan, P., Burnett, A., Straker, L., Davey, P., and Gupta, R. (2009). Discriminating Healthy Controls and Two Clinical Subgroups of Nonspecific Chronic Low Back Pain Patients Using Trunk Muscle Activation and Lumbosacral Kinematics of Postures and Movements. Spine. 34, 1610–1618. doi:10.1097/BRS.0b013e3181aa6175
Dubois, J.-D., Abboud, J., St-Pierre, C., Piché, M., and Descarreaux, M. (2014). Neuromuscular Adaptations Predict Functional Disability Independently of Clinical Pain and Psychological Factors in Patients With Chronic Non-Specific Low Back Pain. J. Electromyogr. Kinesiol. 24, 550–557. doi:10.1016/j.jelekin.2014.04.012
Dupeyron, A., Demattei, C., Kouyoumdjian, P., Missenard, O., Micallef, J. P., and Perrey, S. (2013). Neuromuscular Adaptations After a Rehabilitation Program in Patients With Chronic Low Back Pain: Case Series (Uncontrolled Longitudinal Study). BMC Musculoskelet. Disord. 14, 277. doi:10.1186/1471-2474-14-277
Fairbank, J. C. T., and Pynsent, P. B. (2000). The Oswestry Disability Index. Spine. 25, 2940–2953. doi:10.1097/00007632-200011150-00017
Geisser, M. E., Ranavaya, M., Haig, A. J., Roth, R. S., Zucker, R., Ambroz, C., et al. (2005). A Meta-Analytic Review of Surface Electromyography Among Persons With Low Back Pain and Normal, Healthy Controls. The J. Pain. 6, 711–726. doi:10.1016/j.jpain.2005.06.008
Ghasemi, A., and Zahediasl, S. (2012). Normality Tests for Statistical Analysis: A Guide for Non-Statisticians. Int. J. Endocrinol. Metab. 10, 486–489. doi:10.5812/ijem.3505
Grood, E. S., and Suntay, W. J. (1983). A Joint Coordinate System for the Clinical Description of Three-Dimensional Motions: Application to the Knee. J. Biomech. Eng. 105, 136–144. doi:10.1115/1.3138397
Hidalgo, B., Gilliaux, M., Poncin, W., and Detrembleur, C. (2012). Reliability and Validity of a Kinematic Spine Model During Active Trunk Movement in Healthy Subjects and Patients With Chronic Non-Specific Low Back Pain. J. Rehabil. Med. 44, 756–763. doi:10.2340/16501977-1015
Hosmer, D., Lemeshow, S., and Sturdivant, R. (2013). Applied Logistic Regression. Third edit. New York, NY: Wiley NJ.
Kadaba, M. P., Ramakrishnan, H. K., Wootten, M. E., Gainey, J., Gorton, G., and Cochran, G. V. B. (1989). Repeatability of Kinematic, Kinetic, and Electromyographic Data in normal Adult Gait. J. Orthop. Res. 7, 849–860. doi:10.1002/jor.1100070611
Laird, R. A., Gilbert, J., Kent, P., and Keating, J. L. (2014). Comparing Lumbo-Pelvic Kinematics in People With and Without Back Pain: a Systematic Review and Meta-Analysis. BMC Musculoskelet. Disord. 15, 229. doi:10.1186/1471-2474-15-229
Lima, M., Ferreira, A. S., Reis, F. J. J., Paes, V., and Meziat-Filho, N. (2018). Chronic Low Back Pain and Back Muscle Activity During Functional Tasks. Gait & Posture. 61, 250–256. doi:10.1016/j.gaitpost.2018.01.021
Marras, W. S., Lavender, S. A., Leurgans, S. E., Fathallah, F. A., Ferguson, S. A., Gary Allread, W., et al. (1995). Biomechanical Risk Factors for Occupationally Related Low Back Disorders. Ergonomics. 38, 377–410. doi:10.1080/00140139508925111
Matheve, T., De Baets, L., Bogaerts, K., and Timmermans, A. (2019). Lumbar Range of Motion in Chronic Low Back Pain Is Predicted by Task‐Specific, but Not by General Measures of Pain‐related Fear. Eur. J. Pain. 23, 1171–1184. doi:10.1002/ejp.1384
Mitchell, T., O'Sullivan, P. B., Burnett, A. F., Straker, L., and Smith, A. (2008). Regional Differences in Lumbar Spinal Posture and the Influence of Low Back Pain. BMC Musculoskelet. Disord. 9, 152. doi:10.1186/1471-2474-9-152
Moissenet, F., Rose-Dulcina, K., Armand, S., and Genevay, S. (2021). A Systematic Review of Movement and Muscular Activity Biomarkers to Discriminate Non-Specific Chronic Low Back Pain Patients From an Asymptomatic Population. Sci. Rep. 11, 1–14. doi:10.1038/s41598-021-84034-x
O’Sullivan, P. (2005). Diagnosis and Classification of Chronic Low Back Pain Disorders: Maladaptive Movement and Motor Control Impairments as Underlying Mechanism. Man. Ther. 10, 242–255. doi:10.1016/j.math.2005.07.001
Papi, E., Bull, A. M. J., and McGregor, A. H. (2018). Is There Evidence to Use Kinematic/kinetic Measures Clinically in Low Back Pain Patients? A Systematic Review. Clin. Biomech. 55, 53–64. doi:10.1016/j.clinbiomech.2018.04.006
Portney, L. G., and Watkins, M. P. (2000). Foundations of Clinical Research: Applications to Practice. Upper Saddle River, NJ: Prentice Hall Health.
Sawilowsky, S. S. (2009). New Effect Size Rules of Thumb. J. Mod. App. Stat. Meth. 8, 597–599. doi:10.22237/jmasm/1257035100
Schmid, S., Bangerter, C., Schweinhardt, P., and Meier, M. L. (2021). Identifying Motor Control Strategies and Their Role in Low Back Pain: A Cross-Disciplinary Approach Bridging Neurosciences With Movement Biomechanics. Front. Pain Res. 2, 1–8. doi:10.3389/fpain.2021.715219
Schmid, S., Bruhin, B., Ignasiak, D., Romkes, J., Taylor, W. R., Ferguson, S. J., et al. (2017). Spinal Kinematics During Gait in Healthy Individuals Across Different Age Groups. Hum. Movement Sci. 54, 73–81. doi:10.1016/j.humov.2017.04.001
Seay, J., Selbie, W. S., and Hamill, J. (2008). In Vivolumbo-Sacral Forces and Moments During Constant Speed Running at Different Stride Lengths. J. Sports Sci. 26, 1519–1529. doi:10.1080/02640410802298235
Seerden, S. F. L., Dankaerts, W., Swinnen, T. W., Westhovens, R., de Vlam, K., and Vanwanseele, B. (2019). Multi-Segment Spine and Hip Kinematics in Asymptomatic Individuals During Standardized Return From Forward Bending Versus Functional Box Lifting. J. Electromyogr. Kinesiol. 49, 102352. doi:10.1016/j.jelekin.2019.102352
Shum, G. L. K., Crosbie, J., and Lee, R. Y. W. (2005). Symptomatic and Asymptomatic Movement Coordination of the Lumbar Spine and Hip During an Everyday Activity. Spine. 30, E697–E702. doi:10.1097/01.brs.0000188255.10759.7a
Shum, G. L. K., Crosbie, J., and Lee, R. Y. W. (2007). Movement Coordination of the Lumbar Spine and Hip During a Picking up Activity in Low Back Pain Subjects. Eur. Spine J. 16, 749–758. doi:10.1007/s00586-006-0122-z
Sullivan, M. J. L., Bishop, S. R., and Pivik, J. (1995). The Pain Catastrophizing Scale: Development and Validation. Psychol. Assess. 7, 524–532. doi:10.1037/1040-3590.7.4.524
Vlaeyen, J. W. S., Kole-Snijders, A. M. J., Boeren, R. G. B., and van Eek, H. (1995). Fear of Movement/(re)injury in Chronic Low Back Pain and its Relation to Behavioral Performance. Pain. 62, 363–372. doi:10.1016/0304-3959(94)00279-N
Vogler, D., Paillex, R., Norberg, M., de Goumoëns, P., and Cabri, J. (2008). Validation Transculturelle de l'Oswestry Disability Index en Français. Ann. de Réadaptation de Médecine Physique. 51, 379–385. doi:10.1016/j.annrmp.2008.03.006
von Elm, E., Altman, D. G., Egger, M., Pocock, S. J., Gøtzsche, P. C., and Vandenbroucke, J. P. (2007). The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for Reporting Observational Studies. The Lancet. 370, 1453–1457. doi:10.1016/S0140-6736(07)61602-X
Wernli, K., O'Sullivan, P., Smith, A., Campbell, A., and Kent, P. (2020a). Movement, Posture and Low Back Pain. How Do They Relate? A Replicated Single‐Case Design in 12 People With Persistent, Disabling Low Back Pain. Eur. J. Pain. 24, 1831–1849. doi:10.1002/ejp.1631
Wernli, K., Tan, J.-S., O'Sullivan, P., Smith, A., Campbell, A., and Kent, P. (2020b). Does Movement Change When Low Back Pain Changes? A Systematic Review. J. Orthop. Sports Phys. Ther. 50, 664–670. doi:10.2519/jospt.2020.9635
Keywords: low back pain, motion analysis, lumbar, kinematics, electromyography, angle, angular velocity, muscle activity
Citation: Christe G, Aussems C, Jolles BM and Favre J (2021) Patients With Chronic Low Back Pain Have an Individual Movement Signature: A Comparison of Angular Amplitude, Angular Velocity and Muscle Activity Across Multiple Functional Tasks. Front. Bioeng. Biotechnol. 9:767974. doi: 10.3389/fbioe.2021.767974
Received: 31 August 2021; Accepted: 26 October 2021;
Published: 15 November 2021.
Edited by:
Stefan Schmid, Bern University of Applied Sciences, SwitzerlandReviewed by:
Thomas Matheve, Universiteit Hasselt, BelgiumAndré P. G. Castro, Universidade de Lisboa, Portugal
Copyright © 2021 Christe, Aussems, Jolles and Favre. 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: Guillaume Christe, Guillaume.christe@hesav.ch