Mechanisms of initial endplate failure in the human vertebral body

Endplate failure occurs frequently in osteoporotic vertebral fractures and may be related to the development of high tensile strain. To determine whether the highest tensile strains in the vertebra occur in the endplates, and whether such high tensile strains are associated with the material behavior of the intervertebral disc, we used micro-CT-based finite element analysis to assess tissue-level strains in 22 elderly human vertebrae (81.5±9.6 years) that were compressed through simulated intervertebral discs. In each vertebra, we compared the highest tensile and compressive strains across the different compartments: endplates, cortical shell, and trabecular bone. The influence of Poisson-type expansion of the disc on the results was determined by compressing the vertebrae a second time in which we suppressed the Poisson expansion. We found that the highest tensile strains occurred within the endplates whereas the highest compressive strains occurred within the trabecular bone. The ratio of strain to assumed tissue-level yield strain was the highest for the endplates, indicating that the endplates had the greatest risk of initial failure. Suppressing the Poisson expansion of the disc decreased the amount of highly tensile-strained tissue in the endplates by 79.4±11.3%. These results indicate that the endplates are at the greatest risk of initial failure due to the development of high tensile strains, and that such high tensile strains are associated with the Poisson expansion of the disc. We conclude that initial failure of the vertebra is associated with high tensile strains in the endplates, which in turn are influenced by the material behavior of the disc.     

The relation between intervertebral disc bulging and annular fiber associated strains for simple and complex loading

Mechanical failure of the annulus fibrosus is mostly indicated by tears, fissures, protrusions or disc prolapses. Some of these annulus failures can be caused by a high intradiscal pressure. This has an effect on disc bulging. However, it is not fully understood how disc bulging is related to disc loading. Therefore, the aim of this study was to investigate the annular fiber strains and disc bulging under simple and complex spinal loads. A novel laser scanner was used to image surfaces of six L2-3 segments. Specimens were loaded with 500 N or 7.5 Nm in a spine tester while acquiring surface maps. Loading was applied in the three principal main directions and four combined directions. Disc bulging and tissue surface strains in annulus collagen fiber directions were computed. Two conditions were measured; intact and defect (vertebral body-disc-body units). Axial compression resulted in 2.7% fiber associated strains in intact segments and the defect increased strains up to 6.7%. Disc bulging increased from 0.7 mm to 0.87 mm. Flexion produced 7.2% fiber associated strains and 1.63 mm bulge going up to 17.5% and 2.21 mm after the defect. Highest fiber associated strains were found for the combination of axial rotation plus lateral bending with 24.6% and with a maximal bulging of 1.14 mm. It was found that there is no tight relationship between fiber associated strains and disc bulging. This was especially seen for the load combinations. Highest fiber associated strains were found to be located in small posterolateral regions. Fiber associated strains had a much higher magnitude than previously reported fiber associated strains. The results showed that combined loading is most likely to produce higher associated fiber strains compared to single axis loading.     

Stepwise reduction of functional spinal structures increase disc bulge and surface strains

Previous studies postulated that an axial compression of lumbar intervertebral discs causes a complex strain pattern on the annulus. This pattern is not fully understood, since most studies measured only the uniaxial ultimate tensile strain of the annulus. The aim of this study was to investigate surface strains and their relation to disc bulging. This work was extended to study some defects that are relevant for the intermediate process of finite element modeling. Six specimens (L2-3) with a median age of 51 years were utilized for this in vitro study. Specimens were loaded with pure moments (2.5-7.5Nm) in the principal directions. The anatomy was subsequently reduced in three steps: (1) ligamentous and bony posterior structures, (2) anterior and posterior ligaments and (3) nucleus. Measured were ranges of motion, three-dimensional disc bulging and surface strains of the outer annulus. Lateral bending showed the largest axial strains (9.7%) for intact specimens, which increased to 15.1% after the removal of posterior structures. Disc bulging was largest in flexion with 1.56mm, which increased to 2.06mm after step (1). Defect (2) caused that flexion yielded the largest axial strains with 22.6% and 2.17mm of bulging. We could also determine a constriction effect of these ligaments. Nucleotomy did not essentially increase anterior disc bulging in flexion, but inward disc bulging increased by 0.55mm, in extension. Due to the increase in the complexity of finite element models, it is difficult to obtain data from the literature for validation purposes. This study presents new data, which assist in the development of such models.     

