Lumbar spinal loads vary with body height and weight

Anatomically and mechanically conforming patient-specific spinal fusion cages designed by full-scale topology optimization

Cage subsidence after instrumented lumbar spinal fusion surgery remains a significant cause of treatment failure, specifically for posterior or transforaminal lumbar interbody fusion. Recent advancements in computational techniques and additive manufacturing, have enabled the development of patient-specific implants and implant optimization to specific functional targets. This study aimed to introduce a novel full-scale topology optimization formulation that takes the structural response of the adjacent bone structures into account in the optimization process. The formulation includes maximum and minimum principal strain constraints that lower strain concentrations in the adjacent vertebrae. This optimization approach resulted in anatomically and mechanically conforming spinal fusion cages. Subsidence risk was quantified in a commercial finite element solver for off-the-shelf, anatomically conforming and the optimized cages, in two representative patients. We demonstrated that the anatomically and mechanically conforming cages reduced subsidence risk by 91% compared to an off-the-shelf implant with the same footprint for a patient with normal bone quality and 54% for a patient with osteopenia. Prototypes of the optimized cage were additively manufactured and mechanically tested to evaluate the manufacturability and integrity of the design and to validate the finite element model.

DXA beyond bone mineral density and the REMS technique: new insights for current radiologists practice

Osteoporosis is the most prevalent skeletal disorder, a condition that is associated with significant social and healthcare burden. In the elderly, osteoporosis is commonly associated with sarcopenia, further increasing the risk of fracture. Several imaging techniques are available for a non-invasive evaluation of osteoporosis and sarcopenia. This review focuses on dual-energy X-ray absorptiometry (DXA), as this technique offers the possibility to evaluate bone mineral density and body composition parameters with good precision and accuracy. DXA is also able to evaluate the amount of aortic calcification for cardiovascular risk estimation. Additionally, new DXA-based parameters have been developed in recent years to further refine fracture risk estimation, such as the Trabecular Bone Score and the Bone Strain Index. Finally, we describe the recent advances of a newly developed ultrasound-based technology known as Radiofrequency Echographic Multi-Spectrometry, which represent the latest non-ionizing approach for osteoporosis evaluation at central sites.

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.

Effect of patient-specific factors on regeneration in lumbar spine at healthy disc and total disc replacement. Computer simulation

BACKGROUND AND OBJECTIVE

Degenerative diseases of the spine have a negative impact on the quality of life of patients. This study presents the results of numerical modelling of the mechanical behaviour of the lumbar spine with patient-specific conditions at physiological loads. This paper aims to numerically study the influence of degenerative changes in the spine and the presence of an endoprosthesis on the creation of conditions for tissue regeneration.

METHODS

A numerical model of the mechanical behaviour of lumbar spine at healthy and after total disc replacement under low-energy impacts equivalent to physiological loads is presented. The model is based on the movable cellular automaton method (discrete elements), where the mechanical behaviour of bone tissue is described using the Biot poroelasticity accounting for the presence and transfer of interstitial biological fluid. The nutritional pathways of the intervertebral disc in cases of healthy and osteoporotic bone tissues were predicted based on the analysis of the simulation results according to the mechanobiological principles.

RESULTS

Simulation of total disc replacement showed that osseointegration of the artificial disc plates occurs only in healthy bone tissue. With total disc replacement in a patient with osteoporosis, there is an area of increased risk of bone resorption in the near-contact area, approximately 1 mm wide, around the fixators. Dynamic loads may improve the osseointegration of the implant in pathological conditions of the bone tissue.

CONCLUSIONS

The results obtained in the case of healthy spine and osteoporotic bone tissues correspond to the experimental data on biomechanics and possible methods of IVD regeneration from the position of mechanobiological principles. The results obtained with an artificial disc (with keel-type fixation) showed that the use of this type of endoprosthesis in healthy bone tissues allows to reproduce the function of the natural intervertebral disc and does not contribute to the development of neoplastic processes. In the case of an artificial disc with osteoporosis of bone tissues, there is a zone with increased risk of tissue resorption and development of neoplastic processes in the area near the contact of the implant attachment. This circumstance can be compensated by increasing the loading level.

Comprehensive assessment of global spinal sagittal alignment and related normal spinal loads in a healthy population

Abnormal postoperative global sagittal alignment (GSA) is associated with an increased risk of mechanical complications after spinal surgery. Typical assessment of sagittal alignment relies on a few selected measures, disregarding global complexity and variability of the sagittal curvature. The normative range of spinal loads associated with GSA has not yet been considered in clinical evaluation. The study objectives were to develop a new GSA assessment method that holistically describes the inherent relationships within GSA and to estimate the related spinal loads. Vertebral endplates were annotated on radiographs of 85 non-pathological subjects. A Principal Component Analysis (PCA) was performed to derive a Statistical Shape Model (SSM). Associations between identified GSA variability modes and conventional alignment measures were assessed. Simulations of respective Shape Modes (SMs) were performed using an established musculoskeletal AnyBody model to estimate normal variation in cervico-thoraco-lumbar loads. The first six principal components explained 97.96% of GSA variance. The SSM provides the normative range of GSA and a visual representation of the main variability modes. Normal variation relative to the population mean in identified alignment features was found to influence spinal loads, e.g. the lower bound of the second shape mode (SM2-2σ) corresponds to an increase in L4L5-compression by 378.64 N (67.86%). Six unique alignment features were sufficient to describe GSA almost entirely, demonstrating the value of the proposed method for an objective and comprehensive analysis of GSA. The influence of these features on spinal loads provides a normative biomechanical reference, eventually guiding surgical planning of deformity correction in the future.

Postoperative Sclerotic Modic Changes After Transforaminal Lumbar Interbody Fusion: The Prevalence, Risk Factors, and Impact on Fusion

STUDY DESIGN

A retrospective cohort study.

OBJECTIVE

This study aimed to assess postoperative sclerotic modic changes (MCs) following transforaminal lumbar interbody fusion for lumbar degenerative disc disease, investigating their prevalence, risk factors, and association with clinical outcomes.

SUMMARY OF BACKGROUND DATA

Sclerotic MCs may occur in patients with lumbar degenerative disc disease after lumbar interbody fusion. The incidence and characteristics of postoperative sclerotic MCs, as well as their clinical impact, are unknown.

MATERIALS AND METHODS

The study included 467 patients (510 levels) who underwent single or two-level transforaminal lumbar interbody fusion surgery, divided into a postoperative sclerotic MC group (60 patients, 66 levels) and a non-MC group (407 patients, 444 levels). The time of development and location of postoperative sclerotic MCs, fusion rate, cage subsidence, bilateral process decompression, and cross-link usage were recorded. Preoperative, postoperative, and follow-up visual analogue scale and Oswestry disability index scores were collected. Multivariable logistic regression was used to evaluate factors associated with the development of postoperative sclerotic MCs.

