Statistics in experimental studies on the human spine: Theoretical basics and review of applications

Assessing site-specificity of the biomechanical properties of hamstring aponeuroses using MyotonPRO: A cadaveric study

BACKGROUND

Hamstring muscles are the most frequently reported sites of muscle strain injuries, especially near the bi-articular muscles' myotendinous junction, where aponeurosis provides a connective tissue network linking muscle fibers to the tendon. This study aimed to investigate the reliability and site-specific differences of hamstring aponeuroses under different conditions (formalin and urea) using MyotonPRO.

METHODS

Eight hamstring muscle groups were dissected from four human cadavers (two males and two females) aged 83-93 years. Measurements of the mechanical properties of the aponeuroses from the superficial and deep regions of biceps femoris long head, semitendinosus, and semimembranosus (after formalin solution immersion) were done using MyotonPRO (intra-rater reliability was examined within a 24-h interval), following which the hamstring aponeuroses were measured using a similar procedure after urea solution immersion.

FINDINGS

Test-retest (intra-rater) results revealed that the MyotonPRO measurement of tone, stiffness, relaxation, and creep of cadaveric aponeuroses presented good to excellent reliability (ICC: 0.86 to 0.98). There were no significant differences in tone, stiffness, elasticity, relaxation, and creep among the six sites of hamstring aponeuroses under both formalin and urea conditions. Significant differences between formalin and urea conditions were found in the tone, stiffness, relaxation, and creep of hamstring aponeuroses (P < 0.05).

INTERPRETATION

These results suggested that the biomechanical properties of hamstring aponeuroses showed homogeneity between the sites using MyotonPRO. Urea solution could potentially neutralize the effect of formalin on the biomechanical properties of cadaveric muscle-aponeurosis-tendon units. The present findings might influence the design of subsequent cadaveric studies on hamstring muscle strains.

Biomechanical Assessment of Fracture Loads and Patterns of the Odontoid Process

STUDY DESIGN

Laboratory study.

OBJECTIVE

This study aimed to investigate the biomechanical competence and fracture characteristics of the odontoid process.

SUMMARY OF BACKGROUND DATA

Odontoid fractures of the second cervical vertebra (C2) represent the most common spine fracture type in the elderly. However, very little is known about the underlying biomechanical fracture mechanisms.

MATERIALS AND METHODS

A total of 42 C2 human anatomic specimens were scanned via computed tomography, divided in six groups, and subjected to combined quasistatic loading at -15°, 0°, and 15° in sagittal plane and -50° and 0° in transverse plane until fracturing. Bone mineral density (BMD), height, fusion state of the ossification centers, stiffness, yield load, and ultimate load were assessed.

RESULTS

While lowest values for stiffness, yield load, and ultimate load were observed at load inclination of 15° in sagittal plane, no statistically significant differences were observed between the study groups ( P ≥0.235). BMD correlated positively with yield load ( r2 =0.350, P <0.001) and ultimate load ( r2 =0.955, P <0.001) but not with stiffness ( r2 =0.082, P =0.07). The specimens with clearly distinguishable fusion of the ossification centers revealed less data scattering of the biomechanical outcomes.

CONCLUSION

Load direction plays a subordinate role in traumatic fractures of the odontoid process. BMD was associated with significant correlation to the biomechanical outcomes. Thus, odontoid fractures appear to result from of an interaction between the load magnitude and bone quality.

Motion capture evaluation of sagittal spino-pelvic biomechanics after lumbar spinal fusion

PURPOSE

The spine and pelvis coexist as a dynamic linked system in which spinal and pelvic parameters are correlated. Investigation of this system can inform the understanding and treatment of spinal deformity. Here, we demonstrate the use of motion capture technology to measure spine biomechanical parameters using a novel testing apparatus.

METHODS