Full field strain measurements of collagenous tissue by tracking fiber alignment through vector correlation

Full field strain measurements of biological tissue during loading are often limited to the quantification of fiduciary marker displacements on the tissue surface. These marker measurements can lack the necessary spatial resolution to characterize non-uniform deformation and may not represent the deformation of the load-bearing collagen microstructure. To overcome these potential limitations, a method was developed to track the deformation of the collagen fiber microstructure in ligament tissue. Using quantitative polarized light imaging, fiber alignment maps incorporating both direction and alignment strength at each pixel were generated during facet capsular ligament loading. A grid of virtual markers was superimposed over the tissue in the alignment maps, and the maximization of a vector correlation calculation between fiber alignment maps was used to track marker displacement. Tracking error was quantified through comparisons to the displacements of excised ligament tissue (n=3); separate studies applied uniaxial tension to isolated facet capsular ligament tissue (n=4) to evaluate tracking capabilities during large tissue deformations. The average difference between virtual marker and tissue displacements was 0.07±0.06 pixels. This error in marker location produced principal strain measurements of 1.2±1.6% when markers were spaced 4 pixels apart. During tensile tissue loading, substantial inhomogeneity was detected in the strain field using vector correlation tracking, and the location of maximum strain differed from that produced by standard tracking techniques using coarser meshes. These findings provide a method to directly measure fiber network strains using quantitative fiber alignment data, enabling a better understanding of structure–function relationships in tissues at different length scales.     

Measurement of surface deformation of soft tissue

A method is described for measuring the surface shape and deformations of soft tissue in three dimensions. The method uses close range stereophotogrammetry to record the three-dimensional locations of miniature optical targets applied to the tissue surface. This has been applied to study of human lumbar intervertebral disc. Measurements of the strain along surface annular fibers have been made under varying loads. In this case the maximum expected errors are about 0.15 mm, which corresponds to a strain of less than 1%. Preliminary findings have differed from predictions made in published mathematical models for the disc in that they show very little strain of the annulus in compression loading, but confirm axial torsional loading as liable to produce mechanical disruption of the disc annulus.     

Mechanical damage to the intervertebral disc annulus fibrosus subjected to tensile loading

Damage of the annulus fibrosus is implicated in common spinal pathologies. The objective of this study was to obtain a quantitative relationship between both the number of cycles and the magnitude of tensile strain resulting in damage to the annulus fibrosus. Four rectangular tensile specimens oriented in the circumferential direction were harvested from the outer annulus of 8 bovine caudal discs (n = 32) and subjected to one of four tensile testing protocols: (i) ultimate tensile strain (UTS) test; (ii) baseline cyclic test with 4 series of 400 cycles of baseline cyclic loading (peak strain = 20% UTS); (iii & iv) acute and fatigue damage cyclic tests consisting of 4 x 400 cycles of baseline cyclic loading with intermittent loading to 1 and 100 cycles, respectively, with peak tensile strain of 40%, 60%, and 80% UTS. Normalized peak stress for all mechanically loaded specimens was reduced from 0.89 to 0.11 of the baseline control levels, and depended on the magnitude of damaging strain and number of cycles at that damaging strain. Baseline, acute, and fatigue protocols resulted in permanent deformation of 3.5%, 6.7% and 9.6% elongation, respectively. Damage to the laminate structure of the annulus in the absence of biochemical activity in this study was assessed using histology, transmission electron microscopy, and biochemical measurements and was most likely a result of separation of annulus layers (i.e., delamination). Permanent elongation and stress reduction in the annulus may manifest in the motion segment as sub-catastrophic damage including increased neutral zone, disc bulging, and loss of nucleus pulposus pressure. The preparation of rectangular tensile strip specimens required cutting of collagen fibers and may influence absolute values of results, however, it is not expected to affect the comparisons between loading groups or dose-response reported.     

Finite element modeling of 3D human mesenchymal stem cell-seeded collagen matrices exposed to tensile strain