RESULTS

The prevalence of postoperative sclerotic MCs was 12.8%. The postoperative sclerotic MC group had a higher body mass index (BMI). The postoperative sclerotic MC group demonstrated a fusion rate of 47%, significantly lower than that of the non-MC group (71%) at six months post-operation. At final follow-up, the fusion rate in the postoperative sclerotic MC group was 62%, significantly lower than that of the non-MC group (86%). Postoperative visual analogue scale and Oswestry disability index scores were significantly higher in the group with postoperative sclerotic MCs. BMI and osteoporosis were significantly associated with the development of postoperative sclerotic MCs.

CONCLUSION

Postoperative sclerotic MCs generally appear within the first year after surgery, with a prevalence of 12.8%. The presence of postoperative sclerotic MCs can adversely impact postoperative outcomes. To prevent postoperative sclerotic MCs, the authors postulate extending the immobilization period with external bracing and improving the management of BMI and osteoporosis in the perioperative time window.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

Effect of obesity on spinal loads during load-reaching activities: A subject- and kinematics-specific musculoskeletal modeling approach

Obesity has been associated to increase the risk of low back disorders. Previous musculoskeletal models simulating the effect of body weight on intervertebral joint loads have assumed identical body postures for obese and normal-weight individuals during a given physical activity. Our recent kinematic-measurement studies, however, indicate that obese individuals adapt different body postures (segmental orientations) than normal-weight ones when performing load-reaching activities. The present study, therefore, used a subject- and kinematics-specific musculoskeletal modeling approach to compare spinal loads of nine normal-weight and nine obese individuals each performing twelve static two-handed load-reaching activities at different hand heights, anterior distances, and asymmetry angles (total of 12 tasks × 18 subjects = 216 model simulations). Each model incorporated personalized muscle architectures, body mass distributions, and full-body kinematics for each subject and task. Results indicated that even when accounting for subject-specific body kinematics obese individuals experienced significantly larger (by ∼38% in average) L5-S1 compression (2305 ± 468 N versus 1674 ± 337 N) and shear (508 ± 111 N versus 705 ± 150 N) loads during all reaching activities (p < 0.05 for all hand positions). This average difference of ∼38% was similar to the results obtained from previous modeling investigations that neglected kinematics differences between the two weight groups. Moreover, there was no significant interaction effect between body weight and hand position on the spinal loads; indicating that the effect of body weight on L5-S1 loads was not dependent on the position of hands. Postural differences alone appear, hence, ineffective in compensating the greater spinal loads that obese people experience during reaching activities.

Hybrid musculoskeletal model-based 3D asymmetric lifting prediction and comparison with symmetric lifting

In this study, a 3D asymmetric lifting motion is predicted by using a hybrid predictive model to prevent potential musculoskeletal lower back injuries for asymmetric lifting tasks. The hybrid model has two modules: a skeletal module and an OpenSim musculoskeletal module. The skeletal module consists of a dynamic joint strength based 40 degrees of freedom spatial skeletal model. The skeletal module can predict the lifting motion, ground reaction forces (GRFs), and center of pressure (COP) trajectory using an inverse dynamics-based motion optimization method. The musculoskeletal module consists of a 324-muscle-actuated full-body lumbar spine model. Based on the predicted kinematics, GRFs and COP data from the skeletal module, the musculoskeletal module estimates muscle activations using static optimization and joint reaction forces through the joint reaction analysis tool in OpenSim. The predicted asymmetric motion and GRFs are validated with experimental data. Muscle activation results between the simulated and experimental EMG are also compared to validate the model. Finally, the shear and compression spine loads are compared to NIOSH recommended limits. The differences between asymmetric and symmetric liftings are also compared.

A neuromuscular human body model for lumbar injury risk analysis in a vibration loading environment

BACKGROUND AND OBJECTIVE

Long-term intensive exposure to whole-body vibration substantially increases the risk of low back pain and degenerative diseases in special occupational groups, like motor vehicle drivers, military vehicle occupants, aircraft pilots, etc. This study aims to establish and validate a neuromuscular human body model focusing on improvement of the detailed description of anatomic structures and neural reflex control, for lumbar injury analysis in vibration loading environments.

METHODS

A whole-body musculoskeletal in Opensim codes was first improved by including a detailed anatomic description of spinal ligaments, non-linear intervertebral disc, and lumbar facet joints, and coupling a proprioceptive feedback closed-loop control strategy with GTOs and muscle spindles modeling in Python codes. Then, the established neuromuscular model was multi-levelly validated from sub-segments to the whole model, from regular movements to dynamic responses to vibration loadings. Finally, the neuromuscular model was combined with a dynamic model of an armored vehicle to analyze occupant lumbar injury risk in vibration loadings due to different road conditions and traveling velocities.

RESULT

Based on a series of biomechanical indexes, including lumbar joint rotation angles, the lumbar intervertebral pressures, the displacement of the lumbar segments, and the lumbar muscle activities, the validation results show that the present neuromuscular model is available and feasible in predicting lumbar biomechanical responses in normal daily movement and vibration loading environments. Furthermore, the combined analysis with the armored vehicle model predicted similar lumbar injury risk to the experimental or epidemiologic studies. The preliminary analysis results also showed that road types and travelling velocities have substantial combined effects on lumbar muscle activities, and indicated that intervertebral joint pressure and muscle activity indexes can need to be jointly considered for lumbar injury risk evaluation.

CONCLUSION

In conclusion, the established neuromuscular model is an effective tool to evaluate vibration loading effects on injury risk of the human body and assist vehicle design vibration comfort by directly concerning the human body injury itself.

Subsidence of a partially porous titanium lumbar cage produced by electron beam melting technology

The lumbar intervertebral devices are widely used in the surgical treatment of lumbar diseases. The subsidence represents a serious clinical issue during the healing process, mainly when the interfaces between the implant and the vertebral bodies are not well designed. The aim of this study is the evaluation of subsidence risk for two different devices. The devices have the same shape, but one of them includes a filling micro lattice structure. The effect of the micro lattice structure on the subsidence behavior of the implant was evaluated by means of both experimental tests and finite element analyses. Compressive tests were carried out by using blocks made of grade 15 polyurethane, which simulate the vertebral bone. Non-linear, quasi-static finite element analyses were performed to simulate experimental and physiologic conditions. The experimental tests and the FE analyses showed that the subsidence risk is higher for the device without micro lattice structure, due to the smaller contact surface. Moreover, an overload in the central zone of the contact surface was detected in the same device and it could cause the implant failure. Thus, the micro lattice structure allows a homogenous pressure distribution at the implant-bone interface.

