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

Risk of bony endplate failure during vertebral fracture

The endplate region of the vertebra, which includes the bony endplate (BEP) and underlying subchondral trabecular bone (STB), is critically involved in vertebral fracture (VF). While evidence abounds that failure initiates in the endplate region, the relative risk of failure of the BEP vs. STB has not been established. In this study, micro-finite element models were constructed of L1 vertebrae (n = 21) that were mechanically tested in a prior study and given experimentally matched boundary conditions corresponding to the vertebra's yield point. Volumes of interest (VOIs) were defined corresponding to the BEP and STB; the remainder was defined as the mid-vertebral body (MVB). The proportion of elements within each VOI that yielded was defined as the VOI yield fraction, and this value divided by the yield fraction of the entire model was defined as the normalized yield fraction. While yield fraction did not differ across VOIs (p = 0.179), normalized yield fraction was greater in the BEP than STB and MVB (p < 0.001), indicating a higher risk of yield in the BEP compared to the other two VOIs. None of the yield fractions was correlated with BEP or STB microstructure, and tension (rather than compression) was the dominant mode of tissue level yield. These findings indicate that the BEP, more so than the STB, is likely the site of VF initiation and that current methods of screening for VF risk, because they omit specific analysis of the BEP, are missing the region that matters the most. The endplate region of the vertebra, which includes the bony endplate (BEP) and underlying subchondral bone (SB), is critically involved in vertebral fracture (VF). While evidence abounds that failure initiates in the endplate region, the relative risk of failure of the BEP vs. SB has not been established. In this study, micro-finite element models were constructed of L1 vertebrae (n = 21) that had been mechanically tested in a prior study, and they were given experimentally matched boundary conditions corresponding to the vertebra's yield point. Volumes of interest (VOIs) were defined corresponding to the BEP and SB; the remainder was defined as the mid-vertebral body (MVB). The proportion of yielded elements within each VOI was defined as the VOI yield fraction, and this value divided by the yield fraction of the entire model was defined as the normalized yield fraction. While yield fraction did not differ across VOIs (p = 0.179), normalized yield fraction was greater in the BEP than SB and MVB (p < 0.001), indicating a higher risk of yield in the BEP compared to the other two VOIs. None of the yield fractions was correlated with BEP or SB microstructure, and tension (rather than compression) was the dominant mode of tissue level yield. These findings indicate that the BEP, more so than the SB, is likely the site of VF initiation and that current methods of screening for VF risk, because they omit specific analysis of the BEP, are missing the region that matters the most.

Comparison of predictive value for cage subsidence between MRI-based endplate bone quality and vertebral bone quality scores following transforaminal lumbar interbody fusion: a retrospective propensity-matched study

BACKGROUND CONTEXT

Cage subsidence after lumbar fusion can lead to many adverse outcomes. Low bone mineral density (BMD) is a widely recognized risk factor for cage subsidence. Conventional methods can predict and evaluate BMD, but there are many shortcomings. Recently, MRI-based assessment of bone quality in specific parts of the vertebral body has been proposed, including scores for vertebral bone quality (VBQ) and endplate bone quality (EBQ). However, the predictive accuracy of the two scoring systems for cage subsidence after transforaminal lumbar interbody fusion (TLIF) remains unknown. Therefore, we investigated MRI-based VBQ and EBQ scores for assessing bone quality and compared their predictive value for cage subsidence after TLIF.

PURPOSE

To compare the predictive value between MRI-based VBQ and EBQ scores for cage subsidence after TLIF.

STUDY DESIGN/SETTING

A retrospective case-control study.

PATIENTS SAMPLE

Patients with degenerative lumbar diseases underwent single-level TLIF at our medical center between 2014 and 2020, all of whom had preoperative MRIs available.

OUTCOMES MEASURES

Cage subsidence, disc height, VBQ score, EBQ score, upper and lower vertebral body bone quality (UL-VBQ) score.

METHODS

Data were retrospectively examined for a consecutive sample of 346 patients who underwent TLIF at our medical center between 2014 and 2020. Patients who subsequently experienced cage subsidence or not were matched to each other based on propensity scoring, and the two matched groups (52 patients each) were compared using conditional logistic regression to investigate the association between the potential radiographic factors and cage subsidence. Scores for VBQ and EBQ were assessed for their ability to predict cage subsidence in the matched patients based on the area under the receiver operative characteristic curve (AUC).

RESULTS