The use of human mesenchymal stem cells (hMSCs) in tissue engineering is attractive due to their ability to extensively self-replicate and differentiate into a multitude of cell lineages. It has been experimentally established that hMSCs are influenced by chemical and mechanical signals. However, the combined chemical and mechanical in vitro culture conditions that lead to functional tissue require greater understanding. In this study, finite element models were created to evaluate the local loading conditions on bone marrow-derived hMSCs seeded in three-dimensional collagen matrices exposed to cyclic tensile strain. Mechanical property and geometry data used in the models were obtained experimentally from a previous study in our laboratory and from mechanical testing. Eight finite element models were created to simulate three-dimensional hMSC-seeded collagen matrices exposed to different levels of cyclic tensile strain (10% and 12%), culture media (complete growth and osteogenic differentiating), and durations of culture (7 and 14 days). Through finite element analysis, it was determined that globally applied uniaxial tensile strains of 10% and 12% resulted in local strains up to 18.3% and 21.8%, respectively. Model results were also compared to experimental studies in an attempt to explain observed differences between hMSC response to 10% and 12% cyclic tensile strain.     

A biomechanical model of the lumbosacral joint during lifting activities

A biomechanical model of the lumbosacral region was constructed for the purpose of systematically studying the combined stresses and strains on the local ligaments, muscles and disc tissue during sagittal plane two-handed lifting. The model was validated in two ways. The first validation was a comparison of experimental study results with model predictions. In general predictions compared very reasonably with observed values of several authors with the exception of strain predictions on the articular ligaments. Second, a sensitivity analysis was performed over a wide range of lifting tasks. The predicted stress/strain values followed anticipated patterns and were of reasonable magnitudes. On the basis of the results of the sensitivity analysis it was concluded that typical lifting tasks can lead to excessive disc compressive forces, muscle moment generation requirements, and possibly lumbodorsal fascia strains. Conversely, annulus rupture of a healthy disc due to overstrain appears very unlikely.     

Degeneration affects the fiber reorientation of human annulus fibrosus under tensile load

