The thickness of human vertebral cortical bone and its changes in aging and osteoporosis: a histomorphometric analysis of the complete spinal column from thirty-seven autopsy specimens

Assessment of age-dependent sexual dimorphism in paediatric vertebral size and density using a statistical shape and statistical appearance modelling approach

This work focuses on the growth patterns of the human fourth lumbar vertebra (L4) in a paediatric population, with specific attention to sexual dimorphism. The study aims to understand morphological and density changes in the vertebrae through age-dependent statistical shape and statistical appearance models, which can describe full three-dimensional anatomy. Results show that the main growth patterns are associated with isotropic volumetric vertebral growth, a decrease in the relative size of the vertebral foramen, and an increase in the length of the transverse processes. Moreover, significant sexual dimorphism was demonstrated during puberty. We observe significant age and sex interaction in the anterior vertebral body height (P = 0.005), where females exhibited an earlier increase in rates of vertebral height evolution. Moreover, we also observe an increase in cross-sectional area (CSA) with age (P = 0.020), where the CSA is smaller in females than in males (significant sex effect P = 0.042). Finally, although no significant increase in trabecular bone density with age is observed (P = 0.363), a trend in the statistical appearance model suggests an increase in density with age.

Spine Stereotactic Body Radiation Therapy in Geriatric Patients: Implications of Age and Dose on Iatrogenic Vertebral Compression Fracture Risk

PURPOSE

Stereotactic body radiation therapy (SBRT) is an effective treatment for spinal metastases; however, outcomes specific to a geriatric population have not been described. This study aims to investigate the efficacy and safety of spine SBRT, in particular the rate of iatrogenic vertebral compression fracture (VCF), in patients aged 70 and older.

PATIENTS AND METHODS

From a prospectively maintained single-institutional database of 976 patients and 2407 spinal segments treated with SBRT for vertebral metastases between 2008 and 2021, all patients aged 70 or above were retrospectively reviewed. The primary outcome is the risk of VCF. Secondary outcomes included magnetic resonance imaging-based local failure and overall survival.

RESULTS

A total of 252 consecutive patients with 580 spinal segments treated with spine SBRT were reviewed. The median age was 75.8 (range: 70-90.3) years and the median (interquartile range) follow-up duration was 16.9 (6.4-41.3) months. The median overall survival of the entire cohort was 20.3 months and the 2-year local failure rate was 14.3%. The cumulative risk of VCF at 12 and 24 months were 8.4% and 12.3%, respectively. Significant predictors of VCF on multivariable analyses included greater biologically equivalent dose, baseline fracture, and increasing age. In particular, the 2-year VCF rate and median time to VCF were 30.3% and 3.4 months for those 86 and older, compared with 11.2% and 12.8 months for those younger than 86, respectively (P = .0011).

CONCLUSION AND RELEVANCE

Spine SBRT should be considered in a geriatric population; however, for those 86 and older, we suggest caution due to the significant risk of VCF.

A novel personalized homogenous finite element model to predict the pull-out strength of cancellous bone screws

BACKGROUND

Orthopedic surgeries often involve the insertion of bone screws for various fixation systems. The risk of postoperative screw loosening is usually assessed through experimental or finite element (FE) evaluations of the screw pull-out strength. FE simulations are based on either personalized complex but accurate heterogeneous modeling or non-personalized simple but relatively less accurate homogeneous modeling. This study aimed to develop and validate a novel personalized computed tomography (CT)-based homogeneous FE simulation approach to predict the pull-out force of cancellous bone screws.

METHODS

Twenty FE simulations of L1-L5 vertebral screw pull-out tests were conducted, i.e., 10 heterogeneous and 10 homogenous models. Screws were inserted into the lower-middle region of vertebrae. In our novel homogeneous model, the region around approximately twice the diameter of the screw was used as a bone material reference volume. Subsequently, the overall material property of this region was homogeneously attributed to the entire vertebra, and pull-out simulations were conducted.

RESULTS

The mean error of the predicted pull-out forces by our novel homogenous simulations was ~ 7.9% with respect to our heterogeneous model. When solely the cancellous bone was involved during the pull-out process (i.e., for L1, L2, and L3 vertebral bodies whose cortical bone in the inferior region is thin), the novel homogenous model yielded small mean error of < 6.0%. This error, however, increased to ~ 11% when the screw got involved to the cortical bone (for L4 and L5 vertebrae whose cortical bone in the inferior region is thick).

CONCLUSION

The proposed personalized CT-based homogenous model was highly accurate in estimating the pull-out force especially when only the cancellous bone was involved with the screw.

Biomechanical Analysis of Diffuse Idiopathic Skeletal Hyperostosis and Osteoporosis: Vertebral Fracture Risk Evaluation Using Finite Element Modeling with Clinical and Micro-CT Data in an Elderly Donor

Introduction: Biomechanical analysis of spinal structures is crucial in the evaluation of injuries, the risk of fracture, and age-related changes. Osteoporotic vertebrae are very fragile and therefore constitute a serious risk, especially in the elderly. Methods: At present, clinically relevant decision making in fracture risk assessment is predicated upon finite element analysis (FEA), which utilizes high-resolution computed tomography (CT) scans from clinical practice alongside micro-CT scans from laboratory settings. Of particular interest is the utilization of cortical vertebral body thicknesses, as meticulously measured via micro-CT. The data from a body donation over 80 years old with diffuse idiopathic skeletal hyperostosis (DISH) and osteoporosis (OP) were utilized through FEA to evaluate stresses in cortical and trabecular bone and to predict the stiffness and deformability of the examined vertebral bodies. Results: The investigation revealed a higher density of cortical and cancellous bone in vertebrae affected by DISH. Cortical density was highest in the thoracic section (median 188 µm), while cancellous bone density was 222 HU in the cervical vertebrae. The load on cortical bone increased as we progressed towards the lumbar spine; however, it remained quite constant in cancellous bone. Despite a low bone density, we registered no fractures in vertebrae. Conclusions: The data showed that DISH increased the thickness of the cortical bone and the density of the cancellous bone. The combination of increased cortical and cancellous bone density might reduce the risk of fracture in patients with low bone density. These conclusions emphasize the significance of biomechanical properties in the assessment of fracture risk and have important implications for clinical practice, particularly in relation to the prevention of vertebral fractures in osteoporotic patients with DISH.

Biomechanical Analysis of Lumbar Interbody Fusion Cages With Various Elastic Moduli in Osteoporotic and Non-osteoporotic Lumbar Spine: A Finite Element Analysis

STUDY DESIGN

Finite element analysis (FEA).

OBJECTIVE

This study aimed to explore the effects of cage elastic modulus (Cage-E) on the endplate stress in different bone conditions: osteoporosis (OP) and non-osteoporosis (non-OP). We also explored the correlation between endplate thickness and endplate stress.

METHODS

The FEA models of L4-L5 with lumbar interbody fusion were designed to access the effects of Cage-E on the endplate stress in different bone conditions. Two groups of the Young's moduli of bony structure were assigned to simulate the conditions of OP and non-OP, and the bony endplates were analyzed in 2 kinds of thicknesses: .5 mm and 1.0 mm, with the insertion of cages with different Young's moduli including .5, 1.5, 3, 5, 10, and 20 GPa. After model validation, an axial compressive load of 400 N and a flexion/extension moment of 7.5Nm was performed on the superior surface of L4 vertebral body in order to analyze the distribution of stress.

RESULTS

The maximum Von Mises stress in the endplates increased by up to 100% in the OP model compared with non-OP model under the same condition of cage-E and endplate thickness. In both OP and non-OP models, the maximum endplate stress decreased as the cage-E decreased, but the maximum stress in the lumbar posterior fixation increased as the cage-E decreased. Thinner endplate thickness was associated with increased endplate stress.

CONCLUSION

The endplate stress is higher in osteoporotic bone than non-osteoporotic bone, which explains part of the mechanism of OP-related cage subsidence. It is reasonable to reduce the endplate stress by reducing the cage-E, but we should balance the risk of fixation failure. Endplate thickness is also important when evaluating the cage subsidence risk.

Differences in vertebral bone density between African apes

OBJECTIVES

Low-energy vertebral fractures are a common health concern, especially in elderly people. Interestingly, African apes do not seem to experience as many vertebral fractures and the low-energy ones are even rarer. One potential explanation for this difference is the lower bone density in humans. Yet, only limited research has been done on the vertebral bone density of the great apes and these have mainly included only single vertebrae. Hence the study aim is to expand our understanding of the vertebral microstructure of African apes in multiple spinal segments.

MATERIALS

