Aaron J. Fields

Role of Trabecular Microarchitecture in Whole‐Vertebral Body Biomechanical Behavior

The role of trabecular microarchitecture in whole‐vertebral biomechanical behavior remains unclear, and its influence may be obscured by such factors as overall bone mass, bone geometry, and the presence of the cortical shell. To address this issue, 22 human T9 vertebral bodies (11 female; 11 male; age range: 53–97 yr, 81.5 ± 9.6 yr) were scanned with μCT and analyzed for measures of trabecular microarchitecture, BMC, cross‐sectional area, and cortical thickness. Sixteen of the vertebrae were biomechanically tested to measure compressive strength. To estimate vertebral compressive stiffness with and without the cortical shell for all 22 vertebrae, two high‐resolution finite element models per specimen—one intact model and one with the shell removed—were created from the μCT scans and virtually compressed. Results indicated that BMC and the structural model index (SMI) were the individual parameters most highly associated with strength (R2 = 0.57 each). Adding microarchitecture variables to BMC in a stepwise multiple regression model improved this association (R2 = 0.85). However, the microarchitecture variables in that regression model (degree of anisotropy, bone volume fraction) differed from those when BMC was not included in the model (SMI, mean trabecular thickness), and the association was slightly weaker for the latter (R2 = 0.76). The finite element results indicated that the physical presence of the cortical shell did not alter the relationships between microarchitecture and vertebral stiffness. We conclude that trabecular microarchitecture is associated with whole‐vertebral biomechanical behavior and that the role of microarchitecture is mediated by BMC but not by the cortical shell.

Influence of Vertical Trabeculae on the Compressive Strength of the Human Vertebra

Vertebral strength, a key etiologic factor of osteoporotic fracture, may be affected by the relative amount of vertically oriented trabeculae. To better understand this issue, we performed experimental compression testing, high‐resolution micro–computed tomography (µCT), and micro–finite‐element analysis on 16 elderly human thoracic ninth (T9) whole vertebral bodies (ages 77.5 ± 10.1 years). Individual trabeculae segmentation of the µCT images was used to classify the trabeculae by their orientation. We found that the bone volume fraction (BV/TV) of just the vertical trabeculae accounted for substantially more of the observed variation in measured vertebral strength than did the bone volume fraction of all trabeculae (r2 = 0.83 versus 0.59, p < .005). The bone volume fraction of the oblique or horizontal trabeculae was not associated with vertebral strength. Finite‐element analysis indicated that removal of the cortical shell did not appreciably alter these trends; it also revealed that the major load paths occur through parallel columns of vertically oriented bone. Taken together, these findings suggest that variation in vertebral strength across individuals is due primarily to variations in the bone volume fraction of vertical trabeculae. The vertical tissue fraction, a new bone quality parameter that we introduced to reflect these findings, was both a significant predictor of vertebral strength alone (r2 = 0.81) and after accounting for variations in total bone volume fraction in multiple regression (total R2 = 0.93). We conclude that the vertical tissue fraction is a potentially powerful microarchitectural determinant of vertebral strength.

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.

Innervation of pathologies in the lumbar vertebral end plate and intervertebral disc.

BACKGROUND CONTEXT Magnetic resonance imaging (MRI) has limited diagnostic value for chronic low back pain because of the unclear relationship between any anatomic abnormalities on MRI and pain reported by the patient. Assessing the innervation of end plate and disc pathologies-and determining the relationship between these pathologies and any abnormalities seen on MRI-could clarify the sources of back pain and help identify abnormalities with enhanced diagnostic value. PURPOSE To quantify innervation in the vertebral end plate and intervertebral disc and to relate variation in innervation to the presence of pathologic features observed by histology and conventional MRI. STUDY DESIGN/SETTING A cross-sectional histology and imaging study of vertebral end plates and intervertebral discs harvested from human cadaver spines. METHODS We collected 92 end plates and 46 intervertebral discs from seven cadaver spines (ages 51-67 years). Before dissection, the spines were scanned with MRI to grade for Modic changes and high-intensity zones (HIZ). Standard immunohistochemical techniques were used to localize the general nerve marker protein gene product 9.5. We quantified innervation in the following pathologies: fibrovascular end-plate marrow, fatty end-plate marrow, end-plate defects, and annular tears. RESULTS Nerves were present in the majority of end plates with fibrovascular marrow, fatty marrow, and defects. Nerve density was significantly higher in fibrovascular end-plate marrow than in normal end-plate marrow (p<.001). Of the end plates with fibrovascular and fatty marrow, less than 40% were Modic on MRI. Innervated marrow pathologies collocated with more than 75% of the end plate defects; hence, innervation was significantly higher in end plate defects than in normal end plates (p<.0001). In the disc, nerves were observed in only 35% of the annular tears; in particular, innervation in radial tears tended to be higher than in normal discs (p=.07). Of the discs with radial tears, less than 13% had HIZ on T2 MRI. Innervation was significantly less in radial tears than in fibrovascular end-plate marrow (p=.05) and end-plate defects (p=.02). CONCLUSIONS These findings indicate that vertebral end-plate pathologies are more innervated than intervertebral disc pathologies and that many innervated end-plate pathologies are not detectable on MRI. Taken together, these findings suggest that improved visualization of end-plate pathologies could enhance the diagnostic value of MRI for chronic low back pain.