Among matched patients, those who suffered cage subsidence had significantly higher VBQ score (3.7 vs 3.1, p<.001) and EBQ score (5.0 vs 4.3, p<.001), and regression linked greater risk of subsidence to higher VBQ score (OR 4.557, 95% CI 1.076-19.291, p=.039) and higher EBQ score (OR 5.396, 95% CI 1.158-25.146, p=.032). A cut-off VBQ score of 3.4 predicted the cage subsidence among matched patients with an AUC of 0.799, sensitivity of 84.6%, and specificity of 69.2%. A cut-off EBQ score of 4.7 predicted subsidence with an AUC of 0.829, sensitivity of 76.9%, and specificity of 82.7%.

CONCLUSION

Higher VBQ and EBQ scores are associated with a greater risk of cage subsidence following TLIF, and EBQ may perform better because of greater specificity.

The role of bone metastases on the mechanical competence of human vertebrae

Spine is the most common site for bone metastases. The evaluation of the mechanical competence and failure location in metastatic vertebrae is a biomechanical and clinical challenge. Little is known about the failure behaviour of vertebrae with metastatic lesions. The aim of this study was to use combined micro-Computed Tomography (microCT) and time-lapsed mechanical testing to reveal the failure location in metastatic vertebrae. Fifteen spine segments, each including a metastatic and a radiologically healthy vertebra, were tested in compression up to failure within a microCT. Volumetric strains were measured using Digital Volume Correlation. The images of undeformed and deformed specimens were overlapped to identify the failure location. Vertebrae with lytic metastases experienced the largest average compressive strains (median ± standard deviation: -8506 ± 4748microstrain), followed by the vertebrae with mixed metastases (-7035 ± 15605microstrain), the radiologically healthy vertebrae (-5743 ± 5697microstrain), and the vertebrae with blastic metastases (-3150 ± 4641microstrain). Strain peaks were localised within and nearby the lytic lesions or around the blastic tissue. Failure between the endplate and the metastasis was identified in vertebrae with lytic metastases, whereas failure localised around the metastasis in vertebrae with blastic lesions. This study showed for the first time the role of metastases on the vertebral internal deformations. While lytic lesions lead to failure of the metastatic vertebra, vertebrae with blastic metastases are more likely to induce failure in the adjacent vertebrae. Nevertheless, every metastatic lesion affects the vertebral deformation differently, making it essential to assess how the lesion affects the bone microstructure. These results suggest that the properties of the lesion (type, size, location within the vertebral body) should be considered when developing clinical tools to predict the risk of fracture in patients with metastatic lesions.

Applying a homogeneous pressure distribution to the upper vertebral endplate: Validation of a new loading system, pilot application to human vertebral bodies, and finite element predictions of DIC measured displacements and strains

Image-based personalized Finite Element Models (pFEM) could detect alterations in physiological deformation of human vertebral bodies, but their accuracy has been seldom reported. Meaningful validation experiments should allow vertebral endplate deformability and ensure well-controlled boundary conditions. This study aimed to (i) validate a new loading system to apply a homogeneous pressure on the vertebral endplate during vertebral body compression regardless of endplate deformation; (ii) perform a pilot study on human vertebral bodies measuring surface displacements and strains with Digital Image Correlation (DIC); (iii) determine the accuracy of pFEM of the vertebral bodies. Homogeneous pressure application was achieved by pressurizing a fluid silicone encased in a rubber silicone film acting on the cranial endplate. The loading system was validated by comparing DIC-measured longitudinal strains and lower-end contact pressures, measured on three homogeneous pseudovertebrae of constant transversal section at 2.0 kN, against theoretically calculated values. Longitudinal strains and contact pressures were rather homogeneous, and their mean values close to theoretical calculations (5% underestimation). DIC measurements of surface longitudinal and circumferential displacements and strains were obtained on three human vertebral bodies at 2.0 kN. Complete displacement and strain maps were achieved for anterolateral aspects with random errors ≤0.2 μm and ≤30 μstrain, respectively. Venous plexus and double curvatures limited the completeness and accuracy of DIC data in posterior aspects. pFEM of vertebral bodies, including cortical bone mapping, were built from computed tomography images. In anterolateral aspects, pFEM accuracy of the three vertebrae was: (i) comparable to literature in terms of longitudinal displacements (R2>0.8); (ii) extended to circumferential displacements (pooled data: R2>0.9) and longitudinal strains (zero median error, 95% error: <27%). Circumferential strains were overestimated (median error: 39%). The new methods presented may permit to study how physiological and pathologic conditions influence the ability of vertebral endplates/bodies to sustain loads.