Bone density in the vertebral body of C7, T12, and L3 was measured from 32 Pan troglodytes and 26 Gorilla gorilla using peripheral quantitative computed tomography (pQCT).

RESULTS

There was a clear difference between the three individual vertebrae and consequently the spinal segments in terms of trabecular density and cortical density and thickness. The variation of these bone parameters between the vertebrae differed between the apes but was also different from those reported for humans. The chimpanzees were observed to have overall higher trabecular density, but gorillas had higher cortical density and thickness. Cortical thickness had a relatively strong association with the vertebral size.

DISCUSSION

Despite the similarity in locomotion and posture, the results show slight differences in the bone parameters and their variation between spinal segments in African apes. This variation also differs from humans and appears to indicate a complex influence of locomotion, posture, and body size on the different spinal segments.

Preoperative Assessment of Bone Density Using MRI-Based Vertebral Bone Quality Score Modified for Patients Undergoing Cervical Spine Surgery

STUDY DESIGN

Diagnostic accuracy study.

OBJECTIVES

Previous studies have reported the clinical application of the Vertebral Bone Quality (VBQ) scores for assessing bone density in operative lumbar spine patients. We aim to explore whether the method could be modified and applicable for patients undergoing cervical spine surgery.

METHODS

Adult patients receiving cervical spine surgery for degenerative diseases between September 2020 and March 2022 with non-contrast T1-weighted MRI and DEXA were included. Correlation between cervical VBQ scores and DEXA T-scores was analyzed using Pearson's correlation. Student's t test was used to present the discrepancy between the VBQ of patients with normal bone density (T ≥ -1.0) and patients with osteopenia/osteoporosis (T < -1.0). Statistical significance was set at P < .05.

RESULTS

Eighty-three patients (20 patients with T ≥ -1.0 vs 63 patients with T < -1.0 ) were included. Significant difference was found between the cervical VBQ between groups (2.99 ± .79 vs 3.80 ± .81, P < .001). Interclass correlation coefficient for inter-rater reliability was .82 (95% CI: .70-.93) and .91(95% CI: .84-.97) for intra-rater reliability. The area under the ROC curve was .78 (95% CI: .65-.90). The DEXA T-score of the femoral neck, total hip and the lowest DEXA T-score were found to be significantly correlated with the cervical VBQ score according to Pearson correlation analysis (P < .001).

CONCLUSIONS

This is the first study to apply the VBQ method to assess the bone density in preoperative cervical spine patients. Cervical VBQ scores were significantly correlated with DEXA T-score. With an overall accuracy of .78, the radiation-free and cost-effective method could be a potential tool for screening patients with osteopenia and osteoporosis before surgery.

An additively manufactured model for preclinical testing of cervical devices

Purpose

Composite models have become commonplace for the assessment of fixation and stability of total joint replacements; however, there are no comparable models for the cervical spine to evaluate fixation. The goal of this study was to create the framework for a tunable non-homogeneous model of cervical vertebral body by identifying the relationships between strength, in-fill density, and lattice structure and creating a final architectural framework for specific strengths to be applied to the model.

Methods

The range of material properties for cervical spine were identified from literature. Using additive manufacturing software, rectangular prints with three lattice structures, gyroid, triangle, zig-zag, and a range of in-fill densities were 3D-printed. The compressive and shear strengths for all combinations were calculated in the axial and coronal planes. Eleven unique vertebral regions were selected to represent the distribution of density. Each bone density was converted to strength and subsequently correlated to the lattice structure and in-fill density with the desired material properties. Finally, a complete cervical vertebra model was 3D-printed to ensure sufficient print quality.

Results

Materials testing identified a relationship between in-fill densities and strength for all lattice structures. The axial compressive strength of the gyroid specimens ranged from 1.5 MPa at 10% infill to 31.3 MPa at 100% infill and the triangle structure ranged from 2.7 MPa at 10% infill to 58.4 MPa at 100% infill. Based on these results, a cervical vertebra model was created utilizing cervical cancellous strength values and the corresponding in-fill density and lattice structure combination. This model was then printed with 11 different in-fill densities ranging from 33% gyroid to 84% triangle to ensure successful integration of the non-homogeneous in-fill densities and lattice structures.

Conclusions

The findings from this study introduced a framework for using additive manufacturing to create a tunable, customizable biomimetic model of a cervical vertebra.

Treatment of osteoporosis using a selective androgen receptor modulator ostarine in an orchiectomized rat model

PURPOSE

The selective androgen receptor modulator ostarine has been shown to have advantageous effects on skeletal tissue properties, reducing muscle wasting and improving physical function in males. However, data on effects in male osteoporosis remain limited. In this study, the effects of ostarine on osteoporotic bone were evaluated in a rat model of male osteoporosis and compared with those of testosterone treatments.

METHODS

Eight-month-old male Sprague-Dawley rats were either non-orchiectomized to serve as a healthy control (Non-Orx, Group 1) or orchiectomized (Orx, Groups 2-6) and then grouped (n = 15/group): (1) Non-Orx, (2) Orx, (3) Ostarine Therapy, (4) Testosterone Therapy, (5) Ostarine Prophylaxis and (6) Testosterone Prophylaxis. Prophylaxis treatments started directly after orchiectomy and continued for 18 weeks, whereas Therapy treatments were initiated 12 weeks after Orx. Ostarine and Testosterone were applied orally at daily doses of 0.4 and 50 mg/kg body weight, respectively. The lumbar vertebral bodies and femora were analyzed using biomechanical, micro-CT, ashing, and gene expression analyses.

RESULTS

Ostarine Prophylaxis showed positive effects in preventing osteoporotic changes in cortical and trabecular bone (femoral trabecular density: 26.01 ± 9.1% vs. 20.75 ± 1.2% in Orx and in L4: 16.3 ± 7.3% vs 11.8 ± 2.9% in Orx); biomechanical parameters were not affected; prostate weight was increased (0.62 ± 0.13 g vs 0.18 ± 0.07 g in Orx). Ostarine Therapy increased solely the cortical density of the femur (1.25 ± 0.03 g/cm3 vs. 1.18 ± 0.04 g/cm3 in Orx); other bone parameters remained unaffected. Testosteron Prophylaxis positively influenced cortical density in femur (1.24 ± 0.05 g/cm3 vs. 1.18 ± 0.04 g/cm3 in Orx); Test. Therapy did not change any bony parameters.

CONCLUSION

Ostarine Prophylaxis could be further investigated as a preventative treatment for male osteoporosis, but an androgenic effect on the prostate should be taken into consideration, and combination therapies with other anti-osteoporosis agents could be considered.

Systematic measuring cortical thickness in tibiae for bio-mechanical analysis

BACKGROUND AND OBJECTIVE

Measuring the thickness of cortical bone tissue helps diagnose bone diseases or monitor the progress of different treatments. This type of measurement can be performed visually from CAT images by a radiologist or by semi-automatic algorithms from Hounsfield values. This article proposes a mechanism capable of measuring thickness over the entire bone surface, aligning and orienting all the images in the same direction to have comparable references and reduce human intervention to a minimum. The objective is to batch process large numbers of patients' CAT images obtaining thicknesses profiles of their cortical tissue to be used in many applications.

METHODS

Classical morphological and Deep Learning segmentation is used to extract the area of interest, filtering and interpolation to clean the bones and contour detection and Signed Distance Functions to measure the cortical Thickness. The alignment of the set of bones is achieved by detecting their longitudinal direction, and the orientation is performed by computing their principal component of the center of mass slice.

RESULTS

The method processed in an unattended manner 67% of the patients in the first run and 100% in the second run. The difference in the thickness values between the values provided by the algorithm and the measures done by a radiologist was, on average, 0.25 millimetres with a standard deviation of 0.2.

CONCLUSION

Measuring the cortical thickness of a bone would allow us to prepare accurate traumatological surgeries or study their structural properties. Obtaining thickness profiles of an extensive set of patients opens the way for numerous studies to be carried out to find patterns between bone thickness and the patients' medical, social or demographic variables.

Analysis of stress and stabilization in adolescent with osteoporotic idiopathic scoliosis: finite element method

