Heterogeneity and Spatial Distribution of Intravertebral Trabecular Bone Mineral Density in the Lumbar Spine Is Associated With Prevalent Vertebral Fracture

Anatomically and mechanically conforming patient-specific spinal fusion cages designed by full-scale topology optimization

Cage subsidence after instrumented lumbar spinal fusion surgery remains a significant cause of treatment failure, specifically for posterior or transforaminal lumbar interbody fusion. Recent advancements in computational techniques and additive manufacturing, have enabled the development of patient-specific implants and implant optimization to specific functional targets. This study aimed to introduce a novel full-scale topology optimization formulation that takes the structural response of the adjacent bone structures into account in the optimization process. The formulation includes maximum and minimum principal strain constraints that lower strain concentrations in the adjacent vertebrae. This optimization approach resulted in anatomically and mechanically conforming spinal fusion cages. Subsidence risk was quantified in a commercial finite element solver for off-the-shelf, anatomically conforming and the optimized cages, in two representative patients. We demonstrated that the anatomically and mechanically conforming cages reduced subsidence risk by 91% compared to an off-the-shelf implant with the same footprint for a patient with normal bone quality and 54% for a patient with osteopenia. Prototypes of the optimized cage were additively manufactured and mechanically tested to evaluate the manufacturability and integrity of the design and to validate the finite element model.

Low hounsfield unit values at sagittal section on computed tomography predicts vertebral fracture following short spinal fusion

BACKGROUND

Preoperative identification of osteoporosis during spine surgery is of critical importance. Additionally, the Hounsfield units (HU) measured using computed tomography (CT) have gained considerable attention. This study aimed to propose a more accurate and convenient screening method for predicting vertebral fractures after spinal fusion in elderly patients by analyzing the HU value of different range of interests of thoracolumbar spine.

METHODS

Our sample pool for analysis consisted of 137 elderly female patients aged >70 years who underwent one- or two-level spinal fusion surgery with a diagnosis of adult degenerative lumbar disease. The HU values of the anterior 1/3 of the vertebral bodies based on sagittal plane and those of vertebral bodies based on axial plane at T11-L5 were measured using the perioperative CT. The incidence of postoperative vertebral fractures with respect to the HU value was investigated.

RESULTS

Vertebral fractures were identified in 16 patients during the mean follow-up period of 3.8 years. While no significant association was found between HU value of L1 vertebral body or minimum HU value from axial plane and the incidence of the postoperative vertebral fracture, the minimum vertebral HU value of the anterior 1/3 of vertebral body from sagittal plane was associated with the incidence of the postoperative vertebral fracture. Patients with a minimum anterior 1/3 vertebral HU value of <80 had a higher incidence of postoperative vertebral fractures. The adjacent vertebral fractures occurred at the level of the vertebra with the lowest HU value, with a high probability. The existence of the vertebra with a minimum HU value of <80 within two levels of upper instrumented vertebrae was a risk factor for adjacent vertebral fracture.

CONCLUSION

HU measurement of the anterior 1/3 of vertebral body predicts the risk of vertebral fracture after short spinal fusion surgery.

Sagittal Section Hounsfield Units of the Upper Instrumented Vertebrae as a Predictor of Proximal Junctional Vertebral Fractures Following Adult Spinal Deformity Surgery

STUDY DESIGN

A retrospective observational study.

PURPOSE

This study aimed to determine an accurate and convenient screening method for predicting proximal junctional fractures (PJFr) following surgery for adult spinal deformity (ASD) using computed tomography (CT)-based measurement of Hounsfield units (HUs).

OVERVIEW OF LITERATURE

CT-based measurement of HUs is an alternative tool for assessing bone mineral density. However, the optimal method for predicting adjacent vertebral fractures following spinal fusion using HUs remains unclear.

METHODS

This retrospective observational study included 42 patients who underwent reconstructive surgery for ASD. Elliptical regions of interest (ROIs) on the axial section and rectangular ROIs on the sagittal section were placed at the upper instrumented vertebrae (UIV), UIV+1, and UIV+2. In addition, the HU value of the L2 vertebra was used as the representative.

RESULTS