Vertebral endplate defects are associated with bone mineral density in lumbar degenerative disc disease

PURPOSE

Evidence has shown that lumbar vertebral endplate defects are clinically relevant and closely related to disc degeneration, but the relationship between endplate defects and bone mineral density (BMD) remains unclear. This study aimed to explore the association between endplate defects and BMD-related values in patients with lumbar degenerative disc disease (LDD).

METHODS

Three hundred and twenty-five Chinese adult subjects diagnosed with LDD underwent dual energy X-ray absorptiometry. Endplate defects were classified using lumbar MRI. Groups were subdivided based on the occurrence rates of defect endplates. BMD at the lumbar vertebral and bilateral femur necks was compared between groups, and the association between endplate defects and lumbar BMD-related values was analyzed and adjusted for confounders including age, sex, serum levels of 25-hydroxy vitamin D (25(OH)D), calcium (Ca) and phosphorus (P).

RESULTS

Of 325 patients and 3250 endplates, 59.72% had defects, and 188 patients were divided into the higher defect rate group (occurrence rate > 50%). The higher defect rate group was associated with older age, more common postmenopausal females, higher osteoporosis rates and lower serum Ca and P levels. Lumbar BMD was greater than that at bilateral femur necks and was not equal to osteoporosis diagnosis. Endplate defects were more prevalent in lower segments. The occurrence of endplate defects was positively associated with lumbar BMD-related values in the partial correlation analysis. The association between endplate defects and lumbar BMD varies for subtypes and segments, with a trend of positive association in rim and erosive subtypes after adjusting for confounders.

CONCLUSIONS

The present study demonstrated that the occurrence of endplate defects was associated with greater lumbar BMD values in patients with LDD. This association varies for different defect subtypes and segments. The results indicated that endplate defects should be taken into consideration in osteoporosis treatment to alleviate disc degeneration.

A retrospective study on the efficacy and safety of bone cement in the treatment of endplate fractures

Background

Endplate fractures is an important factor affecting the curative effect of percutaneous kyphoplasty for spinal fracture. The purpose of this study is to investigate the effect of sealing endplate fracture with bone cement on minimally invasive treatment of spinal fracture.

Methods

A total of 98 patients with osteoporotic vertebral fractures combined with endplate fractures treated with bone cement surgery in our hospital were retrospectively analyzed. They were grouped according to whether bone cement was involved in the endplate fractures. Group A: bone cement was not only distributed in the fractured vertebral body, but also dispersed into the endplate fractures. Group B: bone cement was confined to the fractured vertebra but did not diffuse into the cracks of the endplate. The basic information, imaging changes of the fractured vertebral body, VAS score, ODI score, bone cement distribution and postoperative complications of the two groups were analyzed and compared.

Results

The height of the injured vertebra and the kyphotic Cobb angle in the two groups were significantly improved after surgery, but the anterior height of the vertebra in group B was lower than that in group A and the kyphotic Cobb angle was higher than that in group A at the last follow-up (P < 0.05). VAS score and ODI score in 2 groups were significantly improved after operation (P < 0.05), but the VAS score and ODI score in group A were lower than those in group B at the last follow-up (P < 0.05). The incidence of bone cement leakage and adjacent vertebral fracture in group A was higher than that in group B (P < 0.05).

Conclusion

Diffusion of bone cement into the cracks of the endplate may also restore and maintain the height of the injured vertebra, relieve pain and restore lumbar function. However, diffusion of bone cement into the cracks of the endplate can increase the incidence of cement leakage and adjacent vertebral fractures.

Low Hounsfield units on computed tomography are associated with cage subsidence following oblique lumbar interbody fusion (OLIF)

BACKGROUND CONTEXT

Cage subsidence is one of the most common complications following lumbar interbody fusion surgery. Low bone mineral density (BMD) is an important risk factor that contributes to cage subsidence. Hounsfield units (HU) obtained from clinical computed tomography (CT) scans provided a reliable method for determining regional BMD. The association between HU and cage subsidence following oblique lumbar interbody fusion (OLIF) remains unclear.

PURPOSE

The objective of this study is to evaluate the association between vertebral HU value and cage subsidence following OLIF.

STUDY DESIGN/SETTING

A retrospective study.

PATIENT SAMPLE

