Spinal Loading Patterns From Biomechanical Modeling Explain the High Incidence of Vertebral Fractures in the Thoracolumbar Region

Assessment of age-dependent sexual dimorphism in paediatric vertebral size and density using a statistical shape and statistical appearance modelling approach

This work focuses on the growth patterns of the human fourth lumbar vertebra (L4) in a paediatric population, with specific attention to sexual dimorphism. The study aims to understand morphological and density changes in the vertebrae through age-dependent statistical shape and statistical appearance models, which can describe full three-dimensional anatomy. Results show that the main growth patterns are associated with isotropic volumetric vertebral growth, a decrease in the relative size of the vertebral foramen, and an increase in the length of the transverse processes. Moreover, significant sexual dimorphism was demonstrated during puberty. We observe significant age and sex interaction in the anterior vertebral body height (P = 0.005), where females exhibited an earlier increase in rates of vertebral height evolution. Moreover, we also observe an increase in cross-sectional area (CSA) with age (P = 0.020), where the CSA is smaller in females than in males (significant sex effect P = 0.042). Finally, although no significant increase in trabecular bone density with age is observed (P = 0.363), a trend in the statistical appearance model suggests an increase in density with age.

Effect of Different Injury Morphology of the Endplate on Intervertebral Disc Degeneration: Retrospective Cohort Study

OBJECTIVES

To describe a simplified classification scheme for endplate injury morphology based on 3D CT and to examine possible associations between endplate injury morphology and vertebral space and other variables such as type of fracture and disc degeneration in a group of patients with thoracolumbar fractures.

METHODS

This study was a retrospective cohort study. We collected patients with thoracolumbar fractures admitted from January 2015 to August 2020 and divided them into three groups based on the morphology of endplate injury (45 cases of mild endplate injury, 54 cases of moderate endplate injury, and 42 cases of severe endplate injury, SEI). Data of vertebral body and intervertebral space height and angle, the Pfirrmann grade, endplate healing morphology were collected during preoperative, postoperative, and long-term follow-up of patients in each group. One-way analysis of variance (ANOVA), chi-squared test, and repeated measurement ANOVA were used to compare and analyze the influence of endplate injury morphology on patient prognosis.

RESULTS

Most moderate injuries to the endplate (fissure-type injury) and severe injuries (irregular depression-type injury, Schmorl's node-type injury) resulted in significant disc degeneration in the long-term transition. This study also showed significant differences in the height of the anterior margin of the injured spine and the intervertebral space height index during this process.

CONCLUSIONS

The current study suggests that although the region of injury in endplate fissure-type injury is small preoperatively, it may be a major factor in leading to severe disc degeneration, loss of intervertebral height, and Cobb angle in the long term. The results of our study therefore may allow surgeons to predict the prognosis of patients with thoracolumbar fractures and guide their treatment.

Epidemiology of Lumbar Spine Fractures: Twenty-Year Assessment of Nationwide Emergency Department Visit Data

BACKGROUND

Lumbar spine fractures are common injuries associated with substantial morbidity for patients and socioeconomic burden. This study sought to epidemiologically analyze lumbar spine fractures by mechanism of injury and identify temporal trends in patient demographics and disposition, which few studies have previously evaluated.

MATERIALS AND METHODS

A retrospective analysis was done of the US National Electronic Injury Surveillance System (NEISS) database between 2003 and 2022. The sample contained all patients 2 to 101 years old with product-related lumbar fractures presenting to participating institutions' emergency departments. A total of 15,196 unweighted injuries (642,979 weighted injuries) were recorded.

RESULTS

Overall, there was a 20-year incidence rate of 10.14 cases per 100,000 person-years with a 2-fold increase in fracture incidence. Females were more prone to lumbar fracture than males (P=.032). Injuries primarily stemmed from a fall (76.6%). The incidence of lumbar fracture increased most significantly in older patients, with patients 80 years and older showing the greatest annual increase (β=8.771, R2=0.7439, P<.001) and patients 60 to 69 years showing the greatest percent increase with a 3.24-fold increase in incidence. Most (58.9%) of the fractures occurred at home. Females were more often injured at home compared with males (P<.001), who more often sustained lumbar fractures during recreational or athletic activity (P<.001). All patients older than 40 years showed at least a doubling in incidence rate of lumbar fracture between 2003 and 2022.

CONCLUSION

These data demonstrate the pressing need to address poor bone health in the aging population, shown here to have an increasing fracture burden. [Orthopedics. 202x;4x(x):xx-xx.].

Musculoskeletal models determine the effect of a soft active exosuit on muscle activations and forces during lifting and lowering tasks

Exosuits have the potential to mitigate musculoskeletal stress and prevent back injuries during industrial tasks. This study aimed to 1) validate the implementation of a soft active exosuit into a musculoskeletal model of the spine by comparing model predicted muscle activations versus corresponding surface EMG measurements, and 2) evaluate the effect of the exosuit on peak back and hip muscle forces. Fourteen healthy participants performed squat and stoop lift and lower tasks with boxes of 6 and 10 kg, with and without wearing a 2.7 kg soft active exosuit. Participant-specific musculoskeletal models, which included the exosuit, were created in OpenSim. Model validation focused on the back and hip extensors, where temporal agreement between EMG and model estimated muscle activity was generally strong to excellent (average cross-correlation coefficients ranging from 0.84 to 0.98). Root mean square errors of muscle activity (0.05-0.10) were similar with and without the exosuit, and compared well to prior model validation studies without the exosuit (average root mean square errors ranging from 0.05 to 0.19). In terms of performance, the exosuit reduced the estimated peak erector spinae forces during lifting and lowering phases across all lifting tasks but reduced peak hip extensor muscles forces only in a squat lift task of 10 kg. These reductions in total peak muscle forces were approximately 1.7-4.2 times greater than the corresponding exosuit assistance force, which were 146 ± 19 N and 102 ± 14 N at the times of peak erector spinae forces in lifting and lowering, respectively. Overall, the results support the hypothesis that exosuits reduce soft tissue loading, and thereby potentially reduce fatigue and injury risk during manual materials handling tasks. Incorporating exosuits into musculoskeletal models is a valid approach to understand the impact of exosuit assistance on muscle activity and forces.

Optical Marker-Based Motion Capture of the Human Spine: A Scoping Review of Study Design and Outcomes

Biomechanical analysis of the human spine is crucial to understanding injury patterns. Motion capture technology has gained attention due to its non-invasive nature. Nevertheless, traditional motion capture studies consider the spine a single rigid segment, although its alignment changes during movement. Moreover, guidelines that indicate where markers should be placed for a specific exercise do not exist. This study aims to review the methods used to assess spine biomechanics using motion capture systems to determine the marker sets used, the protocols used, the resulting parameters, the analysed activities, and the characteristics of the studied populations. PRISMA guidelines were used to perform a Scoping Review using SCOPUS and Web of Science databases. Fifty-six journal and conference articles from 1997 to 2023 were considered for the analysis. This review showed that Plug-in-Gait is the most used marker set. The lumbar spine is the segment that generates the most interest because of its high mobility and function as a weight supporter. Furthermore, angular position and velocity are the most common outcomes when studying the spine. Walking, standing, and range of movement were the most studied activities compared to sports and work-related activities. Male and female participants were recruited similarly across all included articles. This review presents the motion capture techniques and measurement outcomes of biomechanical studies of the human spine, to help standardize the field. This work also discusses trends in marker sets, study outcomes, studied segments and segmentation approaches.

Modelling the female torso and breast during physical activity: Implications on spinal loading

Methods of modelling the female torso during physical activity often neglect the position and movement of the breast. This novel investigation compares three female torso modelling approaches that differ in complexity (integrated breast, fixed breast, dynamic breast) to determine the effect on spinal joint moments during running and jumping. The commonly used integrated breast model distributed breast mass within the torso, the fixed breast model attached the mass of the breasts to fixed positions on the anterior of the torso, and a new dynamic breast model enabled relative motion between the breasts and anterior torso. Key findings demonstrated minimal differences in lumbar spine moments (<0.05 Nm/kg; 4%) between integrated breast and fixed breast models but greater differences, up to 0.86 Nm/kg (68%) during running and 0.89 Nm/kg (82%) during jumping, when breast motion was included. Thoracic spine moments revealed similar patterns with minimal differences (<0.05 Nm/kg; 11%) between integrated breast and fixed breast models and greater differences, up to 0.48 Nm/kg (92%) during running and 0.63 Nm/kg (66%) during jumping, with the dynamic breast model. Future female musculoskeletal models should consider including breast mass and motion to avoid mis-representing spinal loading in females during running and jumping.

Identification of a lumped-parameter model of the intervertebral joint from experimental data