Assessment of DXA derived bone quality indexes and bone geometry parameters in early breast cancer patients: A single center cross-sectional study

Background

Bone mineral density (BMD) lacks sensitivity in individual fracture risk assessment in early breast cancer (EBC) patients treated with aromatase inhibitors (AIs). New dual-energy X-ray absorptiometry (DXA) based risk factors are needed.

Methods

Trabecular bone score (TBS), bone strain index (BSI) and DXA parameters of bone geometry were evaluated in postmenopausal women diagnosed with EBC. The aim was to explore their association with morphometric vertebral fractures (VFs). Subjects were categorized in 3 groups in order to evaluate the impact of AIs and denosumab on bone geometry: AI-naive, AI-treated minus (AIDen-) or plus (AIDen+) denosumab.

Results

A total of 610 EBC patients entered the study: 305 were AI-naive, 187 AIDen-, and 118 AIDen+. In the AI-naive group, the presence of VFs was associated with lower total hip BMD and T-score and higher femoral BSI. As regards as bone geometry parameters, AI-naive fractured patients reported a significant increase in femoral narrow neck (NN) endocortical width, femoral NN subperiosteal width, intertrochanteric buckling ratio (BR), intertrochanteric endocortical width, femoral shaft (FS) BR and endocortical width, as compared to non-fractured patients. Intertrochanteric BR and intertrochanteric cortical thickness significantly increased in the presence of VFs in AIDen- patients, not in AIDen+ ones. An increase in cross-sectional area and cross-sectional moment of inertia, both intertrochanteric and at FS, significantly correlated with VFs only in AIDen+. No association with VFs was found for either lumbar BSI or TBS in all groups.

Conclusions

Bone geometry parameters are variably associated with VFs in EBC patients, either AI-naive or AI treated in combination with denosumab. These data suggest a tailored choice of fracture risk parameters in the 3 subgroups of EBC patients.

In vivo intervertebral disc mechanical deformation following a treadmill walking "stress test" is inversely related to T1rho relaxation time

OBJECTIVE

To assess the in vivo relationship between the mechanical response of intervertebral discs (IVDs) to dynamic activity and IVD biochemical composition assessed via T1rho relaxation imaging.

DESIGN

Eighteen asymptomatic participants with no history of low back pain (LBP), injury, or surgery underwent magnetic resonance (MR) imaging of their lumbar spine prior to and immediately following a treadmill walking "stress test." Anatomic (SPACE, FLASH) MR images were obtained pre- and post-exercise and utilized to measure IVD mechanical deformation. Quantitative (T1rho) imaging was performed pre-exercise to reflect IVD composition. Pre-exercise anatomic images were also utilized to assess IVD degenerative status based on the modified Pfirrmann scale. To quantify mechanical response, 3D surface models of the L1-L2-L5-S1 IVDs were created from manual segmentations of pre- and post-exercise anatomic images and utilized to assess changes in IVD height. IVD strain (%) was defined as change in IVD height normalized to pre-activity height. Linear mixed models were used to assess the relationships between IVD mechanical deformation (strain), composition (T1rho relaxation time), and degenerative status (Pfirrmann grade).

RESULTS

Increased compressive IVD strain was associated with lower T1rho relaxation times in the nucleus pulposus (NP) of the disc (βT1rho=5.07,CI:[1.52,7.77],Rmarg2=0.52,p=0.005). Thus, an inverse relationship between IVD strain and NP T1rho relaxation time was observed.

CONCLUSION

The in vivo mechanical response of the IVD to the "stress test" was sensitive to differences in NP composition. The results of this study suggest that quantification of in vivo IVD mechanical function and composition may provide insight into IVD health.

Sex-related differences in motor control strategies during isometric contractions in the Sorensen test posture with different external loading

BACKGROUND

There is a sex-related difference in strength and endurance in trunk muscles: males have more strength while females have more endurance.

OBJECTIVE

Investigate sex-related differences in motor control strategies in back muscles during isometric contractions (IC) in the Sorensen test posture.

METHODS

Thirty-six healthy and young volunteers performed different tasks: three maximal voluntary contraction (MVC) tests recorded with a bio-feedback force sensor, and followed by five loaded IC (LIC) tests using bodyweight and loads of 0 to 8 kg with a 2-kg step. Surface electromyography was used to measure the activity of bilateral lumbar paravertebral (LP) and quadratus lumborum (QL) muscles.

RESULTS

Sex-related and load-related differences in EMG amplitudes were highlighted during LIC tests. Females showed significantly (p-value = 0.02) greater EMG amplitude for the highest load (8 kg) than males. Besides, significant differences between low (2 and 4 kg) and high (6 and 8 kg) loads for both LP and QL muscles and for both sexes were observed. Finally, for MVC tests, males produced significantly (p-value = 1.02e-4) greater strength during MVC tests (4.25 ± 1.37 N/kg vs 2.60 ± 0.78 N/kg).

CONCLUSION

The results were clinically relevant to bring attention to load strategy during rehabilitation, particularly in females who seemed to recruit the QL muscle more.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Development of a Computational Model of the Mechanical Behavior of the L4-L5 Lumbar Spine: Application to Disc Degeneration

Degenerative changes in the lumbar spine significantly reduce the quality of life of people. In order to fully understand the biomechanics of the affected spine, it is crucial to consider the biomechanical alterations caused by degeneration of the intervertebral disc (IVD). Therefore, this study is aimed at the development of a discrete element model of the mechanical behavior of the L4-L5 spinal motion segment, which covers all the degeneration grades from healthy IVD to its severe degeneration, and numerical study of the influence of the IVD degeneration on stress state and biomechanics of the spine. In order to analyze the effects of IVD degeneration on spine biomechanics, we simulated physiological loading conditions using compressive forces. The results of modeling showed that at the initial stages of degenerative changes, an increase in the amplitude and area of maximum compressive stresses in the disc is observed. At the late stages of disc degradation, a decrease in the value of intradiscal pressure and a shift in the maximum compressive stresses in the dorsal direction is observed. Such an influence of the degradation of the geometric and mechanical parameters of the tissues of the disc leads to the effect of bulging, which in turn leads to the formation of an intervertebral hernia.

Analysis of Healthcare Push and Pull Task via JACK: Predicted Joint Accuracy during Full-Body Simulation

The posture accuracy of full-body dynamic simulation has been successfully evaluated in JACK Siemens software via analyzing two common push and pull tasks. The difference in joint angles between the actual and predicted human movement directly results in the strength of force exposed on the lumbar spine. In this study, the individual factors, such as body height, body weight, trunk and hip flexion, shoulder movement, and muscle strength between genders, have shown a strong association with the adopted postures and exposed spinal forces during task performance. To provide robust ergonomics analysis, these individual variables should be adequately considered in software design for the long-term goal of injury prevention in diverse occupational workplaces.