Three complete cadaveric spines with skull and pelvis were mounted into a biomechanical testing apparatus. Each lumbar vertebra was monitored by motion capture cameras as the spines underwent maximal anterior and posterior pelvic tilts about two sagittal axes at a controlled speed and applied force. These axes were defined as the sacral axis which passes transversely through the ilium and S1, and the acetabular axis which passes transversely through both acetabula. The experiments were repeated after L4-L5 fusion, and then, after both L4-L5 and T12-S1 fusion with pedicle screw instrumentation. Data were collected for total range of motion and for coupled translation at each functional spinal unit (FSU).

RESULTS

Total range of motion and coupled translation within functional spinal units (FSUs) was decreased after spinal fusion. The displacement of each individual FSU was captured and summarized along with the observed patterns under each experimental condition.

CONCLUSION

Lumbar fusion decreases spinal motion in the sagittal plane in both overall ROM and individual coupled translations of lumbar vertebrae. This was demonstrated using motion capture technology which is useful for quantifying the translations of individual FSUs in a multisegmental spinal model.

Accounting for Biomechanical Measures from Musculoskeletal Simulation of Upright Posture Does Not Enhance the Prediction of Curve Progression in Adolescent Idiopathic Scoliosis

A major clinical challenge in adolescent idiopathic scoliosis (AIS) is the difficulty of predicting curve progression at initial presentation. The early detection of progressive curves can offer the opportunity to better target effective non-operative treatments, reducing the need for surgery and the risks of related complications. Predictive models for the detection of scoliosis progression in subjects before growth spurt have been developed. These models accounted for geometrical parameters of the global spine and local descriptors of the scoliotic curve, but neglected contributions from biomechanical measurements such as trunk muscle activation and intervertebral loading, which could provide advantageous information. The present study exploits a musculoskeletal model of the thoracolumbar spine, developed in AnyBody software and adapted and validated for the subject-specific characterization of mild scoliosis. A dataset of 100 AIS subjects with mild scoliosis and in pre-pubertal age at first examination, and recognized as stable (60) or progressive (40) after at least 6-months follow-up period was exploited. Anthropometrical data and geometrical parameters of the spine at first examination, as well as biomechanical parameters from musculoskeletal simulation replicating relaxed upright posture were accounted for as predictors of the scoliosis progression. Predicted height and weight were used for model scaling because not available in the original dataset. Robust procedure for obtaining such parameters from radiographic images was developed by exploiting a comparable dataset with real values. Six predictive modelling approaches based on different algorithms for the binary classification of stable and progressive cases were compared. The best fitting approaches were exploited to evaluate the effect of accounting for the biomechanical parameters on the prediction of scoliosis progression. The performance of two sets of predictors was compared: accounting for anthropometrical and geometrical parameters only; considering in addition the biomechanical ones. Median accuracy of the best fitting algorithms ranged from 0.76 to 0.78. No differences were found in the classification performance by including or neglecting the biomechanical parameters. Median sensitivity was 0.75, and that of specificity ranged from 0.75 to 0.83. In conclusion, accounting for biomechanical measures did not enhance the prediction of curve progression, thus not supporting a potential clinical application at this stage.

Overdrilling increases the risk of screw perforation in locked plating of complex proximal humeral fractures - A biomechanical cadaveric study

Locked plating of proximal humerus fractures (PHF) is associated with high failure rates (15-37%). Secondary screw perforation is a prominent mode of failure for PHF and typically requires reoperation. The anatomical fracture reduction is an essential factor to prevent fixation failure. However, recent studies indicate that the risk of secondary screw perforation may increase if the articular surface is perforated during predrilling of the screw boreholes (overdrilling). This study aimed to determine whether overdrilling increases the risk of secondary screw perforation in unstable PHF. Nine pairs of human cadaveric proximal humeri were osteotomized to simulate a malreduced and highly unstable 3-part fracture (AO/OTA 11 B1.1), followed by their assignment to two study groups for overdrilling or accurate predrilling in paired design, and fixation with a locking plate. Overdrilling was defined by drilling the calcar screw's boreholes through the articular surface. All humeri were cyclically loaded to screw perforation failure. Number of cycles to initial screw loosening and final perforation failure were analysed. The accurately predrilled group revealed a significantly higher number of cycles to both initial screw loosening (p < 0.01) and final screw perforation failure (p = 0.02), compared to the overdrilled one. This is the first study reporting that drilling to the correct depth significantly increases endurance until screw perforation failure during cyclic loading after locked plating in a highly unstable PHF model. Prevention of overdrilling the boreholes could help reduce failure rates of locked plating. Future work should investigate the prevalence and consequences of overdrilling in clinics.