Pathobiology of Modic changes

PurposeLow back pain (LBP) is the most disabling condition worldwide. Although LBP relates to different spinal pathologies, vertebral bone marrow lesions visualized as Modic changes on MRI have a high specificity for discogenic LBP. This review summarizes the pathobiology of Modic changes and suggests a disease model.MethodsNon-systematic literature review.ResultsChemical and mechanical stimulation of nociceptors adjacent to damaged endplates are likely a source of pain. Modic changes are adjacent to a degenerated intervertebral disc and have three generally interconvertible types suggesting that the different Modic change types represent different stages of the same pathological process, which is characterized by inflammation, high bone turnover, and fibrosis. A disease model is suggested where disc/endplate damage and the persistence of an inflammatory stimulus (i.e., occult discitis or autoimmune response against disc material) create predisposing conditions. The risk to develop Modic changes likely depends on the inflammatory potential of the disc and the capacity of the bone marrow to respond to it. Bone marrow lesions in osteoarthritic knee joints share many characteristics with Modic changes adjacent to degenerated discs and suggest that damage-associated molecular patterns and marrow fat metabolism are important pathogenetic factors. There is no consensus on the ideal therapy. Non-surgical treatment approaches including intradiscal steroid injections, anti-TNF-α antibody, antibiotics, and bisphosphonates have some demonstrated efficacy in mostly non-replicated clinical studies in reducing Modic changes in the short term, but with unknown long-term benefits. New diagnostic tools and animal models are required to improve painful Modic change identification and classification, and to clarify the pathogenesis.ConclusionModic changes are likely to be more than just a coincidental imaging finding in LBP patients and rather represent an underlying pathology that should be a target for therapy.

Trabecular plates and rods determine elastic modulus and yield strength of human trabecular bone

The microstructure of trabecular bone is usually perceived as a collection of plate-like and rod-like trabeculae, which can be determined from the emerging high-resolution skeletal imaging modalities such as micro-computed tomography (μCT) or clinical high-resolution peripheral quantitative CT (HR-pQCT) using the individual trabecula segmentation (ITS) technique. It has been shown that the ITS-based plate and rod parameters are highly correlated with elastic modulus and yield strength of human trabecular bone. In the current study, plate-rod (PR) finite element (FE) models were constructed completely based on ITS-identified individual trabecular plates and rods. We hypothesized that PR FE can accurately and efficiently predict elastic modulus and yield strength of human trabecular bone. Human trabecular bone cores from proximal tibia (PT), femoral neck (FN) and greater trochanter (GT) were scanned by μCT. Specimen-specific ITS-based PR FE models were generated for each μCT image and corresponding voxel-based FE models were also generated in comparison. Both types of specimen-specific models were subjected to nonlinear FE analysis to predict the apparent elastic modulus and yield strength using the same trabecular bone tissue properties. Then, mechanical tests were performed to experimentally measure the apparent modulus and yield strength. Strong linear correlations for both elastic modulus (r(2) = 0.97) and yield strength (r(2) = 0.96) were found between the PR FE model predictions and experimental measures, suggesting that trabecular plate and rod morphology adequately captures three-dimensional (3D) microarchitecture of human trabecular bone. In addition, the PR FE model predictions in both elastic modulus and yield strength were highly correlated with the voxel-based FE models (r(2) = 0.99, r(2) = 0.98, respectively), resulted from the original 3D images without the PR segmentation. In conclusion, the ITS-based PR models predicted accurately both elastic modulus and yield strength determined experimentally across three distinct anatomic sites. Trabecular plates and rods accurately determine elastic modulus and yield strength of human trabecular bone.

Dependence of mechanical properties of trabecular bone on plate–rod microstructure determined by individual trabecula segmentation (ITS)