The Bone Strain Index: An Innovative Dual X-ray Absorptiometry Bone Strength Index and Its Helpfulness in Clinical Medicine

Bone strain Index (BSI) is an innovative index of bone strength that provides information about skeletal resistance to loads not considered by existing indexes (Bone Mineral Density, BMD. Trabecular Bone Score, TBS. Hip Structural Analysis, HSA. Hip Axis Length, HAL), and, thus, improves the predictability of fragility fractures in osteoporotic patients. This improved predictability of fracture facilitates the possibility of timely intervention with appropriate therapies to reduce the risk of fracture. The development of the index was the result of combining clinical, radiographical and construction-engineering skills. In fact, from a physical point of view, primary and secondary osteoporosis, leading to bone fracture, are determined by an impairment of the physical properties of bone strength: density, internal structure, deformation and fatigue. Dual X-ray absorptiometry (DXA) is the gold standard for assessing bone properties, and it allows measurement of the BMD, which is reduced mainly in primary osteoporosis, the structural texture TBS, which can be particularly degraded in secondary osteoporosis, and the bone geometry (HSA, HAL). The authors recently conceived and developed a new bone deformation index named Bone Strain Index (BSI) that assesses the resistance of bone to loads. If the skeletal structure is equated to engineering construction, these three indexes are all considered to determine the load resistance of the construct. In particular, BSI allows clinicians to detect critical information that BMD and TBS cannot explain, and this information is essential for an accurate definition of a patient’s fracture risk. The literature demonstrates that both lumbar and femoral BSI discriminate fractured osteoporotic people, that they predict the first fragility fracture, and further fragility fractures, monitor anabolic treatment efficacy and detect patients affected by secondary osteoporosis. BSI is a new diagnostic tool that offers a unique perspective to clinical medicine to identify patients affected by primary and, specially, secondary osteoporosis. This literature review illustrates BSI’s state of the art and its ratio in clinical medicine.

In Vitro Studies for Investigating Creep of Intervertebral Discs under Axial Compression: A Review of Testing Environment and Results

Creep responses of intervertebral discs (IVDs) are essential for spinal biomechanics clarification. Yet, there still lacks a well-recognized investigation protocol for this phenomenon. Current work aims at providing researchers with an overview of the in vitro creep tests reported by previous studies, specifically specimen species, testing environment, loading regimes and major results, based on which a preliminary consensus that may guide future creep studies is proposed. Specimens used in creep studies can be simplified as a “bone–disc–bone” structure where three mathematical models can be adopted for describing IVDs’ responses. The preload of 10–50 N for 30 min or three cycles followed by 4 h-creep under constant compression is recommended for ex vivo simulation of physiological condition of long-time sitting or lying. It is worth noticing that species of specimens, environment temperature and humidity all have influences on biomechanical behaviors, and thus are summarized and compared through the literature review. All factors should be carefully set according to a guideline before tests are conducted to urge comparable results across studies. To this end, this review also provides a guideline, as mentioned before, and specific steps that might facilitate the community of biomechanics to obtain more repeatable and comparable results from both natural specimens and novel biomaterials.

Patient-Specific Variations in Local Strain Patterns on the Surface of a Trussed Titanium Interbody Cage

Introduction: 3D printed trussed titanium interbody cages may deliver bone stimulating mechanobiological strains to cells attached at their surface. The exact size and distribution of these strains may depend on patient-specific factors, but the influence of these factors remains unknown. Therefore, this study aimed to determine patient-specific variations in local strain patterns on the surface of a trussed titanium interbody fusion cage.

Materials and Methods: Four patients eligible for spinal fusion surgery with the same cage size were selected from a larger database. For these cases, patient-specific finite element models of the lumbar spine including the same trussed titanium cage were made. Functional dynamics of the non-operated lumbar spinal segments, as well as local cage strains and caudal endplate stresses at the operated segment, were evaluated under physiological extension/flexion movement of the lumbar spine.

Results: All patient-specific models revealed physiologically realistic functional dynamics of the operated spine. In all patients, approximately 30% of the total cage surface experienced strain values relevant for preserving bone homeostasis and stimulating bone formation. Mean caudal endplate contact pressures varied up to 10 MPa. Both surface strains and endplate contact pressures varied more between loading conditions than between patients.

Conclusions: This study demonstrates the applicability of patient-specific finite element models to quantify the impact of patient-specific factors such as bone density, degenerative state of the spine, and spinal curvature on interbody cage loading. In the future, the same framework might be further developed in order to establish a pipeline for interbody cage design optimizations.

Bone Strain Index predicts fragility fracture in osteoporotic women: an artificial intelligence-based study

BACKGROUND

We applied an artificial intelligence-based model to predict fragility fractures in postmenopausal women, using different dual-energy x-ray absorptiometry (DXA) parameters.

METHODS

One hundred seventy-four postmenopausal women without vertebral fractures (VFs) at baseline (mean age 66.3 ± 9.8) were retrospectively evaluated. Data has been collected from September 2010 to August 2018. All subjects performed a spine x-ray to assess VFs, together with lumbar and femoral DXA for bone mineral density (BMD) and the bone strain index (BSI) evaluation. Follow-up exams were performed after 3.34 ± 1.91 years. Considering the occurrence of new VFs at follow-up, two groups were created: fractured versus not-fractured. We applied an artificial neural network (ANN) analysis with a predictive tool (TWIST system) to select relevant input data from a list of 13 variables including BMD and BSI. A semantic connectivity map was built to analyse the connections among variables within the groups. For group comparisons, an independent-samples t-test was used; variables were expressed as mean ± standard deviation.

RESULTS

For each patient, we evaluated a total of n = 6 exams. At follow-up, n = 69 (39.6%) women developed a VF. ANNs reached a predictive accuracy of 79.56% within the training testing procedure, with a sensitivity of 80.93% and a specificity of 78.18%. The semantic connectivity map showed that a low BSI at the total femur is connected to the absence of VFs.

CONCLUSION

We found a high performance of ANN analysis in predicting the occurrence of VFs. Femoral BSI appears as a useful DXA index to identify patients at lower risk for lumbar VFs.

Subject-Specific Alignment and Mass Distribution in Musculoskeletal Models of the Lumbar Spine