The angled, lamellar structure of the annulus fibrosus is integral to its load-bearing function. Reorientation of this fiber structure with applied load may contribute to nonlinear mechanical behavior and to large increases in tensile modulus. Fiber reorientation has not yet been quantified for loaded non-degenerated and degenerated annulus fibrosus tissue. The objective of this study was to measure fiber reorientation and mechanical properties (toe- and linear-region modulus, transition strain, and Poisson's ratio) of loaded outer annulus fibrosus tissue using a new application of FFT image processing techniques. This method was validated for quantification of annulus fiber reorientation during loading in this study. We hypothesized that annulus fibrosus fibers would reorient under circumferential tensile load, and that fiber reorientation would be affine. Additionally, we hypothesized that degeneration would affect fiber reorientation, toe-region modulus and Poisson's ratio. Annulus fibrosus fibers were found to reorient toward the loading direction, and degeneration significantly decreased fiber reorientation (the fiber reorientation parameter, m(FFT)=-1.70 degrees /% strain for non-degenerated and -0.95 degrees /% strain for degenerated tissue). Toe-region modulus was significantly correlated with age (r=0.6). Paired t-tests showed no significant difference in the fiber reorientation parameter calculated experimentally with that calculated using an affine prediction. Thus, an affine prediction is a good approximation of fiber reorientation. The findings of this study add to the understanding of overall disc mechanical behavior and degeneration.     

Mechanical properties of the human lumbar anterior longitudinal ligament.

A new technique incorporating a motion analysis system and a materials testing machine was used to investigate regional differences in the tensile mechanical properties of the lumbar spine anterior longitudinal ligament (ALL). Bone-ALL-bone specimens were prepared from young human cadaveric motion segments with no disc or bony pathology. Each specimen was distracted until failure at a constant crosshead displacement rate of 2.5 mm s-1 (approximately 1.0% strain per second). Strains were evaluated from digitized video recordings of markers attached to the ALL at 12 sites along its length and width, including the ligament substance and insertions. The 'overall' strain in the ligament was calculated from the outermost pairs of markers along the ligament length. The average tensile strength, the 'overall' tensile modulus and the 'overall' strain of the ALL at failure were 27.4 MPa (S.D. 5.9), 759 MPa (S.D. 336) and 4.95% (S.D. 1.51), respectively. Large and significant variations in the strains were present along the width and length of the ALL. Peak substance strains were over twofold greater than peak strains at the ligament insertion sites, whereas across the ligament width, peak strains in the outer portion of the ligament were over 40% greater than in the central region. Failure consistently occurred in the ligament mid-substance and ultimate strains at the ligament failure site averaged 12.1% (S.D. 2.3). These results indicate that the strains are highly nonuniform in the normal ALL.     

Methodology to determine failure characteristics of planar soft tissues using a dynamic tensile test

Predicting the injury risk in automotive collisions requires accurate knowledge of human tissues, more particularly their mechanical properties under dynamic loadings. The present methodology aims to determine the failure characteristics of planar soft tissues such as skin, hollow organs and large vessel walls. This consists of a dynamic tensile test, which implies high-testing velocities close to those in automotive collisions. To proceed, I-shaped tissue samples are subjected to dynamic tensile tests using a customized tensile device based on the drop test principle. Data acquisition has especially been adapted to heterogeneous and soft biological tissues given that standard measurement systems (considered to be global) have been completed with a non-contact and full-field strain measurement (considered to be local). This local measurement technique, called the Image Correlation Method (ICM) provides an accurate strain analysis by revealing strain concentrations and avoids damaging the tissue. The methodology has first been applied to human forehead skin and can be further expanded to other planar soft tissues. The failure characteristics for the skin in terms of ultimate stress are 3 MPa +/- 1.5 MPa. The ultimate global longitudinal strains are equal to 9.5%+/-1.9% (Green-Lagrange strain), which contrasts with the ultimate local longitudinal strain values of 24.0%+/-5.3% (Green-Lagrange strain). This difference is a consequence of the tissue heterogeneity, clearly illustrated by the heterogeneous distribution of the local strain field. All data will assist in developing the tissue constitutive law that will be implemented in finite element models.     

Strain transfer in the annulus fibrosus under applied flexion

A detailed understanding of the anatomical and mechanical environment in the intervertebral disc at the scale of the cell is necessary for the design of tissue engineering repair strategies and to elucidate the role of mechanical factors in pathology. The objective of this study was to measure and compare the macroscale to microscale strains in the outer annulus fibrosus in various cellular regions of intact discs over a range of applied flexion. Macroscale strains were measured on the annulus fibrosus surface, and contrasted to in situ microscale strains using novel confocal microscopy techniques for dual labeling of the cell and the extracellular matrix. Fiber oriented surface strains were significantly higher than in situ fiber strains, which implies a mechanism of load redistribution that minimizes strain along the fibers. Non-uniformity of the strains and matrix distortion occurred immediately and most interestingly varied little with increase in flexion (3–16°), suggesting that inter-fiber shear is important in the initial stages of strain redistribution. Fiber oriented intercellular strains were significantly larger and compressive compared to in situ strains in other regions of the extracellular matrix indicating that the mechanical environment in this region may be unique. Further examination of the structural morphology in this pericellular region is needed to fully understand the pathway of strain transfer from the tissue to the cell. This study provides new knowledge on the complex in situ micro-mechanical environment of the annulus fibrosus that is essential to understanding the mechanobiological behavior of this tissue.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

A method to obtain surface strains of soft tissues using a laser scanning device

A three-dimensional laser scanning device was developed allowing surface digitization of musculoskeletal and soft tissue structures under different loads. Image-processing algorithms were formulated for image registration. These were used to determine displacement mapping and then surface strains. Various validation experiments were performed. Accuracy was obtained on a test cylinder after rigid rotation and on a silicon cylinder compressed in four loading steps. The system accuracy (including the scanning and the data evaluation) was +/-0.10% strain in vertical and +/-0.16% strain in shear and circumferential direction for the rigid rotation exhibiting the zero-strain situation. Silicon cylinder compression showed that the accuracy was best for small strains, whereas strains >5% evoked a slight underestimation increasing further with higher strains (error of 0.54% for 7.22% vertical strain). It was possible to increase the accuracy by performing the strain measurements via sub-steps. This had a remaining error of 0.41% for 7.22% vertical strain. A further experiment was carried out in order to acquire the surface strain of a human lumbar intervertebral disc while it was forced to flexion and extension. This study introduced a laser-based scanning method to obtain soft tissue surface strains. It is important to know the strain distribution of musculoskeletal structures and soft tissues. This could help to better understand the mechanical loading of biological structures e.g. the processes in fracture healing. These data could also be used to assist in the validation process for finite-element models.     

Determination of the tensile mechanical properties of the segmented mitral valve annulus

The mitral valve annulus is a complex and irregular component of the mitral valve apparatus, serving both a structural and sphincteric role. We have sought to determine the mechanical properties of the mitral valve annulus segmentally. Twenty porcine hearts were dissected to isolate the annulus. The annulus was segmented into four sections: anterior, posterior, and left and right commissural sections. Ten of these were tensile tested to failure as control samples. The remaining ten were digested in order to fully isolate the annulus from the myocardium, and subsequently tensile tested to failure. Histological samples of each segment were analysed to determine collagen/annular content. Whole segments of muscular annulus were tensile tested to failure; the stress and strain at failure and location of failure were determined in these larger specimens. Our results demonstrated that the anterior annulus is stiffer than the posterior segment by a factor of approximately 27 at a 2% strain level, and approximately 13 at a 6% strain. There is a trend in the results that identifies that the muscular annulus is stiffest at the right commissural segment, while the posterior segment tends to be the least stiff. The stiffness of the samples can be correlated with the area associated with the dense collagen annulus using histological analysis. Finally, the weakest section of the mitral valve annulus was identified as the intersection of the right commissural segment and the posterior segment.     

Mechanisms for mechanical damage in the intervertebral disc annulus fibrosus

Intervertebral disc degeneration results in disorganization of the laminate structure of the annulus that may arise from mechanical microfailure. Failure mechanisms in the annulus were investigated using composite lamination theory and other analyses to calculate stresses in annulus layers, interlaminar shear stress, and the region of stress concentration around a fiber break. Scanning electron microscopy (SEM) was used to evaluate failure patterns in the annulus and evaluate novel structural features of the disc tissue. Stress concentrations in the annulus due to an isolated fiber break were localized to approximately 5 microm away from the break, and only considered a likely cause of annulus fibrosus failure (i.e., radial tears in the annulus) under extreme loading conditions or when collagen damage occurs over a relatively large region. Interlaminar shear stresses were calculated to be relatively large, to increase with layer thickness (as reported with degeneration), and were considered to be associated with propagation of circumferential tears in the annulus. SEM analysis of intervertebral disc annulus fibrosus tissue demonstrated a clear laminate structure, delamination, matrix cracking, and fiber failure. Novel structural features noted with SEM also included the presence of small tubules that appear to run along the length of collagen fibers in the annulus and a distinct collagenous structure representative of a pericellular matrix in the nucleus region.     

Live free or die: stretch-induced apoptosis occurs when adaptive reorientation of annulus fibrosus cells is restricted

High matrix strains in the intervertebral disc occur during physiological motions and are amplified around structural defects in the annulus fibrosus (AF). It remains unknown if large matrix strains in the human AF result in localized cell death. This study investigated strain amplitudes and substrate conditions where AF cells were vulnerable to stretch-induced apoptosis. Human degenerated AF cells were subjected to 1 Hz-cyclic tensile strains for 24h on uniformly collagen coated substrates and on substrates with 40 μm stripes of collagen that restricted cellular reorientation. AF cells were capable of responding to stretch (stress fibers and focal adhesions aligned perpendicular to the direction of stretch), but were vulnerable to stretch-induced apoptosis when cytoskeletal reorientation was restricted, as could occur in degenerated states due to fibrosis and crosslink accumulation and at areas where high strains occur (around structural defects, delaminations, and herniations).     

Application of a fiber-reinforced continuum theory to multiple deformations of the annulus fibrosus.

Accurate tissue stress predictions for the annulus fibrosus are essential for understanding the factors that cause or contribute to disc degeneration and mechanical failure. Current computational models used to predict in vivo disc stresses utilize material laws for annular tissue that are not rigorously validated against experimental data. Consequently, predictions of disc stress resulting from physical activities may be inaccurate and therefore unreliable as a basis for defining mechanical-biologic injury criteria. To address this need we present a model for the annulus as an isotropic ground substance reinforced with two families of collagen fibers, and an approach for determining the material constants by simultaneous consideration of multiple experimental data sets. Two strain energy functions for the annulus are proposed and used in the theory to derive the constitutive equations relating the stress to pure stretch deformations. These equations are applied to four distinct experimental protocols and the material constants are determined from a simultaneous, nonlinear regression analysis. Good agreement between theory and experiment is achieved when the invariants are included within multiple, separate exponentials in the strain energy function.