Objective: To explore the effect of osteoporosis on the stress, stability, and lumbar intervertebral disc of AIS lumbar vertebrae by finite element method. Better understand the biomechanical characteristics of osteoporotic scoliosis.Methods: Based on the CT images of normal lumbar vertebrae and lumbar vertebrae with AIS, the finite element models were established to simulate the estimated osteoporosis by changing the Young's modulus of cortical bone, cancellous bone, and endplate. Four finite element models of normal lumbar, osteoporotic lumbar, normal AIS lumbar and osteoporotic AIS lumbar were established, and the same load and boundary conditions were applied respectively. The displacement, stress, and intervertebral disc strain of the four models were compared to explore the effect of osteoporosis on the stability and injury risk of AIS.Results: After suffering from osteoporosis, under the same load, the displacement of lumbar spine increases, the stability decreases, and the stability of AIS lumbar spine decrease more obviously, especially under extension load. Suffering from osteoporosis will increase the stress of lumbar spine, AIS lumbar spine increases more obviously, and the stress is more concentrated, Osteoporotic lumbar spine only affects the strain of intervertebral disc when AIS lumbar spine bends on the concave side, resulting in greater strain behind the concave side of intervertebral disc.Conclusions: AIS patients with OP have lower lumbar stability, a higher risk of fracture of lumbar vertebrae, and spinal nerves are more likely to be compressed by intervertebral discs. OP can aggravate the scoliosis of lumbar vertebrae.

The Effect of Concave-Side Intertransverse Ligament Laxity on the Stress of AIS Lumbar Spine Based on Finite Element Method

(1) Background: Scoliosis has the mechanical characteristic of asymmetric stress distribution, which is one of the reasons for the aggravation of scoliosis. Bracing therapy is the best treatment for AIS, but it is difficult and costly to operate. Is it possible to reduce pressure in the concave side by relaxing the ITL in the concave side of scoliosis, so as to improve the abnormal stress distribution of scoliosis? In this paper, a finite element method was used to simulate the effect of the relaxation of concave-side ITL on the stress of a lumbar spine with scoliosis, which provides some guidance for the treatment of scoliosis. (2) Methods: Using CT images of a patient with scoliosis whose Cobb Angle was 43° and Lordosis Angle was 45, a scoliosis lumbar was established, and Young's modulus of the ITL of the concave-side lumbar spine was reduced by 95% to simulate ligament relaxation. By comparing the stress condition of the model vertebral body with no ligament relaxation, the effect of concave-side ITL relaxation on the mechanical characteristics of scoliosis lumbar spine was explored. (3) Results: An effective and complete model of the lumbar spine was established. The concave ITL relaxed, which only had a great impact on the bending loads. After the ligament was relaxed, the stability of the spine was reduced. Stress concentration on the concave side of vertebrae and the IVD was aggravated. Under loads on the convex side, the maximum stress on the vertebral body and the IVD increased significantly, making lumbar vertebrae more vulnerable to injury. (4) Conclusions: Laxity of the ITL on the concave side of the AIS lumbar only affects the bending load. Laxity of the concave-side ligament will reduce the stability of the lumbar, aggravate the uneven stress distribution of scoliotic lumbar vertebrae, increase the risk of IVD injury, and be unfavorable for the scoliotic lumbar spine. Relaxation of the concave ITL alone is not an effective way to treat scoliosis.

Patient-specific numerical investigation of the correction of cervical kyphotic deformity based on a retrospective clinical case

Little research has been reported on evaluating the safety of the fixation construct in cervical kyphosis correction. In this study, we proposed a principal-strain criterion to evaluate the safety of the fixation construct and validated the modeling method against a retrospective case of anterior cervical discectomy fusion (ACDF). From C2 to T2 vertebra bodies, fixation instruments were reconstructed and positioned as per postoperative computed tomography (CT) scans. Head weight (HW) and various moments estimated from isometric strength data were imposed onto the C2. The postoperative stability of non-surgical segments, deformations surrounding the screw trajectories, and contact slipping on zygapophysial joints were analyzed. The model was validated against the reality that the patient had a good fusion and deformity correction. The ACDF restricted the range of motions (ROMs) of cervical segments and lent stability to vertebra fusion, no failure was found in the finite element (FE) model of cervical vertebrae. The deformation surrounding the screw trajectories were concentrated to the lateral sides of trajectories, recommending that the shape of the anterior cervical plate conforming to the curvature of the vertebra and screws fully inserted into vertebrae reduced the deformation concentration around the screw trajectories.

Subsidence following cervical discectomy and implant-to-bone ratio

Background

Implant subsidence is an undesirable effect after anterior cervical discectomy and fusion (ACDF). We investigated the relation between the rate of implant subsidence and the ratio of the implant surface area to the surface area of the adjacent bone.

Methods

We operated 170 disc spaces in a group of 104 patients. Two types of implants were used: 1) PEEK (polyetheretherketone) cages and 2) titanium-coated (TC) PEEK cages. Patients were randomised to receive a specific implant using a randomisation table. All implants had a surface area of 1.61 cm². Based on computed tomography images, bone surface areas were calculated for vertebral bodies immediately adjacent to the interbody implants. The implant-to-bone surface ratio was then calculated for each disc space. Implant subsidence was assessed over 12 months of follow-up, and associations between implant subsidence, the type of implant, and the implant-to-bone surface ratio were investigated.

Results

Twelve months after the surgery, computed tomography was performed on 86 patients (144 disc spaces). Furthermore, in 166 disc spaces and 102 patients, conventional radiographs were obtained.

Subsidence was observed in 21% of the examined intervertebral spaces, and it was more frequently associated with higher values of bone surface area and lower values of the implant-to-bone surface ratio. The type of implant (PEEK vs TC-PEEK cages) did not significantly influence the rate of implant subsidence.

Conclusions

Implant subsidence was significantly related to the value of a coefficient representing the ratio of the implant's surface area to the bone surface area of the adjacent vertebral bodies, with subsidence occurring significantly more rarely for coefficient values ≥ 0.37.

The effects of pre-stressed rods contoured by different bending techniques on posterior instrumentation of L4-L5 lumbar spine segment: A finite element study

Posterior pedicle screw instrumentation (PPSI) is a well-known method in lumbar spine surgery. Understanding how PPSI affects the biomechanics of the lumbar spine is an important issue. In particular, during PPSI operations, surgeons bend rods according to the patients’ spinal curvatures based on their own experiences. As a result, residual stresses develop on the rods due to this bending. Although many finite element-based biomechanical studies have been performed for PPSI, studies comparing the effects of residual stresses arising from rod contouring on the construct stresses are lacking. Thus, this study aimed to investigate the effects of residual stress in PPSI using rods contoured with a French bender and an in-situ bender, as well as comparing the differences in stress increment with straight and contoured rods for different physiological motions. Accordingly, a finite element (FE) model of the L4-L5 lumbar spine segment was developed for PPSI and the effects of residual stresses on rods were investigated by using different bending methods. In the simulations, it was found that rods contoured with a French bender with residual stress resulted in significantly more increased stress in PPSI compared to those contoured with an in-situ bender. The results of this study emphasize that increased stress in PPSI due to the residual stresses for different physiological motions may increase the risk of PPSI failures. Additionally, the finite element modeling approach employed here could be used as a fundamental tool for further investigations of topics such as PPSI fatigue life and failure studies.

Does Graft Position Affect Subsidence After Anterior Cervical Discectomy and Fusion?

STUDY DESIGN

Retrospective cohort study.

OBJECTIVE

Implant subsidence is an important prognostic factor after anterior cervical discectomy and fusion (ACDF). Our purpose in this study was to investigate whether graft position affects subsidence after ACDF and to determine if there is a difference in clinical results based on allograft subsidence and position.

METHODS

We reviewed 92 patients who underwent single-level ACDF with allograft and plate between January 2012 and October 2018. Treatment levels were divided based on allograft position within 2 mm of the posterior margin of the augmented plate (Anterior group) or at greater than 2 mm (Center group). Subsidence was defined as segmental vertebral body height decrease of 2 mm or more at 1 year compared to 1 week after surgery.

RESULTS

Overall subsidence prevalence was 15 (16%) cases. Subsidence was 11% in the Anterior group (8/73) and 39% in the Center group (7/19; P = .012). The subsidence group showed smaller graft footprint size (graft/endplate ratio) compared with the nonsubsidence group, and pseudarthrosis occurred frequently in the subsidence. There was no significant difference in clinical results according to graft position. These findings indicate that anterior graft position reduces risk of subsidence (95% confidence interval = 0.085-0.949).

CONCLUSION

To minimize risk of subsidence, the graft should be positioned on the anterior position of the surgical-level endplate. In addition, using a graft appropriate for endplate size will further reduce subsidence.

Measuring the thickness of vertebral endplate and shell using digital tomosynthesis