Musculoskeletal modeling is a well-established method in spine biomechanics and generally employed for investigations concerning both the healthy and the pathological spine. It commonly involves inverse kinematics and optimization of muscle activity and provides detailed insight into joint loading. The aim of the present work was to develop and validate a procedure for the automatized generation of semi-subject-specific multi-rigid body models with an articulated lumbar spine. Individualization of the models was achieved with a novel approach incorporating information from annotated EOS images. The size and alignment of bony structures, as well as specific body weight distribution along the spine segments, were accurately reproduced in the 3D models. To ensure the pipeline’s robustness, models based on 145 EOS images of subjects with various weight distributions and spinopelvic parameters were generated. For validation, we performed kinematics-dependent and segment-dependent comparisons of the average joint loads obtained for our cohort with the outcome of various published in vivo and in situ studies. Overall, our results agreed well with literature data. The here described method is a promising tool for studying a variety of clinical questions, ranging from the evaluation of the effects of alignment variation on joint loading to the assessment of possible pathomechanisms involved in adjacent segment disease.

Increasing BMI increases lumbar intervertebral disc deformation following a treadmill walking stress test

High body mass index (BMI) and obesity have been implicated as risk factors for lumbar degenerative disc disease and low back pain. Despite this, there is limited in vivo data to quantify how obesity influences the mechanical function of intervertebral discs (IVD) in response to activities of daily living. Recently, our lab has developed methodologies to non-invasively measure in vivo IVD deformation resulting from activities of daily living using magnetic resonance (MR) imaging and solid modeling techniques. This pilot study expands on these methodologies to assess how BMI influences IVD deformation following treadmill walking in eight asymptomatic individuals. Ordinary least squares regression analyses revealed a statistically significant relationship between BMI and compressive deformation (strain (%)) in the L5-S1 IVD (R² = 0.61, p < 0.05). This relationship was weaker in the L3-L4 (R² = 0.28, p > 0.05) and L4-L5 IVDs (R² = 0.28, p > 0.05). Importantly, no relationship between pre-exercise disc height and BMI was identified (p > 0.05). Therefore, the results of this study suggest that BMI may alter the mechanical response of lumbar spine IVDs, particularly at the L5-S1 level. Furthermore, the observed relationship between increased BMI and IVD compressive deformation, in the absence of a detected relationship between pre-exercise disc height and BMI, suggests that changes in IVD mechanical function may be more sensitive to alterations in disc health than static clinical imaging alone. This finding highlights the importance of quantifying disc mechanical function when examining the relationship between BMI and IVD degeneration.

Beyond Bone Mineral Density: A New Dual X-Ray Absorptiometry Index of Bone Strength to Predict Fragility Fractures, the Bone Strain Index

For a proper assessment of osteoporotic fragility fracture prediction, all aspects regarding bone mineral density, bone texture, geometry and information about strength are necessary, particularly in endocrinological and rheumatological diseases, where bone quality impairment is relevant. Data regarding bone quantity (density) and, partially, bone quality (structure and geometry) are obtained by the gold standard method of dual X-ray absorptiometry (DXA). Data about bone strength are not yet readily available. To evaluate bone resistance to strain, a new DXA-derived index based on the Finite Element Analysis (FEA) of a greyscale of density distribution measured on spine and femoral scan, namely Bone Strain Index (BSI), has recently been developed. Bone Strain Index includes local information on density distribution, bone geometry and loadings and it differs from bone mineral density (BMD) and other variables of bone quality like trabecular bone score (TBS), which are all based on the quantification of bone mass and distribution averaged over the scanned region. This state of the art review illustrates the methodology of BSI calculation, the findings of its in reproducibility and the preliminary data about its capability to predict fragility fracture and to monitor the follow up of the pharmacological treatment for osteoporosis.

Bone strain index as a predictor of further vertebral fracture in osteoporotic women: An artificial intelligence-based analysis

BACKGROUND

Osteoporosis is an asymptomatic disease of high prevalence and incidence, leading to bone fractures burdened by high mortality and disability, mainly when several subsequent fractures occur. A fragility fracture predictive model, Artificial Intelligence-based, to identify dual X-ray absorptiometry (DXA) variables able to characterise those patients who are prone to further fractures called Bone Strain Index, was evaluated in this study.

METHODS

In a prospective, longitudinal, multicentric study 172 female outpatients with at least one vertebral fracture at the first observation were enrolled. They performed a spine X-ray to calculate spine deformity index (SDI) and a lumbar and femoral DXA scan to assess bone mineral density (BMD) and bone strain index (BSI) at baseline and after a follow-up period of 3 years in average. At the end of the follow-up, 93 women developed a further vertebral fracture. The further vertebral fracture was considered as one unit increase of SDI. We assessed the predictive capacity of supervised Artificial Neural Networks (ANNs) to distinguish women who developed a further fracture from those without it, and to detect those variables providing the maximal amount of relevant information to discriminate the two groups. ANNs choose appropriate input data automatically (TWIST-system, Training With Input Selection and Testing). Moreover, we built a semantic connectivity map usingthe Auto Contractive Map to provide further insights about the convoluted connections between the osteoporotic variables under consideration and the two scenarios (further fracture vs no further fracture).

RESULTS

TWIST system selected 5 out of 13 available variables: age, menopause age, BMI, FTot BMC, FTot BSI. With training testing procedure, ANNs reached predictive accuracy of 79.36%, with a sensitivity of 75% and a specificity of 83.72%. The semantic connectivity map highlighted the role of BSI in predicting the risk of a further fracture.

CONCLUSIONS

Artificial Intelligence is a useful method to analyse a complex system like that regarding osteoporosis, able to identify patients prone to a further fragility fracture. BSI appears to be a useful DXA index in identifying those patients who are at risk of further vertebral fractures.

Comparison of different lifting analysis tools in estimating lower spinal loads - Evaluation of NIOSH criterion

Excessive loads on the human spine is recognized as a risk factor for back injuries/pain. Various lifting analysis tools such as musculoskeletal models, regression equations and NIOSH (National Institute for Occupational Safety and Health) lifting equation (NLE) have been proposed to evaluate and mitigate associated risks during manual material handling activities. Present study aims to compare predicted spinal loads from 5 different lifting analysis tools as well as to critically evaluate the NIOSH recommended weight limit (RWL). Spinal loads were estimated under different symmetric/asymmetric lifting tasks in which hand-load mass at each task was set based on RWL from NLE. Estimated intradiscal pressures (IDPs) of various tools were also compared with in vivo measurements. We compared RWL by NLE versus our estimations of RWL calculated from our regression equations using biomechanical criteria (compression <3400 N with/without shear <1000, 1250 or 1500 N). Our regression equations followed by OpenSim, AnyBody, simple polynomial and 3DSSPP satisfactorily predicted L4-L5 IDPs. Lifting analysis tools estimated comparable spinal compression forces (mean Pearson's r = 0.80; standard deviation of relative difference = 26%) while in shear, differences were greater (mean Pearson's r = 0.68; standard deviation of relative difference = 56%). NLE estimations of RWL were conservative in comparison with our estimations for lean individuals (BMI < 25 kg/m²) when compression <3400 N and shear <1250 N were considered as the biomechanical criteria. For heavier individuals, however, NLE estimations of RWL generated spinal compression >3400 N (NIOSH biomechanical safety threshold) as well as shear >1000 N. Although RWLs estimated by NLE was body weight independent, body weight substantially altered RWLs estimated from our regression equations. For improved estimation of the risk of injury, more accurate failure criteria for spinal segments are essential.