The Effect of Degeneration on Internal Strains and the Mechanism of Failure in Human Intervertebral Discs Analyzed Using Digital Volume Correlation (DVC) and Ultra-High Field MRI

The intervertebral disc (IVD) plays a main role in absorbing and transmitting loads within the spinal column. Degeneration alters the structural integrity of the IVDs and causes pain, especially in the lumbar region. The objective of this study was to investigate non-invasively the effect of degeneration on human 3D lumbar IVD strains (n = 8) and the mechanism of spinal failure (n = 10) under pure axial compression using digital volume correlation (DVC) and 9.4 Tesla magnetic resonance imaging (MRI). Degenerate IVDs had higher (p < 0.05) axial strains (58% higher), maximum 3D compressive strains (43% higher), and maximum 3D shear strains (41% higher), in comparison to the non-degenerate IVDs, particularly in the lateral and posterior annulus. In both degenerate and non-degenerate IVDs, peak tensile and shear strains were observed close to the endplates. Inward bulging of the inner annulus was observed in all degenerate IVDs causing an increase in the AF compressive, tensile, and shear strains at the site of inward bulge, which may predispose it to circumferential tears (delamination). The endplate is the spine's “weak link” in pure axial compression, and the mechanism of human vertebral fracture is associated with disc degeneration. In non-degenerate IVDs the locations of failure were close to the endplate centroid, whereas in degenerate IVDs they were in peripheral regions. These findings advance the state of knowledge on mechanical changes during degeneration of the IVD, which help reduce the risk of injury, optimize treatments, and improve spinal implant designs. Additionally, these new data can be used to validate computational models.

Assessment of trunk muscle activation and intervertebral load in adolescent idiopathic scoliosis by musculoskeletal modelling approach

Adolescent idiopathic scoliosis (AIS) is a three-dimensional deformity of the spine, the aetiology and pathogenesis of which are poorly understood. Unfortunately, biomechanical data describing trunk muscle activation and intervertebral load, which can contribute to understanding the pathomechanics of the AIS spine, cannot be measured in vivo due to the invasiveness of the procedures. The present study provides the biomechanical characterization of the spinal loads in scoliotic subjects by exploiting musculoskeletal modelling approach, allowing for calculating biomechanical measures in an assigned posture. A spine model with articulated ribcage previously developed in AnyBody software was applied. The predicted outcomes were evaluated in the upright posture, depending on scoliosis severity and curve type, in a population of 132 scoliotic subjects with mild, moderate, and severe scoliosis. Radiographic-based three dimensional reconstruction of vertebral orientations and scaling of body segments and trunk muscle cross-section area guaranteed geometrical subject-specificity. Validation analysis supporting the application of the model was performed. Trunk muscles were found more activated in the convex side of the scoliotic curve, in agreement with reference in vivo measurements, with progressive increase with scoliosis severity. The intervertebral lateral shear was found positively correlated with the severity of the scoliosis, demonstrating that the transferred load is not a priori orthogonal to vertebral endplate in the frontal plane, and thus questioning the assumption of the 'follower load' approach in case of experimental or computational study on the scoliotic spine. The study opens the way for the subject-specific characterization of scoliosis in assigned loading and motion conditions.