The vertebral endplate and cortical shell play an important structural role and contribute to the overall strength of the vertebral body, are at highest risk of initial failure, and are involved in degenerative disease of the spine. The ability to accurately measure the thickness of these structures is therefore important, even if difficult due to relatively low resolution clinical imaging. We posit that digital tomosynthesis (DTS) may be a suitable imaging modality for measurement of endplate and cortical shell thickness owing to the ability to reconstruct multiplanar images with good spatial resolution at low radiation dose. In this study, for 25 cadaveric L1 vertebrae, average and standard deviation of endplate and cortical shell thickness were measured using images from DTS and microcomputed tomography (μCT). For endplate thickness measurements, significant correlations between DTS and μCT were found for all variables when comparing thicknesses measured in both the overall endplate volume (R² = 0.25-0.54) and when measurements were limited to a central range of coronal or sagittal slices (R² = 0.24-0.62). When compared to reference values from the overall shell volume, DTS thickness measurements were generally nonsignificant. However, when measurement of cortical shell thickness was limited to a range of central slices, DTS outcomes were significantly correlated with reference values for both sagittal and coronal central regions (R² = 0.21-0.49). DTS may therefore offer a means for measurement of endplate thickness and, within a limited sagittal or coronal measurement volume, for measurement of cortical shell thickness.

Comparison of L1 CT-attenuation and cortical thickness in predicting osteoporosis by opportunistic CT

BACKGROUND

In vertebrae, the amount of cortical bone has been estimated at 30-60%, but 45-75% of axial load on a vertebral body is borne by cortical bone.

OBJECTIVE

To compare the role of L1 CT-attenuation and cortical thickness in predicting osteoporosis by opportunistic CT and explore cortical thickness value in osteoporosis.

METHODS

We collected data of 94 patients who underwent DXA and thoracic and/or abdominal CT to demonstrate an entire L1 for other indications in routine practice. Patients were divided into three groups according to T-score: osteoporosis, osteopenia, or normal. CT-attenuation value and cortical thickness of L1 were measured. ANOVA analysis was utilized to analyze CT-attenuation and cortical thickness among the three groups. Sensitivity, specificity, and area under the curve (AUC) predicting low BMD were determined using ROC. Pearson correlations were employed to describe relationship between L1 BMD and CT-attenuation value, BMD, as well as cortical thickness.

RESULTS

The mean cortical thickness was 0.83±0.11, 0.72±0.10, and 0.64±0.09 mm for normal, osteopenia, and osteoporotic subgroups, respectively. A statistically significant difference was observed in cortical thickness and CT-attenuation value among these three subgroups. A mean CT-attenuation value threshold of > 148.7 yielded 73.0% sensitivity and 86.0% specificity for distinguishing low BMD from normal with an AUC = 0.83. Pearson correlation analysis indicated that BMD was positively correlated with CT-attenuation (r = 0.666, P < 0.001) and cortical thickness (r = 0.604, P < 0.001).

CONCLUSIONS

L1 CT-attenuation and cortical thickness measured on opportunistic CT can help predict osteoporosis. Compared with cortical thickness, CT-attenuation is a more sensitive and accurate index for distinguishing low BMD from normal.

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to ^89,90Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of ⁹⁰Sr and ⁸⁹Sr intake to humans.

Influence of posterior pedicle screw fixation at L4-L5 level on biomechanics of the lumbar spine with and without fusion: a finite element method

Background

Posterior pedicle screw (PS) fixation, a common treatment method for widespread low-back pain problems, has many uncertain aspects including stress concentration levels, effects on adjacent segments, and relationships with physiological motions. A better understanding of how posterior PS fixation affects the biomechanics of the lumbar spine is needed. For this purpose, a finite element (FE) model of a lumbar spine with posterior PS fixation at the L4–L5 segment level was developed by partially removing facet joints (FJs) to imitate an actual surgical procedure. This FE study aimed to investigate the influence of the posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion by determining which physiological motions have the most increase in posterior instrumentation (PI) stresses and FJ loading.

Results

It was determined that posterior PS fixation increased FJ loading by approximately 35% and 23% at the L3–L4 adjacent level with extension and lateral bending motion, respectively. This increase in FJ loading at the adjacent level could point to the possibility that adjacent segment disease has developed or progressed after posterior lumbar interbody fusion. Furthermore, analyses of peak von Mises stresses on PI showed that the maximum PI stresses of 272.1 MPa and 263.7 MPa occurred in lateral bending and flexion motion before fusion, respectively.

Conclusions

The effects of a posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion were investigated for all physiological motions. This model could be used as a fundamental tool for further studies, providing a better understanding of the effects of posterior PS fixation by clearing up uncertain aspects.

Topological design and biomechanical evaluation for 3D printed multi-segment artificial vertebral implants

Customized spinal implants fabricated by additive manufacturing have been increasingly used clinically to restore the physiological functions. However, the mechanisms and methods about the design for the spinal implants are not clear, especially for the reconstruction of multi-segment vertebral. This study aims to develop a novel multi-objective optimization methodology based on various normal spinal activities, to design the artificial vertebral implant (AVI) with lightweight, high-strength and high-stability. The biomechanical performance for two types of AVI was analyzed and compared under different loading conditions by finite element method. These implants were manufactured via selective laser melting technology and evaluated via compressive testing. Results showed the maximum Mises stress of the optimized implant under various load cases were about 41.5% of that of the trussed implant, and below fatigue strength of 3D printed titanium materials. The optimized implant was about 2 times to trussed implant in term of the maximum compression load and compression stiffness to per unit mass, which indicated the optimized implant can meet the safety requirement. Finally, the optimized implant has been used in clinical practice and good short-term clinical outcomes were achieved. Therefore, the novel developed method provides a favorable guarantee for the design of 3D printed multi-segment artificial vertebral implants.

Maintaining Bone Health in the Lumbar Spine: Routine Activities Alone Are Not Enough

Public health organisations typically recommend a minimum amount of moderate intensity activities such as walking or cycling for two and a half hours a week, combined with some more demanding physical activity on at least 2 days a week to maintain a healthy musculoskeletal condition. For populations at risk of bone loss in the lumbar spine, these guidelines are particularly relevant. However, an understanding of how these different activities are influential in maintaining vertebral bone health is lacking. A predictive structural finite element modelling approach using a strain-driven algorithm was developed to study mechanical stimulus and bone adaptation in the lumbar spine under various physiological loading conditions. These loading conditions were obtained with a previously developed full-body musculoskeletal model for a range of daily living activities representative of a healthy lifestyle. Activities of interest for the simulations include moderate intensity activities involving limited spine movements in all directions such as, walking, stair ascent and descent, sitting down and standing up, and more demanding activities with large spine movements during reaching and lifting tasks. For a combination of moderate and more demanding activities, the finite element model predicted a trabecular and cortical bone architecture representative of a healthy vertebra. When more demanding activities were removed from the simulations, areas at risk of bone degradation were observed at all lumbar levels in the anterior part of the vertebral body, the transverse processes and the spinous process. Moderate intensity activities alone were found to be insufficient in providing a mechanical stimulus to prevent bone degradation. More demanding physical activities are essential to maintain bone health in the lumbar spine.

Evaluation of myostatin as a possible regulator and marker of skeletal muscle-cortical bone interaction in adults

INTRODUCTION

Bone mass was recently reported to be related to skeletal muscle mass in humans, and a decrease in cortical bone is a risk factor for osteoporosis. Because circulating myostatin is a factor that primarily controls muscle metabolism, this study examined the role of myostatin in bone mass-skeletal muscle mass interactions.

METHODS

The subjects were 375 middle-aged community residents with no history of osteoporosis or sarcopenia who participated in a health check-up. Cortical bone thickness and cancellous bone density were measured by ultrasonic bone densitometry in a health check-up survey. The subjects were divided into those with low cortical bone thickness (LCT) or low cancellous bone density (LBD) and those with normal values (NCT/NBD). Bone metabolism markers (TRACP-5b, etc.), skeletal muscle mass, serum myostatin levels, and lifestyle were then compared between the groups.

RESULTS

The percentage of diabetic participants, TRACP-5b, and myostatin levels were significantly higher, and the frequency of physical activity, skeletal muscle mass, grip strength, and leg strength were significantly lower in the LCT group than in the NCT group. The odds ratio (OR) of high myostatin levels in the LCT group compared with the NCT group was significant (OR 2.17) even after adjusting for related factors. Between the low cancellous bone density (LBD) and normal cancellous bone density (NBD) groups, significant differences were observed in the same items as between the LCT and NCT groups, but no significant differences were observed in skeletal muscle mass and blood myostatin levels. The myostatin level was significantly negatively correlated with cortical bone thickness and skeletal muscle mass.

CONCLUSIONS