Through predictive simulations, multibody models can aid the treatment of spinal pathologies by identifying optimal surgical procedures. Critical to achieving accurate predictions is the definition of the intervertebral joint. The joint pose is often defined by virtual palpation. Intervertebral joint stiffnesses are either derived from literature, or specimen-specific stiffnesses are calculated with optimisation methods. This study tested the feasibility of an optimisation method for determining the specimen-specific stiffnesses and investigated the influence of the assigned joint pose on the subject-specific estimated stiffness. Furthermore, the influence of the joint pose and the stiffness on the accuracy of the predicted motion was investigated. A computed tomography based model of a lumbar spine segment was created. Joints were defined from virtually palpated landmarks sampled with a Latin Hypercube technique from a possible Cartesian space. An optimisation method was used to determine specimen-specific stiffnesses for 500 models. A two-factor analysis was performed by running forward dynamic simulations for ten different stiffnesses for each successfully optimised model. The optimisations calculated a large range of stiffnesses, indicating the optimised specimen-specific stiffnesses were highly sensitive to the assigned joint pose and related uncertainties. A limited number of combinations of optimised joint stiffnesses and joint poses could accurately predict the kinematics. The two-factor analysis indicated that, for the ranges explored, the joint pose definition was more important than the stiffness. To obtain kinematic prediction errors below 1 mm and 1° and suitable specimen-specific stiffnesses the precision of virtually palpated landmarks for joint definition should be better than 2.9 mm.

Differences in vertebral bone density between African apes

OBJECTIVES

Low-energy vertebral fractures are a common health concern, especially in elderly people. Interestingly, African apes do not seem to experience as many vertebral fractures and the low-energy ones are even rarer. One potential explanation for this difference is the lower bone density in humans. Yet, only limited research has been done on the vertebral bone density of the great apes and these have mainly included only single vertebrae. Hence the study aim is to expand our understanding of the vertebral microstructure of African apes in multiple spinal segments.

MATERIALS

Bone density in the vertebral body of C7, T12, and L3 was measured from 32 Pan troglodytes and 26 Gorilla gorilla using peripheral quantitative computed tomography (pQCT).

RESULTS

There was a clear difference between the three individual vertebrae and consequently the spinal segments in terms of trabecular density and cortical density and thickness. The variation of these bone parameters between the vertebrae differed between the apes but was also different from those reported for humans. The chimpanzees were observed to have overall higher trabecular density, but gorillas had higher cortical density and thickness. Cortical thickness had a relatively strong association with the vertebral size.

DISCUSSION

Despite the similarity in locomotion and posture, the results show slight differences in the bone parameters and their variation between spinal segments in African apes. This variation also differs from humans and appears to indicate a complex influence of locomotion, posture, and body size on the different spinal segments.

Spinal joint moment prediction following simulated breast surgery using a female whole-body musculoskeletal model

This study aimed to use a musculoskeletal model to predict changes in spinal moments following simulated breast surgery. A female full body musculoskeletal model with a fully articulated thoracolumbar spine and independent moveable breast segments was customised for this study. Key findings suggest that the simulated removal of breast tissue (750 g to 1501 g) can reduce the magnitude of lumbar spine extensor moments by >0.05 Nm/kg during walking and jogging. A customised female whole-body musculoskeletal model is capable of providing a first approximation of changes in spinal loading following simulated breast surgery.

Musculoskeletal spine modeling in large patient cohorts: how morphological individualization affects lumbar load estimation

Introduction: Achieving an adequate level of detail is a crucial part of any modeling process. Thus, oversimplification of complex systems can lead to overestimation, underestimation, and general bias of effects, while elaborate models run the risk of losing validity due to the uncontrolled interaction of multiple influencing factors and error propagation. Methods: We used a validated pipeline for the automated generation of multi-body models of the trunk to create 279 models based on CT data from 93 patients to investigate how different degrees of individualization affect the observed effects of different morphological characteristics on lumbar loads. Specifically, individual parameters related to spinal morphology (thoracic kyphosis (TK), lumbar lordosis (LL), and torso height (TH)), as well as torso weight (TW) and distribution, were fully or partly considered in the respective models according to their degree of individualization, and the effect strengths of these parameters on spinal loading were compared between semi- and highly individualized models. T-distributed stochastic neighbor embedding (T-SNE) analysis was performed for overarching pattern recognition and multiple regression analyses to evaluate changes in occurring effects and significance. Results: We were able to identify significant effects (p < 0.05) of various morphological parameters on lumbar loads in models with different degrees of individualization. Torso weight and lumbar lordosis showed the strongest effects on compression (β ≈ 0.9) and anterior-posterior shear forces (β ≈ 0.7), respectively. We could further show that the effect strength of individual parameters tended to decrease if more individual characteristics were included in the models. Discussion: The induced variability due to model individualization could only partly be explained by simple morphological parameters. Our study shows that model simplification can lead to an emphasis on individual effects, which needs to be critically assessed with regard to in vivo complexity. At the same time, we demonstrated that individualized models representing a population-based cohort are still able to identify relevant influences on spinal loading while considering a variety of influencing factors and their interactions.

Efficacy and Safety of Conservative Treatment Compared With Surgical Treatment for Thoracolumbar Fracture With Score 4 Thoracolumbar Injury Classification and Severity (TLICS): A Systematic Review and Meta-analysis

STUDY DESIGN

This was a systematic review and meta-analysis.

OBJECTIVE

The clinical outcomes, radiologic outcome, and complications were compared between surgical treatment and conservative treatment of thoracolumbar fractures with a Thoracolumbar Injury Classification and Severity (TLICS) score of 4.

SUMMARY OF BACKGROUND DATA

The thoracolumbar fracture is the main reason leading to the spinal cord injury. Some studies suggested that the treatment of TLICS=4 is a "gray zone." Hence, the efficacy and safety of surgical treatment and conservative treatment of thoracolumbar fractures with scores 4 TLICS was still debated.

MATERIALS AND METHODS

A comprehensive search of PubMed, Embase, and the Cochrane Library, Chinese National Knowledge Infrastructure (CNKI), Chongqing VIP Database (VIP), and Wan Fang Database was performed up to October 2021. Relevant studies were identified using specific eligibility criteria and data was extracted and analyzed based on primary and secondary outcomes.

RESULTS

A total of 10 studies involving 555 patients were included (3 randomized controlled trials and 7 retrospective studies). There was no significant difference of hospital time (standardized mean difference=0.24, 95% CI: -1.50 to 1.97, P =0.79) and Oswestry Disability Index (mean difference=2.97, 95% CI: -1.07 to 7.01, P =0.15) between surgery and nonsurgery. The length of returning to work was shorter in surgical treatment (standardized mean difference=1.27, 95% CI: 0.07-2.46, P =0.04). Visual Analog Scale in surgical treatment was lower at 1, 3, and 6 months (respectively, P <0.00001, P =0.003, and P =0.02). However, there existed no significant difference between surgical treatment and nonsurgical treatment at 12 and >24 months (respectively, P =0.18 and 0.17). Cobb angle was lower in surgical treatment at postoperative at 6, 12, and >24 months (respectively, P =0.005, P <0.00001, P =0.002, and P =0.0002). Finally, the surgical treatment had a lower incidence of complications (odds ratio=3.89, 95% CI: 1.90-7.94, P =0.0002).

CONCLUSIONS

Current evidence recommended that surgical treatment is superior to conservative treatment of TLICS score of 4 at the early follow-up. Surgical treatment had lower Cobb angle, Visual Analog Scale scores, and complications compared with a nonsurgical TLICS score of 4. However, these findings needed to be verified further by multicenter, double-blind, and large-sample randomized controlled trials.

Biomechanical analysis of complications following T10-Pelvis spinal fusion: A population based computational study

Proximal junctional kyphosis (PJK) and proximal junctional failure (PJF) are challenging complications of long fusion constructs for the treatment of adult spinal deformity. The objective of this study is to understand the biomechanical stresses proximal to the upper instrumentation of a T10-pelvis fusion in a large patient cohort. The pre-fusion models were subject-specific thoracolumbar spine models that incorporate the height, weight, spine curvature, and muscle morphology of 250 individuals from the Framingham Heart Study Multidetector CT Study. To create post-fusion models, the subject-specific models were further modified to eliminate motion between the intervertebral joints from T10 to the pelvis. OpenSim analysis tools were used to calculate the medial lateral shear force, anterior posterior shear force, and compressive force on the T9 vertebra during the static postures. Differences between pre-fusion and post-fusion T9 biomechanics were consistent between increased segmental mobility and unchanged segmental mobility conditions. For all static postures, compression decreased (p < 0. 0005). Anterior-posterior shear force significantly increased (p < 0. 0005) during axial twist and significantly increased (p < 0. 0005) during trunk flexion. Medial lateral shear force significantly increased (p < 0. 0005) during axial twist. This computational study provided the first use of subject-specific models to investigate the biomechanics of long spinal fusions. Patients undergoing T10-Pelvis fusion were predicted to have increased shear forces and decreased compressive force at the T9 vertebra, independent of change in segmental mobility. The computational model shows potential for the investigation of spinal fusion biomechanics to reduce the risk of PJK or PJF.

The effect of a soft active back support exosuit on trunk motion and thoracolumbar spine loading during squat and stoop lifts