PJFr occurred in 28.6% of patients within 2 years following surgery. The HU values obtained from the axial sections of L2, UIV, UIV+1, and UIV+2 were not significantly associated with the incidence of PJFr within 2 years, except for the ROI set in the lower region of the L2 vertebra. However, the HU value of the anterior third of the UIV in the sagittal section was significantly lower in the PJFr group than in the nonPJFr group (87.0 vs. 160.3, p =0.001). A UIV HU value of <100 was associated with a higher incidence of PJFr than an HU vaue of >100 (p <0.05).

CONCLUSIONS

Measurements of HU in the anterior one-third of the UIV in the sagittal section demonstrated predictive ability for PJFr following ASD surgery. A UIV HU value of <100 emerged as a risk factor for PJFr.

Low Hounsfield unit values on computed tomography as a potential predictor of vertebral fracture in patients with rheumatoid arthritis: The KURAMA cohort study

OBJECTIVE

Hounsfield units (HU) measured using computed tomography (CT) have gained considerable attention for the detection of osteoporosis. This study aimed to investigate whether opportunistic CT could predict vertebral fractures in patients with rheumatoid arthritis (RA).

METHODS

A total of 233 patients with RA who underwent chest CT were included in this study. The HU values of the anterior 1/3 of the vertebral bodies based on the sagittal plane at T11-L2 after reconstruction were measured. The incidence of vertebral fractures was investigated with respect to the HU value.

RESULTS

Vertebral fractures were identified in 32 patients during a mean follow-up period of 3.8 years. In patients who experienced vertebral fractures within 2 years of CT imaging, the HU values of the vertebral bodies (T11-L2) were lower than those in patients who did not experience fractures. Receiver operating characteristic curve analysis identified that a T11 HU value of <125 was a risk factor for vertebral fracture within 2 years. Multivariate analysis showed that a T11 HU value of <125 and the existence of prevalent vertebral fractures were significant risk factors for fracture.

CONCLUSION

HU measurements of the anterior 1/3 of the vertebral body are a potential predictor for vertebral fractures in patients with RA.

Vertebral Fracture Risk Thresholds from Phantom-Less Quantitative Computed Tomography-Based Finite Element Modeling Correlate to Phantom-Based Outcomes

INTRODUCTION

Osteoporosis indicates weakened bones and heightened fracture susceptibility due to diminished bone quality. Dual-energy x-ray absorptiometry is unable to assess bone strength. Volumetric bone mineral density (vBMD) from quantitative computed tomography (QCT) has been used to establish guidelines as equivalent measurements for osteoporosis. QCT-based finite element analysis (FEA) has been implemented using calibration phantoms to establish bone strength thresholds based on the established vBMD. The primary aim was to validate vertebral failure load thresholds using a phantom-less approach with previously established thresholds, advancing a phantom-free approach for fracture risk prediction.

METHODOLOGY

A controlled cohort of 108 subjects (68 females) was used to validate sex-specific vertebral fracture load thresholds for normal, osteopenic, and osteoporotic subjects, obtained using a QCT/FEA-based phantom-less calibration approach and two material equations.

RESULTS

There were strong prediction correlations between the phantom-less and phantom-based methods (R2: 0.95 and 0.97 for males, and R2: 0.96 and 0.98 for females) based on the two equations. Bland Altman plots and paired t-tests showed no significant differences between methods. Predictions for bone strengths and thresholds using the phantom-less method matched those obtained using the phantom calibration and those previously established, with ≤4500 N (fragile) and ≥6000 N (normal) bone strength in females, and ≤6500 N (fragile) and ≥8500 N (normal) bone strength in males.

CONCLUSION

Phantom-less QCT-based FEA can allow for prospective and retrospective studies evaluating incidental vertebral fracture risk along the spine and their association with spine curvature and/or fracture etiology. The findings of this study further supported the application of phantom-less QCT-based FEA modeling to predict vertebral strength, aiding in identifying individuals prone to fractures. This reinforces the rationale for adopting this method as a comprehensive approach in predicting and managing fracture risk.

New Insights in the Pathophysiology, Epidemiology, and Response to Treatment of Osteoporotic Vertebral Fractures

CONTEXT

Vertebral fractures (VFs) make up an important but challenging group of fractures often caused by osteoporosis. Osteoporotic fractures pose unique diagnostic challenges in generally requiring imaging for diagnosis. The objective of this narrative mini-review is to provide an overview of these recent advances in our knowledge of VF pathophysiology and epidemiology with particular focus on endocrine diseases, prevention, and treatment.

EVIDENCE ACQUISITION