A decrease in cortical bone thickness was associated with a decrease in skeletal muscle mass accompanied by an increase in the blood myostatin level. Blood myostatin may regulate the bone-skeletal muscle relationship and serve as a surrogate marker of bone metabolism, potentially linking muscle mass to bone structure.

The Influence of Axial Compression on the Cellular and Mechanical Function of Spinal Tissues; Emphasis on the Nucleus Pulposus and Annulus Fibrosus: A Review

In light of the correlation between chronic back pain and intervertebral disc (IVD) degeneration, this literature review seeks to illustrate the importance of the hydraulic response across the nucleus pulposus (NP)-annulus fibrosus (AF) interface, by synthesizing current information regarding injurious biomechanics of the spine, stemming from axial compression. Damage to vertebrae, endplates (EPs), the NP, and the AF, can all arise from axial compression, depending on the segment's posture, the manner in which it is loaded, and the physiological state of tissue. Therefore, this movement pattern was selected to illustrate the importance of the bracing effect of a pressurized NP on the AF, and how injuries interrupting support to the AF may contribute to IVD degeneration.

Mechanical properties and fluid permeability of gyroid and diamond lattice structures for intervertebral devices: functional requirements and comparative analysis

Current intervertebral fusion devices present multiple complication risks such as a lack of fixation, device migration and subsidence. An emerging solution to these problems is the use of additively manufactured lattice structures that are mechanically compliant and permeable to fluids, thus promoting osseointegration and reducing complication risks. Strut-based diamond and sheet-based gyroid lattice configurations having a pore diameter of 750 µm and levels of porosity of 60, 70 and 80% are designed and manufactured from Ti-6Al-4V alloy using laser powder bed fusion. The resulting structures are CT-scanned, compression tested and subjected to fluid permeability evaluation. The stiffness of both structures (1.9-4.8 GPa) is comparable to that of bone, while their mechanical resistance (52-160 MPa) is greater than that of vertebrae (3-6 MPa), thus decreasing the risks of wither bone or implant failure. The fluid permeability (5-57 × 10^-9 m²) and surface-to-volume ratios (~3) of both lattice structures are close to those of vertebrae. This study shows that both types of lattice structures can be produced to suit the application specifications within certain limits imposed by physical and equipment-related constraints, providing potential solutions for reducing the complication rate of spinal devices by offering a better fixation through osseointegration.

Human Thoracolumbar Spine Tolerance to Injury and Mechanisms From Caudo-Cephalad Loading: A Parametric Modeling Study

The aims of this investigation were to delineate the internal biomechanics of the spine under vertical impact vector and assess the probability of injury. Male and female whole-body human finite element models were used. The restrained occupants were positioned on the seat, and caudo-cephalad impacts were applied to the base. Different acceleration-time profiles (50-200 ms pulse durations, 11-46 g peak accelerations) were used as inputs in both models. The resulting stress-strain profiles in the cortical and cancellous bones were evaluated at different vertebral levels. Using the peak transmitted forces at the thoracolumbar disc level as the response variable, the probability of injury for the male spine was obtained from experimental risk curves for the various pulses. Results showed that the shorter pulse durations and rise times impart greater loading on the thoracolumbar spine. The analysis of von Mises stress and strain distributions showed that the compression-related fractures are multifaceted with contributions from both the cortical and cancellous bony components of the body. Profiles are provided in the paper. The intervertebral disc may be involved in the fracture mechanism, because it acts as a medium of load transfer between adjacent vertebrae. Injury risks for the shortest pulse was 63%, and for the widest pulse it was close to zero, and injury probabilities for other pulses are given. The present modeling study provides insights into the mechanisms of internal load transfer and describes injury risk levels from caudal to cephalad impacts.

Pedicle Screws Challenged: Lumbar Cortical Density and Thickness Are Greater in the Posterior Elements Than in the Pedicles

STUDY DESIGN

Controlled laboratory study.

OBJECTIVE

To measure the total bone mineral density (BMD), cortical volume, and cortical thickness in seven different anatomical regions of the lumbar spine.

METHODS

Using computed tomography (CT) images, 3 cadaveric spines were digitally isolated by applying filters for cortical and cancellous bone. Each spine model was separated into 5 lumbar vertebrae, followed by segmentation of each vertebra into 7 anatomical regions of interest using 3-dimensional software modeling. The average Hounsfield units (HU) was determined for each region and converted to BMD with calibration phantoms of known BMD. These BMD measurements were further analyzed by the total volume, cortical volume, and cancellous volume. The cortical thickness was also measured. A similar analysis was performed by vertebral segment. St Mary's Medical Center's Institutional Review Board approved this study. No external funding was received for this work.

RESULTS

The lamina and inferior articular process contained the highest total BMD, thickest cortical shell, and largest percent volumes of cortical bone. The vertebral body demonstrated the lowest BMD. The BMDs of the L4 and L5 segments were lower; however, there were no statistically significant differences in BMD between the L1-L5 vertebral segments.

CONCLUSION

Extrapedicular regions of the lumbar vertebrae, including the lamina and inferior articular process, contain denser bone than the pedicles. Since screw pullout strength relies greatly on bone density, the lamina and inferior articular processes may offer stronger fixation of the lumbar spine.

Biomechanical characterization of vertebral body replacement in situ: Effects of different fixation strategies

BACKGROUND AND OBJECTIVE

Artificial vertebral implant with a lateral or posterior screw-rod fixation system are usually employed in lumbar reconstruction surgery to rebuild the lumbar spine after partial resection due to a tumor or trauma. However, few studies have investigated the effect of the various fixation systems on the biomechanics of the reconstructed lumbar system. This study aims to evaluate the influence of different surgical fixation strategies on the biomechanical performance of a reconstructed lumbar spine system in terms of the strength and long-term stability.

METHODS

Two typical lumbar spine reconstruction case models that correspond to lateral or posterior fixation systems were built based on the clinical data. Finite element analyses were performed, and comparisons were made between the two models based on the predicted stress distribution of the reconstructed lumbar spine model, bone-growth area of the endplate, and the range of motion under various normal daily activities.

RESULTS

The load from the upper vertebral body was found to be effectively transmitted onto the lower vertebral body by a vertebral implant with the lateral fixation system; this was favorable for bone growth after surgery. However, significantly high stresses were concentrated around the interaction region between the screws and bone, owing to the uneven lateral fixation structure; this may increase the risk of bone fractures and screw loosening in the long term. For the posterior fixation case, stably posterior fixation structure was favorable to maintain stability for the reconstructed lumbar spine. However, the load was mainly transmitted via the fixation rod rather than the vertebral implant, owing to the stress shielding effect. Therefore, the predicted strain on the endplate were insufficient for bone ingrowth under most of the spinal activates, which could cause bone loss and prosthesis loosening.

CONCLUSIONS

In this study, the comparisons of the reconstructed lumbar spine system with lateral and posterior fixation strategies were conducted. The Pros and Cons of these two fixation strategies was deeply discussed and the associated clinical issues were provided. The results of this study will have a clear impact in understanding the biomechanics of the lumbar spine with different fixation strategies and providing necessary instructions to the design and application of the lumbar spinal fixation system.

Trajectory correction based on the gradient impulse response function improves high-resolution UTE imaging of the musculoskeletal system

PURPOSE

UTE sequences typically acquire data during the ramping up of the gradient fields, which makes UTE imaging prone to eddy current and system delay effects. The purpose of this work was to use a simple gradient impulse response function (GIRF) measurement to estimate the real readout gradient waveform and to demonstrate that precise knowledge of the gradient waveform is important in the context of high-resolution UTE musculoskeletal imaging.

METHODS

The GIRF was measured using the standard hardware of a 3 Tesla scanner and applied on 3D radial UTE data (TE: 0.14 ms). Experiments were performed on a phantom, in vivo on a healthy knee, and in vivo on patients with spine fractures. UTE images were reconstructed twice, first using the GIRF-corrected gradient waveforms and second using nominal-corrected waveforms, correcting for the low-pass filter characteristic of the gradient chain.

RESULTS

Images reconstructed with the nominal-corrected gradient waveforms exhibited blurring and showed edge artifacts. The blurring and the edge artifacts were reduced when the GIRF-corrected gradient waveforms were used, as shown in single-UTE phantom scans and in vivo dual-UTE gradient-echo scans in the knee. Further, the importance of the GIRF-based correction was indicated in UTE images of the lumbar spine, where thin bone structures disappeared when the nominal correction was employed.

CONCLUSION

The presented GIRF-based trajectory correction method using standard scanner hardware can improve the quality of high-resolution UTE musculoskeletal imaging.