Back support exosuits aim to reduce tissue demands and thereby risk of injury and pain. However, biomechanical analyses of soft active exosuit designs have been limited. The objective of this study was to evaluate the effect of a soft active back support exosuit on trunk motion and thoracolumbar spine loading in participants performing stoop and squat lifts of 6 and 10 kg crates, using participant-specific musculoskeletal models. The exosuit did not change overall trunk motion but affected lumbo-pelvic motion slightly, and reduced peak compressive and shear vertebral loads at some levels, although shear increased slightly at others. This study indicates that soft active exosuits have limited kinematic effects during lifting, and can reduce spinal loading depending on the vertebral level. These results support the hypothesis that a soft exosuit can assist without limiting trunk movement or negatively impacting skeletal loading and have implications for future design and ergonomic intervention efforts.

Percutaneous monoplanar screws versus hybrid fixed axial and polyaxial screws in intermediate screw fixation for traumatic thoracolumbar burst fractures: a case-control study

BACKGROUND

To compare the clinical and radiological outcomes of monoplanar screws (MSs) versus hybrid fixed axial and polyaxial screws (HSs) in percutaneous short-segment intermediate screw fixation (PSISF) for traumatic thoracolumbar burst fractures (TTBFs) in patients without neurologic impairment.

METHODS

A consecutive series of 100 patients with single-segment TTBFs and no neurologic impairment who underwent PSISF with 6 monoplanar screws (MS group) or correct were retrospectively enrolled. The demographic data, radiologic evaluation indicators, perioperative indicators and clinical assessment indicators were analysed between the MS group and HS group.

RESULTS

The demographic data and perioperative indicators were not significantly different in the two groups (P > 0.05). The postoperative anterior vertebral height ratio (AVHR), kyphosis Cobb angle (KCA), vertebral wedge angle (VWA) and spinal canal encroachment rate (SCER) were significantly improved in both groups (*P < 0.05). The MS group obtained better correction than the HS group in terms of improvement in the AVHR, KCA and VWA after surgery (*P < 0.05). At the last follow-up, the MS group had less correction loss of AVHR, KCA and VWA (*P < 0.05). The MS group presented greater improvement in the SCER at the last follow-up (*P < 0.05). The visual analogue scale (VAS) score and Oswestry Disability Index (ODI) score of all patients were significantly better postoperatively than those preoperatively (*P < 0.05), and the scores collected at each follow-up visit did not differ significantly between the two groups (P > 0.05). In the MS group, no internal fixation failure was observed during the follow-up period, but, in the HS group, two cases of internal fixation failure were observed at the last follow-up (one case of rod loosening and one case of screw breakage).

CONCLUSIONS

Both MSs and HSs fixation are effective treatments for TTBFs and have comparable clinical outcomes. In contrast, MSs fixation can improve the correction effect, better improve the SCER, and further reduce correction loss as well as reduce the incidence of instrumentation failure. Therefore, MSs fixation might be a better option for treating TTBFs in patients without neurological deficits.

Using inertial measurement units to estimate spine joint kinematics and kinetics during walking and running

Optical motion capture (OMC) is considered the best available method for measuring spine kinematics, yet inertial measurement units (IMU) have the potential to collect data outside the laboratory. When combined with musculoskeletal modeling, IMU technology may be used to estimate spinal loads in real-world settings. To date, IMUs have not been validated for estimates of spinal movement and loading during both walking and running. Using OpenSim Thoracolumbar Spine and Ribcage models, we compare IMU and OMC estimates of lumbosacral (L5/S1) and thoracolumbar (T12/L1) joint angles, moments, and reaction forces during gait across six speeds for five participants. For comparisons, time series are ensemble averaged over strides. Comparisons between IMU and OMC ensemble averages have low normalized root mean squared errors (< 0.3 for 81% of comparisons) and high, positive cross-correlations (> 0.5 for 91% of comparisons), suggesting signals are similar in magnitude and trend. As expected, joint moments and reaction forces are higher during running than walking for IMU and OMC. Relative to OMC, IMU overestimates joint moments and underestimates joint reaction forces by 20.9% and 15.7%, respectively. The results suggest using a combination of IMU technology and musculoskeletal modeling is a valid means for estimating spinal movement and loading.

Characterizing Bone Phenotypes Related to Skeletal Fragility Using Advanced Medical Imaging

PURPOSE OF REVIEW

Summarize the recent literature that investigates how advanced medical imaging has contributed to our understanding of skeletal phenotypes and fracture risk across the lifespan.

RECENT FINDINGS

Characterization of bone phenotypes on the macro-scale using advanced imaging has shown that while wide bones are generally stronger than narrow bones, they may be more susceptible to age-related declines in bone strength. On the micro-scale, HR-pQCT has been used to identify bone microarchitecture phenotypes that improve stratification of fracture risk based on phenotype-specific risk factors. Adolescence is a key phase for bone development, with distinct sex-specific growth patterns and significant within-sex bone property variability. However, longitudinal studies are needed to evaluate how early skeletal growth impacts adult bone phenotypes and fracture risk. Metabolic and rare bone diseases amplify fracture risk, but the interplay between bone phenotypes and disease remains unclear. Although bone phenotyping is a promising approach to improve fracture risk assessment, the clinical availability of advanced imaging is still limited. Consequently, alternative strategies for assessing and managing fracture risk include vertebral fracture assessment from clinically available medical imaging modalities/techniques or from fracture risk assessment tools based on clinical risk factors. Bone fragility is not solely determined by its density but by a combination of bone geometry, distribution of bone mass, microarchitecture, and the intrinsic material properties of bone tissue. As such, different individuals can exhibit distinct bone phenotypes, which may predispose them to be more vulnerable or resilient to certain perturbations that influence bone strength.

Biomimetic Inspired Hydrogels for Regenerative Vertebral Body Stenting

PURPOSE OF REVIEW

This review aims to explore the potential of biomimetic hydrogels as an alternative to bone cement in vertebral body stenting (VBS), a minimally invasive treatment for vertebral compression fractures.

RECENT FINDINGS

The use of bone cement in VBS procedures can lead to complications such as incomplete fracture reduction and cement leakage. Biomimetic hydrogels have gained significant attention as potential biomaterial alternatives for VBS due to their unique properties, including tuneable therapeutic and mechanical properties. Over the past decade, there has been significant advancements in the development of biomimetic hydrogels for bone regeneration, employing a wide range of approaches to enhance the structural and functional properties of hydrogels. Biomimetic hydrogels hold significant promise as safer and reparative alternatives to bone cement for VBS procedures. However, further research and development in this field are necessary to explore the full potential of hydrogel-based systems for vertebral bone repair.

Advances in the Development of Nano-Engineered Mechanically Robust Hydrogels for Minimally Invasive Treatment of Bone Defects

Injectable hydrogels were discovered as attractive materials for bone tissue engineering applications given their outstanding biocompatibility, high water content, and versatile fabrication platforms into materials with different physiochemical properties. However, traditional hydrogels suffer from weak mechanical strength, limiting their use in heavy load-bearing areas. Thus, the fabrication of mechanically robust injectable hydrogels that are suitable for load-bearing environments is of great interest. Successful material design for bone tissue engineering requires an understanding of the composition and structure of the material chosen, as well as the appropriate selection of biomimetic natural or synthetic materials. This review focuses on recent advancements in materials-design considerations and approaches to prepare mechanically robust injectable hydrogels for bone tissue engineering applications. We outline the materials-design approaches through a selection of materials and fabrication methods. Finally, we discuss unmet needs and current challenges in the development of ideal materials for bone tissue regeneration and highlight emerging strategies in the field.

EMG Validation of a Subject-Specific Thoracolumbar Spine Musculoskeletal Model During Dynamic Activities in Older Adults

Musculoskeletal models can uniquely estimate in vivo demands and injury risk. In this study, we aimed to compare muscle activations from subject-specific thoracolumbar spine OpenSim models with recorded muscle activity from electromyography (EMG) during five dynamic tasks. Specifically, 11 older adults (mean = 65 years, SD = 9) lifted a crate weighted to 10% of their body mass in axial rotation, 2-handed sagittal lift, 1-handed sagittal lift, and lateral bending, and simulated a window opening task. EMG measurements of back and abdominal muscles were directly compared to equivalent model-predicted activity for temporal similarity via maximum absolute normalized cross-correlation (MANCC) coefficients and for magnitude differences via root-mean-square errors (RMSE), across all combinations of participants, dynamic tasks, and muscle groups. We found that across most of the tasks the model reasonably predicted temporal behavior of back extensor muscles (median MANCC = 0.92 ± 0.07) but moderate temporal similarity was observed for abdominal muscles (median MANCC = 0.60 ± 0.20). Activation magnitude was comparable to previous modeling studies, and median RMSE was 0.18 ± 0.08 for back extensor muscles. Overall, these results indicate that our thoracolumbar spine model can be used to estimate subject-specific in vivo muscular activations for these dynamic lifting tasks.

Comparison of thoracolumbar versus non-thoracolumbar osteoporotic vertebral compression fractures in risk factors, vertebral compression degree and pre-hospital back pain

BACKGROUND

Thoracolumbar spine is at high risk of osteoporotic vertebral compression fractures (OVCF). This study aimed to identify the differences in risk factors, vertebral compression degree and back pain characteristics of thoracolumbar OVCF (TL-OVCF) and non-thoracolumbar OVCF (nTL-OVCF).

METHODS