Adults with degenerative spinal conditions underwent single-level OLIF at our institution from October 2017 and August 2020 OUTCOME MEASURES: Cage subsidence, disc height, vertebral body global HU value, upper and lower instrumented vertebrae HU value, endplate HU value, fusion rate.

METHODS

This retrospective study was conducted on patients who underwent single-level OLIF at one institution between October 2017 and August 2020. Cage subsidence was measured using the CT scan postoperatively based on the cage protrusion through the vertebral endplates. The HU values were measured from preoperative CT according to previously reported methods.

RESULTS

A total of 70 patients with a mean follow-up of 15.4 months were included in the analysis. The subsidence rate was 25.7% (n=18/70). The average cage subsidence was 2.2 mm, with a range of 0-7.7 mm. No significant difference was found in age, sex, or body mass index (BMI) between the two groups. The mean global HU value of the lumbar vertebral body (L1-5) was 142.7±30.1 in nonsubsidence and 103.7±11.5 in subsidence (p=.004). The upper instrumented vertebrae (UIV) HU value was 141.4±29.7 in the nonsubsidence and 101.1±10.2 in subsidence, (p=.005). The lower instrumented vertebrae (LIV) HU value was 147.4±34.9 in nonsubsidence and 108.1±13.7 in subsidence, (p<.001). The AUC of the UIV HU value was 0.917 (95% CI: 0.853-0.981), and the most appropriate threshold of the HU value was 115 (sensitivity: 84.6%, specificity: 100%). The AUC of the LIV HU value was 0.893 (95%CI: 0.819-0.966), and the most appropriate threshold of the HU value was 125 (sensitivity: 76.9%, specificity: 100%). The mean upper endplate HU value was 235.4±50.9, and the mean lower endplate HU value was 193.4±40.3. No significant difference (upper endplate p=.314, lower endplate p=.189) was observed between the two groups.

CONSLUSIONS

Lower preoperative vertebral body HU values were associated with cage subsidence after single-level OLIF. However, the endplate HU values were not associated with cage subsidence. Preoperative HU measurement is useful in the prediction of the cage subsidence.

Structure-function relationships of the human vertebral endplate

BACKGROUND

Although deformation and fracture of the vertebral endplate have been implicated in spinal conditions such as vertebral fracture and disc degeneration, few biomechanical studies of this structure are available. The goal of this study was to quantify the mechanical behavior of the vertebral endplate.

METHODS

Eight-five rectangular specimens were dissected from the superior and/or inferior central endplates of human lumbar spine segments L1 to L4. Micro-computed tomography (μCT) imaging, four-point-bend testing, and ashing were performed to quantify the apparent elastic modulus and yield stress (modulus and yield stress, respectively, of the porous vertebral endplate), tissue yield stress (yield stress of the tissue of the vertebral endplate, excluding pores), ultimate strain, fracture strain, bone volume fraction (BV/TV), bone mineral density (BMD), and various measures of tissue density and composition (tissue mineral density, ash fraction, and ash density). Regression was used to assess the dependence of mechanical properties on density and composition.

RESULTS

Wide variations in elastic and failure properties, and in density and tissue composition, were observed. BMD and BV/TV were good predictors of many of the apparent-level mechanical properties, including modulus, yield stress, and in the case of the inferior vertebral endplate, failure strains. Similar values of the mechanical properties were noted between superior and inferior vertebral endplates. In contrast to the dependence of apparent stiffness and strength on BMD and BV/TV, none of the mechanical properties depended on any of the tissue-level density measurements.

CONCLUSION

The dependence of many of the mechanical properties of the vertebral endplate on BV/TV and BMD suggests possibilities for noninvasive assessment of how this region of the spine behaves during habitual and injurious loading. Further study of the nonmineral components of the endplate tissue is required to understand how the composition of this tissue may influence the overall mechanical behavior of the vertebral endplate.

Finite element analysis of wedge and biconcave deformity in four different height restoration after augmentation of osteoporotic vertebral compression fractures

Purpose

Biomechanical comparison of wedge and biconcave deformity of different height restoration after augmentation of osteoporotic vertebral compression fractures was analyzed by three-dimensional finite element analysis (FEA).

Methods

Three-dimensional finite element model (FEM) of T11-L2 segment was constructed from CT scan of elderly osteoporosis patient. The von Mises stresses of vertebrae, intervertebral disc, facet joints, displacement, and range of motion (ROM) of wedge and biconcave deformity were compared at four different heights (Genant 0–3 grade) after T12 vertebral augmentation.