Selection of suitable pedicle screw for degenerated cortical and cancellous bone of human lumbar spine: A finite element study

Pedicular arthrodesis is the traditional procedure in terms of increase in the biomechanical stability with higher fixation rate. The current work aims to identify the effect of three spinal pedicle screws considering cortical and cancellous degeneracy condition. Lumbar section L2-L3 is utilized and various load and moment conditions were applied to depict the various biomechanical parameters for selection of suitable screw. Three dimensional model is considered in finite element analysis to identify the various responses of pedicle screw at bone screw juncture. Computed tomography (CT) images of a healthy male were considered to generate the finite element vertebral model. Generated intact model was further utilized to develop the other implanted models of degenerated cortical and cancellous bone models. The three fused instrumented models with different cortical and cancellous degeneracy conditions were analyzed in finite element analysis. The results were obtained as stress pattern at bone screw boundary and intervertebral disc stress. FE simulated results represents significant changes in the von Mises stress due to various load and moment conditions on degenerated bones during different body movement conditions. Results have shown that among all pedicle screws, the 6.0 mm diameter screw reflects very less stress values at the juncture. Multiple results on biomechanical aspects obtained during the FE study can be considered to design a new stabilization device and may be helpful to plan surgery of critical sections.

ANALYSIS OF MUSCLE STRENGTH EFFECTS ON EXERCISE PERFORMANCE USING DYNAMIC STABILIZATION EXERCISE DEVICE

Muscle strength may vary depending on the pathological issues and static life habits. These conditions lead to abnormal spinal loads and change muscle strength as well as activation patterns, thereby causing spinal disorders. In this study, the effects of muscle strength on the spine stabilization exercise were analyzed using a whole-body tilt device. Musculoskeletal modeling was performed and the results were validated through a comparison with the electromyography (EMG) analysis results. Based on the validated basic model, modeling was performed for the whole-body tilt device. To examine the exercise effect and muscle activation while the maximum muscle force capacity (MFC) was varied from 30[Formula: see text]N/cm2 to 60[Formula: see text]N/cm2 and 90[Formula: see text]N/cm2, the muscle force was predicted through inverse dynamics analysis. When MFC was 30[Formula: see text]N/cm2, the posterior direction of the tilt could not be analyzed (no solution found). When MFC was 60[Formula: see text]N/cm2, it could be analyzed, but the muscle force was predicted to be higher compared to when MFC was 90[Formula: see text]N/cm2. It was confirmed that muscle strength is a very important element for maintaining postural activities and performing exercise. Therefore, for rehabilitation patients and elderly people with weak muscle strength, hard or extreme exercise may cause musculoskeletal injuries.

Nonlinear stress-dependent recovery behavior of the intervertebral disc

The intervertebral disc exhibits complex mechanics due to its heterogeneous structure, inherent viscoelasticity, and interstitial fluid-matrix interactions. Sufficient fluid flow into the disc during low loading periods is important for maintaining mechanics and nutrient transport. However, there is a lack of knowledge on the effect of loading magnitude on time-dependent recovery behavior and the relative contribution of multiple recovery mechanisms during recovery. In most studies that have evaluated disc recovery behavior, a single load condition has been considered, making it difficult to compare findings across studies. Hence, the objective of this study was to quantify unloaded disc recovery behavior after compressive creep loading under a wide range of physiologically relevant stresses (0.2-2 MPa). First, the repeatability of disc recovery behavior was assessed. Once repeatable recovery behavior was confirmed, each motion segment was subject to three cycles of creep-recovery loading, where each cycle consisted of a 24-h creep at a pre-assigned load (100, 200, 300, 600, 900, or 1200 N), followed by an 18-h recovery period at a nominal load (10 N). Results showed that disc recovery behavior was strongly influenced by the magnitude of loading. The magnitude of instantaneous and time-dependent recovery deformations increased nonlinearly with an increase in compressive stress during creep. In conclusion, this study highlights that elastic deformation, intrinsic viscoelasticity, and poroelasticity all have substantial contributions to disc height recovery during low loading periods. However, their relative contributions to disc height recovery largely depend on the magnitude of loading. While loading history does not influence the contribution of the short-term recovery, the contribution of long-term recovery is highly sensitive to loading magnitude.

Bone strain index in the prediction of vertebral fragility refracture

Dual-energy x-ray absorptiometry (DXA) can provide quantitative (bone mineral density, BMD) and qualitative (trabecular bone score, TBS) indexes of bone status, able to predict fragility fractures in most osteoporotic patients. A new qualitative index of bone strength, based on finite element analysis and named bone strain index (BSI), has been recently developed from lumbar DXA scan. We present the preliminary results about the BSI ability to predict a refracture in patients with fragility fractures. A total of 143 consecutive fractured patients with primary osteoporosis (121 females) performed a spine x-ray examination for the calculation of spine deformity index (SDI) and a DXA densitometry for BMD, TBS, and BSI at basal time and in the follow-up. A refracture was considered as a one-unit increase in SDI. For each unit increase of the investigated indexes, the hazard ratio of refracture, 95% confidence interval, p value, and proportionality test p value were for BSI 1.201, 0.982−1.468, 0.074, and 0.218; for lumbar BMD 0.231, 0.028−1.877, 0.170, and 0.305; and for TBS 0.034, 0.001−2.579, 0.126, and 0.518, respectively. BSI was the index predictive of refracture nearest to statistical significance. If confirmed, it may be used for a better risk assessment of osteoporotic patients.

Artificial neural network analysis of bone quality DXA parameters response to teriparatide in fractured osteoporotic patients

