Abstract
The central nervous system (CNS) dynamically employs a sophisticated weighting strategy of sensory input, including vision, vestibular and proprioception signals, towards attaining optimal postural control during different conditions. Non-specific low back pain (NSLBP) patients frequently demonstrate postural control deficiencies which are generally attributed to challenges in proprioceptive reweighting, where they often rely on an ankle strategy regardless of postural conditions. Such impairment could lead to potential loss of balance, increased risk of falling, and Low back pain recurrence. In this study, linear and non-linear indicators were extracted from center-of-pressure (COP) and trunk sagittal angle data based on 4 conditions of vibration positioning (vibration on the back, ankle, none or both), 2 surface conditions (foam or rigid), and 2 different groups (healthy and non-specific low back pain patients). Linear discriminant analysis (LDA) was performed on linear and non-linear indicators to identify the best sensory condition towards accurate distinction of non-specific low back pain patients from healthy controls. Two indicators: Phase Plane Portrait ML and Entropy ML with foam surface condition and both ankle and back vibration on, were able to completely differentiate the non-specific low back pain groups. The proposed methodology can help clinicians quantitatively assess the sensory status of non-specific low back pain patients at the initial phase of diagnosis and throughout treatment. Although the results demonstrated the potential effectiveness of our approach in Low back pain patient distinction, a larger and more diverse population is required for comprehensive validation.
1 Introduction
Likely driven by ageing and population increase, the Global Burden of disease Report reveals a significant increase in both the number of Low back pain (LBP) patients and the prevalence of LBP in all age groups worldwide in the last 2 decades (1990–2019). The same report also points out that although the prevalence of LBP increases with increasing age, currently, the 50–54 years age group has the highest LBP prevalence (Vos et al., 2021). Expectedly, such trend in an active population is associated with grave health and economic consequences. (Airaksinen et al., 2006; Balagué et al., 2012). Globally, LBP continues to be the main cause of years lived with disability or YLDs. Health economists have estimated that the caring costs for 15% of people with low back pain are equivalent to 85% of those associated with the remaining population (Hashemi et al., 1997; Hashemi et al., 1998; Filiz et al., 2005). In the US alone, health economists estimate that up to 624 billion dollars are spent annually on direct and indirect expenses associated with LBP (Dagenais et al., 2008). Importantly, it has been shown that there is no particular identifiable pathoanatomical etiology for LBP in 85–95 percent of patients who visit primary care doctors, hence the common diagnosis of non-specific low back pain (NSLBP) (Airaksinen et al., 2006; Hancock et al., 2011). The exact mechanisms and causes of NSLBP, which remain elusive, are therefore highly relevant in low back pain investigations (Costa et al., 2013).
Optimal postural control is task-dependent and dynamically variable during the performance of different daily activities. The human brain performs diverse targeted control adjustments on the ankle, knee, hip and spine in order to maintain postural stability (Allum et al., 1998). The central nervous system (CNS), in turn, adapts and selects the appropriate sensory input signals according to changes in postural conditions. The CNS, subsequently, reweights (dynamically adjusts the weight assigned to a particular signal) visual, vestibular, and proprioceptive input to generate the appropriate muscle forces required for the particular task. This allows effective control of the center of mass, resulting in proper equilibrium of the body (Brumagne et al., 2004; Carver et al., 2006). Previous investigations explored different models and methods for preserving equilibrium (Horak and Nashner, 1986; Runge et al., 1999). Several hypothesized postural control strategies can be interpreted by the infamous inverted pendulum model. The ankle strategy suggests that stability is maintained by balancing the body initially around the ankle joint (Horak and Nashner, 1986). This strategy is adequate for simple conditions, such as standing on a rigid surface, but is less successful in more complex conditions, including foam and other unstable surfaces. According to the inverted model theory, sufficient motion has to be generated at the trunk and the hip joint (hip strategy) for optimal postural control during more complex conditions (Horak and Nashner, 1986). On the other hand, multi segmental postural control models delineate that postural control is provided not only by specific modifications at one joint, but also include various modifications at different joints performed by the CNS (Morasso and Schieppati, 1999; Schieppati et al., 2002; Kiemel et al., 2008).
The literature demonstrates controversy regarding the causal association between postural control and LBP. Several studies concluded that postural control Impairment can be considered as one of the causal factors associated with non-specific low back pain (NSLBP) (Peterka, 2002; Horak, 2006; Ruhe et al., 2011), while others maintained that that there is no consensus regarding the link between posture and low back pain (Swain et al., 2020). Although NSLBP Patients have demonstrated impaired postural control during complicated conditions (standing on unstable support surfaces, eyes closed, etc.), only a few prospective investigations looked into the cause-and-effect relationship between the development of LBP and impaired proprioceptive postural control. One study demonstrated that improper lumbar repositioning and increased posterior pelvic tilt while sitting enhance the likelihood of developing LBP in nursing students (Mitchell et al., 2010). Another study indicated that LBP in college athletes was exacerbated by delayed trunk muscular responses during sitting (Cholewicki et al., 2005). Proprioceptive impairment was implicated as the underlying mechanism in both prospective studies, however direct assessment of the proprioceptive system was missing from both investigations. Due to the lack of evaluation and agreement on the association between posture and low back pain, no clear causation can be construed regarding the role of postural control impairment in the development of LBP. On the other hand, it has been demonstrated that compared to healthy individuals, NSLBP patients exhibit altered proprioceptive postural control, in which the patients adopt a body and trunk stiffening strategy and rely more on ankle proprioception to control their posture during quiet upright standing (Brumagne et al., 2008). This highlights the potential of using the different sensory reweighting strategy as means of distinguishing NSLBP patients from healthy individuals.
Many previous studies focused on the motor output of postural control in NSLBP patients (Ruhe et al., 2011). Others advocated the importance of considering the sensory inputs, particularly in terms of weighting the various proprioceptive signals, a critical consideration in optimal postural control (Gandevia et al., 1992; Lackner and DiZio, 2005; Carver et al., 2006). It is indeed well established that musculoskeletal injuries (Haghighat et al., 2021a; Haghighat et al., 2021b; Haghighat et al., 2021c; Haghighat et al., 2021d; Davoudi et al., 2022a), neuromusculoskeletal disease, such as stroke, Parkinson’s, Multiple Sclerosis, Chronic obstructive pulmonary disease (Davoudi et al., 2022b) and other psychological factors including anxiety, are associated with disrupting the function of the sensory systems (Jamali et al., 2019; Dehmiyani et al., 2022). Change in sensory input function has also been observed among patients with NSLBP (Newcomer et al., 2000). These changes, typically impact several physiological functions including reduced sensory acuity (Sharma and Pai, 1997), altered muscle recruitment patterns (van Dieën et al., 2003; Hodges et al., 2009), and reorganization of the somatosensory regions of the brain cortex (Flor et al., 1997). Peripheral simulation of muscle spindles has revealed that LBP patients rely on their ankle muscle proprioception more than the optimal weighted ankle and back muscle proprioception strategy common in people free of back pain (Brumagne et al., 2008; Claeys et al., 2011). It is believed that in this patient population, the CNS increases ankle proprioception gain, reflecting impaired proprioceptive weighing capacity. Indeed, in challenging postural circumstances, such as standing on unstable surfaces, when the ankle strategy is no longer an appropriate option for the CNS, less reliable postural control is detected in LBP patients (Kiers et al., 2012), and sensory input integration is less optimal due to the decrease in proprioceptive signals from back muscles leading to impaired postural control (della Volpe et al., 2006; Brumagne et al., 2008; Claeys et al., 2011; Mok et al., 2011). Abnormal spinal loading, pain and recurrence of NSLBP are believed to result from such decrease in the proprioceptive function and the consequent impairment of postural control (Claeys, 2013).
In the past several decades, many relevant studies investigated the biomechanics of the low back by using various approaches, tools and technologies to distinguish NSLBP patients from healthy individuals. Marras et al., (1993) and Marras et al., (1995) found significant differences in angular velocity and acceleration discriminating low back pain patients from healthy controls by using a custom-designed three axial goniometer. Classification techniques were also utilized based on coefficients assigned to individuals in eleven subgroups (patients with 10 different levels of low back pain and one healthy cohort) to quantify the presence and intensity of LBP. The combined spinal movement in LBP patients, measured with a three dimensional electromagnetic tracking system, demonstrated reduced movement as compared with healthy individuals (Barrett et al., 1999). Some research groups used Ultrasonic measurement systems (Vogt et al., 2001) to study the impact of non-specific low back pain on 3-D kinematics, while others focused on dynamic pelvic and thoracic oscillations, which were observed to increase as a function of NSLBP. Several studies (Ferreira et al., 2004; van der Hulst et al., 2010; Bervis et al., 2020) focused on abdominal and lumbar muscle activity in chronic low back pain patients, demonstrating increased muscle activity of the erector spinae and rectus abdominis in CLBP patients during walking and moving a flexi bar. Laird et al. (2014) measured the lumbo-pelvic movement in people with and without low back pain, demonstrating that compared to those without LBP, LBP patients showed decreased lumbar ROM, decreased proprioception, and moved more slowly. Scholtes et al. (2009) investigated the lumbopelvic movement in LBP patients and healthy control cohorts measured by a three dimensional motion capture system. They discovered that those with low back pain had an earlier lumbopelvic rotation and a larger maximum lumbopelvic rotation angle compared to healthy individuals during knee flexion and hip lateral rotation. Ashouri et al. (2017) analyzed trunk motion of LBP and healthy controls during several tasks measured by IMU sensors. A Support Vector Machine classifier methodology was utilized for data classification, which enabled the distinction between the two groups with high accuracy and sensitivity. More recently, smart wearables and sensor-based methodologies have yielded encouraging results for the quantification and classification of LBP (Davoudi et al., 2020). These novel technologies, along with big data management tools (AI, ML, SVM, etc.) promise a paradigm shift towards low-cost portable continuous monitoring personalized precision medicine.
Although the literature, including the above-mentioned studies, is rich with quantitative investigations exploring the various dimensions of NSLBP, the multifactorial etiology of the disease and its negative impact on today’s world remain unresolved challenges. Elusive issues include understanding the role of motor control in NSLBP patients and the distinctive motor control strategies adopted by these patients for preserving postural stability. To augment the findings in literature, this investigation was designed to scrutinize the best sensory condition to distinguish NSLBP patients from healthy controls leveraging different surface, and vibration conditions, as well as, linear and non-linear mathematical tools. We hypothesize that the sensory condition of standing on a foam surface with induced vibration, in conjunction with non-linear mathematical tools, can provide a suitable framework to differentiate LBP patients from healthy individuals.
2 Materials and methods
2.1 Ethics statement
The University Internal Ethics Board evaluated and authorized the study involving human subjects (IRB Shahid Beheshti University of Medical Sciences, No: IR.SBMU.RETECH.REC.1396.1392). Prior to beginning the trial, participants were fully informed and given their consent.
2.2 Participants
The experimental set up in this study is similar to our prior investigation on chronic low back pain patients (Shokouhyan et al., 2020) including male participants who suffered LBP at least three times in the recent year or one time in a 3 months duration. The participants were divided evenly between the two groups: an NSLBP group and a healthy control group based on two questionnaires. The ODI (Oswestry Disability Index) (Fairbank and Pynsent, 2000) and the NPRS (quantification of pain) questionnaires were completed by each individual for the assessment of LBP and to rate the pain on a numerical scale, respectively (Joos et al., 1991). Patients who reported an ODI >6 or NPRS >0 were categorized as part of the “NSLBP” group. In this study, however, all males in the “healthy” group reported zero on both the NPRS and ODI questionnaires. The experimental protocol dictated that the experiment would be postponed if any participant reported pain, although none of participants reported any pain during the current experimental tests. The NSLBP patients had to be free of any neurological, respiratory, or cardiovascular disease, as well as any spinal, neck, chest, or lumbar surgeries, in order to meet the inclusion criteria. 5 subjects were excluded due to not meeting the inclusion criteria and the experimental tests were performed on 20 Chronic Non-Specific Low Back Pain (CNSLBP) patients. Demographic data, including age, height, weight, and BMI index are presented in Table 1.
TABLE 1
| Variables | Healthy | CNSLBP | p-value |
|---|---|---|---|
| N (Gender) | 20 (Male) | 20 (Male) | |
| Age | 0.2 | ||
| Height (cm) | 0.316 | ||
| Weight (kg) | 0.26 | ||
| BMI | 0.97 |
Participants’ demographic information for both healthy and low back patients.
2.3 Muscle proprioception
According to prior investigations, muscle vibration is the most common method for altering proprioception inputs (Goodwin et al., 1972; Roll and Vedel, 1982). In order to bias the soleus and lumbar muscles’ proprioception in this investigation, we developed an in-house vibrating apparatus with four brushless DC motors. The triceps surae, which is located in the calf of the lower legs, and the longissimus and multifidus muscles, which span the lumbar vertebrae L3 to L5, were also targeted for placement of the device. Based on previous studies, the vibration frequency of our device was adjusted to 70 Hz, with an amplitude of around 0.5 mm, to provide the best altered proprioceptive data (Goodwin et al., 1972; Roll and Vedel, 1982; Cordo and Gurfinkel, 2004).
2.4 Procedure
A force plate was used for collecting the center-of-pressure (COP) body fluctuation data for each subject (Bertec, Columbus, OH, United States). The kinematic data was captured using a Vicon optical motion capture system with reflective markers synced to the force plate. According to the literature, the markers were attached to the C7, T12, lower sternum (xiphoid process), clavicle (Incisura jugularis), right scapula, right and left sides of the PSIS (posterior superior iliac spine), and ASIS (anterior superior iliac spine). On both devices, the sampling frequency was set to 100 Hz. The motor straps were fastened to the bilateral multifidus muscles and the muscular spindle of the triceps surae on each leg. Each subject completed 8 randomized different trials while having their eyesight obscured (by wearing an eye mask) as follows: 1) standing on a motionless rigid surface (without any vibrator-induced movement); 2) standing on a rigid surface with the activation of the triceps vibrators; 3) standing on a rigid surface with the activation of the multifidus vibrators; 4) standing on a rigid surface with the activation of both the triceps and multifidus vibrators; 5) standing on a motionless foam surface; 6) standing on a foam surface with the activation of the triceps vibrators; 7) standing on a foam surface with the activation of the multifidus vibrators; and 8) standing on a foam surface with the activation of both the triceps and multifidus vibrators. The trunk angles were recorded in the three anatomical planes for each trial, and the COP data was obtained in the anterior posterior (AP) and medial lateral (ML) directions. Each trial lasted 30 s: for the first 10 s of each trial, the subject stood on the force platform without experiencing any vibration (the balancing phase); for the next 20 s, the motors were turned on at a frequency of 70 Hz (the vibration phase).
2.5 Linear indicators of COP time series
According to the literature, a second-order Butterworth non-linear filter was used to filter the data (Ghomashchi et al., 2011) using a 5 Hz cutoff frequency. The mean total velocity, the standard deviation of displacement, and the phase plane portrait for both anterior-posterior (AP) and medial-lateral (ML) directions were obtained in accordance with Table 2 in order to analyze the center-of-pressure data. is the average of the balance time series, represents each point of the vibration time series, and denotes the length of the time series.
TABLE 2
| Parameter | Formula |
|---|---|
| SD of Amplitude (mm) | |
| AP Eq. 1 | |
| ML (Eq. 2) | |
| SD of Amplitude () | |
| AP (Eq. 3) | |
| ML (Eq. 4) | |
| Mean total velocity () | |
| (Eq. 5) | |
| Phase plane portrait (Arbitrary unit) | |
| AP (Eq. 6) | |
| ML (Eq. 7) | |
| AP-ML (Eq. 8) | |
Linear parameters of COP analysis (Salavati et al., 2009).
Unit of measures are as follow: and are in mm, , and are in , , and are in arbitary unit.
2.6 RQA and Lyapunov non-linear indicators of COP time series and trunk sagittal angle
The phase space of the COP and trunk sagittal angle data was reconstructed in MATLAB (MathWorks, Inc., Natick, Massachusetts, United States) using average Mutual Information (AMI) and False Nearest Neighbors (FNN) methods, which are the conventional techniques for finding the time-delay parameter and the embedding dimension parameter, respectively (Horak et al., 2003). The AP and ML directions of the COP, and trunk sagittal angle were the three signals that were used to rebuild the phase space for each subject. The COP and trunk angle’s space embedding dimensions were typically 3 in most cases. The first minimum relative for each participant was considered as the time delay.
By using Recurrence Quantification Analysis (RQA) approach, the features of complexity and amount of recurrence can be calculated as dynamic properties of an observed time series. Riley et al. (1999) extracted numerical criteria based on diagonal lines in recurrence plot (RP). In this study, RQA quantitative measurements were computed by the RQA software (Webber Jr, 2009), developed by Webber Jr and Zbilut (2005). According to Riley et al. (1999), the neighborhood radius was calculated as 2.5 percent of the mean distance using the Euclidean norm.
Upon phase space reconstruction for the COP and trunk angle time series, short and long terms Lyapunov were calculated as stability indicators (more negative indicates more stable) of the dynamic path in reconstructed phase domain according to Rosenstein et al. (1993). Based on the initial slope of the curve during the first few sample intervals, the short-term time () scale was obtained. Similarly, the slope of the function following the rising interval was used to get the long-term Lyapunov () exponent.
2.7 Linear discriminant analysis
Normal distribution was verified using a MATLAB Lillie test ((MathWorks, Inc., Natick, MA, United States) for all the calculated parameters (COP and trunk angle for 2 groups, 4 vibration areas and 2 placement conditions). All the parameters were normally distributed, and hence, analysis of variance (ANOVA) and multiple analysis of variance (MANOVA) tests were used to compare linear and non-linear indicators of the COP and trunk data in order to investigate statistically significant differences. The group (healthy or NSLBP), vibration location (triceps, multifidus, none, and both), and foot positioning condition were taken into consideration as the independent variables in this study (rigid or foam) resulting in a experimental design configuration. At a level of , the results were considered significant (Field, 2013). A scheme was developed (Figure 1) to identify the best sensory conditions and performance parameters to separate the groups. Linear Discriminant Analysis (LDA) was performed on linear, RQA and Lyapunov parameters separately. Representatives of each analysis variable group were compared with each other, and the best dependent variables with the best sensory condition were determined. This procedure can be seen in Figure 1.
FIGURE 1

Flow chart of statistical analysis.
3 Results
The results of the ODI and NPRS questionnaires demonstrate significant differences between the healthy participants and the NSLBP group, as shown in Table 3. The ANOVA results are depicted in Supplementary Tables SA1–SA3. The linear discriminant analysis demonstrated that and Entropy ML were the best parameters to distinguish NSLBP and Healthy groups in this study, with a performance classification percentage of 100% and 70%, respectively, while the foam surface in conjunction with vibration presented the best sensory condition (Figure 2). As shown in Figure 2, NSLBP and Healthy individuals of this study can be successfully discriminated with these two parameters. The two groups can be separated by a line, although the line formula is not clear due to a possible overlap in a bigger membership group.
TABLE 3
| Questioners | Healthy (SD) | Patient (SD) | Significant difference |
|---|---|---|---|
| ODI-2 (0–100) | 0 | 12.3 (3.6) | Yes |
| NPRS (0–10) | 0 | 2.5 (1.2) | Yes |
Oswestry disabity inventory questionnaire and pain scale results from participants.
FIGURE 2

Results of the discriminant analysis.
The correlation Matrix of the significant relationship amongst the eight chosen performance parameters and selected primary condition are indicated in Table 4. It is not surprising that many of these candidates are highly correlated. As the table reveals, the final model that completely separated the two groups kept one measure from the linear and one from the RQA analyses representatives ( and Entropy ML), in conjunction with the foam surface and vibration of both the ankle and back.
TABLE 4
| Discriminant parameters | ||||||||
|---|---|---|---|---|---|---|---|---|
| 1 | 0.94a | 0.77a | -0.1 | 0.47a | -0.71a | 0.53a | 0.6a | |
| 1 | 0.63a | -0.11 | 0.42a | -0.64a | 0.41a | 0.65a | ||
| 1 | 0.3 | 0.57a | -0.78a | 0.77a | 0.6a | |||
| 1 | 0.53a | -0.37b | 0.57a | 0.36b | ||||
| 1 | -0.78a | 0.71a | 0.69a | |||||
| 1 | -0.73a | -0.75a | ||||||
| 1 | 0.74a | |||||||
| 1 |
Correlation matrix for LDA analysis representatives.
Correlation is significant at the 0.01 level.
Correlation is significant at the 0.05 level.
Expresses phase plane portrait in ML direction for trial num#8 (foam + both vibrations).
4 Discussion
This study aimed to identify the best performance indicators and sensory conditions to differentiate NSLBP patients from healthy controls, incorporating various surface conditions (rigid or foam), muscle vibration conditions (triceps, multifidus, none, or both), and leveraging multiple linear (LDA) and non-linear (Lyapunov and RQA)) mathematical tools. The findings indicated that there were significant differences in the majority of the parameters used in this experiment (linear parameters, RQA, and Lyapunov components) between the NSLBP and healthy cohorts. The identification of the best indicators and sensory conditions could enable clinicians to distinguish individuals with proprioception impairment. Our results show that using linear discriminant analysis, only two parameters ( , and Entropy ML), as well as the sensory condition of foam surface with vibration on both the ankle and back, are sufficient to separate two cohorts of Healthy and NSLBP individuals. Although the use of foam surface has been reported exclusively as a sensory condition to discriminate the LBP patients from controls (Salavati et al., 2009), our investigation showed that the addition of vibration has a powerful impact confirming the hypothesis (Figure 2). Importantly, this implies that the current results could provide a direction for future studies aiming to investigate a diagnostic framework for people with proprioception disorders. In our previous work (Davoudi et al., 2020) we classified Chronic Non-specific Low Back pain patients (CNLBP) into three subgroups of low, medium and high risk by linear discriminant analysis based on special sensory conditions and IMU sensors. Our current results suggest that Low Back Pain patients could be potentially distinguished from healthy individuals by employing a sensory condition and linear indicators, although future work is required for further validation.
Psychological variables have been shown to impact mobility and postural balance, in addition to LBP (Pincus et al., 2002). One meta-analysis study found a significant, albeit modest, link between negative emotional characteristics (such as anxiety, heightened fear of movement or pain) and more stiff and protective spinal movement (Christe et al., 2021a). The study suggested further investigation of “specific and personalized assessment of psychological variables, pain severity, and spinal motor behavior” (Christe et al., 2021b). Another study Wernli K et al. (2022) demonstrated that each participant had an individual recovery journey from conscious and non-conscious protection to conscious and non-conscious non-protection. Prior to and following cognitive functional therapy intervention for 12 patients with disabling LBP, pre- and post-quantitative measures of movement, posture, psychological factors, pain, and activity limitation integrated well with the findings from qualitative interviews. Targeting psychological variables may thus be crucial for regaining normal posture and should be further explored in similar studies.
Koch et al. (2020) used COM and surface electromyography in NSLBP and healthy individuals during standing on rigid or foam surface with open and occluded eyes. Although the COM value was higher in low back pain patients, which is consistent with further use of the ankle strategy, the COM and EMG parameters did not show any significant differences between low back pain patients and healthy individuals. By contrast, in this study, most indicators revealed significant differences between the NSLBP patients and healthy groups. It seems therefore that that applying vibration on ankle and back muscles had a important impact on biasing the motor control impairment in non-specific low back pain patients, which enabled the use of simple indicators instead of EMG for group distinction. Nogueira et al. (2020) found multiple discriminant variables based on COP data analysis, although their study did not investigate the effects of unstable surface and/or the effect of vibration on motor control impairment in low back pain patients. Posture control changes radically on an unstable surface (foam) or by biasing the proprioception signals (for example by muscle vibration). Ito et al. (2018) showed that changes in ankle or hip strategy proprioceptive signals lead to changes in NSLBP patients at different frequencies. On the other hand, there are multiple studies which could distinguish the LBP patients form healthy individuals based on trunk motion analysis (McDowell et al., 2018). These studies, however, do not provide much insight on associated motor control impairment. In this work, we employed the two main factors (surface and muscle vibration) with major impact on determining the motor control strategy for preserving the body stability. This approach enables us to compare and differentiate the motor control strategy in low back pain patients with that in healthy individuals towards diagnosing proprioceptive impairment and making informed treatment decisions. In contrast with our hypothesis, linear indicators could be of more value in future studies for group discrimination (100%), while the ML direction seems also to play a key role in the distinction of the two cohorts.
There were several noteworthy limitations in this study. Firstly, linear discriminant analysis requires more input data to correctly examine the differences of indicators in both cohorts of low back pain patients and healthy individuals. This study should therefore be expanded to include more healthy and low back pain patients. Secondly, this investigation was conducted only males due to logistics and hence looking at females and examining gender differences is an important future direction. Another limitation of this study is that only linear discriminant analysis was performed. Machine learning and/or non-linear discriminant methods can be utilized for better accuracy and sensitivity.
5 Conclusion
The main contribution in this study is the development of a relatively simple methodology leveraging linear and non-linear mathematical tools to investigate the main factors (surface and muscle vibration) impacting motor control strategy towards the quantitative discrimination between NSLBP patients with proprioceptive disorders from healthy individuals. LDA analysis identified two indicators (PPPML and EntropyML) on a foam surface, with both ankle and back vibration, that could discriminate the two groups successfully. The current study sheds light on the need to provide a standardized quantitative platform to distinguish NSLBP populations from healthy groups and offers a starting point for future studies to help clinicians assess the sensory status of patients towards preventive medicine and precise patient-centric rehabilitation. Further work is needed to increase the size and diversity of the population studied in order to validate the present results while cautiously recommending it for clinical applications.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
SS, MD, MA, MT, and SiB, contributed to the experimental design and protocol. SS and MA contributed to the analysis of the data. KK, MP, and SiB contributed to the problem statement, introduction and discussion.
Funding
The Sharif University of Technology supported this study. The Biomechanics Laboratory of the Faculty of Rehabilitation at Shahid Beheshti University of Medical Sciences, hosted all of the laboratory experiments used in this study.
Acknowledgments
We appreciate Mr. Hamzeh Asadi and Ms. Zeynab Najafi for their significant contributions to the manufacturing of the vibrating gadget. We also would like to express our gratitude to Dr. Dehghan for her help with the data acquisition.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.1078805/full#supplementary-material
References
1
AiraksinenO.BroxJ. I.CedraschiC.HildebrandtJ.Klaber-MoffettJ.KovacsF.et al (2006). European guidelines for the management of chronic nonspecific low back pain. Eur. spine J.15 (2), s192. 10.1007/s00586-006-1072-1
2
AllumJ.BloemB.CarpenterM.HulligerM.Hadders-AlgraM. (1998). Proprioceptive control of posture: A review of new concepts. Gait posture8 (3), 214–242. 10.1016/s0966-6362(98)00027-7
3
AshouriS.AbediM.AbdollahiM.ManshadiF. D.ParnianpourM.KhalafK. (2017). A novel approach to spinal 3-D kinematic assessment using inertial sensors: Towards effective quantitative evaluation of low back pain in clinical settings. Comput. Biol. Med.89, 144–149. 10.1016/j.compbiomed.2017.08.002
4
BalaguéF.MannionA. F.PelliséF.CedraschiC. (2012). Non-specific low back pain–Authors' reply. Lancet379 (9829), 1874–1875. 10.1016/s0140-6736(12)60804-6
5
BarrettC.SingerK.DayR. (1999). Assessment of combined movements of the lumbar spine in asymptomatic and low back pain subjects using a three-dimensional electromagnetic tracking system. Man. Ther.4 (2), 94–99. 10.1054/math.1999.0175
6
BervisS.KahriziS.ParnianpourM.AmirmoezziY.ShokouhyanS.MoteallehA. (2020). Amplitude of electromyographic activity of trunk and lower extremity muscles during oscillatory forces of flexi-bar on stable and unstable surfaces in people with nonspecific low back pain. J. Biomed. Phys. Eng.12, 521–534. 10.31661/jbpe.v0i0.1912-1022
7
BrumagneS.CordoP.VerschuerenS. (2004). Proprioceptive weighting changes in persons with low back pain and elderly persons during upright standing. Neurosci. Lett.366 (1), 63–66. 10.1016/j.neulet.2004.05.013
8
BrumagneS.JanssensL.KnapenS.ClaeysK.Suuden-JohansonE. (2008). Persons with recurrent low back pain exhibit a rigid postural control strategy. Eur. Spine J.17 (9), 1177–1184. 10.1007/s00586-008-0709-7
9
CarverS.KiemelT.JekaJ. J. (2006). Modeling the dynamics of sensory reweighting. Biol. Cybern.95 (2), 123–134. 10.1007/s00422-006-0069-5
10
CholewickiJ.SilfiesS. P.ShahR. A.GreeneH. S.ReevesN. P.AlviK.et al (2005). Delayed trunk muscle reflex responses increase the risk of low back injuries. Spine30 (23), 2614–2620. 10.1097/01.brs.0000188273.27463.bc
11
ChristeG.CrombezG.EddS.OpsommerE.JollesB. M.FavreJ. (2021a). Relationship between psychological factors and spinal motor behaviour in low back pain: A systematic review and meta-analysis. Pain162 (3), 672–686. 10.1097/j.pain.0000000000002065
12
ChristeG.NzambaJ.DesarzensL.LeubaA.DarlowB.PichonnazC. (2021b). Physiotherapists’ attitudes and beliefs about low back pain influence their clinical decisions and advice. Musculoskelet. Sci. Pract.53, 102382. 10.1016/j.msksp.2021.102382
13
ClaeysK.BrumagneS.DankaertsW.KiersH.JanssensL. (2011). Decreased variability in postural control strategies in young people with non-specific low back pain is associated with altered proprioceptive reweighting. Eur. J. Appl. Physiol.111 (1), 115–123. 10.1007/s00421-010-1637-x
14
ClaeysK. (2013). The role of proprioceptive postural control in the development and maintenance of low back pain: A cross-sectional and a prospective study. Doctoral dissertation. Leuven: Procopia.
15
CordoP. J.GurfinkelV. S. (2004). “Motor coordination can be fully understood only by studying complex movements,” in Progress in brain research (Amsterdam, Netherlands: Elsevier), 29–38.
16
CostaL. d. C. M.KoesB. W.PranskyG.BorkanJ.MaherC. G.SmeetsR. J. (2013). Primary care research priorities in low back pain: An update. Spine38 (2), 148–156. 10.1097/brs.0b013e318267a92f
17
DagenaisS.CaroJ.HaldemanS. (2008). A systematic review of low back pain cost of illness studies in the United States and internationally. spine J.8 (1), 8–20. 10.1016/j.spinee.2007.10.005
18
DavoudiM.FarsaniM. K.BabaeeT.RanjbarH.ShokouhyanS. M.GhaznaviA.et al (2022a). The effects of combining high-top shoes with twister wrap orthoses on balance parameters of children with spastic diplegic cerebral palsy. J. Biomed. Phys. Eng.12 (1), 91–100. 10.31661/jbpe.v0i0.2106-1358
19
DavoudiM.ShokouhyanS. M.AbediM.MeftahiN.RahimiA.RashediE.et al (2020). A practical sensor-based methodology for the quantitative assessment and classification of chronic non specific low back patients (NSLBP) in clinical settings. Sensors20 (10), 2902. 10.3390/s20102902
20
DavoudiM.ShokouhyanS. M.AbediM.RezaeiM.ParnianpourM.KhalafK.et al (2022b). Falls incidence and associated risk factors among people with chronic obstructive pulmonary disease (COPD).
21
DehmiyaniA.MehdizadehH.AzadA.CheraghifardM.JamaliS.DavoudiM.et al (2022). Apathy exacerbates postural control impairments in stroke survivors: The potential effects of cognitive dual-task for improving postural control. Neuropsychologia174, 108344. 10.1016/j.neuropsychologia.2022.108344
22
della VolpeR.PopaT.GinanneschiF.SpidalieriR.MazzocchioR.RossiA. (2006). Changes in coordination of postural control during dynamic stance in chronic low back pain patients. Gait posture24 (3), 349–355. 10.1016/j.gaitpost.2005.10.009
23
FairbankJ. C.PynsentP. B. (2000). The Oswestry disability index. Spine25 (22), 2940–2953. 10.1097/00007632-200011150-00017
24
FerreiraP. H.FerreiraM. L.HodgesP. W. (2004). Changes in recruitment of the abdominal muscles in people with low back pain: Ultrasound measurement of muscle activity. Spine29 (22), 2560–2566. 10.1097/01.brs.0000144410.89182.f9
25
FieldA. (2013). Discovering statistics using IBM SPSS statistics. California, United States: sage.
26
FilizM.CakmakA.OzcanE. (2005). The effectiveness of exercise programmes after lumbar disc surgery: A randomized controlled study. Clin. Rehabil.19 (1), 4–11. 10.1191/0269215505cr836oa
27
FlorH.BraunC.ElbertT.BirbaumerN. (1997). Extensive reorganization of primary somatosensory cortex in chronic back pain patients. Neurosci. Lett.224 (1), 5–8. 10.1016/s0304-3940(97)13441-3
28
GandeviaS.McCloskeyD.BurkeD. (1992). Kinaesthetic signals and muscle contraction. Trends Neurosci.15 (2), 62–65. 10.1016/0166-2236(92)90028-7
29
GhomashchiH.EstekiA.NasrabadiA. M.SprottJ. C.BahrPeymaF. (2011). Dynamic patterns of postural fluctuations during quiet standing: A recurrence quantification approach. Int. J. Bifurc. Chaos21 (04), 1163–1172. 10.1142/s021812741102891x
30
GoodwinG. M.McCloskeyD. I.MatthewsP. B. (1972). Proprioceptive illusions induced by muscle vibration: Contribution by muscle spindles to perception?Science175 (4028), 1382–1384. 10.1126/science.175.4028.1382
31
HaghighatF.EbrahimiS.RezaieM.ShafieeE.ShokouhyanS. M.MoteallehA.et al (2021a). Trunk, pelvis, and knee kinematics during running in females with and without patellofemoral pain. Gait Posture89, 80–85. 10.1016/j.gaitpost.2021.06.023
32
HaghighatF.RezaieM.EbrahimiS.ShokouhyanS. M.MoteallehA.ParnianpourM. (2021b). Coordination variability during walking and running in individuals with and without patellofemoral pain Part 1: Lower limb intersegmental coordination variability. J. Med. Biol. Eng.41 (3), 295–304. 10.1007/s40846-021-00603-5
33
HaghighatF.RezaieM.EbrahimiS.ShokouhyanS. M.MoteallehA.ParnianpourM. (2021c). Coordination variability during walking and running in individuals with and without patellofemoral pain Part 2: Proximal segments coordination variability. J. Med. Biol. Eng.41 (3), 305–313. 10.1007/s40846-021-00598-z
34
HaghighatF.RezaieM.EbrahimiS.ShokuhianS. M.MoteallehA.SalehiR.et al (2021d). The correlation between intersegmental coordination variability and frontal plane hip kinematics during running in persons with patellofemoral pain. J. Biomed. Phys. Eng.
35
HancockM. J.MaherC. G.LaslettM.HayE.KoesB. (2011). Discussion paper: What happened to the ‘bio’in the bio-psycho-social model of low back pain?Eur. Spine J.20 (12), 2105–2110. 10.1007/s00586-011-1886-3
36
HashemiL.WebsterB. S.ClancyE. A. (1998). Trends in disability duration and cost of workers' compensation low back pain claims (1988-1996). J. Occup. Environ. Med.40 (12), 1110–1119. 10.1097/00043764-199812000-00011
37
HashemiL.WebsterB. S.ClancyE. A.VolinnE. (1997). Length of disability and cost of workers' compensation low back pain claims. J. Occup. Environ. Med.39 (10), 937–945. 10.1097/00043764-199710000-00005
38
HodgesP.van den HoornW.DawsonA.CholewickiJ. (2009). Changes in the mechanical properties of the trunk in low back pain may be associated with recurrence. J. biomechanics42 (1), 61–66. 10.1016/j.jbiomech.2008.10.001
39
HorakF. B.NashnerL. M. (1986). Central programming of postural movements: Adaptation to altered support-surface configurations. J. neurophysiology55 (6), 1369–1381. 10.1152/jn.1986.55.6.1369
40
HorakF. B. (2006). Postural orientation and equilibrium: What do we need to know about neural control of balance to prevent falls?Age ageing35 (2), ii7–ii11. 10.1093/ageing/afl077
41
HorakJ.LadislavK.RaidlA. (2003). Deterministic chaos and its physical application. San Francisco: Academia.
42
ItoT.SakaiY.MoritaY.YamazakiK.IgarashiK.NishioR.et al (2018). Proprioceptive weighting ratio for balance control in static standing is reduced in elderly patients with non–specific low back pain. Spine43 (24), 1704–1709. 10.1097/brs.0000000000002817
43
JamaliS.AzadA.MehdizadehH.DoostdarA.HoseinpourF.MehdizadehM.et al (2019). Time-course investigation of postural sway variability: Does anxiety exacerbate the sensory reweighting impairment in chronic stroke survivors?Neuropsychologia127, 185–194. 10.1016/j.neuropsychologia.2019.01.023
44
JoosE.PeretzA.BeguinS.FamaeyJ. P. (1991). Reliability and reproducibility of visual analogue scale and numeric rating scale for therapeutic evaluation of pain in rheumatic patients. J. Rheumatol.18 (8), 1269–1270.
45
KiemelT.ElahiA. J.JekaJ. J. (2008). Identification of the plant for upright stance in humans: Multiple movement patterns from a single neural strategy. J. neurophysiology100 (6), 3394–3406. 10.1152/jn.01272.2007
46
KiersH.BrumagneS.Van DieenJ.van der WeesP.VanheesL. (2012). Ankle proprioception is not targeted by exercises on an unstable surface. Eur. J. Appl. Physiol.112 (4), 1577–1585. 10.1007/s00421-011-2124-8
47
KochC.Garcia-AugundezA.GöbelS.HänselF. (2020). A case control study to investigate differences in motor control between individuals with and without non-specific low back pain during standing. Plos one15 (7), e0234858. 10.1371/journal.pone.0234858
48
LacknerJ. R.DiZioP. (2005). Vestibular, proprioceptive, and haptic contributions to spatial orientation. Annu. Rev. Psychol.56, 115–147. 10.1146/annurev.psych.55.090902.142023
49
LairdR. A.GilbertJ.KentP.KeatingJ. L. (2014). Comparing lumbo-pelvic kinematics in people with and without back pain: A systematic review and meta-analysis. BMC Musculoskelet. Disord.15 (1), 229. 10.1186/1471-2474-15-229
50
MarrasW. S.FergusonS. A.KimJ.-Y.CrowellR. R.BoseS.SimonS. R.et al (1995). The classification of anatomic-and symptom-based low back disorders using motion measure models. Spine20 (23), 2531–2546. 10.1097/00007632-199512000-00013
51
MarrasW. S.LavenderS. A.LeurgansS. E.RajuluS. L.AllreadS. W. G.FathallahF. A.et al (1993). The role of dynamic three-dimensional trunk motion in occupationally-related low back disorders. Spine18 (5), 617–628. 10.1097/00007632-199304000-00015
52
McDowellC.SmykM.SungP. S. (2018). Compensatory strategy between trunk-hip kinematics and reaction time following slip perturbation between subjects with and without chronic low back pain. J. Electromyogr. Kinesiol.43, 68–74. 10.1016/j.jelekin.2018.09.005
53
MitchellT.O'SullivanP. B.BurnettA.StrakerL.SmithA.ThorntonJ.et al (2010). Identification of modifiable personal factors that predict new-onset low back pain: A prospective study of female nursing students. Clin. J. pain26 (4), 275–283. 10.1097/ajp.0b013e3181cd16e1
54
MokN. W.BrauerS. G.HodgesP. W. (2011). Changes in lumbar movement in people with low back pain are related to compromised balance. Spine36 (1), E45–E52. 10.1097/BRS.0b013e3181dfce83
55
MorassoP. G.SchieppatiM. (1999). Can muscle stiffness alone stabilize upright standing?J. neurophysiology82 (3), 1622–1626. 10.1152/jn.1999.82.3.1622
56
NewcomerK. L.LaskowskiE. R.YuB.JohnsonJ. C.AnK.-N. (2000). Differences in repositioning error among patients with low back pain compared with control subjects. Spine25 (19), 2488–2493. 10.1097/00007632-200010010-00011
57
NogueiraJ. F.CarrascoA. C.PelegrinelliA. R. M.GuenkaL. C.SilvaM. F.BelaL. F. D.et al (2020). Posturography comparison and discriminant analysis between individuals with and without chronic low back pain. J. Manip. Physiological Ther.43 (5), 469–475. 10.1016/j.jmpt.2019.01.002
58
PeterkaR. (2002). Sensorimotor integration in human postural control. J. neurophysiology88 (3), 1097–1118. 10.1152/jn.2002.88.3.1097
59
PincusT.BurtonA. K.VogelS.FieldA. P. (2002). A systematic review of psychological factors as predictors of chronicity/disability in prospective cohorts of low back pain. Spine27 (5), E109–E120. 10.1097/00007632-200203010-00017
60
RileyM.BalasubramaniamR.TurveyM. (1999). Recurrence quantification analysis of postural fluctuations. Gait posture9 (1), 65–78. 10.1016/s0966-6362(98)00044-7
61
RollJ.VedelJ. (1982). Kinaesthetic role of muscle afferents in man, studied by tendon vibration and microneurography. Exp. Brain Res.47 (2), 177–190. 10.1007/BF00239377
62
RosensteinM. T.CollinsJ. J.De LucaC. J. (1993). A practical method for calculating largest Lyapunov exponents from small data sets. Phys. D. Nonlinear Phenom.65 (1-2), 117–134. 10.1016/0167-2789(93)90009-p
63
RuheA.FejerR.WalkerB. (2011). Center of pressure excursion as a measure of balance performance in patients with non-specific low back pain compared to healthy controls: A systematic review of the literature. Eur. Spine J.20 (3), 358–368. 10.1007/s00586-010-1543-2
64
RungeC.ShupertC.HorakF.ZajacF. (1999). Ankle and hip postural strategies defined by joint torques. Gait posture10 (2), 161–170. 10.1016/s0966-6362(99)00032-6
65
SalavatiM.HadianM. R.MazaheriM.NegahbanH.EbrahimiI.TalebianS.et al (2009). Test–retest reliabty of center of pressure measures of postural stability during quiet standing in a group with musculoskeletal disorders consisting of low back pain, anterior cruciate ligament injury and functional ankle instability. Gait posture29 (3), 460–464. 10.1016/j.gaitpost.2008.11.016
66
SchieppatiM.GiordanoA.NardoneA. (2002). Variability in a dynamic postural task attests ample flexibility in balance control mechanisms. Exp. Brain Res.144 (2), 200–210. 10.1007/s00221-002-1028-6
67
ScholtesS. A.GombattoS. P.Van DillenL. R. (2009). Differences in lumbopelvic motion between people with and people without low back pain during two lower limb movement tests. Clin. Biomech.24 (1), 7–12. 10.1016/j.clinbiomech.2008.09.008
68
SharmaL.PaiY.-C. (1997). Impaired proprioception and osteoarthritis. Curr. Opin. rheumatology9 (3), 253–258. 10.1097/00002281-199705000-00013
69
ShokouhyanS. M.DavoudiM.HoviattalabM.AbediM.BervisS.ParnianpourM.et al (2020). Linear and non-linear dynamic methods toward investigating proprioception impairment in non-specific low back pain patients. Front. Bioeng. Biotechnol.8, 584952
70
SwainC. T.PanF.OwenP. J.SchmidtH.BelavyD. L. (2020). No consensus on causality of spine postures or physical exposure and low back pain: A systematic review of systematic reviews. J. biomechanics102, 109312. 10.1016/j.jbiomech.2019.08.006
71
van der HulstM.Vollenbroek-HuttenM. M.RietmanJ. S.HermensH. J. (2010). Lumbar and abdominal muscle activity during walking in subjects with chronic low back pain: Support of the “guarding” hypothesis?J. Electromyogr. Kinesiol.20 (1), 31–38. 10.1016/j.jelekin.2009.03.009
72
van DieënJ. H.SelenL. P.CholewickiJ. (2003). Trunk muscle activation in low-back pain patients, an analysis of the literature. J. Electromyogr. Kinesiol.13 (4), 333–351. 10.1016/s1050-6411(03)00041-5
73
VogtL.PfeiferK.PortscherM.BanzerW. (2001). Influences of nonspecific low back pain on three-dimensional lumbar spine kinematics in locomotion. Spine26 (17), 1910–1919. 10.1097/00007632-200109010-00019
74
VosT.LimS. S.AbbafatiC.AbbasK. M.AbbasiM.AbbasifardM. (2020). Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet396 (10258), 1204–1222.
75
WebberC.JrMarwanN.FacchiniA.GiulianiA. (2009). Simpler methods do it better: Success of Recurrence Quantification Analysis as a general purpose data analysis tool. Physics Letters A373 (41), 3753–6.
76
WebberC. L.JrZbilutJ. P. (2005). Recurrence quantification analysis of nonlinear dynamical systems. Tutorials Contemp. nonlinear methods Behav. Sci.94, 26–94.
77
WernliK.SmithA.CollF.CampbellA.KentP.O'SullivanP. (2022). From protection to non‐protection: A mixed methods study investigating movement, posture and recovery from disabling low back pain. Eur. J. Pain26, 2097–2119. 10.1002/ejp.2022
Summary
Keywords
non-specific low back pain, linear discriminant analysis, COP, posture control, proprioceptive impairment
Citation
Shokouhyan SM, Davoudi M, Hoviattalab M, Abedi M, Bervis S, Parnianpour M, Brumagne S and Khalaf K (2022) Distinction of non-specific low back pain patients with proprioceptive disorders from healthy individuals by linear discriminant analysis. Front. Bioeng. Biotechnol. 10:1078805. doi: 10.3389/fbioe.2022.1078805
Received
24 October 2022
Accepted
28 November 2022
Published
13 December 2022
Volume
10 - 2022
Edited by
Stefan Schmid, Bern University of Applied Sciences, Switzerland
Reviewed by
Michael Lukas Meier, Balgrist University Hospital, Switzerland
Scott Tagliaferri, Deakin University, Australia
Updates
Copyright
© 2022 Shokouhyan, Davoudi, Hoviattalab, Abedi, Bervis, Parnianpour, Brumagne and Khalaf.
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: Kinda Khalaf, kinda.khalaf@ku.ac.ae
This article was submitted to Biomechanics, a section of the journal Frontiers in Bioengineering and Biotechnology
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.