OVCF patients hospitalized in a spine center between June 2016 and October 2020 were retrospectively studied. Demographics, comorbidity, spine trauma, bone mineral density, duration of pre-hospital back pain, extent of vertebral marrow edema, and degree of vertebral compression of patients with nTL-OVCF were summarized and compared to those with TL-OVCF.

RESULTS

A total of 944 patients with acute single-segment OVCF were included. There were 708 (75.0%) TL-OVCF located in T11-L2 and 236 (25.0%) nTL-OVCF in lower lumbar (L3-L5) and middle thoracic (T5-T10) spine. The female-male ratio was 4.1 in nTL-OVCF and differed not significantly from TL-OVCF. The middle thoracic OVCF were older and had higher comorbidity of coronary heart disease (21.3%) and cerebral infarction (36.3%) than TL-OVCF (12.1% and 20.6%). In nTL-OVCF the ratio of apparent spine trauma (44.9%) and pre-hospital back pain ≤ 1 week (47.5%) was lower than in TL-OVCF (66.9% and 62.6%). The T-score value of lumbar spine was - 2.99 ± 1.11, - 3.24 ± 1.14, - 3.05 ± 1.40 in < 70, 70-80, > 80 years old TL-OVCF and differed not significantly from nTL-OVCF. The lower lumbar OVCF had more cranial type of vertebral marrow edema (21.8%) and fewer concurrent lumbodorsal fasciitis (30.8%) than TL-OVCF (16.8% and 43.4%). In TL-OVCF the anterior-posterior vertebral height ratio was lower with back pain for > 4 weeks than for ≤ 1, 1-2, and 2-4 weeks. In nTL-OVCF the degree of vertebral compression differed not significantly with pre-hospital back pain for ≤ 1, 1-2, 2-4, and > 4 weeks.

CONCLUSIONS

Thoracolumbar spine has 2-folds higher risk of OVCF than non-thoracolumbar spine. Non-thoracolumbar OVCF are not associated with female gender, apparent spine trauma or poor bone mineral density, but tend to maintain the degree of vertebral compression and cause longer duration of pre-hospital back pain.

Smartphone use affects gait performance, spinal kinematics and causes spinal musculoskeletal discomfort in young adults

BACKGROUND

Smartphone use may lead to alterations in spinal kinematics and musculoskeletal discomfort.

OBJECTIVES

The aim of this study was to evaluate the effect of smartphone use on spinal kinematics, and to examine the relationship between smartphone addiction, spinal discomfort, and gait parameters.

DESIGN

Cross-Sectional Study.

METHODS

The study included 42 healthy adults aged 18-30 years. A photographic method was used for spinal kinematic evaluation in sitting, standing and at the end of a 3-min walk. GAITRite electronic walkway was used for spatiotemporal gait parameters. Smartphone addiction was evaluated with the Smartphone Addiction Scale - Short Version (SAS-SV). The Cornell Musculoskeletal System Discomfort Questionnaire (CMDQ) was used to evaluate feelings of discomfort and pain.

RESULTS

There was an increase in head, cervical, and thoracic flexion angles while sitting, standing, and at the end of a 3-min walk. Similarly, an increase in thoracolumbar and lumbar flexion angles was observed only in the sitting position (p < 0.05). While using a smartphone during walking, cadence, walking speed, step length decreased, while step duration and double support duration increased (p < 0.05). A statistically significant correlation was determined between the SAS-SV and CMDQ scores (p < 0.05).

CONCLUSION

The study showed that smartphone use has an impact on spinal kinematics during sitting, standing and at the end of a 3 min-walk, as well as on the spatiotemporal parameters of gait. This study suggest that smartphone addiction should be taken into consideration due to its potential to cause musculoskeletal discomfort and there may be a need to raise public awareness on this matter.

Effects of Bone Cross-Link Bridging on Fracture Mechanism and Surgical Outcomes in Elderly Patients with Spine Fractures

STUDY DESIGN

This study adopted a cross-sectional study design.

PURPOSE

This study was designed to investigate the effects of bone cross-link bridging on fracture mechanism and surgical outcomes in vertebral fractures using the maximum number of vertebral bodies with bony bridges between adjacent vertebrae without interruption (maxVB).

OVERVIEW OF LITERATURE

The complex interplay of bone density and bone bridging in the elderly can complicate vertebral fractures, necessitating a better understanding of fracture mechanics.

METHODS

We examined 242 patients (age >60 years) who underwent surgery for thoracic to lumbar spine fractures from 2010 to 2020. Subsequently, the maxVB was classified into three groups: maxVB (0), maxVB (2-8), and maxVB (9-18), and parameters, including fracture morphology (new Association of Osteosynthesis classification), fracture level, and neurological deficits were compared. In a sub-analysis, 146 patients with thoracolumbar spine fractures were classified into the three aforementioned groups based on the maxVB and compared to determine the optimal operative technique and evaluate surgical outcomes.

RESULTS

Regarding the fracture morphology, the maxVB (0) group had more A3 and A4 fractures, whereas the maxVB (2-8) group had less A4 and more B1 and B2 fractures. The maxVB (9-18) group exhibited an increased frequency of B3 and C fractures. Regarding the fracture level, the maxVB (0) group tended to have more fractures in the thoracolumbar transition region. Furthermore, the maxVB (2-8) group had a higher fracture frequency in the lumbar spine area, whereas the maxVB (9-18) group had a higher fracture frequency in the thoracic spine area than the maxVB (0) group. The maxVB (9-18) group had fewer preoperative neurological deficits but a higher reoperation rate and postoperative mortality than the other groups.

CONCLUSIONS

The maxVB was identified as a factor influencing fracture level, fracture type, and preoperative neurological deficits. Thus, understanding the maxVB could help elucidate fracture mechanics and assist in perioperative patient management.

Dependence of trunk muscle size and position on age, height, and weight in a multi-ethnic cohort of middle-aged and older men and women

Trunk muscle size and location relative to the spine are key factors affecting their capacity to assist in trunk movement, strength, and function. There remains limited information on how age, weight and height affect these measurements across multiple spinal levels, and prior studies had limited samples in terms of size and ethnicity. In this study, we measured trunk muscles in coronal plane slices at T4 - L4 of CT scans acquired in 507 participants, aged 40-90 years, from the community-based Framingham Heart Study. Mixed-effects linear regressions, stratified by sex, determined the contributions of age, height and weight, to muscle cross-sectional area (CSA), the distance from the vertebral body centroid (CD), and the in-plane angle of the line between the vertebral body and the muscle centroids (CA). Muscle CSA decreased with higher age by an average of -0.8% per year, but weight (average 0.8% per kg) and height (average -0.05% per cm) had mixed results, with both positive and negative effects depending on muscle group and level. Muscle CD increased with weight by an average of 0.3% per kg, but had mixed effects for age (average 0.8% per year) and height (average 0.1% per cm). Muscle CA had mixed associations with age (average 0.05% per year), weight (average 0.01% per kg) and height (average -0.05% per cm). A prediction program created with these results provides a simple approach for estimating probable values for trunk muscle size and position in the absence of medical imaging.

Comparison of acute single versus multiple osteoporotic vertebral compression fractures in radiographic characteristic and bone fragility

BACKGROUND

Osteoporotic vertebral compression fractures (OVCF) are common in aged population with bone fragility. This study aimed to identify the radiographic and bone fragility characteristic of acute single and multiple OVCF.

METHODS

OVCF patients hospitalized in a spine center between June 2016 and October 2020 were retrospectively studied. Demographics, comorbidity, bone mineral density, spine trauma, duration of pre-hospital back pain, anatomical location and distribution pattern of OVCF, extent of vertebral marrow edema, and degree of vertebral compression of patients with multi-segment vertebral fractures (MSVF) were summarized and compared to those with single segment vertebral fractures (SSVF).

RESULTS

A total of 1182 patients with 1530 acute fractured vertebrae were included. There were 944 SSVF (79.9%) and 238 MSVF (20.1%) simultaneously involving two (MSVF-2) or three and more vertebra (MSVF-3/m). The Female-Male ratio was 4.4 and differed not significantly between SSVF and MSVF. Females in SSVF were younger than males while MSVF-2 tended to occur in older females. L1, T12, and L2 were the three most frequently fractured vertebra and MSVF involved more vertebra in thoracic and lumbar spine. 31.1% in MSVF-2 and 83.1% in MSVF-3/m had at least two vertebral fractures in adjacent. The fractured thoracolumbar vertebra in MSVF was less compressed than that in SSVF. Apparent spine trauma was reported by 61.4% of SSVF, 44.1% of MSVF-2, and 36.3% of MSVF-3/m, while early hospitalization with pre-hospital back pain ≤ 1 week was 58.9% in SSVF, 45.3% in MSVF-2, and 25.9% in MSVF-3/m. Only females aged 70-80 years old in MSVF-3/m showed lower baseline bone mineral density than in MSVF-2 and SSVF. MSVF were not associated with increased comorbidity of hypertension, diabetes, coronary heart disease, cerebral infarction, and chronic pulmonary disease.

CONCLUSIONS

20% of acute OVCF can involve multiple vertebra without significant spine trauma or lower baseline bone mineral density. Multiple OVCF tend to occur in adjacent vertebra with less thoracolumbar vertebral compression but longer duration of pre-hospital back pain.