We searched PubMed on May 23, 2022, for studies of VFs in humans. Results were limited to papers available as full-text publications in English, published from 2020 and onward. This yielded 3457 citations. This was supplemented by earlier publications selected to add context to the recent findings.

EVIDENCE SYNTHESIS

Studies addressed VF risk in hyperthyreosis, hyperparathyroidism, acromegaly, Cushing syndrome, primary aldosteronism, and diabetes. For pharmaceutical treatment, new studies or analyses were identified for romosozumab and for weekly teriparatide. Several studies, including studies in the immediate pipeline, were intervention studies with vertebroplasty or kyphoplasty, including combination with stem cells or pharmaceuticals.

CONCLUSIONS

Endocrinologists should be aware of the high likelihood of osteoporotic VFs in patients with endocrine diseases. Though licensed treatments are able to substantially reduce the occurrence of VFs in patients with osteoporosis, the vast majority of recent or ongoing randomized controlled trials in the VF area focus on advanced invasive therapy of the fracture itself.

Differential Cortical Volumetric Bone Mineral Density within the Human Rib

INTRODUCTION

The human rib provides a vital role in the protection of thoracic contents. Rib fractures are linked to injuries and health complications that can be fatal. Current clinical methods to assess fracture risk and bone quality are insufficient to quantify intra-element differences in bone mineral density (BMD) or to identify at-risk populations. Utilizing quantitative computed tomography (QCT) provides accurate measures of volumetric BMD (vBMD) along the length of the rib which can help delineate factors influencing differential fracture risk.

METHODOLOGY

One mid-level rib was obtained from 54 post-mortem human subjects (PMHS) and scanned using QCT. Volumes of interest (VOIs) were created for sites at 30%, 50%, and 75% of rib total curve length. Mean Hounsfield units (HU) from each VOI were converted to vBMD using a scan-specific cortical phantom calibration curve. Additionally, rib and lumbar areal BMD (aBMD) were obtained from a sub-sample of 33 PMHS.

RESULTS

Significant differences in vBMD were found between all sites within the rib (p<0.01). When analyzed by sex, vBMD between the 30% and 50% site were no longer different in either males or females (p>0.05). Separating the sample into discrete age groups demonstrated the relative differences in vBMD between sites diminished with age. Further, age as a continuous variable significantly predicted rib vBMD at all sites (p<0.05), but with little practical or clinical utility (R², 14.7- 22.8%). Similarly, only small amounts of variation in rib vBMD were explained from DXA lumbar and rib aBMD (R² , 1.1-21.8%).

CONCLUSIONS

vBMD significantly decreased from the posterior (30%) site to the anterior (75%) site within the rib which may represent adaptation to localized mechanical loading. These differences could result in differential fracture risk across the rib. As thoracic injury can be fatal, using comprehensive assessments of bone quality that accounts for variation within the rib may provide more accurate identification of at-risk populations.

Cortex or cancellous-which is early for the decrease of bone content for vertebral body in health?

OBJECTIVE

To obtain the cortex and cancellous parameter of the vertebral bone of healthy subjects using QCT. To explore which is earlier or faster for bone loss with age.

MATERIALS AND METHODS

733 physical examiners underwent chest low-dose CT examination were recruited, from April 1, 2021 to October 1, 2021. QCT sequence was used to obtain the bone mineral density of T12-L2 vertebral body without additional radiation. The mass and area of vertebral cortex and cancellous at the central level of L2 vertebral body were measured. The age -related characteristics of vertebral cortex and cancellous between male and female was analyzed and compared.

RESULTS

The vBMD of T12-L2 vertebral body decreased with age. Significant differences were found in volumetric bone mineral density (vBMD) of T12-L2 vertebral body. For female, significant differences were found in bone content involving cortical mass, cancellous mass, cortical area, cancellous area, cortical mass/cancellous mass and cortical area/cancellous area in different age groups, respectively. The cortical mass decreased with age in female. The cancellous mass of female increased and then decreased with peak at 31-40 y. The cortical area of female decreased gradually before 71 y. The cancellous area of female increased and then decreased with peak at 51-60 y. The values of mass ratio and area ratio in female showed a slowly downward trend with age. Significant differences of bone content between non-menopausal and menopausal women were found except the cancellous mass. For male, no significant differences were found in all parameters of bone content.

CONCLUSION