Material property calibration is more important than element size and number of different materials on the finite element modelling of vertebral bodies: A Taguchi study

Finite element (FE) modelling of a vertebral body (VB) is considered challenging due to the many parameters involved such as element size and type, and material properties. Previous studies have reported how these parameters affect the mechanical behaviour of a VB model; however, most studies just compared results without any specific statistical tool to quantify their influence. The Taguchi Method (TM) has been successfully used in manufacturing and biomechanics to evaluate process parameters and to determine optimum set-up conditions. This study aimed to evaluate the influence of the main finite element modelling parameters on the mechanical behaviour of a VB model using the Taguchi Method. A FE model was developed based on a C2 juvenile porcine vertebral body and three of the most commonly used modelling parameters were evaluated using TM in terms of the change in the predicted stiffness in comparison to experimental values: element size, number of different material properties for VB (based on grey-scale bins) and calibration factor for grey-scale to density to Young's Modulus equation. The influence of the combined factors was also assessed. The Taguchi analysis showed that the three factors are independent. The calibration factor is the main contributor, accounting for 97% of the predicted stiffness, with the value of 0.03 most closely aligning the numerical and experimental results. Element size accounted for 2% of the predicted stiffness, with 0.75 mm being the optimal, while the number of grey-scale bins influenced the results by less than 1%. Our findings indicate that the calibration factor is the main modelling parameter, with the element size and number of bins accounting for less than 3% of the predicted stiffness. Therefore, calibration of material properties should be done based on a large number of samples to ensure reliable results.

Biomechanical Analysis of Allograft Spacer Failure as a Function of Cortical-Cancellous Ratio in Anterior Cervical Discectomy/Fusion: Allograft Spacer Alone Model

The design and ratio of the cortico-cancellous composition of allograft spacers are associated with graft-related problems, including subsidence and allograft spacer failure. Methods: The study analyzed stress distribution and risk of subsidence according to three types (cortical only, cortical cancellous, cortical lateral walls with a cancellous center bone) and three lengths (11, 12, 14 mm) of allograft spacers under the condition of hybrid motion control, including flexion, extension, axial rotation, and lateral bending,. A detailed finite element model of a previously validated, three-dimensional, intact C3–7 segment, with C5–6 segmental fusion using allograft spacers without fixation, was used in the present study. Findings: Among the three types of cervical allograft spacers evaluated, cortical lateral walls with a cancellous center bone exhibited the highest stress on the cortical bone of spacers, as well as the endplate around the posterior margin of the spacers. The likelihood of allograft spacer failure was highest for 14 mm spacers composed of cortical lateral walls with a cancellous center bone upon flexion (PVMS, 270.0 MPa; 250.2%) and extension (PVMS: 371.40 MPa, 344.2%). The likelihood of allograft spacer subsidence was also highest for the same spacers upon flexion (PVMS, 4.58 MPa; 28.1%) and extension (PVMS: 12.71 MPa, 78.0%). Conclusion: Cervical spacers with a smaller cortical component and of longer length can be risk factors for allograft spacer failure and subsidence, especially in flexion and extension. However, further study of additional fixation methods, such as anterior plates/screws and posterior screws, in an actual clinical setting is necessary.

Photoacoustic imaging for surgical guidance: Principles, applications, and outlook

Minimally invasive surgeries often require complicated maneuvers and delicate hand-eye coordination and ideally would incorporate "x-ray vision" to see beyond tool tips and underneath tissues prior to making incisions. Photoacoustic imaging has the potential to offer this feature but not with ionizing x-rays. Instead, optical fibers and acoustic receivers enable photoacoustic sensing of major structures-such as blood vessels and nerves-that are otherwise hidden from view. This imaging process is initiated by transmitting laser pulses that illuminate regions of interest, causing thermal expansion and the generation of sound waves that are detectable with conventional ultrasound transducers. The recorded signals are then converted to images through the beamforming process. Photoacoustic imaging may be implemented to both target and avoid blood-rich surgical contents (and in some cases simultaneously or independently visualize optical fiber tips or metallic surgical tool tips) in order to prevent accidental injury and assist device operators during minimally invasive surgeries and interventional procedures. Novel light delivery systems, counterintuitive findings, and robotic integration methods introduced by the Photoacoustic & Ultrasonic Systems Engineering Lab are summarized in this invited Perspective, setting the foundation and rationale for the subsequent discussion of the author's views on possible future directions for this exciting frontier known as photoacoustic-guided surgery.

Assessment of femoral neck bone metabolism using 18F-sodium fluoride PET/CT imaging

BACKGROUND

Standard of care metabolic bone disease assessment relies on changes to bone quantity, which can only be detected after structural changes occur.

PURPOSE

To investigate the usefulness of Bone Metabolism Score (BMS), derived from fluorine 18 labeled sodium fluoride (¹⁸F-NaF) PET/CT imaging as a biomarker of localized metabolic changes at the femoral neck.

METHODS

In this retrospective study, 139 participants (68 females and 71 males, ages 21-75 years) that had undergone ¹⁸F-NaF PET/CT were included. BMS was calculated as the ratio of standard uptake value (SUV) in the bone region to that of the total region. Correlations and linear regressions of BMS with age, CT-derived bone mineral density (BMD), body mass index (BMI), height, and weight were conducted. Differences in BMS between women younger and older than the age of 50 years were assessed. Inter- and intra-operator reproducibility was evaluated by coefficient of variation (CV) and intra-class correlation coefficient (ICC).

RESULTS

Among females, age was negatively correlated with left and right whole BMS (5.61% and 4.90% drop in BMS per decade of life) and left and right cortical BMS (10.50% and 10.09% drop in BMS per decade of life). BMS of women older than 50 years was lower than BMS of women younger than 50 years (P < .0001). Among males, age was negatively correlated with left and right whole BMS (4.29% and 4.25% drop in BMS per decade of life) and left and right cortical BMS (9.13% and 10.30% drop in BMS per decade of life). BMD was positively correlated with whole (r = 0.80, P < .0001) and cortical (r = 0.92, P < .0001) BMS.

CONCLUSIONS

BMS could provide functional insight regarding bone metabolism in the femoral neck to complement bone health status assessed through conventional structural imaging. The methodology described herein could be potentially useful for assessing hip fracture risk in individuals when BMD tests provide borderline determination of bone disease.

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.

Effect of osteoporosis on internal fixation after spinal osteotomy: A finite element analysis

BACKGROUND

Severe kyphotic deformity can affect the quality of life of the elderly and is commonly treated by an osteotomy. Considering that the elderly often suffer from osteoporosis, the safety and efficacy of internal fixation are particularly important. The aim of this study was to analyse the effect of osteoporosis on internal fixation after spinal osteotomy.

METHODS

One patient with a thoracolumbar kyphotic deformity who underwent spinal osteotomy was included. The CT images of the entire spine were used to construct a finite element model of the spine internal fixation after osteotomy. Material parameters were assigned to osteoporosis and normal bone groups, and the loads were used to simulate different working conditions, including axial compression, flexion, extension and lateral bending.

FINDINGS

Compared with normal bone mass, the pressure on osteotomized vertebrae was reduced by 8.32%, 1.92%, 36.79% and 79.80% in mild osteoporosis model during axial compression, flexion, extension and lateral bending, respectively. The pressure on screws and rods was increased in an osteoporosis model under axial compression. During flexion and lateral bending, the pressure on screws was increased but was decreased on rods. The opposite result was found during extension. With the degree of osteoporosis increases, the change of stress is more obvious.

INTERPRETATION

Under different bone mass conditions, the distribution patterns of stress in vertebrae, screws and rods were relatively similar. Collectively, the stress levels of vertebral bone were decreased and the stress levels of the screw/rod system were increased in an osteoporosis model compared to a normal bone model. Hence, osteoporosis may increase the risk of fracture and internal fixation failure.

Differences in Trabecular Bone, Cortical Shell, and Endplate Microstructure Across the Lumbar Spine

Background

The quality of the vertebral body structures such as endplate, cortex, and trabecular bone is important for understanding the performance of implants, particularly at the bone-implant interface. Although vertebral body structures have been analyzed separately in the literature, there is no comprehensive study to assess these anatomical measurements along with their interrelationships in the lumbar spine. Therefore, the purpose of this study was to assess variations in trabecular microstructure, vertebral endplate thickness and concavity, and vertebral body cortex thickness within the lumbar spine.

Methods

A total of 80 lumbar vertebrae (L1-L5) were dissected from 16 human cadaver specimens and imaged with microcomputed tomography to determine trabecular microstructure, vertebral cortex thickness, endplate thickness, and maximum endplate concavity depth. A paired t test and regression analysis were used to determine significant differences (P < .05) between different vertebral levels and correlations between the analyzed anatomical parameters.