Labour's pain: strenuous subsistence work, mechanical wear-and-tear and musculoskeletal pain in a non-industrialized population

Musculoskeletal pain is the most debilitating human health condition. Neurophysiological pain mechanisms are highly conserved and promote somatic maintenance and learning to avoid future harm. However, some chronic pain might be more common owing to mismatches between modern lifestyles and traits that originally evolved under distinct premodern conditions. To inform assumptions about factors affecting chronic pain vulnerability prior to industrialization, we assess pain prevalence, perceived causes, and predictors among Tsimane forager-horticulturalists. Habitual subsistence work is the primary reported cause of pain throughout life for both sexes, and pain is more common with age, especially in the back, and for those with more musculoskeletal problems. Sex differences in pain are relatively weak, and we find no association between women's reproductive history and pain, contrary to the hypothesis that reproduction causes women's greater pain susceptibility. Age-standardized current pain prevalence is 1.7-8.2 times higher for Tsimane than other select populations, and Tsimane chronic pain prevalence is within the range of variation observed elsewhere. Chronic low back pain is not a 'mismatch disease' limited to post-industrialized populations. Hominin musculoskeletal changes supporting bipedalism probably imposed health costs, which, after millions of years of evolution, remain an epidemiological burden that may be exacerbated by modern conditions.

Applying a homogeneous pressure distribution to the upper vertebral endplate: Validation of a new loading system, pilot application to human vertebral bodies, and finite element predictions of DIC measured displacements and strains

Image-based personalized Finite Element Models (pFEM) could detect alterations in physiological deformation of human vertebral bodies, but their accuracy has been seldom reported. Meaningful validation experiments should allow vertebral endplate deformability and ensure well-controlled boundary conditions. This study aimed to (i) validate a new loading system to apply a homogeneous pressure on the vertebral endplate during vertebral body compression regardless of endplate deformation; (ii) perform a pilot study on human vertebral bodies measuring surface displacements and strains with Digital Image Correlation (DIC); (iii) determine the accuracy of pFEM of the vertebral bodies. Homogeneous pressure application was achieved by pressurizing a fluid silicone encased in a rubber silicone film acting on the cranial endplate. The loading system was validated by comparing DIC-measured longitudinal strains and lower-end contact pressures, measured on three homogeneous pseudovertebrae of constant transversal section at 2.0 kN, against theoretically calculated values. Longitudinal strains and contact pressures were rather homogeneous, and their mean values close to theoretical calculations (5% underestimation). DIC measurements of surface longitudinal and circumferential displacements and strains were obtained on three human vertebral bodies at 2.0 kN. Complete displacement and strain maps were achieved for anterolateral aspects with random errors ≤0.2 μm and ≤30 μstrain, respectively. Venous plexus and double curvatures limited the completeness and accuracy of DIC data in posterior aspects. pFEM of vertebral bodies, including cortical bone mapping, were built from computed tomography images. In anterolateral aspects, pFEM accuracy of the three vertebrae was: (i) comparable to literature in terms of longitudinal displacements (R2>0.8); (ii) extended to circumferential displacements (pooled data: R2>0.9) and longitudinal strains (zero median error, 95% error: <27%). Circumferential strains were overestimated (median error: 39%). The new methods presented may permit to study how physiological and pathologic conditions influence the ability of vertebral endplates/bodies to sustain loads.

Stiffness and Strain Properties Derived From Digital Tomosynthesis-Based Digital Volume Correlation Predict Vertebral Strength Independently From Bone Mineral Density

Vertebral fractures are the most common osteoporotic fractures, but their prediction using standard bone mineral density (BMD) measurements from dual energy X-ray absorptiometry (DXA) is limited in accuracy. Stiffness, displacement, and strain distribution properties derived from digital tomosynthesis-based digital volume correlation (DTS-DVC) have been suggested as clinically measurable metrics of vertebral bone quality. However, the extent to which these properties correlate to vertebral strength is unknown. To establish this relationship, two independent experiments, one examining isolated T11 and the other examining L3 vertebrae within the L2-L4 segments from cadaveric donors were utilized. Following DXA and DTS imaging, the specimens were uniaxially compressed to fracture. BMD, bone mineral content (BMC), and bone area were recorded for the anteroposterior and lateromedial views from DXA, stiffness, endplate to endplate displacement and distribution statistics of intravertebral strains were calculated from DTS-DVC and vertebral strength was measured from mechanical tests. Regression models were used to examine the relationships of strength with the other variables. Correlations of BMD with vertebral strength varied between experimental groups (R2adj = 0.19-0.78). DTS-DVC derived properties contributed to vertebral strength independently from BMD measures (increasing R2adj to 0.64-0.95). DTS-DVC derived stiffness was the best single predictor (R2adj = 0.66, p < 0.0001) and added the most to BMD in models of vertebral strength for pooled T11 and L3 specimens (R2adj = 0.95, p < 0.0001). These findings provide biomechanical relevance to DTS-DVC calculated properties of vertebral bone and encourage further efforts in the development of the DTS-DVC approach as a clinical tool.

Vitamin D status, including serum levels and sun exposure are associated or correlated with bone mass measurements diagnosis, and bone density of the spine

BACKGROUND

Osteoporosis is a health complication worldwide, especially in developing countries. The prevalence was reported to be 18.3% globally. While the effect of biochemical factors on fracture risk/odds has been documented, the association/correlation between serum 25(OH)D levels, vitamin D dietary intake, and sun exposure with bone mineral density (BMD) remains controversial. This study aimed to evaluate the association and correlation between vitamin D status, including serum levels, dietary intakes, and sun exposure with BMD. We hypothesized that vitamin D-related factors would have different correlations/associations with BMD, which would help better evaluate future studies' results.

METHODS

A total of 186 individuals were included in this study (winter 2020). BMD was measured by Dual-energy X-ray absorptiometry. Blood serum levels of 25(OH)D, phosphorus, calcium, parathyroid hormone (PTH), and calcitonin were tested using standard lab tests. Valid and reliable questionnaires were used for sun exposure assessment and dietary intakes.

RESULTS

There was a significant protective association between spine BMD (classifications, two groups) (OR = 0.69, 95%CI: 0.50-0.94; p-value = 0.023), BMD diagnosis (classifications, two groups) (OR = 0. 69, 95%CI: 0.49-0.87; p-value = 0.036) and sun exposure. There was a significant and moderate correlation between Spine measurements (Spine BMD: Pearson correlation coefficient = 0.302, p-value = 0.046; Spine T-score: Pearson correlation coefficient = 0.322, p-value = 0.033, Spine Z-score: Pearson correlation coefficient = 0.328, p-value = 0.030) and serum 25(OH)D. In addition, participants with osteopenia and osteoporosis significantly consume a higher amount of soluble fiber than the normal BMD group. There was no significant correlation between vitamin D intake and BMD.

CONCLUSION

In conclusion, serum 25(OH)D levels and sun exposure are correlated and associated with BMD. However, prospective studies are needed to investigate the association between dietary vitamin D intake and BMD.

Spinal Fractures during Touristic Motorboat Sea Cruises: An Underestimated and Avoidable Phenomenon

PURPOSE

Each summer, many vacationers enjoy the Mediterranean Sea shores. Among the recreational nautical activities, motorboat cruise is a popular choice that leads to a significant number of thoracolumbar spine fractures at our clinic. This phenomenon seems to be underreported, and its injury mechanism remains unclear. Here, we aim to describe the fracture pattern and propose a possible mechanism of injury.

METHODS

We retrospectively reviewed the clinical, radiological, and contextual parameters of all motorboat-related spinal fracture cases during a 14-year period (2006-2020) in three French neurosurgical level I centers bordering the Mediterranean Sea. Fractures were classified according to the AOSpine thoracolumbar classification system.

RESULTS

A total of 79 patients presented 90 fractures altogether. Women presented more commonly than men (61/18). Most of the lesions occurred at the thoracolumbar transition region between T10 and L2 (88.9% of the levels fractured). Compression A type fractures were seen in all cases (100%). Only one case of posterior spinal element injury was observed. The occurrence of neurological deficit was rare (7.6%). The most commonly encountered context was a patient sitting at the boat's bow, without anticipating the trauma, when the ship's bow suddenly elevated while crossing another wave, resulting in a "deck-slap" mechanism hitting and propelling the patient in the air.

CONCLUSIONS

Thoracolumbar compression fractures are a frequent finding in nautical tourism. Passengers seated at the boat's bow are the typical victims. Some specific biomechanical patterns are involved with the boat's deck suddenly elevating across the waves. More data with biomechanical studies are necessary to understand the phenomenon. Prevention and safety recommendations should be given before motorboat use to fight against these avoidable fractures.

Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges

Numerical models of the musculoskeletal system as investigative tools are an integral part of biomechanical and clinical research. While finite element modeling is primarily suitable for the examination of deformation states and internal stresses in flexible bodies, multibody modeling is based on the assumption of rigid bodies, that are connected via joints and flexible elements. This simplification allows the consideration of biomechanical systems from a holistic perspective and thus takes into account multiple influencing factors of mechanical loads. Being the source of major health issues worldwide, the human spine is subject to a variety of studies using these models to investigate and understand healthy and pathological biomechanics of the upper body. In this review, we summarize the current state-of-the-art literature on multibody models of the thoracolumbar spine and identify limitations and challenges related to current modeling approaches.