Teriparatide is a bone-forming therapy for osteoporosis that increases bone quantity and texture, with uncertain action on bone geometry. No data are available regarding its influence on bone strain. To investigate teriparatide action on parameters of bone quantity and quality and on Bone Strain Index (BSI), also derived from DXA lumbar scan, based on the mathematical model finite element method. Forty osteoporotic patients with fractures were studied before and after two years of daily subcutaneous 20 mcg of teriparatide with dual X-ray photon absorptiometry to assess bone mineral density (BMD), hip structural analysis (HSA), trabecular bone score (TBS), BSI. Spine deformity index (SDI) was calculated from spine X-ray. Shapiro-Wilks, Wilcoxon and Student's t test were used for classical statistical analysis. Auto Contractive Map was used for Artificial Neural Network Analysis (ANNs). In the entire population, the ameliorations after therapy regarded BSI (-13.9%), TBS (5.08%), BMD (8.36%). HSA parameters of femoral shaft showed a worsening. Dividing patients into responders (BMD increase >10%) and non-responders, the first presented TBS and BSI ameliorations (11.87% and -25.46%, respectively). Non-responders presented an amelioration of BSI only, but less than in the other subgroup (-6.57%). ANNs maps reflect the mentioned bone quality improvements. Teriparatide appears to ameliorate not only BMD and TBS, but also BSI, suggesting an increase of bone strength that may explain the known reduction in fracture risk, not simply justified by BMD increase. BSI appears to be a sensitive index of TPD effect. ANNs appears to be a valid tool to investigate complex clinical systems.

A new finite element based parameter to predict bone fracture

Dual Energy X-Ray Absorptiometry (DXA) is currently the most widely adopted non-invasive clinical technique to assess bone mineral density and bone mineral content in human research and represents the primary tool for the diagnosis of osteoporosis. DXA measures areal bone mineral density, BMD, which does not account for the three-dimensional structure of the vertebrae and for the distribution of bone mass. The result is that longitudinal DXA can only predict about 70% of vertebral fractures. This study proposes a complementary tool, based on Finite Element (FE) models, to improve the DXA accuracy. Bone is simulated as elastic and inhomogeneous material, with stiffness distribution derived from DXA greyscale images of density. The numerical procedure simulates a compressive load on each vertebra to evaluate the local minimum principal strain values. From these values, both the local average and the maximum strains are computed over the cross sections and along the height of the analysed bone region, to provide a parameter, named Strain Index of Bone (SIB), which could be considered as a bone fragility index. The procedure is initially validated on 33 cylindrical trabecular bone samples obtained from porcine lumbar vertebrae, experimentally tested under static compressive loading. Comparing the experimental mechanical parameters with the SIB, we could find a higher correlation of the ultimate stress, σULT, with the SIB values (R2adj = 0.63) than that observed with the conventional DXA-based clinical parameters, i.e. Bone Mineral Density, BMD (R2adj = 0.34) and Trabecular Bone Score, TBS (R2adj = -0.03). The paper finally presents a few case studies of numerical simulations carried out on human lumbar vertebrae. If our results are confirmed in prospective studies, SIB could be used-together with BMD and TBS-to improve the fracture risk assessment and support the clinical decision to assume specific drugs for metabolic bone diseases.

Subject-specific finite element analysis of a lumbar cage produced by electron beam melting

The aim of this study was the analysis of the mechanical behaviour of a partially porous lumbar custom-made cage by means of a subject-specific finite element analysis (FEA). The cage, made of Ti6Al4V ELI alloy, was produced via electron beam melting (EBM) process and surgically implanted in a female subject, 50 years old. The novelty of this study was the customized design of the cage and of its internal structure, which is impossible to obtain with the traditional production techniques. The 3D model of the spine was obtained from the computed tomography (CT) of the patient. Moreover, high-resolution industrial CT was also used to reconstruct a 3D model of the cage, with its real (as-produced) features, such as superficial roughness, morphology of the bulk and of the porous structure. The workflow was divided in several steps: the main finite element analyses were non-linear and quasi-static regarding: the rhombic dodecahedron (RD) unit cell of the porous structure; the device; the whole L4-L5 motion segment with the implanted cage. Stress distribution was calculated under compression load for all models. For the RD unit cell, the maximum stress appeared at the connected cross nodes, where notch effect was present. For the cage subjected to a load of 1 kN, the porous structure did not present any functional failure. For the whole biomechanical system subjected to a physiological load of 360 N, the calculated stress in the bone was smaller than its yield strength value. On the axial view, a zone with higher compressive stresses was present on the L5 vertebral body. This was due to the contact stress between the cage and the vertebra. From the comparison between FE results and the CT images of the spine, bone remodelling was supposed, with the formation of new bone. Graphical abstract Workflow showing the phases of the research.

Biomechanical Assessment of the NIOSH Lifting Equation in Asymmetric Load-Handling Activities Using a Detailed Musculoskeletal Model

OBJECTIVE

To assess adequacy of the National Institute for Occupational Safety and Health (NIOSH) Lifting Equation (NLE) in controlling lumbar spine loads below their recommended action limits during asymmetric load-handling activities using a detailed musculoskeletal model, that is, the AnyBody Modeling System.

BACKGROUND

The NIOSH committee employed simplistic biomechanical models for the calculation of the spine compressive loads with no estimates of the shear loads. It is therefore unknown whether the NLE would adequately control lumbar compression and shear loads below their recommended action limits during asymmetric load-handling activities.

METHOD

Twenty-four static stoop lifting tasks at different load asymmetry angles, heights, and horizontal distances were performed by one normal-weight (70 kg) and one obese (93 kg) individual. For each task, the recommended weight limit computed by the NLE and body segment angles measured by a video-camera system (VICON) were prescribed in the participant-specific models developed in the AnyBody Modeling System that estimated spinal loads.

RESULTS

For both individuals, the NLE adequately controlled L5-S1 loads below their recommended action limits for all activities performed in upright postures. Both individuals, however, experienced compressive and/or shear L5-S1 loads beyond the recommended action limits when lifting was performed near the floor with large load asymmetry.

CONCLUSION

The NLE failed to control spinal loads below the recommended limits during asymmetric lifting tasks performed near the floor.

APPLICATION

The NLE should be used with caution for extreme tasks involving load handling near the floor with large load asymmetry.

Computational and Experimental Fatigue Analysis of Contoured Spinal Rods

Posterior fixation with contoured rods is an established methodology for the treatment of spinal deformities. Both uniform industrial preforming and intraoperative contouring introduce tensile and compressive plastic deformations, respectively, at the concave and at the convex sides of the rod. The purpose of this study is to develop a validated numerical framework capable of predicting how the fatigue behavior of contoured spinal rods is affected by residual stresses when loaded in lordotic and kyphotic configurations. Established finite element models (FEM) describing static contouring were implemented as a preliminary simulation step and were followed by subsequent cyclical loading steps. The equivalent Sines stress distribution predicted in each configuration was compared to that in straight rods (SR) and related to the corresponding experimental number of cycles to failure. In the straight configuration, the maximum equivalent stress (441 MPa) exceeds the limit curve, as confirmed by experimental rod breakage after around 1.9 × 105 loading cycles. The stresses further increased in the lordotic configuration, where failure was reached within 2.4 × 104 cycles. The maximum equivalent stress was below the limit curve for the kyphotic configuration (640 MPa), for which a run-out of 106 cycles was reached. Microscopy inspection confirmed agreement between numerical predictions and experimental fatigue crack location. The contouring technique (uniform contouring (UC) or French bender (FB)) was not related to any statistically significant difference. Our study demonstrates the key role of residual stresses in altering the mean stress component, superposing to the tensile cyclic load, potentially explaining the higher failure rate of lordotic rods compared to kyphotic ones.