The changes of vertebral BMD, bone content of cortex and cancellous have different characteristics in different age. The change of cortex in female maybe earlier and faster than that of cancellous, especially in menopausal women.

A Three-Dimensional Cement Quantification Method for Decision Prediction of Vertebral Recompression after Vertebroplasty

Objective

Proposing parameters to quantify cement distribution and increasing accuracy for decision prediction of vertebroplasty postoperative complication.

Methods

Finite element analysis was used to biomechanically assess vertebral mechanics (n = 51) after percutaneous vertebroplasty (PVP) or kyphoplasty (PKP). The vertebral space was divided into 27 portions. The numbers of cement occupied portions and numbers of cement-endplate contact portions were defined as overall distribution number (oDN) and overall endplate contact number (oEP), respectively. And cement distribution was parametrized by oDN and oEP. The determination coefficients of vertebral mechanics and parameters (R 2) can validate the correlation of proposed parameters with vertebral mechanics.

Results

oDN and oEP were mainly correlated with failure load (R 2 = 0.729) and stiffness (R 2 = 0.684), respectively. oDN, oEP, failure load, and stiffness had obvious difference between the PVP group and the PKP group (P < 0.05). The regional endplate contact number in the front column is most correlated with vertebral stiffness (R 2 = 0.59) among all regional parameters. Cement volume and volume fraction are not dominant factors of vertebral augmentation, and they are not suitable for postoperative fracture risk prediction.

Conclusions

Proposed parameters with high correlation on vertebral mechanics are promising for clinical utility. The oDN and oEP can strongly affect augmented vertebral mechanics thus is suitable for postoperative fracture risk prediction. The parameters are beneficial for decision-making process of revision surgery necessity. Parametrized methods are also favorable for surgeon's preoperative planning. The methods can be inspirational for clinical image recognition development and auxiliary diagnosis.

Focal osteoporosis defect is associated with vertebral compression fracture prevalence in a bone mineral density-independent manner

Introduction

Focal osteoporosis defect has shown a high association with the bone fragility and osteoporotic fracture prevalence. However, no routine computed tomography (CT)-based vertebral focal osteoporosis defect measurement and its association with vertebral compression fracture (VCF) were discussed yet. This study aimed to develop a routine CT-based measurement method for focal osteoporosis defect quantification, and to assess its association with the VCF prevalence.

Materials and Methods

A total of 205 cases who underwent routine CT scanning, were retrospectively reviewed and enrolled into either the VCF or the control group. The focal bone mineral content loss (focal BMC loss), measured as the cumulated demineralization within bone void space, was proposed for focal osteoporosis defect quantification. Its scan-rescan reproducibility and its correlation with trabecular bone mineral density (BMD) and apparent microarchitecture parameters were evaluated. The association between focal BMC loss and the prevalence of VCF was studied by logistic regression.

Results

The measurement of focal BMC loss showed high reproducibility (RMSSD = 0.011 mm, LSC = 0.030 mm, ICC = 0.97), and good correlation with focal bone volume fraction (r = 0.79, P < 0.001), trabecular bone separation (r = 0.76, P < 0.001), but poor correlation with trabecular BMD (r = 0.37, P < 0.001). The focal BMC loss was significantly higher in the fracture group than the control (1.03 ± 0.13 vs. 0.93 ± 0.11 mm; P < 0.001), and was associated with prevalent VCF (1.87, 95% CI = 1.31-2.65, P < 0.001) independent of trabecular BMD level.

Discussion

As a surrogate measure of focal osteoporosis defect, focal BMC Loss independently associated with the VCF prevalence. It suggests that focal osteoporosis defect is a common manifestation that positively contributed to compression fracture risk and can be quantified with routine CT using focal BMC Loss.

3D assessment of intervertebral disc degeneration in zebrafish identifies changes in bone density that prime disc disease