Association between sagittal alignment and loads at the adjacent segment in the fused spine: a combined clinical and musculoskeletal modeling study of 205 patients with adult spinal deformity

PURPOSE

Sagittal malalignment is a risk factor for mechanical complications after surgery for adult spinal deformity (ASD). Spinal loads, modulated by sagittal alignment, may explain this relationship. The aims of this study were to investigate the relationships between: (1) postoperative changes in loads at the proximal segment and realignment, and (2) absolute postoperative loads and postoperative alignment measures.

METHODS

A previously validated musculoskeletal model of the whole spine was applied to study a clinical sample of 205 patients with ASD. Based on clinical and radiographic data, pre-and postoperative patient-specific alignments were simulated to predict loads at the proximal segment adjacent to the spinal fusion.

RESULTS

Weak-to-moderate associations were found between pre-to-postop changes in lumbar lordosis, LL (r =  - 0.23, r =  - 0.43; p < 0.001), global tilt, GT (r = 0.26, r = 0.38; p < 0.001) and the Global Alignment and Proportion score, GAP (r = 0.26, r = 0.37; p < 0.001), and changes in compressive and shear forces at the proximal segment. GAP score parameters, thoracic kyphosis measurements and the slope of upper instrumented vertebra were associated with changes in shear. In patients with T10-pelvis fusion, moderate-to-strong associations were found between postoperative sagittal alignment measures and compressive and shear loads, with GT showing the strongest correlations (r = 0.75, r = 0.73, p < 0.001).

CONCLUSIONS

Spinal loads were estimated for patient-specific full spinal alignment profiles in a large cohort of patients with ASD pre-and postoperatively. Loads on the proximal segments were greater in association with sagittal malalignment and malorientation of proximal vertebra. Future work should explore whether they provide a causative mechanism explaining the associated risk of proximal junction complications.

Validation of an EMG submaximal method to calibrate a novel dynamic EMG-driven musculoskeletal model of the trunk: Effects on model estimates

BACKGROUND

Multijoint EMG-assisted optimization models are reliable tools to predict muscle forces as they account for inter- and intra-individual variations in activation. However, the conventional method of normalizing EMG signals using maximum voluntary contractions (MVCs) is problematic and introduces major limitations. The sub-maximal voluntary contraction (SVC) approaches have been proposed as a remedy, but their performance against the MVC approach needs further validation particularly during dynamic tasks.

METHODS

To compare model outcomes between MVC and SVC approaches, nineteen healthy subjects performed a dynamic lifting task with two loading conditions.

RESULTS

Results demonstrated that these two approaches produced highly correlated results with relatively small absolute and relative differences (<10 %) when considering highly-aggregated model outcomes (e.g. compression forces, stability indices). Larger differences were, however, observed in estimated muscle forces. Although some model outcomes, e.g. force of abdominal muscles, were statistically different, their effect sizes remained mostly small (η_G² ≤ 0.13) and in a few cases moderate (η_G² ≤ 0.165).

CONCLUSION

The findings highlight that the MVC calibration approach can reliably be replaced by the SVC approach when the true MVC exertion is not accessible due to pain, kinesiophobia and/or the lack of proper training.

Risk factors and strategies for recovery quality, postoperative pain, and recurrent fractures between percutaneous kyphoplasty and percutaneous vertebroplasty in elderly patients with thoracolumbar compression fractures: a retrospective comparative cohort study

Background

With the increase of clinical cases and the improvement of operation, we found that recurrent fracture of the adjacent vertebral body is a common long-term complication of percutaneous kyphoplasty (PKP). However, the mechanism of re-fracture of adjacent vertebrae after PKP has not been unified. Therefore, through retrospective study, this paper discussed the risk factors and countermeasures affecting the quality of rehabilitation, postoperative pain and recurrent fracture in elderly PKP patients.

Methods

From December 2019 to May 2021, 313 patients with osteoporotic spinal fractures were analyzed retrospectively. Cases were allocated to percutaneous vertebroplasty (PVP; n=130) and PKP (n=183) groups according to the modes of operation. Visual analogue scale (VAS), Cobb angle, and Oswestry disability index (ODI) were evaluated. Based on the occurrence of new fractures, the PKP cohort (n=15) and control cohort (n=32) were classified. Questionnaires analyzed the postoperative re-fractures of people with different characteristics, and the influencing factors of postoperative re-fracture were measured by multivariate logistic regression analysis.

Results

The postoperative VAS scores were significantly lower in the PKP group. The ODI scores in the PKP group were considerably lower than those in the PVP group after surgery. Univariate analysis indicated that age, number of injured vertebrae, history of complicated fracture, number of operative vertebrae, and bone mineral density (BMD) were remarkably correlated with recurrent fracture after PKP. Logistic regression analysis indicated that age, operative vertebral body, BMD, and the number of injured vertebrae were independent risk factors for recurrent fracture after PKP. BMI, BMD, low back soft tissue injury, postoperative vertebral height recovery rate, sagittal Cobb angle improvement rate, total diffusion coefficient of bone cement, short-term complications, non-union, and recurrent fracture were the main risk factors of residual low back pain after PKP.

Conclusions

The clinical efficacy of PKP in elderly patients with thoracolumbar vertebral compression fracture is superior to that of PVP. Clinical attention should be paid to identifying high-risk factors for complications after PKP, and preventive measures should be implemented to help reduce the occurrence of recurrent fractures and postoperative residual pain.

The association between vertical laminar fracture and recurrent kyphosis after implant removal of Thoracolumbar burst fracture: a retrospective study

BACKGROUND

Surgeons often encounter recurrent kyphosis of Cobb angle following thoracolumbar burst fracture surgery. Some factors affecting postoperative correction loss have been studied in previous studies, but few have examined the relationship between laminar fractures and postoperative loss of correction.

METHODS

The clinical data of 86 patients with thoracolumbar burst fracture who met the inclusion criteria and were admitted to our Department of Spine Surgery between 2013 and 2020 was retrospectively analyzed. To examine the association between laminar fracturs and postoperative correction loss, demographic and radiographic characteristics of the two groups were analyzed.

RESULTS

The presence or absence of laminar fractures was statistically different between the two groups (P < 0.05). Binary logistic regression analysis showed that laminar fractures and preoperative Cobb were statistically significant in the two groups. There were statistically significant differences in the degree of injury of laminar fractures in the coronal plane between the two groups (P < 0.05).

CONCLUSION

This study investigated that the presence or absence of laminar fractures and preoperative Cobb contribute to loss of correction after thoracolumbar burst fracture surgery. There was a statistically significant difference between full-length and partial-length laminar fractures on the loss of postoperative correction of thoracolumbar burst fractures with laminar fractures.

Effects of geometric individualisation of a human spine model on load sharing: neuro-musculoskeletal simulation reveals significant differences in ligament and muscle contribution

In spine research, two possibilities to generate models exist: generic (population-based) models representing the average human and subject-specific representations of individuals. Despite the increasing interest in subject specificity, individualisation of spine models remains challenging. Neuro-musculoskeletal (NMS) models enable the analysis and prediction of dynamic motions by incorporating active muscles attaching to bones that are connected using articulating joints under the assumption of rigid body dynamics. In this study, we used forward-dynamic simulations to compare a generic NMS multibody model of the thoracolumbar spine including fully articulated vertebrae, detailed musculature, passive ligaments and linear intervertebral disc (IVD) models with an individualised model to assess the contribution of individual biological structures. Individualisation was achieved by integrating skeletal geometry from computed tomography and custom-selected muscle and ligament paths. Both models underwent a gravitational settling process and a forward flexion-to-extension movement. The model-specific load distribution in an equilibrated upright position and local stiffness in the L4/5 functional spinal unit (FSU) is compared. Load sharing between occurring internal forces generated by individual biological structures and their contribution to the FSU stiffness was computed. The main finding of our simulations is an apparent shift in load sharing with individualisation from an equally distributed element contribution of IVD, ligaments and muscles in the generic spine model to a predominant muscle contribution in the individualised model depending on the analysed spine level.

The Role of Multifidus in the Biomechanics of Lumbar Spine: A Musculoskeletal Modeling Study

BACKGROUND

The role of multifidus in the biomechanics of lumbar spine remained unclear.

PURPOSE

This study aimed to investigate the role of multifidus in the modeling of lumbar spine and the influence of asymmetric multifidus atrophy on the biomechanics of lumbar spine.

METHODS

This study considered five different multifidus conditions in the trunk musculoskeletal models: group 1 (with entire multifidus), group 2 (without multifidus), group 3 (multifidus with half of maximum isometric force), group 4 (asymmetric multifidus atrophy on L5/S1 level), and group 5 (asymmetric multifidus atrophy on L4/L5 level). In order to test how different multifidus situations would affect the lumbar spine, four trunk flexional angles (0°, 30°, 60°, and 90°) were simulated. The calculation of muscle activation and muscle force was done using static optimization function in OpenSim. Then, joint reaction forces of L5/S1 and L4/L5 levels were calculated and compared among the groups.

