Regional variations in trabecular architecture of the lumbar vertebra: associations with age, disc degeneration and disc space narrowing

Characterization of mechanical stiffness using additive manufacturing and finite element analysis: potential tool for bone health assessment

BACKGROUND

Bone health and fracture risk are known to be correlated with stiffness. Both micro-finite element analysis (μFEA) and mechanical testing of additive manufactured phantoms are useful approaches for estimating mechanical properties of trabecular bone-like structures. However, it is unclear if measurements from the two approaches are consistent. The purpose of this work is to evaluate the agreement between stiffness measurements obtained from mechanical testing of additive manufactured trabecular bone phantoms and μFEA modeling. Agreement between the two methods would suggest 3D printing is a viable method for validation of μFEA modeling.

METHODS

A set of 20 lumbar vertebrae regions of interests were segmented and the corresponding trabecular bone phantoms were produced using selective laser sintering. The phantoms were mechanically tested in uniaxial compression to derive their stiffness values. The stiffness values were also derived from in silico simulation, where linear elastic μFEA was applied to simulate the same compression and boundary conditions. Bland-Altman analysis was used to evaluate agreement between the mechanical testing and μFEA simulation values. Additionally, we evaluated the fidelity of the 3D printed phantoms as well as the repeatability of the 3D printing and mechanical testing process.

RESULTS