Results

L1 vertebra had significantly better (P < .02) trabecular bone microstructure (eg, trabecular bone volume fraction) than all other lumbar vertebrae. However, L1 vertebra also had significantly thinner (P ≤ .02) anterior, left, and right cortices compared to all other vertebral levels. Within L3-L5 intervertebral disc spaces, cranial endplates had significantly greater (P ≤ .03) thickness and maximum concavity depth compared to their respective caudal endplates. No strong correlations were observed between trabecular bone microstructure, maximum endplate concavity depth, vertebral cortex, and endplate thickness parameters.

Conclusions

Detailed reference data of these anatomical parameters for each lumbar vertebral body can aid in improved understanding of bone quality, particularly when assessing different implant designs and fixation approaches. Moreover, such anatomical knowledge may help clinicians with optimal implant design selection and surgical placement of these devices into their respective locations.

Can low-frequency guided waves at the tibia paired with machine learning differentiate between healthy and osteopenic/osteoporotic subjects? A pilot study

PURPOSE

Axial transmission quantitative acoustics (ax-QA) has shown to be a promising tool for assessing bone health and properties in a safe, inexpensive, and portable manner. This study investigated the efficacy of low-frequency ax-QA measured at the tibia, paired with a support vector machine (SVM) approach for combining multiple acoustic indicators, to diagnose osteoporosis as defined by bone mineral density.

METHODS

This pilot study measured 41 female subjects using ax-QA (flexural mode, 3 kHz) at the tibia and using dual X-ray absorptiometry (DXA) at the lumbar spine, femoral neck, and distal radius. For each location, a threshold classifier and SVM were trained to differentiate between healthy and non-healthy subjects based on the phase velocity at different frequencies. Receiver Operating Characteristics and area under curve values (AUC) were used to assess the classifiers' performances for various thresholds and class-weights.

RESULTS

The SVM outperformed the threshold classifier for all three bone locations at low false positive rates. While differentiation between healthy and non-healthy bone states was poor for the spine (AUC: 0.56 ± 0.04), good to moderate performances were observed for the radius (AUC: 0.83 ± 0.03) and hip (AUC: 0.71 ± 0.04).

CONCLUSIONS

Low-frequency ax-QA has demonstrated potential for complementing DXA in screening for osteoporosis at the radius and hip. Through further addition of acoustic indicators ax-QA could provide a diagnostic alternative in third-world countries, and bring bone health screening and monitoring into the hands of clinicians and general health practitioners everywhere.

Evaluation of ostarine as a selective androgen receptor modulator in a rat model of postmenopausal osteoporosis

Selective androgen receptor modulators (SARMs) have shown beneficial effects on muscle wasting, general physical function and bone properties in male mammals. However, data on the effects of SARMs in postmenopausal osteoporotic bone are scarce. We evaluated the effects of the SARM drug ostarine on postmenopausal osteoporotic bone in a rat osteoporosis model. Ovariectomy was performed on 46 of 56 3-month-old female Sprague-Dawley rats. Eight weeks after ovariectomy, ostarine was orally administered daily for 5 weeks in dosages of 0.04 (low, OVX + Ost. 0.04), 0.4 (intermediate, OVX + Ost. 0.4), and 4 mg/kg (high, OVX + Ost. 4) body weight. Another ovariectomized group received no ostarine. Lumbar vertebrae and femora were removed for biomechanical, gene expression, ashing, and computer tomography analyses. Low dose showed no effects. The effects of intermediate and high doses were comparable overall. Improvements were mainly seen in structural properties such as bone mineral density and bone volume density. However, the effects in femora were superior to effects in vertebrae. Ostarine treatment for 5 weeks did not improve significantly biomechanical properties. mRNA expression of the receptor activator of NF-κB ligand decreased after treatment, and uterine weight increased. Serum levels of phosphorus increased following ostarine treatment in intermediate and high-dose groups. Short-term treatment of osteoporotic bone with ostarine leads to improvement of several microstructural bone indices. While we did not observe changes in biomechanics, it is conceivable that longer treatment may also improve biomechanical properties. Further studies are needed to characterize longer time effects and side effects of ostarine in osteoporosis.

Improved accuracy in the assessment of vertebral cortical thickness by quantitative computed tomography using the Iterative Convolution OptimizatioN (ICON) method

Vertebral whole bone strength is substantially affected by cortical bone properties. Disease and therapy may affect cancellous and cortical bone differently. Unlike Dual X-ray Absorptiometry (DXA), Quantitative Computed Tomography (QCT) permits selective assessment of cortical and cancellous bone, but image quality limits the accuracy. We present an image processing method specifically adopted to thin cortices that substantially improves accuracy. Ten human vertebrae embedded in epoxy resin were imaged using clinical QCT and High-Resolution QCT (HR-QCT) protocols, both acquired on a clinical whole body CT scanner, whereas high resolution peripheral QCT (HR-pQCT) was used as gold standard. Microstructural variables and BMD were calculated using in-house software StructuralInsight for QCT and HR-QCT and the manufacturer's μCT evaluation software for HR-pQCT. An adjusted measure, a deconvolved cortical thickness (dcCt.Th), corrected for partial volume effects, was derived applying the new Iterative Convolution OptimizatioN (ICON) method. Direct measurements of cortical thickness (Ct.Th) showed substantial overestimation with mean ± standard deviation of 1.8 ± 0.5 mm for QCT and 1.5 ± 0.3 mm for HR-QCT compared to 0.37 ± 0.07 mm using HR-pQCT. Correlations of both QCT (r² = 0.05, p > 0.5.) and HR-QCT (r² = 0.38, p = 0.060) with the gold standard HR-pQCT were not significant. Also QCT-based BMD and BMC as well as HR-QCT-based BMD did not show a significant correlation with the gold standard approach. Only HR-QCT-based BMC showed a modest correlation (r² = 0.59, p = 0.01) After applying ICON corrections, dcCt.Th resulted in 0.52 ± 0.09 mm for QCT and 0.43 ± 0.07 mm for HR-QCT, both significantly correlated to HR-pQCT (r² = 0.75, p = 0.0012 and r² = 0.93, p < 0.0001, respectively). The average overestimation bias of Ct.Th was reduced from (402 ± 157)% to (45 ± 17)% for QCT and from (330 ± 69)% to (19 ± 8)% for HR-QCT. Due to inaccurate segmentation uncorrected QCT-based Ct.Th measures as well as BMD and BMC showed no correlation to HR-pQCT and thus such bias cortical data can be misleading. The application of ICON reduced random overestimation bias to about 50 μm and 20 μm for QCT and HR-QCT, respectively, leading to a recovery of a significant correlation with the reference data of HR-pQCT. This reveals the potential for fairly accurate assessment of cortical thickness, allowing to better characterize cortical mechanical competence. These results warrant testing of the performance in vivo.

Symposium: Imaging modalities for drug-related osteonecrosis of the jaw (6), assessment of mandibular metabolism due to long-term administration of an anti-resorptive agent by bone scintigraphy (secondary publication)

It is not yet known why anti-resorptive agent-related osteonecrosis specifically affects the jaw. Here we assessed changes in the bone metabolism of the mandible in response to long-term bisphosphonate (BP) therapy, and we compared the bone metabolism changes of the mandible with those of other bone sites using a quantitative analysis by bone scintigraphy. The region of interest was selected by identifying without an abnormal accumulation of the mandible, humerus, second and fourth lumbar vertebrae, iliac crest, intertrochanteric femur and diaphysis. Bone scintigraphy images were quantified using a value we termed the ‘bone uptake value (BUV)’. In the low-dose bisphosphonate (LBP) group (n = 21), the patients were undergoing osteoporosis treatment with low-dose BP. The high-dose BP (HBP) group consisted of 12 bone metastasis patients undergoing high-dose BP treatment. The Control group was 47 subjects with oral disease who had never been treated with an anti-resorptive agent. Our analyses demonstrated that with long-term BP administration, the bone metabolism of the iliac crest and intertrochanteric femur was suppressed but that of the mandible was enhanced. There was no significant difference in bone metabolism with either the low-dose BP or high-dose BP treatment. The effects of the long-term administration of BP were site-specific.

Effect of intracortical bone properties on the phase velocity and cut-off frequency of low-frequency guided wave modes (20-85 kHz)