RESULTS

The models without multifidus had the highest normalized compressive forces on the L4/L5 level in trunk flexion tasks. In extreme cases produced by group 2 models, the normalized compressive forces on L4/L5 level were 444% (30° flexion), 568% (60° flexion), and 576% (90° flexion) of upper body weight, which were 1.82 times, 1.63 times, and 1.13 times as large as the values computed by the corresponding models in group 1. In 90° flexion, the success rate of simulation in group 2 was 49.6%, followed by group 3 (84.4%), group 4 (89.6%), group 5 (92.8%), and group 1 (92.8%).

CONCLUSIONS

The results demonstrate that incorporating multifidus in the musculoskeletal model is important for increasing the success rate of simulation and decreasing the incidence of overestimation of compressive load on the lumbar spine. Asymmetric multifidus atrophy has negligible effect on the lower lumbar spine in the trunk flexion posture. The results highlighted the fine-tuning ability of multifidus in equilibrating the loads on the lower back and the necessity of incorporating multifidus in trunk musculoskeletal modeling.

Using static postures to estimate spinal loading during dynamic lifts with participant-specific thoracolumbar musculoskeletal models

Static biomechanical simulations are sometimes used to estimate in vivo kinetic demands because they can be solved efficiently, but this ignores any potential inertial effects. To date, comparisons between static and dynamic analyses of spinal demands have been limited to lumbar joint differences in young males performing sagittal lifts. Here we compare static and dynamic vertebral compressive and shear force estimates during axial, lateral, and sagittal lifting tasks across all thoracic and lumbar vertebrae in older men and women. Participant-specific thoracolumbar full-body musculoskeletal models estimated vertebral forces from recorded kinematics both with and without consideration of dynamic effects, at an identified frame of peak vertebral loading. Static analyses under-predicted dynamic compressive and resultant shear forces, by an average of about 16% for all three lifts across the thoracic and lumbar spine but were highly correlated with dynamic forces (average r² > .95). The study outcomes have the potential to enable standard clinical and occupational estimates using static analyses.

Towards Single Camera Human 3D-Kinematics

Markerless estimation of 3D Kinematics has the great potential to clinically diagnose and monitor movement disorders without referrals to expensive motion capture labs; however, current approaches are limited by performing multiple de-coupled steps to estimate the kinematics of a person from videos. Most current techniques work in a multi-step approach by first detecting the pose of the body and then fitting a musculoskeletal model to the data for accurate kinematic estimation. Errors in training data of the pose detection algorithms, model scaling, as well the requirement of multiple cameras limit the use of these techniques in a clinical setting. Our goal is to pave the way toward fast, easily applicable and accurate 3D kinematic estimation. To this end, we propose a novel approach for direct 3D human kinematic estimation D3KE from videos using deep neural networks. Our experiments demonstrate that the proposed end-to-end training is robust and outperforms 2D and 3D markerless motion capture based kinematic estimation pipelines in terms of joint angles error by a large margin (35% from 5.44 to 3.54 degrees). We show that D3KE is superior to the multi-step approach and can run at video framerate speeds. This technology shows the potential for clinical analysis from mobile devices in the future.

Relationship between Intervertebral Disc Compression Force and Sagittal Spinopelvic Lower Limb Alignment in Elderly Women in Standing Position with Patient-Specific Whole Body Musculoskeletal Model

The intervertebral disc loading based on compensated standing posture in patients with adult spinal deformity remains unclear. We analyzed the relationship between sagittal alignment and disc compression force (Fm). In 14 elderly women, the alignment of the sagittal spinopelvic and lower extremities was measured. Fm was calculated using the Anybody Modeling System. Patients were divided into low sagittal vertical axis (SVA) and high SVA groups. Comparisons between the two groups were performed and the relationship between the Fm and each parameter was examined using Spearman's correlation coefficient (r). The mean lumbar Fm in the high SVA group was 67.6%; significantly higher than that in the low SVA group (p = 0.046). There was a negative correlation between cervical Fm with T1 slope (r = -0.589, p = 0.034) and lumbar Fm with lumbar lordosis (r = -0.566, p = 0.035). Lumbar Fm was positively correlated with center of gravity-SVA (r = 0.615, p = 0.029), T1 slope (r = 0.613, p = 0.026), and SVA (r = 0.612, p = 0.020). The results suggested sagittal malalignment increased the load on the thoracolumbar and lower lumbar discs and was associated with cervical disc loading.

How to implement guidelines and models of care

In subjects older than 50 years, the presence of clinical risk factors (CRFs) for fractures or a recent fracture is the cornerstone for case finding. In patients who are clinically at high short- and long-term risk of fractures (those with a recent clinical fracture or with multiple CRFs), further assessment with bone mineral density (BMD) measurement using dual-energy absorptiometry (DXA), imaging of the spine, fall risk evaluation and laboratory examination contributes to treatment decisions according to the height and modifiability of fracture risk. Treatment is available with anti-resorptive and anabolic drugs, and from the start of treatment a lifelong strategy is needed to decide about continuous, intermittent, and sequential therapy. Implementation of guidelines requires further initiatives for improving case finding, public awareness about osteoporosis and national policies on reimbursement of assessment and therapy.

Evaluation of Load-To-Strength Ratios in Metastatic Vertebrae and Comparison With Age- and Sex-Matched Healthy Individuals

Vertebrae containing osteolytic and osteosclerotic bone metastases undergo pathologic vertebral fracture (PVF) when the lesioned vertebrae fail to carry daily loads. We hypothesize that task-specific spinal loading patterns amplify the risk of PVF, with a higher degree of risk in osteolytic than in osteosclerotic vertebrae. To test this hypothesis, we obtained clinical CT images of 11 cadaveric spines with bone metastases, estimated the individual vertebral strength from the CT data, and created spine-specific musculoskeletal models from the CT data. We established a musculoskeletal model for each spine to compute vertebral loading for natural standing, natural standing + weights, forward flexion + weights, and lateral bending + weights and derived the individual vertebral load-to-strength ratio (LSR). For each activity, we compared the metastatic spines' predicted LSRs with the normative LSRs generated from a population-based sample of 250 men and women of comparable ages. Bone metastases classification significantly affected the CT-estimated vertebral strength (Kruskal-Wallis, p < 0.0001). Post-test analysis showed that the estimated vertebral strength of osteosclerotic and mixed metastases vertebrae was significantly higher than that of osteolytic vertebrae (p = 0.0016 and p = 0.0003) or vertebrae without radiographic evidence of bone metastasis (p = 0.0010 and p = 0.0003). Compared with the median (50%) LSRs of the normative dataset, osteolytic vertebrae had higher median (50%) LSRs under natural standing (p = 0.0375), natural standing + weights (p = 0.0118), and lateral bending + weights (p = 0.0111). Surprisingly, vertebrae showing minimal radiographic evidence of bone metastasis presented significantly higher median (50%) LSRs under natural standing (p < 0.0001) and lateral bending + weights (p = 0.0009) than the normative dataset. Osteosclerotic vertebrae had lower median (50%) LSRs under natural standing (p < 0.0001), natural standing + weights (p = 0.0005), forward flexion + weights (p < 0.0001), and lateral bending + weights (p = 0.0002), a trend shared by vertebrae with mixed lesions. This study is the first to apply musculoskeletal modeling to estimate individual vertebral loading in pathologic spines and highlights the role of task-specific loading in augmenting PVF risk associated with specific bone metastatic types. Our finding of high LSRs in vertebrae without radiologically observed bone metastasis highlights that patients with metastatic spine disease could be at an increased risk of vertebral fractures even at levels where lesions have not been identified radiologically.

Determination of the 3D Human Spine Posture from Wearable Inertial Sensors and a Multibody Model of the Spine

Determination of spine posture is of great interest for the effective prevention, evaluation, treatment and evolution monitoring of spinal disorders. Limitations of traditional imaging systems, including cost, radiation exposure (for X-ray based systems), projection volume issues and subject positioning requirements, etc., make non-invasive motion assessment tools effective alternatives for clinical and non-clinical use. In this work, a procedure was developed to obtain a subject-specific multibody model of the spine using either inertial or optical sensors and, based on this multibody model, to estimate the locations and orientations of the 17 vertebrae constituting the thoracolumbar spine. The number and calibration of the sensors, angular offsets, scaling difficulties and gender differences were addressed to achieve an accurate 3D-representation of the spine. The approach was validated by comparing the estimated positions of the sensors on 14 healthy subjects with those provided by an optical motion capture system. A mean position error of lower than 12 mm was obtained, thus showing that the proposed method can offer an effective non-invasive tool for the assessment of spine posture.

Skin marker-based subject-specific spinal alignment modeling: A feasibility study