Individual trabecula segmentation (ITS) technique can decompose the trabecular bone network into individual trabecular plates and rods and is capable of quantifying the plate/rod-related microstructural characteristics of trabecular bone. This novel technique has been shown to be able to provide in-depth insights into micromechanics and failure mechanisms of human trabecular bone, as well as to distinguish the fracture status independent of area bone mineral density in clinical applications. However, the plate/rod microstructural parameters from ITS have never been correlated to experimentally determined mechanical properties of human trabecular bone. In this study, on-axis cylindrical trabecular bone samples from human proximal tibia (n=22), vertebral body (n=10), and proximal femur (n=21) were harvested, prepared, scanned using micro computed-tomography (µCT), analyzed with ITS and mechanically tested. Regression analyses showed that the plate bone volume fraction (pBV/TV) and axial bone volume fraction (aBV/TV) calculated by ITS analysis correlated the best with elastic modulus (R(2)=0.96-0.97) and yield strength (R(2)=0.95-0.96). Trabecular plate-related microstructural parameters correlated highly with elastic modulus and yield strength, while most rod-related parameters were found inversely and only moderately correlated with the mechanical properties. In addition, ITS analysis also identified that trabecular bone at human femoral neck had the highest trabecular plate-related parameters while the other sites were similar with each other in terms of plate-rod microstructure.

Vertebral fragility and structural redundancy

The mechanisms of age‐related vertebral fragility remain unclear, but may be related to the degree of “structural redundancy” of the vertebra; ie, its ability to safely redistribute stress internally after local trabecular failure from an isolated mechanical overload. To better understand this issue, we performed biomechanical testing and nonlinear micro‐CT–based finite element analysis on 12 elderly human thoracic ninth vertebral bodies (age 76.9 ± 10.8 years). After experimentally overloading the vertebrae to measure strength, we used nonlinear finite element analysis to estimate the amount of failed tissue and understand the failure mechanisms. We found that the amount of failed tissue per unit bone mass decreased with decreasing bone volume fraction (r2 = 0.66, p < 0.01). Thus, for the weak vertebrae with low bone volume fraction, overall failure of the vertebra occurred after failure of just a tiny proportion of the bone tissue (<5%). This small proportion of failed tissue had two sources: the existence of fewer vertically oriented load paths to which load could be redistributed from failed trabeculae; and the vulnerability of the trabeculae in these few load paths to undergo bending‐type failure mechanisms, which further weaken the bone. Taken together, these characteristics suggest that diminished structural redundancy may be an important aspect of age‐related vertebral fragility: vertebrae with low bone volume fraction are highly susceptible to collapse because so few trabeculae are available for load redistribution if the external loads cause any trabeculae to fail.

Alterations in intervertebral disc composition, matrix homeostasis and biomechanical behavior in the UCD-T2DM rat model of type 2 diabetes.

Type 2 diabetes (T2D) adversely affects many tissues, and the greater incidence of discogenic low back pain among diabetic patients suggests that the intervertebral disc is affected too. Using a rat model of polygenic obese T2D, we demonstrate that diabetes compromises several aspects of disc composition, matrix homeostasis, and biomechanical behavior. Coccygeal motion segments were harvested from 6‐month‐old lean Sprague‐Dawley rats, obese Sprague‐Dawley rats, and diabetic obese UCD‐T2DM rats (diabetic for 69 ± 7 days). Findings indicated that diabetes but not obesity reduced disc glycosaminoglycan and water contents, and these degenerative changes correlated with increased vertebral endplate thickness and decreased endplate porosity, and with higher levels of the advanced glycation end‐product (AGE) pentosidine. Consistent with their diminished glycosaminoglycan and water contents and their higher AGE levels, discs from diabetic rats were stiffer and exhibited less creep when compressed. At the matrix level, elevated expression of hypoxia‐inducible genes and catabolic markers in the discs from diabetic rats coincided with increased oxidative stress and greater interactions between AGEs and one of their receptors (RAGE). Taken together, these findings indicate that endplate sclerosis, increased oxidative stress, and AGE/RAGE‐mediated interactions could be important factors for explaining the greater incidence of disc pathology in T2D.

Multi-scale modeling of the human vertebral body: comparison of micro-CT based high-resolution and continuum-level models.

The overall goal of this study was to assess the mechanistic fidelity of continuum-level finite element models of the vertebral body, which represent a promising tool for understanding and predicting clinical fracture risk. Two finite element (FE) models were generated from micro-CT scans of each of 13 T9 vertebral bodies--a micro-FE model at 60-micron resolution and a coarsened, continuum-level model at 0.96-mm resolution. Two previously-reported continuum-level modulus-density relationships for human vertebral bone were parametrically varied to investigate their effects on model fidelity using the micro-CT models as a gold standard. We found that the modulus-density relation, particularly that assigned to the peripheral bone, substantially altered the regression coefficients, but not the degree of correlation between continuum and micro-FE predictions of whole-vertebral stiffness. The major load paths through the vertebrae compared well between the continuum-level and micro-FE models (von-Mises distribution), but the distributions of minimum principal strain were notably different. We conclude that continuum-level models provide robust measures of whole-vertebral behavior, describe well the load transfer paths through the vertebra, but provide strain distributions that are markedly different than the volume-averaged micro-scale strains. Appreciation of these multi-scale differences should improve interpretation of results from these sorts of continuum models and may improve their clinical utility.