We observed good agreement between the mechanically tested stiffness and μFEA stiffness, with R2 of 0.84 and normalized root mean square deviation of 8.1%. We demonstrate that the overall trabecular bone structures are printed in high fidelity (Dice score of 0.97 (95% CI, [0.96,0.98]) and that mechanical testing is repeatable (coefficient of variation less than 5% for stiffness values from testing of duplicated phantoms). However, we noticed some defects in the resin microstructure of the 3D printed phantoms, which may account for the discrepancy between the stiffness values from simulation and mechanical testing.

CONCLUSION

Overall, the level of agreement achieved between the mechanical stiffness and μFEA indicates that our μFEA methods may be acceptable for assessing bone mechanics of complex trabecular structures as part of an analysis of overall bone health.

Potential Role for Stem Cell Regenerative Therapy as a Treatment for Degenerative Disc Disease and Low Back Pain: A Systematic Review

Back pain is the single leading cause of disability worldwide. Despite the prevalence and morbidity of lower back pain, we still lack a gold-standard treatment that restores the physiological function of degenerated intervertebral discs. Recently, stem cells have emerged as a promising strategy for regenerative therapy for degenerative disc disease. In this study, we review the etiology, pathogenesis, and developing treatment strategies for disc degeneration in low back pain with a focus on regenerative stem cell therapies. A systematic search of PubMed/MEDLINE/Embase/Clinical Trials.gov databases was conducted for all human subject abstracts or studies. There was a total of 10 abstracts and 11 clinical studies (1 RCT) that met the inclusion criteria. The molecular mechanism, approach, and progress of the different stem cell strategies in all studies are discussed, including allogenic bone marrow, allogenic discogenic cells, autologous bone marrow, adipose mesenchymal stem cells (MSCs), human umbilical cord MSC, adult juvenile chondrocytes, autologous disc derived chondrocytes, and withdrawn studies. Clinical success with animal model studies is promising; however, the clinical outcomes of stem cell regenerative therapy remain poorly understood. In this systematic review, we found no evidence to support its use in humans. Further studies on efficacy, safety, and optimal patient selection will establish whether this becomes a viable, non-invasive therapeutic option for back pain.

Biomechanical evaluation of position and bicortical fixation of anterior lateral vertebral screws in a porcine model

Although an anterior approach with anterior lateral screw fixation has been developed for stabilizing the thoracolumbar spine clinically, screw loosening still occurs. In this novel in vitro study, we attempted to elucidate the optimal screw position in the lateral lumbar vertebra and the effect of bicortical fixation. A total of 72 fresh-frozen lumbar vertebrae from L1-6 were harvested from 12 mature pigs and randomly assigned to two modalities: bicortical fixation (n = 36) and unicortical fixation (n = 36). Six groups of screw positions in the lateral vertebral body in each modality were designated as central-anterior, central-middle, central-posterior, lower-anterior, lower-middle, and lower- posterior; 6 specimens were used in each group. The correlations between screw fixation modalities, screw positions and axial pullout strength were analyzed. An appropriate screw trajectory and insertional depth were confirmed using axial and sagittal X-ray imaging prior to pullout testing. In both bicortical and unicortical fixation modalities, the screw pullout force was significantly higher in the posterior or middle position than in the anterior position (p < 0.05), and there was no significant differences between the central and lower positions. The maximal pullout forces from the same screw positions in unicortical fixation modalities were all significantly lower, decreases that ranged from 32.7 to 74%, than those in bicortical fixation modalities. Our study using porcine vertebrae showed that screws in the middle or posterior position of the lateral vertebral body had a higher pullout performance than those in the anterior position. Posteriorly positioned lateral vertebral screws with unicortical fixation provided better stability than anteriorly positioned screws with bicortical fixation.

Polg mtDNA mutator mice reveal limited involvement of vertebral bone loss in premature aging-related thoracolumbar hyperkyphosis

Background

Age-related hyperkyphosis is multifactorial and involves alterations of vertebral bone, intervertebral discs (IVD) and paraspinal muscles. The relative contribution of these tissues remains unclear. Here, we compared differences in vertebral bone microarchitecture and IVD thickness between prematurely aging mice with spinal hyperkyphosis and wild type littermates.

Methods

Thoracolumbar vertebral columns were dissected from homozygous Polg ^D257A and age-matched wild type littermates. Micro-computed tomography was performed to quantify cortical and trabecular bone parameters at anterior and posterior portions of T8-L4 vertebrae. In addition, vertebral shape, transaxial facet joint orientation and IVD thickness were quantified. Differences in anterior/posterior ratios between genotypes were compared by Student's t-test and association between vertebral bone and IVD parameters was investigated using Pearson correlation analysis.

Results

Hyperkyphotic homozygous mice displayed generalized osteopenia that was more pronounced at the posterior compared with anterior portion of thoracolumbar vertebrae. An increase in the anterior/posterior ratio of trabecular bone parameters was revealed at the thoracolumbar junction (T13-L1). Polg ^D257A displayed diffuse loss of cortical bone thickness, yet anterior/posterior ratios were unchanged. Despite generalized and regional bone loss, vertebral shape was unaffected. PolG ^D257A mice showed a 10-20 % reduction of IVD thickness at both thoracic and lumbar levels, with only minimal histopathological changes. IVD thickness was negatively correlated with anterior/posterior ratios of trabecular bone parameters, as well as with more coronally oriented facet joints, but negatively correlated with the anterior/posterior ratio of cortical bone thickness.

Conclusions

Aging-induced regional changes of vertebral trabecular and cortical bone did not lead to altered vertebral shape in Polg ^D257A mice but may indirectly cause hyperkyphosis through reduction of IVD thickness. These findings suggest a limited role for aging-induced bone loss in spinal hyperkyphosis and warrants further research on the involvement of paraspinal muscle degeneration.

Fundamentals of Intervertebral Disc Degeneration

Lumbar disc degeneration is one of the leading causes of chronic low back pain. The degenerative cascade is often initiated by an imbalance between catabolic and anabolic processes in the intervertebral discs. As a consequence of extracellular matrix degradation, neoinnervation and neovascularization take place. Ultimately, this degenerative process results in disc bulging and loss of nucleus pulposus and water content and subsequent loss of disc height. Most patients respond to conservative management and surgical interventions well initially, yet a significant number of patients continue to suffer from chronic low back pain. Because of the high prevalence of long-term discogenic pain, regenerative biological therapies, including gene therapies, growth factors, cellular-based injections, and tissue-engineered constructs, have attracted significant attention in light of their potential to directly address the degenerative process. Understanding the pathophysiology of degenerative disc disease is important in both refining existing technologies and developing innovative techniques to reverse the degenerative processes in the discs. In this review, we aimed to cover the underlying pathophysiology of degenerative disc disease as well as its associated risk factors and give a comprehensive summary about the developmental, structural, radiological, and biomechanical properties of human intervertebral discs.

Hydrostatic integrity of the intervertebral disc assessed by MRI

Hydrostatic integrity of the intervertebral disc (IVD) is lost during the process of degeneration. Invasive pressure profilometry (IPP) can quantify it, however, is not applicable for clinical use. We aimed to investigate correlations between IPP and MRI findings to assess non-invasive MRI based methods for prediction of hydrostatic integrity of the intervertebral disc. The pressure profiles of 39 lumbar spinal segments originating from 22 human cadavers were recorded during axial compression in the neutral, the flexed and the extended positions. Disc pressure profiles were measured and mathematically transformed to a novel metric that quantifies pressure profile heterogeneity across the disc. The relationship between pressure profile inhomogeneity ("pressure score") and clinically established magnetic resonance-based classifications systems and demographic parameters was then tested using Spearman correlation tests. Pressure profile inhomogeneities were correlated with IVD degeneration (according to Pfirrmann, rho = 0.43, p = 0.006), endplate defects (according to Rajasekaran, rho = 0.39, p = 0.013), segmental degeneration (according to Farshad, rho = 0.41, p = 0.009) and age (rho = 0.32, p = 0.049). Modic changes per se did not affect the pressure profiles significantly (p = 0.23) and pressure scores did not correlate with BMI (rho = -0.21, p = 0.2). Heterogeneity of segmental IVD pressure profiles is a unique measure of disc function. We demonstrate that established clinical methods for MRI characterization of the intervertebral disc, the endplate and overall segmental degeneration all correlate with the hydrostatic integrity of the IVD and can be used for its assessment.

Trabecular Architecture and Mechanical Heterogeneity Effects on Vertebral Body Strength

PURPOSE OF REVIEW

We aimed to synthesize the recent work on the intra-vertebral heterogeneity in density, trabecular architecture and mechanical properties, its implications for fracture risk, its association with degeneration of the intervertebral discs, and its implications for implant design.

RECENT FINDINGS

As compared to the peripheral regions of the centrum, the central region of the vertebral body exhibits lower density and more sparse microstructure. As compared to the anterior region, the posterior region shows higher density. These variations are more pronounced in vertebrae from older persons and in those adjacent to degenerated discs. Mixed results have been reported in regard to variation along the superior-inferior axis and to relationships between the heterogeneity in density and vertebral strength and fracture risk. These discrepancies highlight that, first, despite the large amount of study of the intra-vertebral heterogeneity in microstructure, direct study of that in mechanical properties has lagged, and second, more measurements of vertebral loading are needed to understand how the heterogeneity affects distributions of stress and strain in the vertebra. These future areas of study are relevant not only to the question of spine fractures but also to the design and selection of implants for spine fusion and disc replacement. The intra-vertebral heterogeneity in microstructure and mechanical properties may be a product of mechanical adaptation as well as a key determinant of the ability of the vertebral body to withstand a given type of loading.

Association of vertebral endplate microstructure with bone strength in men and women

Epidemiological and biomechanical evidence indicates that the risk of vertebral fracture differs between men and women, and that vertebral fracture frequently involves failure of the endplate region. The goal of this study was to compare the bone microstructure of the endplate region-defined as the (bony) vertebral endplate and underlying subchondral trabecular bone-between sexes and to determine whether any such sex differences are associated with vertebral strength. The bone density (volume fraction, apparent density and tissue mineral density) of the superior-most 2 mm of the vertebra, and the bone density and trabecular architecture of the next 5 mm were quantified using micro-computed tomography in human T8 (12 female, 16 male) and L1 (13 female, 12 male) vertebrae. Average density of the vertebra (integral bone mineral density (BMD)) was determined by quantitative computed tomography and compressive strength by mechanical testing. Few differences were found between male and female vertebrae in the density of the endplate region; none were found in trabecular architecture. However, whereas endplate volume fraction was positively correlated with integral BMD in male vertebrae (r = 0.654, p < .001), no correlation was found in the female vertebrae (r = 0.157, p = .455). Accounting for the density of the endplate region improved predictions of vertebral strength (p < .034) and eliminated sex-specificity in the strength prediction that was based on integral BMD alone. These results suggest that the density of the endplate region influences vertebral fracture and that non-invasive assessment of this region's density can contribute to predictions of vertebral strength in men and women.

Bisphosphonate treatment changes regional distribution of trabecular microstructure in human lumbar vertebrae

BACKGROUND

In osteoporosis patients, antiresorptive treatments such as alendronate reduce the resorption of trabecular bone and thus minimize vertebral fracture risk. However, fracture risk reduction efficacy of antiresorptive drugs varies between skeletal sites and is highest for vertebral bone. In human vertebrae, cancellous bone is distributed heterogeneously between regions. This microstructural heterogeneity is changing with patient age and is likely to play a major role in vertebral failure mechanisms and fracture susceptibility. Whether antiresorptive treatment affects the heterogeneity of vertebral microstructure in osteoporosis has not been unraveled.

METHODS

Our aim was to assess whether antiresorptive treatment would have a region-dependent influence on vertebral trabecular bone. Therefore, we used high-resolution peripheral quantitative computed tomography (HR-pQCT), microcomputed tomography (microCT) and uniaxial compression testing to determine the structure and mechanical properties of trabecular bone cores from anterior and posterior regions of 22 lumbar vertebrae from elderly osteoporotic women. We analyzed age-matched ex vivo bone samples from bisphosphonate-treated female osteoporosis patients (age: 82 ± 7y, bisphosphonate treatment period: 4 ± 2 years) along treatment-naïve female controls (82 ± 7y).

RESULTS

MicroCT analysis showed a significantly lower bone volume fraction (p = 0.006) and lower trabecular number (p = 0.003) for the anterior bone cores compared to posterior bone cores in the treatment-naïve group. The bisphosphonate-treated group had a more homogeneous bone volume distribution and did not show significant regional differences in bone volume, it however also displayed significantly different trabecular numbers (p = 0.016). In bone cores of the bisphosphonate-treated group, trabeculae were thicker in comparison to treatment-naïve controls (p = 0.011). Differences in bone volume further resulted in different maximum forces during compression testing between the samples. In addition, the percental difference between BV/TV_μCT in anterior and posterior bone cores was lower in bisphosphonate-treated vertebrae when vertebrae with directly adjacent fractures (n = 3) were excluded.

CONCLUSION

In conclusion, regional trabecular bone microstructure in lumbar vertebrae of bisphosphonate-treated women was more homogeneous compared to treatment-naïve controls. Bisphosphonate treatment, which specifically targets resorption surfaces common in anterior vertebral bone, might have resulted in a region-specific preservation of vertebral microstructure and loading capacity. This could have positive implications for the reduction of wedge fracture risk and add to the explanation of the higher efficacy of fracture risk reduction in vertebrae in comparison to other fracture regions.

Subregional areal bone mineral density (aBMD) is a better predictor of heterogeneity in trabecular microstructure of vertebrae in young and aged women than subregional trabecular bone score (TBS)

BACKGROUND

Currently, bone densitometry fails to identify nearly half of those elderly patients at immediate fracture risk. To improve clinical assessment of vertebral fracture risk, we aimed to determine how the DXA-based 2D parameter Trabecular Bone Score (TBS) relates to subregional variability in 3D trabecular microstructure in young and elderly women compared to aBMD.

METHODS

T12 vertebrae from 29 women (11 young: 32 ± 6 years, 18 aged: 71 ± 5 years) were DXA-scanned ex vivo in anterior-posterior (AP) and lateral projection providing vertebral aBMD and TBS. Additionally, aBMD and TBS were measured for three horizontal (superior, mid-horizontal, inferior) and three vertical subregions (anterior, mid-vertical, posterior) and related to 3D microstructure indices, i.e. bone volume per tissue volume (BV/TV), trabecular number and thickness (Tb.N, Tb.Th), based on HRpQCT.

RESULTS

Subregional high-resolution tomography showed significant differences in trabecular parameters for both age groups: In horizontal subregions, BV/TV was lowest superiorly, Tb.Th was highest mid-horizontally, and Tb.N was lowest mid-horizontally and highest inferiorly. Correspondingly, aBMD varied between horizontal subregions, with differences depending on projection direction. TBS varied only in lateral projections of the aged group, with lower values for the mid-horizontal subregion. In vertical subregions, BV/TV, Tb.N, and aBMD were highest posteriorly for both groups. TBS did not differ between vertical subregions. Regression analysis showed aBMD as a predictor explained more of the variance in subregional 3D microstructure compared to TBS. Stepwise multi-regression analysis revealed only three combinations of subregion, projection, and group where aBMD and TBS were both significant predictors.

CONCLUSIONS

Subregional aBMD reflects variations in trabecular bone microstructure better than subregional TBS for trisected regions. Specifically, lateral aBMD identifies microstructural heterogeneities independent of age and may improve prediction of vertebral strength and susceptibility to specific fracture types.

Disc degeneration promotes regional inhomogeneity in the trabecular morphology of loaded rat tail vertebrae

Background

There is a close relationship between the vertebral trabecular morphology and the condition of the associated disc.

Objective

The relationship between disc degeneration and vertebral trabecular inhomogeneity is unclear. This study aimed to analyse the regional changes of vertebral trabecular morphology after disc degeneration.

Methods

Thirty male Sprague–Dawley rats were randomly assigned to five groups. Group 1 served as an experimental group for the assessment of disc degeneration induced by needle puncture. Group 2 served as a sham group for trabecular morphology analysis. In Group 3, rats had their tail bent between the eighth and tenth coccygeal vertebrae. In Group 4, the tail of rats was bent with a compression load of 4.5 N. In Group 5, rats first underwent disc degeneration induced by a needle puncture before their tail was bent with a compressive load of 4.5 N. Magnetic resonance imaging was performed on all groups, and histological examination was performed on rodents from Group 1. The ninth coccygeal vertebrae of rats from Groups 2–5 were scanned by Micro-computed tomography. Trabecular morphologic changes were assessed in the concave and convex regions by bone volume fraction, trabecular number, trabecular thickness and trabecular separation.

Results

Vertebral trabecular morphology in the concave region improved significantly, whereas the convex region was of significantly lower trabecular morphologic parameters with disc degeneration. The difference in trabecular morphologic parameters between the convex and concave regions increased significantly after disc degeneration.

Conclusion

Disc degeneration promotes regional inhomogeneity in the vertebral trabecular morphology, with the convex region of the vertebrae having the worse trabecular bone morphology than the concave region.

The translational potential of this article

Our study indicates that disc degeneration promotes regional inhomogeneity in the vertebral trabecular morphology. Regional variations in trabecular microarchitecture are helpful to predict vertebral fragility. This may help to elucidate the mechanisms by which disc degeneration contributes to vertebral fracture.

Correspondence between bone mineral density and intervertebral disc degeneration across age and sex

The distribution of bone tissue within the vertebra can modulate vertebral strength independently of average density and may change with age and disc degeneration. Our results show that the age-associated decrease in bone density is spatially non-uniform and associated with disc health, suggesting a mechanistic interplay between disc and vertebra.

PURPOSE

While the decline of bone mineral density (BMD) in the aging spine is well established, the extent to which age influences BMD distribution within the vertebra is less clear. Measures of regional BMD (rBMD) may improve predictions of vertebral strength and suggest how vertebrae might adapt with intervertebral disc degeneration. Thus, we aimed to assess how rBMD values were associated with age, sex, and disc height loss (DHL).

METHODS

We measured rBMD in the L3 vertebra of 377 participants from the Framingham Heart Study (41-83 years, 181 M/196 F). Integral (Int.BMD) and trabecular BMD (Tb.BMD) were measured from QCT images. rBMD ratios (anterior/posterior, superior/mid-transverse, inferior/mid-transverse, and central/outer) were calculated from the centrum. A radiologist assigned a DHL severity score to adjacent intervertebral discs (L2-L3 and L3-L4).

RESULTS

Int.BMD and Tb.BMD were both associated with age, though the decrease across age was greater in women (Int.BMD, - 2.6 mg/cm3 per year; Tb.BMD, - 2.6 mg/cm3 per year) than men (Int.BMD, - 0.5 mg/cm3 per year; Tb.BMD, - 1.2 mg/cm3 per year). The central/outer (- 0.027/decade) and superior/mid-transverse (- 0.018/decade) rBMD ratios were negatively associated with age, with similar trends in men and women. Higher Int.BMD or Tb.BMD was associated with increased odds of DHL after adjusting for age and sex. Low central/outer ratio and high anterior/poster and superior/mid-transverse ratios were also associated with increased odds of DHL.

CONCLUSIONS

Our results indicate that the distribution of bone within the L3 vertebra is different across age, but not between sexes, and is associated with disc degeneration.

Dynamic adaptation of vertebral endplate and trabecular bone following annular injury in a rat model of degenerative disc disease

BACKGROUND CONTEXT

Degenerative disc disease (DDD) is associated with longitudinal remodeling of paravertebral tissues. Although chronic vertebral changes in advanced stages of DDD are well-studied, very little data exists on acute vertebral bone remodeling at the onset and progression of DDD.

PURPOSE

To longitudinally characterize bony remodeling in a rodent model of disc injury-induced DDD.

STUDY DESIGN

In vivo animal study involving a rat annulus fibrosus injury model of DDD.

METHODS

Eight female Lewis rats were assigned to intervertebral disc (IVD) injury (Puncture) or sham surgery (Sham). All rats underwent anterior, transperitoneal approach to the lumbar spine, and Puncture rats underwent annulus fibrosus injury at the L3-L4 and L5-L6 IVDs (n = 8 per group). Live micro computed tomography imaging (10-μm voxel size) was performed 1 week before surgery and postoperatively at 2-week intervals up to a 12-week endpoint. Bone morphology and densitometry of the cranial vertebral body and bony endplate were analyzed and reported with respect to the preoperative baseline scan. Sagittal Safranin-O/Fast-Green and Toluidine Blue histology evaluated using the Rutges IVD score and a custom vertebral endplate score.

RESULTS

Vertebral trabecular tissue mineral density (TMD), vertebral trabecular spacing, endplate TMD, and endplate apparent bone mineral density were all significantly greater in Puncture compared with Sham at 4 weeks and each subsequent timepoint. Puncture rats exhibited marginally lower endplate total volume. Anterior endplate osteophyte formation and central physeal ossification were observed in Puncture rats. Endpoint histological analysis demonstrated moderate evidence of IVD degeneration, indicating that vertebral bone adaptation occurs in the acute phases of DDD onset and progression.

CONCLUSIONS

Annulus injury-induced DDD leads to acute and progressive changes to the morphology and densitometry of bone in the adjacent vertebral bodies and endplates.

Structure-function relationships at the human spinal disc-vertebra interface

Damage at the intervertebral disc-vertebra interface associates with back pain and disc herniation. However, the structural and biomechanical properties of the disc-vertebra interface remain underexplored. We sought to measure mechanical properties and failure mechanisms, quantify architectural features, and assess structure-function relationships at this vulnerable location. Vertebra-disc-vertebra specimens from human cadaver thoracic spines were scanned with micro-computed tomography (μCT), surface speckle-coated, and loaded to failure in uniaxial tension. Digital image correlation (DIC) was used to calculate local surface strains. Failure surfaces were scanned using scanning electron microscopy (SEM), and adjacent sagittal slices were analyzed with histology and SEM. Seventy-one percent of specimens failed initially at the cartilage endplate-bone interface of the inner annulus region. Histology and SEM both indicated a lack of structural integration between the cartilage endplate (CEP) and bone. The interface failure strength was increased in samples with higher trabecular bone volume fraction in the vertebral endplates. Furthermore, failure strength decreased with degeneration, and in discs with thicker CEPs. Our findings indicate that poor structural connectivity between the CEP and vertebra may explain the structural weakness at this region, and provide insight into structural features that may contribute to risk for disc-vertebra interface injury. The disc-vertebra interface is the site of failure in the majority of herniation injuries. Here we show new structure-function relationships at this interface that may motivate the development of diagnostics, prevention strategies, and treatments to improve the prognosis for many low back pain patients with disc-vertebra interface injuries. © 2017 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 36:192-201, 2018.

Is bone density associated with intervertebral disc pressure in healthy and degenerated discs?

The coupling of the intervertebral disc (IVD) and vertebra as a biomechanical unit suggests that changes in the distribution of pressure within the IVD (intradiscal pressure, IDP) as a result of disc degeneration can influence the distribution of bone density within the vertebra, and vice versa. The goal of this study was to assess the correspondence between IDP and bone density in the adjacent vertebrae, with emphasis on how this correspondence differs between healthy and degenerated IVDs. Bone density of the endplates and subchondral bone in regions adjacent to the anterior and posterior annulus fibrosus (aAF and pAF, respectively) and nucleus pulposus (NP) was measured via quantitative computed tomography (QCT) in 61 spine segments (T7-9, T9-11, T10-12; 71±14years). IDP was measured in the aAF, NP, and pAF regions in 26 of the spine segments (68±16years) while they were tested in flexed (5°) or erect postures. Disc degeneration was assessed by multiple grading schemes. No correlation was found between bone density and IDP in either posture (p>0.104). Regional variations in IDP and, to a greater extent bone density, were found to change with advancing degeneration: both IDP (p=0.045) and bone density (p=0.024) decreased in the NP region relative to the aAF region. The finding of only a modest correspondence between degeneration-associated changes in IDP and bone density may arise from complexity in how IDP relates to mechanical force transmission through the endplate and from limitations of the available IVD grading schemes in estimating the mechanical behavior of the IVD.

Influence of the shape of the micro-finite element model on the mechanical properties calculated from micro-finite element analysis

Assessing the biomechanical properties of trabecular bone is of major biological and clinical significance for the research of bone diseases, fractures and their treatments. Micro-finite element (µFE) models are becoming increasingly popular for investigating the biomechanical properties of trabecular bone. The shapes of µFE models typically include cube and cylinder. Whether there are differences between cubic and cylindrical µFE models has not yet been studied. In the present study, cubic and cylindrical µFE models of human vertebral trabecular bone were constructed. A 1% strain was prescribed to the model along the superior-inferior direction. E values were calculated from these models, and paired t-tests were performed to determine whether these were any differences between E values obtained from cubic and cylindrical models. The results demonstrated that there were no statistically significant differences in the E values between cubic and cylindrical models, and there were no significant differences in Von Mises stress distributions between the two models. These findings indicated that, to construct µFE models of vertebral trabecular bone, cubic or cylindrical models were both feasible. Choosing between the cubic or cylindrical µFE model is dependent upon the specific study design.

Cranio-caudal asymmetries in trabecular architecture reflect vertebral fracture patterns

Clinically, vertebral fractures often occur in the upper lumbar spine and involve the superior endplate of a vertebra (which is immediately caudal to a disc). Knowledge that the cranial endplate of a disc is thicker and has greater bone mineral density (BMD) than the corresponding caudal endplate helps to explain this phenomenon. In this study, we investigated structural differences in vertebral trabeculae on either side of a lumbar disc to provide further insight into vertebral fracture risk. As the focus is trabecular difference within a spinal motion segment, we define cranial and caudal vertebral trabeculae relative to the disc. Ninety-two spinal motion segments from 46 cadaveric lumbar spines (males, mean age 50years, range 21-63years) were studied. Disc narrowing on radiography and spread of barium sulfate (BaSO₄) on discography were measured to indicate disc degeneration. Micro-computed tomography (μCT) images were obtained at a resolution of 82μm for each vertebra and processed to include only vertebral trabeculae. Using image processing, the vertebral trabeculae were divided into superior and inferior halves, and then into central and peripheral regions which were approximately opposite to the disc pulposus and annulus, and further into anterior and posterior sub-regions. Microarchitecture measurements for each vertebral region were obtained to determine the differences between the cranial and caudal trabeculae (relative to disc) and their associations with age and disc degeneration within each spinal motion segment. Data from the upper (L1/2-L3/4) and lower (L4/5) lumbar segments were analyzed separately. In the upper lumbar region, the trabeculae cranial to a disc on average had 5.3% greater BMD and trabecular bone volume, 3.6% greater trabecular number, 9.7% greater connectivity density, and 3.7% less trabecular separation than the corresponding caudal trabeculae (P<0.05 for all). Similar trends were observed in peripheral, anterior and posterior regions, but not in central region. No structural difference was observed in the trabeculae of L4/5 segment. Structural asymmetries of vertebral trabeculae were not associated with age, disc degeneration, or disc narrowing. Vertebral trabecular parameters cranial to the disc were greater than caudally in the upper but not in the lower lumbar region. Findings further explain why vertebral fractures are more common in the upper lumbar region and more frequently involve the endplate caudal to a disc.

Relationship between spinal osteoarthritis and vertebral fractures in men older than 50 years: data from the Camargo Cohort Study

Spinal osteoarthritis has been suggested as a risk factor for vertebral fractures. However, results are conflicting: most of the data are focused on the lumbar region, and referred to postmenopausal women, whereas data for men are scarce. The aim of this study is to assess the relationship between spinal osteoarthritis and vertebral fractures in men over 50 years of age. We conducted a cross-sectional study, nested in a prospective population-based cohort, including 507 community-dwelling men, 93 of them with at least one vertebral fracture. Vertebral fractures, osteophytosis, and disc space narrowing (DSN) were assessed by lateral thoracic and lumbar radiographs. Anthropometric, clinical, and densitometric variables were also analyzed. A multiple logistic regression model was performed. Eighty-five percent of vertebral fractures were located at the thoracic spine. Osteophytosis and DSN showed a bimodal distribution, with major frequency peaks at mid- and distal lumbar spine. The three distributions overlapped around the T9 vertebra. We did not find any relationship between lumbar osteoarthritis and vertebral fractures. Nevertheless, thoracic osteophytosis (OR, 1.84; 95 % CI, 1.05-3.17; p = 0.03) and DSN (OR, 2.52; 95 % CI, 1.43-4.46; p = 0.001) were found to be independently associated with prevalent vertebral fractures, after adjusting for confounders. Our results suggest a positive relationship between radiologic osteoarthritic changes at the thoracic spine and prevalent vertebral fractures in men more than 50 years of age. Osteoarthritis may act as a local risk factor, in addition to other mechanical factors, resulting in a greater propensity to fracture, especially at the mid-thoracic region.

Quantitative, 3D Visualization of the Initiation and Progression of Vertebral Fractures Under Compression and Anterior Flexion

The biomechanical mechanisms leading to vertebral fractures are not well understood. Clinical and laboratory evidence suggests that the vertebral endplate plays a key role in failure of the vertebra as a whole, but how this role differs for different types of vertebral loading is not known. Mechanical testing of human thoracic spine segments, in conjunction with time-lapsed micro-computed tomography, enabled quantitative assessment of deformations occurring throughout the entire vertebral body under axial compression combined with anterior flexion ("combined loading") and under axial compression only ("compression loading"). The resulting deformation maps indicated that endplate deflection was a principal feature of vertebral failure for both loading modes. Specifically, the onset of endplate deflection was temporally coincident with a pronounced drop in the vertebra's ability to support loads. The location of endplate deflection, and also vertebral strength, were associated with the porosity of the endplate and the microstructure of the underlying trabecular bone. However, the location of endplate deflection and the involvement of the cortex differed between the two types of loading. Under the combined loading, deflection initiated, and remained the largest, at the anterior central endplate or the anterior ring apophysis, depending in part on health of the adjacent intervertebral disc. This deflection was accompanied by outward bulging of the anterior cortex. In contrast, the location of endplate deflection was more varied in compression loading. For both loading types, the earliest progression to a mild fracture according to a quantitative morphometric criterion occurred only after much of the failure process had occurred. The outcomes of this work indicate that for two physiological loading modes, the vertebral endplate and underlying trabecular bone are critically involved in vertebral fracture. These outcomes provide a strong biomechanical rationale for clinical methods, such as algorithm-based qualitative (ABQ) assessment, that diagnose vertebral fracture on the basis of endplate depression. © 2015 American Society for Bone and Mineral Research.

A Predictive Model of Vertebral Trabecular Anisotropy From Ex Vivo Micro-CT

Spine-related disorders are amongst the most frequently encountered problems in clinical medicine. For several applications such as 1) to improve the assessment of the strength of the spine, as well as 2) to optimize the personalization of spinal interventions, image-based biomechanical modeling of the vertebrae is expected to play an important predictive role. However, this requires the construction of computational models that are subject-specific and comprehensive. In particular, they need to incorporate information about the vertebral anisotropic micro-architecture, which plays a central role in the biomechanical function of the vertebrae. In practice, however, accurate personalization of the vertebral trabeculae has proven to be difficult as its imaging in vivo is currently infeasible. Consequently, this paper presents a statistical approach for accurate prediction of the vertebral fabric tensors based on a training sample of ex vivo micro-CT images. To the best of our knowledge, this is the first predictive model proposed and validated for vertebral datasets. The method combines features selection and partial least squares regression in order to derive optimal latent variables for the prediction of the fabric tensors based on the more easily extracted shape and density information. Detailed validation with 20 ex vivo T12 vertebrae demonstrates the accuracy and consistency of the approach for the personalization of trabecular anisotropy.