Results

In wedge deformity, the stress of T12 decreased as the vertebral height in neutral position, flexion, extension, and left axial rotation, whereas increased sharply in bending at Genant 0; L1 and L2 decreased in all positions excluding flexion of L2, and T11 increased in neutral position, flexion, extension, and right axial rotation at Genant 0. No significant changes in biconcave deformity. The stress of T11-T12, T12-L1, and L1-L2 intervertebral disc gradually increased or decreased under other positions in wedge fracture, whereas L1-L2 no significant change in biconcave fracture. The utmost overall facet joint stress is at Genant 3, whereas there is no significant change under the same position in biconcave fracture. The displacement and ROM of the wedge fracture had ups and downs, while a decline in all positions excluding extension in biconcave fracture.

Conclusions

The vertebral restoration height after augmentation to Genant 0 affects the von Mises stress, displacement, and ROM in wedge deformity, which may increase the risk of fracture, whereas restored or not in biconcave deformity.

Influence of endplate size and implant positioning of vertebral body replacements on biomechanics and outcome

BACKGROUND

Spinal stabilization by an anterior vertebral body replacement is frequently used in patients suffering from destroyed vertebral bodies. The aim of this study was to analyse (i) the choice of endplate size and positioning of vertebral body replacements in daily patient care and (ii) if these factors have an influence on clinical and radiological outcomes.

METHOD

Patients' outcomes were analysed three years after vertebral body replacement implantation using the visual analogue scale spine score. Safe zones on the vertebral body endplates were defined. Overall endplate coverage and implant subsidence were evaluated by CT and X-ray. Compression tests were performed on 22 lumbar vertebral bodies to analyse endplates sizes' influence on subsidence.

FINDING

Mean coverage of the vertebral body's superior and inferior endplates by the vertebral body replacement was 27.8% and 30.8%, respectively. Mean overlap of the safe zone by the implant was 49.8% and 40.6%. Mean subsidence was 1.1 ± 1.2 mm, but it did not have any effect on the outcome. In the compression tests, no significant difference (p = 0.468) was found between the two endplate sizes.

INTERPRETATION

Coverage of vertebral body endplates and positioning of implants in the safe zone did not entirely comply with the given recommendations. The amount of endplate coverage had no influence on subsidence or long-term outcomes in daily patient care. On the other hand, correct positioning of the implant may influence its subsidence.

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.

The Association Between Endplate Changes and Risk for Early Severe Cage Subsidence Among Standalone Lateral Lumbar Interbody Fusion Patients

STUDY DESIGN

Retrospective case series.

OBJECTIVE

The aim of this study was to investigate the association of Modic type endplate changes with the risk of severe subsidence after standalone lateral lumbar interbody fusion (SA-LLIF).

SUMMARY OF BACKGROUND DATA

It has been reported that certain endplate radiolographic features are associated with higher regional bone mineral density (BMD) in the adjacent vertebrae in the lumbar spine. It remains unclear whether these changes have protective effects against osteoporotic complications such as cage subsidence after lumbar surgery.

METHODS

We reviewed patients undergoing SA-LLIF from 2007 to 2016 with a follow-up >6 months. Cage subsidence was assessed utilizing the grading system by Marchi et al. As potential contributing factors for cage subsidence, we measured the endplate volumetric BMD (EP-vBMD) and the standard trabecular volumetric BMD measurement in the vertebral body. Modic changes (MC) on magnetic resonance imaging were measured as a qualitative factor for endplate condition. Univariate analysis and multivariate logistic regression analyses with a generalized mixed model were conducted.

RESULTS

Two hundred six levels in 97 patients were included in the final analysis. Mean age (± SD) was 66.7 ± 10.7. Sisty-sdpercent of the patients were female. Severe subsidence was observed in 66 levels (32.0%). After adjusting for age, bone morphogenetic protein (BMP) use, and number of levels fused, the presence of MC type 2 was significantly associated with lower risk of severe subsidence (OR = 0.28 [0.09-0.88], P = 0.029). Whereas, EP-vBMD did not demonstrate a statistical significance (p = 0.600).

CONCLUSION

The presence of a Modic type 2 change was significantly associated with lower odds of severe subsidence after SA-LLIF. Nonetheless, this significant association was independent from regional EP-vBMD values. This finding suggests that microstructural and/or material property changes associated with Modic type 2 changes might have a protective effect in this patient population.

LEVEL OF EVIDENCE

4.