Musculoskeletal models have the potential to improve diagnosis and optimize clinical treatment by predicting accurate outcomes on an individual basis. However, the subject-specific modeling of spinal alignment is often strongly simplified or is based on radiographic assessments, exposing subjects to unnecessary radiation. We therefore developed and introduced a novel skin marker-based approach for modeling subject-specific spinal alignment and evaluated its feasibility by comparing the predicted L1/L2 spinal loads during various functional activities with the loads predicted by the generically scaled models as well as with in vivo measured data obtained from the OrthoLoad database. Spinal loading simulations resulted in considerably higher compressive forces for both scaling approaches over all simulated activities, and AP shear forces that were closer or similar to the in vivo data for the subject-specific approach during upright standing activities and for the generic approach during activities that involved large flexions. These results underline the feasibility of the proposed method and associated workflow for inter- and intra-subject investigations using musculoskeletal simulations. When implemented into standard model scaling workflows, it is expected to improve the accuracy of muscle activity and joint loading simulations, which is crucial for investigations of treatment effects or pathology-dependent deviations.

Risk factors affecting vertebral collapse and kyphotic progression in postmenopausal osteoporotic vertebral fractures

INTRODUCTION

We aimed to investigate the risk factors that affect vertebral deformity 6 months after osteoporotic vertebral fractures (OVFs) at the time of injury.

MATERIALS AND METHODS

From May 2017 to May 2020, 70 postmenopausal women with OVFs were evaluated for age; body mass index; number of previous OVFs; total 25-hydroxy vitamin D [25(OH)D] levels; posterior wall injury on computed tomography; cross-sectional area (CSA) of the psoas major, erector spinae, and multifidus; fat infiltration; vertebral instability (VI) upon admission; collapse rate (CR); and kyphotic angle (KA) at 6 months after injury. A multiple regression analysis was conducted to identify the risk factors for the CR and KA.

RESULTS

The CR was correlated with posterior wall injury (r = 0.295, p = 0.022), 25(OH)D levels (r =  - 0.367, p = 0.002), and VI (r = 0.307, p = 0.010). In the multiple regression analysis, the 25(OH)D levels (p = 0.032) and VI (p = 0.035) were significant risk factors for the CR at the 6-month follow-up. The KA was correlated with the 25(OH)D levels (r =  - 0.262, p = 0.031) and VI (r = 0.298, p = 0.012). In the multiple regression analysis, the CSA of the psoas major (p = 0.011) and VI (p < 0.001) were significant risk factors for the KA at the 6-month follow-up.

CONCLUSION

In cases with large VI at the time of injury, the CR and KA were significantly higher at 6 months after injury. Moreover, the CR was affected by the 25(OH)D level, while the KA was affected by the CSA of the psoas major upon admission.

Using Statistical Shape and Appearance Modelling to characterise the 3D shape and material properties of human lumbar vertebrae: A proof of concept study

Patient variation affects the outcomes of a range of spinal interventions, from disc replacement to vertebral fixation and vertebroplasty. Statistical Shape and Appearance Modelling (SSAM) can be used to describe anatomical variation and pathological differences within the population. To better understand how bone density and shape variation affect load transfer with respect to surgical treatments, Finite Element (FE) models can be generated from a SSAM. The aim for this study is to understand whether geometric and density variation as well as multiple vertebral levels can be incorporated into a single SSAM and whether this can be used to investigate the relationships between, and effects of, the various modes of variation. FE models of 14 human lumbar vertebrae that had been μCT imaged and validated through experimental testing were used as input specimens for a SSAM. The validity of the SSAM was evaluated by using principal component analysis to identify the primary modes of geometric and bone density variation and comparing to those in the input set. FE models were generated from the SSAM to examine the response to loading. The mean error between the input set and generated models for volume, mean density and FE compressive stiffness were 10%, 3% and 10% respectively. Principal Component (PC) 1 captured the majority of the bone density variation. The remaining PCs described specific geometric variation. The FE models generated from the SSAM showed the variations in vertebral stiffness as a result of complex relationships between bone density and shape. The SSAM created has limited data for its input set, however, it acts as a proof of concept for the novel combination of material and shape variation into a single shape model. This approach and the tools developed can be applied to wider patient groups and treatment scenarios to improve patient stratification and to optimise treatments.

Effect of Yoga on Health-Related Outcomes in People at Risk of Fractures: A Systematic Review

We summarized the effects of yoga on health-related outcomes and adverse events in men and postmenopausal women ≥50 years-old at increased risk of fracture, to inform the updated Osteoporosis Canada clinical practice guidelines. Six databases were searched for observational studies, randomized controlled trials and case series. Certainty of evidence was assessed using the Grading of Recommendations, Assessment, Development and Evaluation handbook. Nine studies were included and reported using narrative syntheses due to the limited available evidence. Overall, the available evidence was of very low certainty. There was no effect of yoga on health-related quality of life in randomized trials. Effects on other health-related outcomes were mixed or not available in the literature. Five studies reported no adverse events directly related to the study intervention, and two studies did not report whether adverse events occurred. However, two case series reported vertebral fractures related to yoga participation, possibly due to excessive spinal flexion. Due to the limited and very low certainty evidence, guideline developers will need to draw indirect evidence from yoga studies among middle aged or older adults that are not at fracture risk. (PROSPERO: CRD42019124898) NOVELTY BULLETS: • Evidence in general was of very low certainty. • Yoga had no effect on health-related quality of life in randomized trials. Evidence was mixed or unavailable for other outcomes. • Case studies reported yoga poses involving spinal flexion coincided with incidents of vertebral compression fracture among older adults with increased fracture risk.

Walking Biomechanics and Spine Loading in Patients With Symptomatic Lumbar Spinal Stenosis

Symptomatic lumbar spinal stenosis is a leading cause of pain and mobility limitation in older adults. It is clinically believed that patients with lumbar spinal stenosis adopt a flexed trunk posture or bend forward and alter their gait pattern to improve tolerance for walking. However, a biomechanical assessment of spine posture and motion during walking is broadly lacking in these patients. The purpose of this study was to evaluate lumbar spine and pelvic sagittal angles and lumbar spine compressive loads in standing and walking and to determine the effect of pain and neurogenic claudication symptoms in patients with symptomatic lumbar spinal stenosis. Seven participants with symptomatic lumbar spinal stenosis, aged 44–82, underwent a 3D opto-electronic motion analysis during standing and walking trials in asymptomatic and symptomatic states. Passive reflective marker clusters (four markers each) were attached to participants at T1, L1, and S2 levels of the spine, with additional reflective markers at other spinal levels, as well as the head, pelvis, and extremities. Whole-body motion data was collected during standing and walking trials in asymptomatic and symptomatic states. The results showed that the spine was slightly flexed during walking, but this was not affected by symptoms. Pelvic tilt was not different when symptoms were present, but suggests a possible effect of more forward tilt in both standing (p = 0.052) and walking (p = 0.075). Lumbar spine loading during symptomatic walking was increased by an average of 7% over asymptomatic walking (p = 0.001). Our results did not show increased spine flexion (adopting a trunk-flexed posture) and only indicate a trend for a small forward shift of the pelvis during both symptomatic walking and standing. This suggests that provocation of symptoms in these patients does not markedly affect their normal gait kinematics. The finding of increased spine loading with provocation of symptoms supports our hypothesis that spine loading plays a role in limiting walking function in patients with lumbar spinal stenosis, but additional work is needed to understand the biomechanical cause of this increase.

Type, size, and position of metastatic lesions explain the deformation of the vertebrae under complex loading conditions

BACKGROUND

Bone metastases may lead to spine instability and increase the risk of fracture. Scoring systems are available to assess critical metastases, but they lack specificity, and provide uncertain indications over a wide range, where most cases fall. The aim of this work was to use a novel biomechanical approach to evaluate the effect of lesion type, size, and location on the deformation of the metastatic vertebra.

METHOD

Vertebrae with metastases were identified from 16 human spines from a donation programme. The size and position of the metastases, and the Spine Instability Neoplastic Score (SINS) were evaluated from clinical Quantitative Computed Tomography images. Thirty-five spine segments consisting of metastatic vertebrae and adjacent healthy controls were biomechanically tested in four different loading conditions. The strain distribution over the entire vertebral bodies was measured with Digital Image Correlation. Correlations between the features of the metastasis (type, size, position and SINS) and the deformation of the metastatic vertebrae were statistically explored.

RESULTS

The metastatic type (lytic, blastic, mixed) characterizes the vertebral behaviour (Kruskal-Wallis, p = 0.04). In fact, the lytic metastases showed more critical deformation compared to the control vertebrae (average: 2-fold increase, with peaks of 14-fold increase). By contrast, the vertebrae with mixed or blastic metastases did not show a clear trend, with deformations similar or lower than the controls. Once the position of the lytic lesion with respect to the loading direction was taken into account, the size of the lesion was significantly correlated with the perturbation to the strain distribution (r² = 0.72, p < 0.001). Conversely, the SINS poorly correlated with the mechanical evidence, and only in case of lytic lesions (r² = 0.25, p < 0.0001).

CONCLUSION

These results highlight the relevance of the size and location of the lytic lesion, which are marginally considered in the current clinical scoring systems, in driving the spinal biomechanical instability. The strong correlation with the biomechanical evidence indicates that these parameters are representative of the mechanical competence of the vertebra. The improved explanatory power compared to the SINS suggests including them in future guidelines for the clinical practice.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.