Biomechanical analysis of manual material handling movement in healthy weight and obese workers

The risk of back injury during work remains high today for manual materials handler. The purpose of this study is to identify the potential presence of compensatory strategies in obese and non-obese handlers and evaluate the impact these strategies have on trunk kinematics and kinetics. The biomechanical and ergonomic impacts in 17 obese and 20 healthy-weight handlers were evaluated. The task studied consisted in moving boxes from a conveyor to a hand trolley and back. The results show that the anthropometric characteristics of obese handlers are linked to a significant increase in peak lumbar loading during lifting and lowering of boxes. Few postural differences between the two groups were observed. These results suggest that the excess weight of an obese worker has a significant added effect on the musculoskeletal structures of the back, which exposes obese handlers to a higher risk of developing a musculoskeletal disorder during load handling.

Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion

Spinal fusion is a standard surgical treatment for patients suffering from low back pain attributed to disc degeneration. However, results are somewhat variable and unpredictable. With fusion the kinematic behaviour of the spine is altered. Fusion and/or stabilizing implants carrying considerable load and prevent rotation of the fused segments. Associated with these changes, a risk for accelerated disc degeneration at the adjacent levels to fusion has been demonstrated. However, there is yet no method to predict the effect of fusion surgery on the adjacent tissue levels, i.e. bone and disc. The aim of this study was to develop a coupled and patient-specific mechanoregulated model to predict disc generation and changes in bone density after spinal fusion and to validate the results relative to patient follow-up data. To do so, a multiscale disc mechanoregulation adaptation framework was developed and coupled with a previously developed bone remodelling algorithm. This made it possible to determine extra cellular matrix changes in the intervertebral disc and bone density changes simultaneously based on changes in loading due to fusion surgery. It was shown that for 10 cases the predicted change in bone density and degeneration grade conforms reasonable well to clinical follow-up data. This approach helps us to understand the effect of surgical intervention on the adjacent tissue remodelling. Thereby, providing the first insight for a spine surgeon as to which patient could potentially be treated successfully by spinal fusion and in which patient has a high risk for adjacent tissue changes.

Relationship between anthropometric measures and sagittal spinal curvatures in adult male handball players

Purpose. Increasing anthropometric measures bring considerable spinal loads during sports practice, which inversely affects the adaptation abilities of the spinal structures; this in turn influences the spinal curvatures. The study was conducted to explore the relationship between anthropometric measures and sagittal spinal curvatures in handball players. Method. The total of 83 male handball players were divided into 2 groups, depending on their body height: group 1 (age, 23.62 ± 2.07 years) consisted of 40 handball players with height above average, group 2 (age, 24.63 ± 2.58 years) consisted of 43 handball players with height below average. The thoracic and lumbar curvatures and trunk height were measured with the Formetric III 4D spine and posture analysis system. Results. The thoracic kyphosis of group 1 was significantly higher than that of group 2 (p = 0.038), without a significant difference in lumbar lordosis (p = 0.312), and significant difference in the coefficient of compensation between thoracic kyphosis and lumbar lordosis (p = 0.026). Group 1 showed strong positive correlation between body height and kyphotic angle (r = 0.897), and moderate positive correlation with lordosis angle (r = 0.496). In group 2, there was weak positive correlation with kyphotic angle (r = 0.381), and weak negative correlation with lumbar lordosis angle (r = -0.355). Conclusions. Increasing body height of handball players is associated with bigger kyphotic and lordotic angles. Owing to frequent sagittal asymmetric overloading of the spine during handball training, exercises that help maintain good posture and correct the thoracic kyphosis are required, especially for taller players.

Thoracolumbar spine loading associated with kinematics of the young and the elderly during activities of daily living

Excessive mechanical loading of the spine is a critical factor in vertebral fracture initiation. Most vertebral fractures develop spontaneously or due to mild trauma, as physiological loads during activities of daily living might exceed the failure load of osteoporotic vertebra. Spinal loading patterns are affected by vertebral kinematics, which differ between elderly and young individuals. In this study, the effects of age-related changes in spine kinematics on thoracolumbar spinal segmental loading during dynamic activities of daily living were investigated using combined experimental and modeling approach. Forty-four healthy volunteers were recruited into two age groups: young (N = 23, age = 27.1 ± 3.8) and elderly (N = 21, age = 70.1 ± 3.9). The spinal curvature was assessed with a skin-surface device and the kinematics of the spine and lower extremities were recorded during daily living tasks (flexion-extension and stand-sit-stand) with a motion capture system. The obtained data were used as input for a musculoskeletal model with a detailed thoracolumbar spine representation. To isolate the effect of kinematics on predicted loads, other model properties were kept constant. Inverse dynamics simulations were performed in the AnyBody Modeling System to estimate corresponding spinal loads. The maximum compressive loads predicted for the elderly motion patterns were lower than those of the young for L2/L3 and L3/L4 lumbar levels during flexion and for upper thoracic levels during stand-to-sit (T1/T2-T8/T9) and sit-to-stand (T3/T4-T6/T7). However, the maximum loads predicted for the lower thoracic levels (T9/T10-L1/L2), a common site of vertebral fractures, were similar compared to the young. Nevertheless, these loads acting on the vertebrae of reduced bone quality might contribute to a higher fracture risk for the elderly.

High stature and body mass might affect the occurrence of Schmorl’s nodes

Schmorl’s nodes are vertical herniation of intervertebral discs into the body of neighbouring vertebral endplate. Notwithstanding extensive studies, no consensus has been reached in the subject of their possible etiology. It is hypothesized that physical stress, trauma and high axial loading are the key factors in the occurrence of this pathology. The main objective of the current work is to reevaluate the relationship between stature and body mass and Schmorl’s nodes. For this purpose, skeletal samples from Lithuania (44 males and 19 females) and Poland (97 males and 60 females) were used. The study confirmed that Schmorl’s nodes are age-independent, and more frequent in males (12.63% on the superior and 19.32% on the inferior surface of vertebrae) than in females (6.23% and 12.29% respectively). Obtained results also suggest that high stature (e.g. Spearmann correlation for superior: R=0.20 p=0.017, and inferior: R=0.31 p=0.000 surface of vertebrae) and body mass (R=0.25, p=0.002 and R=0.32, p<0.001, respectively) are factors that increase the risk of Schmorl’s nodes. Authors hypothesize that the afore-mentioned body size traits alter loadings acting on intervertebral discs, and rigidity of the spine.