The assessment of intracortical bone properties is of interest since early-stage osteoporosis is associated with resorption in the endosteal region. However, understanding the interaction between ultrasonic guided waves and the cortical bone structure remains challenging. The purpose of this work is to investigate the effect of intracortical bone properties on the ultrasonic response obtained at low-frequency (<100 kHz) using an axial transmission configuration. The semi-analytical finite element method was used to simulate the propagation of guided waves in a waveguide with realistic geometry and material properties. An array of 20 receivers was used to calculate the phase velocity and cut-off frequency of the excited modes using the two-dimensional Fourier transform. The results show that the position of the emitter around the circumference of the bone is an important parameter to control since it can lead to variations of up to 10 dB in the amplitude of the transmitted modes. The cut-off frequency of the high order modes was, however, only slightly affected by the circumferential position of the emitter, and was sensitive mainly to the axial shear modulus. The phase velocity and cut-off frequency in the 20-85 kHz range are promising parameters for the assessment of intracortical properties.

Diffuse reflectance spectroscopy, a potential optical sensing technology for the detection of cortical breaches during spinal screw placement

Safe and accurate placement of screws remains a critical issue in open and minimally invasive spine surgery. We propose to use diffuse reflectance (DR) spectroscopy as a sensing technology at the tip of a surgical instrument to ensure a safe path of the instrument through the cancellous bone of the vertebrae. This approach could potentially reduce the rate of cortical bone breaches, thereby resulting in fewer neural and vascular injuries during spinal fusion surgery. In our study, DR spectra in the wavelength ranges of 400 to 1600 nm were acquired from cancellous and cortical bone from three human cadavers. First, it was investigated whether these spectra can be used to distinguish between the two bone types based on fat, water, and blood content along with photon scattering. Subsequently, the penetration of the bone by an optical probe was simulated using the Monte-Carlo (MC) method, to study if the changes in fat content along the probe path would still enable distinction between the bone types. Finally, the simulation findings were validated via an experimental insertion of an optical screw probe into the vertebra aided by x-ray image guidance. The DR spectra indicate that the amount of fat, blood, and photon scattering is significantly higher in cancellous bone than in cortical bone (p  <  0.01), which allows distinction between the bone types. The MC simulations showed a change in fat content more than 1 mm before the optical probe came in contact with the cortical bone. The experimental insertion of the optical screw probe gave similar results. This study shows that spectral tissue sensing, based on DR spectroscopy at the instrument tip, is a promising technology to identify the transition zone from cancellous to cortical vertebral bone. The technology therefore has the potential to improve the safety and accuracy of spinal screw placement procedures.

3-D Subject-Specific Shape and Density Estimation of the Lumbar Spine From a Single Anteroposterior DXA Image Including Assessment of Cortical and Trabecular Bone

Dual Energy X-ray Absorptiometry (DXA) is the standard exam for osteoporosis diagnosis and fracture risk evaluation at the spine. However, numerous patients with bone fragility are not diagnosed as such. In fact, standard analysis of DXA images does not differentiate between trabecular and cortical bone; neither specifically assess of the bone density in the vertebral body, which is where most of the osteoporotic fractures occur. Quantitative computed tomography (QCT) is an alternative technique that overcomes limitations of DXA-based diagnosis. However, due to the high cost and radiation dose, QCT is not used for osteoporosis management. We propose a method that provides a 3-D subject-specific shape and density estimation of the lumbar spine from a single anteroposterior (AP) DXA image. A 3-D statistical shape and density model is built, using a training set of QCT scans, and registered onto the AP DXA image so that its projection matches it. Cortical and trabecular bone compartments are segmented using a model-based algorithm. Clinical measurements are performed at different bone compartments. Accuracy was evaluated by comparing DXA-derived to QCT-derived 3-D measurements for a validation set of 180 subjects. The shape accuracy was 1.51 mm at the total vertebra and 0.66 mm at the vertebral body. Correlation coefficients between DXA and QCT-derived measurements ranged from 0.81 to 0.97. The method proposed offers an insightful 3-D analysis of the lumbar spine, which could potentially improve osteoporosis and fracture risk assessment in patients who had an AP DXA scan of the lumbar spine without any additional examination.

A finite element study on intra-operative corrective forces and evaluation of screw density in scoliosis surgeries

Scoliosis is an abnormal sideways curvature of the spine and rib cage, which may need surgical treatments. Most of the corrective maneuvers in scoliosis surgeries are based on surgeon’s experience; hence, there is great interest of understanding how the correction ratio can be influenced by the magnitude of forces and moments. Therefore, the objective of this study was to develop and validate a detailed finite element model of the thoracolumbar which can be used to simulate the scoliosis surgeries based on patient-specific clinical images. The validated models of five patients were carefully developed, and the surgery procedures were simulated and the corrective forces were estimated using inverse finite element analysis during the surgery. Furthermore, parametric studies including the influences of the corrective force magnitude and screw density were evaluated. The results showed that the maximum estimated correction force and moment were 173 (±55.43) N and 10.67 (±2.02) N m, respectively, which were aligned with measured clinical observations. The sensitivity analysis on the magnitude of applied force to the screws showed that correction ratio was slightly increased in level 1 (i.e. FB = 1.3 ×  F) but decreased in level 2 (i.e. FB = 1.6 ×  F). In addition, the parametric study on increasing the number of pedicle screws showed that there was no significant difference between lower and higher screw density. However, the stress distribution was significantly greater using higher screw density during correction maneuvers. In conclusion, this study shows a direct relationship between the applied force/moment and screw density and the correction ratio up to a border line which should be defined accurately. This detailed computational modeling can be used in clinic in hope of achieving the optimum outcome of scoliosis surgery using individual patient-specific characterization.

Photoacoustic imaging of a human vertebra: implications for guiding spinal fusion surgeries

It is well known that there are structural differences between cortical and cancellous bone. However, spinal surgeons currently have no reliable method to non-invasively determine these differences in real-time when choosing the optimal starting point and trajectory to insert pedicle screws and avoid surgical complications associated with breached or weakened bone. This paper explores 3D photoacoustic imaging of a human vertebra to noninvasively differentiate cortical from cancellous bone for this surgical task. We observed that signals from the cortical bone tend to appear as compact, high-amplitude signals, while signals from the cancellous bone have lower amplitudes and are more diffuse. In addition, we discovered that the location of the light source for photoacoustic imaging is a critical parameter that can be adjusted to non-invasively determine the optimal entry point into the pedicle. Once inside the pedicle, statistically significant differences in the contrast and SNR of signals originating from the cancellous core of the pedicle (when compared to signals originating from the surrounding cortical bone) were obtained with laser energies of 0.23-2.08 mJ (p  <  0.05). Similar quantitative differences were observed with an energy of 1.57 mJ at distances  ⩾6 mm from the cortical bone of the pedicle. These quantifiable differences between cortical and cancellous bone (when imaging with an ultrasound probe in direct contact with each bone type) can potentially be used to ensure an optimal trajectory during surgery. Our results are promising for the introduction and development of photoacoustic imaging systems to overcome a wide range of longstanding challenges with spinal surgeries, including challenges with the occurrence of bone breaches due to misplaced pedicle screws.

Low Dose of Bisphosphonate Enhances Sclerostin Antibody-Induced Trabecular Bone Mass Gains in Brtl/+ Osteogenesis Imperfecta Mouse Model

Osteogenesis imperfecta (OI) is a genetic disorder characterized by altered bone quality and imbalanced bone remodeling, leading to skeletal fractures that are most prominent during childhood. Treatments for OI have focused on restoring pediatric bone density and architecture to recover functional strength and consequently reduce fragility. Though antiresorptive agents like bisphosphonates (BPs) are currently the most common intervention for the treatment of OI, a number of studies have shown efficacy of sclerostin antibody (SclAb) in inducing gains in bone mass and reducing fragility in OI mouse models. In this study, the effects of the concurrent use of BP and SclAb were evaluated during bone growth in a mouse harboring an OI-causing Gly→Cys mutation on col1a1. A single dose of antiresorptive BP facilitated the anabolic action of SclAb by increasing availability of surfaces for new bone formation via retention of primary trabeculae that would otherwise be remodeled. Chronic effects of concurrent administration of BP and SclAb revealed that accumulating cycles conferred synergistic gains in trabecular mass and vertebral stiffness, suggesting a distinct advantage of both therapies combined. Cortical gains in mass and strength occurred through SclAb alone, independent of presence of BP. In conclusion, these preclinical results support the scientific hypothesis that minimal antiresorptive treatment can amplify the effects of SclAb during early stages of skeletal growth to further improve bone structure and rigidity, a beneficial outcome for children with OI. © 2018 American Society for Bone and Mineral Research.