Back pain is a common condition with a high social impact and represents a global health burden. Intervertebral disc disease (IVDD) is one of the major causes of back pain; no therapeutics are currently available to reverse this disease. The impact of bone mineral density (BMD) on IVDD has been controversial, with some studies suggesting osteoporosis as causative for IVDD and others suggesting it as protective for IVDD. Functional studies to evaluate the influence of genetic components of BMD in IVDD could highlight opportunities for drug development and repurposing. By taking a holistic 3D approach, we established an aging zebrafish model for spontaneous IVDD. Increased BMD in aging, detected by automated computational analysis, is caused by bone deformities at the endplates. However, aged zebrafish spines showed changes in bone morphology, microstructure, mineral heterogeneity, and increased fragility that resembled osteoporosis. Elements of the discs recapitulated IVDD symptoms found in humans: the intervertebral ligament (equivalent to the annulus fibrosus) showed disorganized collagen fibers and herniation, while the disc center (nucleus pulposus equivalent) showed dehydration and cellular abnormalities. We manipulated BMD in young zebrafish by mutating sp7 and cathepsin K, leading to low and high BMD, respectively. Remarkably, we detected IVDD in both groups, demonstrating that low BMD does not protect against IVDD, and we found a strong correlation between high BMD and IVDD. Deep learning was applied to high-resolution synchrotron µCT image data to analyze osteocyte 3D lacunar distribution and morphology, revealing a role of sp7 in controlling the osteocyte lacunar 3D profile. Our findings suggest potential avenues through which bone quality can be targeted to identify beneficial therapeutics for IVDD.

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 severity of hepatic disorder is related to vertebral microstructure deterioration in cadaveric donors with liver cirrhosis

Patients with liver cirrhosis (LC) commonly suffer from osteoporosis and vertebral fracture, but data about their vertebral micro-architectural changes are still limited. This study aimed to evaluate the potential differences in trabecular micro-architecture of lumbar vertebrae between male LC patients and healthy controls, in relation to etiology and pathohistological scoring of the liver disorder. After pathohistological examination of liver tissue, micro-computed tomography was performed on the vertebral samples included into: alcoholic liver cirrhosis group (ALC; n = 16; age: 59 ± 8 years), non-alcoholic liver cirrhosis group (non-ALC; n = 15; age: 69 ± 10 years) and control group (n = 16; age: 58 ± 6 years). Our data showed significant impairment of the trabecular microstructure in the lumbar vertebrae from LC donors, regardless of the alcoholic/non-alcoholic origin of liver disorder, as illustrated by lower BV/TV, Tb.Th, and Tb.N compared with controls (p < .05). Moreover, depredation in trabecular micro-architecture was inversely associated with pathohistological scores (p < .05), indicating that severity of liver disorder could be an important predictor of reduced vertebral strength in LC. We noticed significant micro-architectural deterioration in the trabecular compartment of the lumbar vertebrae of male individuals with alcoholic and non-alcoholic LC, which was associated with the severity of the liver disease. Thus, clinical assessment of fracture risk should be advised for all LC patients, regardless of the alcoholic origin of liver cirrhosis. Additionally, adequate and timely treatment of liver disorder may decelerate the progression of bone impairment in LC patients.

The micro-structural analysis of lumbar vertebrae in alcoholic liver cirrhosis

Although vertebral fracture is more common among alcoholic liver cirrhosis patients when compared to general population, current data on three-dimensional micro-architecture are scarce. Our study showed significant trabecular deterioration in lumbar vertebrae obtained from alcoholic liver cirrhosis donors, suggesting that they should be advised to undergo early-stage screening for osteoporosis.

PURPOSE

Recent studies showed an increased incidence of vertebral fractures in alcoholic liver cirrhosis (ALC) patients, while data about vertebral micro-structure are still limited. The aim of this study was to compare trabecular and cortical micro-architecture of lumbar vertebrae between ALC patients and healthy age- and sex-matched controls.

METHODS

Our study included lumbar vertebral samples of male cadaveric donors, divided into ALC (n = 20, age: 59 ± 6 years) and control group (n = 20, age: 59 ± 8 years). Following pathohistological verification of liver cirrhosis, trabecular and cortical bone micro-architecture was analyzed by micro-computed tomography (micro-CT).

RESULTS

Micro-CT evaluation of the trabecular bone in lumbar vertebrae showed a significant decrease in bone volume fraction, trabecular thickness, trabecular number, and connectivity (p < 0.01). In contrast to trabecular deterioration, prominent alteration in cortical parameters was not observed in lumbar vertebrae of ALC patients (p > 0.05).

CONCLUSIONS

Our data indicate that susceptibility to non-traumatic fractures in ALC patients could be explained by alterations in trabecular bone micro-architecture. Thus, we genuinely recommend osteological screening of the lumbar spine for all ALC patients in order to evaluate individual fracture risk. Graphical abstract.