Imaging the femoral cortex: thickness, density and mass from clinical CT

Individual thigh muscle and proximal femoral features predict displacement in femoral neck Fractures: An AI-driven CT analysis

INTRODUCTION

Hip fractures, particularly among the elderly, impose a significant public health burden due to increased morbidity and mortality. Femoral neck fractures, commonly resulting from low-energy falls, can lead to severe complications such as avascular necrosis, and often necessitate total hip arthroplasty. This study harnesses AI to enhance musculoskeletal assessments by performing automatic muscle segmentation on whole thigh CT scans and detailed cortical measurements using the StradView program. The primary aim is to improve the prediction and prevention of severe femoral neck fractures, ultimately supporting more effective rehabilitation and treatment strategies.

METHODS

This study measured anatomical features from whole thigh CT scans of 60 femoral neck fracture patients. An AI-driven individual muscle segmentation model (a dice score of 0.84) segmented 27 muscles in the thigh region, to calculate muscle volumes. Proximal femoral bone parameters were measured using StradView, including average cortical thickness, inner density and FWHM at four regions. Correlation analysis evaluated relationships between muscle features, cortical parameters, and fracture displacement. Machine learning models (Random Forest, SVM and Multi-layer Perceptron) predicted displacement using these variables.

RESULTS

Correlation analysis showed significant associations between femoral neck displacement and trabecular density at the femoral neck/intertrochanter, as well as volumes of specific thigh muscles such as the Tensor fasciae latae. Machine learning models using a combined feature set of thigh muscle volumes and proximal femoral parameters performed best in predicting displacement, with the Random Forest model achieving an F1 score of 0.91 and SVM model 0.93.

CONCLUSION

Decreased volumes of the Tensor fasciae latae, Rectus femoris, and Semimembranosus muscles, coupled with reduced trabecular density at the femoral neck and intertrochanter, were significantly associated with increased fracture displacement. Notably, our SVM model-integrating both muscle and femoral features-achieved the highest predictive performance. These findings underscore the critical importance of muscle strength and bone density in rehabilitation planning and highlight the potential of AI-driven predictive models for improving clinical outcomes in femoral neck fractures.

Cranial vault thickness, its internal organization, and its relationship with endocranial shape in Neanderthals and modern humans

Compared to the more elongated crania of Neanderthals, modern humans have a rounder, more globular cranial vault. The factors contributing to this globular cranial and endocranial morphology remain poorly understood. Cranial vault thickness (CVT) plays a role in shaping the braincase. It has been proposed that CVT variation in hominins reflects distinct stimuli influencing the cranial vault layers to different degrees. We aim to determine to what extent CVT differences could explain the well-documented endocranial shape differences between modern humans and Neanderthals. Additionally, we quantify the contributions of each cranial vault layer-the inner and outer tables and the diploë-to elucidate the processes driving CVT variation. We analyzed crania from seven Neanderthals and 75 modern humans using a ray casting method to measure cortical, diploic, and overall CVT. We generated morphometric maps of thickness distributions for each cranial vault layer. We then quantified the endocranial shape using geometric morphometrics and analyzed its relationship with CVT using two-block partial least squares analysis and regression models. Our findings reveal that Neanderthals generally have greater CVT and thicker layers of the cranial vault than modern humans. However, their ranges overlap with those of modern humans. The thicker vaults of Neanderthals are primarily driven by the diploic layer, with notable differences in thickness distribution patterns. Additionally, the inner and outer cortical tables exhibit distinct thickness distribution patterns between the two groups. Structural differences between the parietal bones of Neanderthals and modern humans are particularly pronounced. Furthermore, we observed a unique correspondence between CVT and endocranial shape in modern humans, which is not shared by Neanderthals.

AFSegNet: few-shot 3D ankle-foot bone segmentation via hierarchical feature distillation and multi-scale attention and fusion

Accurate segmentation of ankle and foot bones from CT scans is essential for morphological analysis. Ankle and foot bone segmentation challenges due to the blurred bone boundaries, narrow inter-bone gaps, gaps in the cortical shell, and uneven spongy bone textures. Our study endeavors to create a deep learning framework that harnesses advantages of 3D deep learning and tackles the hurdles in accurately segmenting ankle and foot bones from clinical CT scans. A few-shot framework AFSegNet is proposed considering the computational cost, which comprises three 3D deep-learning networks adhering to the principles of progressing from simple to complex tasks and network structures. Specifically, a shallow network first over-segments the foreground, and along with the foreground ground truth are used to supervise a subsequent network to detect the over-segmented regions, which are overwhelmingly inter-bone gaps. The foreground and inter-bone gap probability map are then input into a network with multi-scale attentions and feature fusion, a loss function combining region-, boundary-, and topology-based terms to get the fine-level bone segmentation. AFSegNet is applied to the 16-class segmentation task utilizing 123 in-house CT scans, which only requires a GPU with 24 GB memory since the three sub-networks can be successively and individually trained. AFSegNet achieves a Dice of 0.953 and average surface distance of 0.207. The ablation study and comparison with two basic state-of-the-art networks indicates the effectiveness of the progressively distilled features, attention and feature fusion modules, and hybrid loss functions, with the mean surface distance error decreased up to 50 %.

In-silico study of the biomechanical effects of proximal-fibular osteotomy on knee joint contact pressure in varus-valgus misalignment

The aim of this work is to investigate in-silico the biomechanical effects of a proximal fibular osteotomy (PFO) on a knee joint with different varus/valgus deformities on the progression of knee osteoarthritis (KOA). A finite element analysis (FEA) of a human lower extremity consisting of the femoral, tibial and fibular bones and the cartilage connecting them was designed. The FEA was performed in a static standing primitive position to determine the contact pressure (CP) distribution and the location of the center of pressure (CoP). The analysis examined the relationship between these factors and the degree of deformation of the hip-knee angle in the baseline condition. The results suggested that PFO could be a simple and effective surgical treatment for patients with associated genu varum. This work also reported that a possible CP homogenization and a CoP correction can be achieved for medial varus deformities after PFO. However, it reduced its effectiveness for tibial origin valgus misalignment and worsened in cases of femoral valgus misalignment.

Anteromedial Cortical Support in Reduction of Trochanteric Hip Fractures: From Definition to Application

➤ The concept of anteromedial cortical support (AMCS) serves as valuable guidance for the intraoperative reduction of trochanteric hip fractures.➤ Positive medial cortical support (MCS) and positive or neutral anterior cortical support (ACS) are desirable. Some evidence has suggested that positive MCS is potentially superior to neutral MCS.➤ Experimental studies underscore the vital importance of the anteromedial wall and reveal why positive MCS potentially outperforms neutral MCS.➤ Incorporating the AMCS concept, the Chang reduction quality criteria (CRQC) are a reliable alternative approach to evaluate the reduction quality of trochanteric hip fractures.

How effective is proximal fibular osteotomy in redistributing joint pressures? Insights from an HTO comparative in-silico study

BACKGROUND

Knee osteoarthritis (KOA) represents a widespread degenerative condition among adults that significantly affects quality of life. This study aims to elucidate the biomechanical implications of proximal fibular osteotomy (PFO), a proposed cost-effective and straightforward intervention for KOA, comparing its effects against traditional high tibial osteotomy (HTO) through in-silico analysis.

METHODS

Using medical imaging and finite element analysis (FEA), this research quantitatively evaluates the biomechanical outcomes of a simulated PFO procedure in patients with severe medial compartment genu-varum, who have undergone surgical correction with HTO. The study focused on evaluating changes in knee joint contact pressures, stress distribution, and anatomical positioning of the center of pressure (CoP). Three models are generated for each of the five patients investigated in this study, a preoperative original condition model, an in-silico PFO based on the same original condition data, and a reversed-engineered HTO in-silico model.

RESULTS

The novel contribution of this investigation is the quantitative analysis of the impact of PFO on the biomechanics of the knee joint. The results provide mechanical evidence that PFO can effectively redistribute and homogenize joint stresses, while also repositioning the CoP towards the center of the knee, similar to what is observed post HTO. The findings propose PFO as a potentially viable and simpler alternative to conventional surgical methods for managing severe KOA, specifically in patients with medial compartment genu-varum.

CONCLUSION

This research also marks the first application of FEA that may support one of the underlying biomechanical theories of PFO, providing a foundation for future clinical and in-silico studies.

Experimental study on determining the degree of bone healing by wall thickness ratio analysis

To verify the reliability and accuracy of wall thickness ratio analysis to determine the degree of bone healing, fracture models were established with 6 beagles. X-ray, micro-CT, and CT scans were performed at 24 weeks. The healthy side and the affected side were used to simulate the three-dimensional geometric model after internal fixation, and the mesh was divided. The mean and median CT wall thickness values were obtained through the wall thickness analysis. X-ray, CT, micro-CT, and gross appearance were used to determine the degree of bone healing, which was compared with wall thickness analysis. There was a positive correlation between the average CT value and the median wall thickness. The correlation coefficient analysis of the median wall thickness ratio (R2) and healing index ratio (R3) showed a positive correlation. The results of the wall thickness ratio (R2) and the healing index ratio (R3) were used to determine bone healing, and the results were consistent with the results of the actual mechanical test and image analysis. The results of wall thickness ratio analysis were significantly correlated with the degree of bone healing. This method is simple, rapid, and practical to analyze and judge the degree of bone healing.

A Correspondence-Based Network Approach for Groupwise Analysis of Patient-Specific Spatiotemporal Data

Methods for statistically analyzing patient-specific data that vary both spatially and over time are currently either limited to summary statistics or require elaborate surface registration. We propose a new method, called correspondence-based network analysis, which leverages particle-based shape modeling to establish correspondence across a population and preserve patient-specific measurements and predictions through statistical analysis. Herein, we evaluated this method using three published datasets of the hip describing cortical bone thickness of the proximal femur, cartilage contact stress, and dynamic joint space between control and patient cohorts to evaluate activity- and group-based differences, as applicable, using traditional statistical parametric mapping (SPM) and our proposed spatially considerate correspondence-based network analysis approach. The network approach was insensitive to correspondence density, while the traditional application of SPM showed decreasing area of the region of significance with increasing correspondence density. In comparison to SPM, the network approach identified broader and more connected regions of significance for all three datasets. The correspondence-based network analysis approach identified differences between groups and activities without loss of subject and spatial specificity which could improve clinical interpretation of results.

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.

Micron-resolution Imaging of Cortical Bone under 14 T Ultrahigh Magnetic Field

Compact, mineralized cortical bone tissues are often concealed on magnetic resonance (MR) images. Recent development of MR instruments and pulse techniques has yielded significant advances in acquiring anatomical and physiological information from cortical bone despite its poor 1 H signals. This work demonstrates the first MR research on cortical bones under an ultrahigh magnetic field of 14 T. The 1 H signals of different mammalian species exhibit multi-exponential decays of three characteristic T2 or T2 * values: 0.1-0.5 ms, 1-4 ms, and 4-8 ms. Systematic sample comparisons attribute these T2 /T2 * value ranges to collagen-bound water, pore water, and lipids, respectively. Ultrashort echo time (UTE) imaging under 14 T yielded spatial resolutions of 20-80 microns, which resolves the 3D anatomy of the Haversian canals. The T2 * relaxation characteristics further allow spatial classifications of collagen, pore water and lipids in human specimens. The study achieves a record of the spatial resolution for MR imaging in bone and shows that ultrahigh-field MR has the unique ability to differentiate the soft and organic compartments in bone tissues.

Positive medial cortical support versus anatomical reduction for trochanteric hip fractures: Finite element analysis and biomechanical testing

BACKGROUND AND OBJECTIVES

The anatomical reduction (AR) is usually considered the best option for fractures. Nevertheless, in unstable trochanteric hip fractures (UTHF), previous clinical reports found that the positive medial cortical support (PMCS, an over-reduction technique) attained higher mechanical stability, but this challenging clinical finding still needs experimental validation.

METHODS

This study constructed in-silico and biomechanical PMCS and AR models, with the use of the most clinically-representative geometry design of fracture models, the multi-directional design in FE analysis, and the subject-specific (osteoporotic) bone material properties, to make the models better mimic the actual condition in clinical settings. Then multiple performance variables (von-Mises stress, strain, integral axial stiffness, displacement, structural changes, etc.) were assessed to uncover details of integral and regional stability.

RESULTS

Among in-silico comparison, PMCS models showed significantly lower maximum displacement than AR models, and the maximum von Mises stress of implants (MVMS-I) was significantly lower in PMCS models than in AR models (highest MVMS-I in -30°-A3-AR of 1055.80 ± 93.37 MPa). Besides, PMCS models had significantly lower maximum von Mises stress along fracture surfaces (MVMS-F) (highest MVMS-F in 30°-A2-AR of 416.40 ± 38.01 MPa). Among biomechanical testing comparison, PMCS models showed significantly lower axial displacement. Significantly lower change of neck-shaft angle (CNSA) was observed in A2-PMCS models. A fair amount of AR models converted into the obvious negative medial cortical support (NMCS) condition, whereas all PMCS models kept the PMCS condition. The results were also validated through comparison to previous clinical data.

CONCLUSIONS

The PMCS is superior to the AR in the UTHF surgery. The current study opens up the second thought of the role of over-reduction technique in bone surgery.

3D bioprinted autologous bone particle scaffolds for cranioplasty promote bone regeneration with both implanted and native BMSCs

Although autologous bone (AB) grafting is considered to be the gold standard for cranioplasty, unresolved problems remain, such as surgical-site infections and bone flap absorption. In this study, an AB scaffold was constructed via three-dimensional (3D) bedside-bioprinting technology and used for cranioplasty. To simulate the skull structure, a polycaprolactone shell was designed as an external lamina, and 3D-printed AB and a bone marrow-derived mesenchymal stem cell (BMSC) hydrogel was used to mimic cancellous bone for bone regeneration. Ourin vitroresults showed that the scaffold exhibited excellent cellular affinity and promoted osteogenic differentiation of BMSCs in both two-dimensional and 3D culture systems. The scaffold was implanted in beagle dog cranial defects for up to 9 months, and the scaffold promoted new bone and osteoid formation. Furtherin vivostudies indicated that transplanted BMSCs differentiated into vascular endothelium, cartilage, and bone tissues, whereas native BMSCs were recruited into the defect. The results of this study provide a method for bedside bioprinting of a cranioplasty scaffold for bone regeneration, which opens up another window for clinical applications of 3D printing in the future.

Advancing Osteoporosis Evaluation Procedures: Detailed Computational Analysis of Regional Structural Vulnerabilities in Osteoporotic Bone

Osteoporotic fractures of the femur are associated with poor healing, disability, reduced quality of life, and high mortality rates within 1 year. Moreover, osteoporotic fractures of the femur are still considered to be an unsolved problem in orthopedic surgery. In order to more effectively identify osteoporosis-related fracture risk and develop advanced treatment approaches for femur fractures, it is necessary to acquire a greater understanding of how osteoporosis alters the diaphyseal structure and biomechanical characteristics. The current investigation uses computational analyses to comprehensively examine how femur structure and its associated properties differ between healthy and osteoporotic bones. The results indicate statistically significant differences in multiple geometric properties between healthy femurs and osteoporotic femurs. Additionally, localized disparities in the geometric properties are evident. Overall, this approach will be beneficial in the development of new diagnostic procedures for highly detailed patient-specific detection of fracture risk, for establishing novel injury prevention treatments, and for informing advanced surgical solutions.

On the influence of computed tomography's slice thickness on computer tomography based finite element analyses results

BACKGROUND

Patient-specific autonomous finite element analyses of femurs, based on clinical computed tomography scans may be used to monitor the progression of bone-related diseases. Some CT scan protocols provide lower resolution (slice thickness of 3 mm) that affects the accuracy. To investigate the impact of low-resolution scans on the CT-based finite element analyses results, identical CT raw data were reconstructed twice to generate a 1 mm ("gold standard") and a 3 mm slice thickness scans.

METHODS

CT-based finite element analyses of twenty-four femurs (twelve patients) under stance and sideways fall loads were performed based on 1 and 3 mm slice thickness scans. Bone volume, load direction, and strains were extracted at different locations along the femurs and differences were evaluated.

FINDINGS

Average differences in bone volume were 1.0 ± 1.5%. The largest average difference in strains in stance position was in the neck region (11.0 ± 13.4%), whereas in other regions these were much smaller. For sidewise fall loading, the average differences were at most 9.2 ± 16.0%.

INTERPRETATION

Whole-body low dose CT scans (3 mm-slice thickness) are suboptimal for monitoring strain changes in patient's femurs but may allow longitudinal studies if larger than 5% in all areas and larger than 12% in the upper neck. CT-based finite element analyses with slice thickness of 3 mm may be used in clinical practice for patients with smoldering myeloma to associate changes in strains with progression to active myeloma if above ∼10%.

Cross-sectional geometry of the femoral diaphyseal cortical bones: analysis of central mass distribution

A detailed analysis of differences in skeletal shape among many individuals is expected to reveal the mechanical significance behind various morphological features. To confirm the distribution of the cortical bone region in cross sections, the relative position of the central mass distribution (CMD) of the cortical bone region to the CMD of the entire cross section was examined. A total of 90 right human femoral skeletons were examined using clinical multi-slice computed tomography. For nine cross sections of each femur, we determined the CMD of the whole area, including both cortical bone and medullary areas, as CMD-W, and that of the cortical bone region in the same cross section as CMD-C, and they were compared. The medial and anterior portion of the cortex was relatively thick just below the lesser trochanter. The posterior cortical bone tended to be relatively thick in the region from the center to the distal part of the diaphysis. Females had a significantly more medially deviated CMD than males throughout the entire diaphysis. These results suggest that femurs with advanced cortical bone thinning tend to have a concentration of cortical bone in their medial portion. CMD-C was located farther from the diaphysis axis as the degree of medial bending increased. Conversely, the greater the lateral bending of the diaphysis, the closer CMD-C was to the diaphysis axis. As the amount of bone decreases with age, self-adjustment could occur so that the cortical bone's critical area remains to prevent a decrease in mechanical strength.

The West China Hospital radiographic classification for fibrous dysplasia in femur and adjacent bones: A retrospective analysis of 205 patients

Objective

This study aims to investigate the reliability and clinical outcome of a newly developed classification system for patients with fibrous dysplasia (FD) of the femur and adjacent bones, optimizing its evaluation and management.

Methods

A total of 205 patients (121 female and 84 male) with FD in the femur and adjacent bones were included in this retrospective study. All affected femurs were measured and treated based on this classification at our institution between 2009 and 2019. Based on previous studies and extensive clinical follow‐up, we cautiously proposed the West China Hospital radiographic classification for FD in femur and adjacent bones following corresponding treatment options. There are five types with five radiographic features, including proximal femur bone loss, coxa vara, femoral shaft deformity, genu valgum, and hip arthritis. The intraobserver and interobserver reproducibility of this classification was assessed by four observers using the Cohen kappa statistic. The clinical outcome was evaluated using the criteria of Guille.

Results

At a median follow‐up of 60 months (range 6–120), 205 patients (median 34.8 years old, range 18–73 years old) were categorized into the following five types: Type I 31.7%, Type II 30.2%, Type III 20.5%, Type IV 10.2%, and Type V 7.3%. The mean interobserver and intraobserver kappa scores were 0.85 (range 0.77–0.89) and 0.85 (range 0.79–0.92), respectively. For clinical outcomes, there was no significant difference in the postoperative Guille score for Type I patients (mean 9.01 ± 1.22). There was a significant increase in the postoperative Guille score in Type II, III, IV, and V, compared to the preoperative values (P < 0.01). For complications, two Type III patients reported pain, and one Type III patient had mild‐to‐moderate Trendelenburg gait. One Type IV patient had a mild Trendelenburg gait. And two Type V patients still had mild limping.

Conclusion

This classification is reproducible and serves as a tool for evaluating and treating FD in the femur and adjacent bones. Therefore, we recommend this classification for the diagnosis and treatment of FD‐related deformities in the femur and adjacent bones.

Cetaceans Humerus Radiodensity by CT: A Useful Technique Differentiating between Species, Ecophysiology, and Age

Cetaceans are mammals that underwent a series of evolutionary adaptations to live in the aquatic environment, including morphological modifications of various anatomical structures of the skeleton and their bone mineral density (BMD); there are few studies on the latter. BMD is related to the radiodensity measured through computed tomography (CT) in Hounsfield units (HU). This work aimed to test and validate the usefulness of studying humeral bone radiodensity by CT of two cetacean species (the Atlantic spotted dolphin and the pygmy sperm whale) with different swimming and diving habits. The radiodensity was analysed at certain levels following a new protocol based on a review of previous studies. Humeral radiodensity values were related to four aspects: species, diving behaviour, swimming activity level, and age. We observed that the consistent differences in the radiodensity of the cortical bone of the distal epiphysis between animals of different life-history categories suggest that this bone portion could be particularly useful for future ontogenetic studies. Hence, this technique may be helpful in studying and comparing species with different ecophysiologies, particularly distinguishing between swimming and diving habits.

A deep learning-based approach to automatic proximal femur segmentation in quantitative CT images

Automatic CT segmentation of proximal femur has a great potential for use in orthopedic diseases, especially in the imaging-based assessments of hip fracture risk. In this study, we proposed an approach based on deep learning for the fast and automatic extraction of the periosteal and endosteal contours of proximal femur in order to differentiate cortical and trabecular bone compartments. A three-dimensional (3D) end-to-end fully convolutional neural network (CNN), which can better combine the information among neighbor slices and get more accurate segmentation results by 3D CNN, was developed for our segmentation task. The separation of cortical and trabecular bones derived from the QCT software MIAF-Femur was used as the segmentation reference. Two models with the same network structures were trained, and they achieved a dice similarity coefficient (DSC) of 97.82% and 96.53% for the periosteal and endosteal contours, respectively. Compared with MIAF-Femur, it takes half an hour to segment a case, and our CNN model takes a few minutes. To verify the excellent performance of our model for proximal femoral segmentation, we measured the volumes of different parts of the proximal femur and compared it with the ground truth, and the relative errors of femur volume between predicted result and ground truth are all less than 5%. This approach will be expected helpful to measure the bone mineral densities of cortical and trabecular bones, and to evaluate the bone strength based on FEA.

Can neck fractures in proximal humeri be predicted by CT-based FEA?

BACKGROUND

Proximal humeri fractures at anatomical and surgical neck (∼5% and ∼50% incidence respectively) are frequent in elderly population. Yet, neither in-vitro experiments nor CT-based finite element analyses (CTFEA) have investigated these in depth. Herein we enhance (Dahan et al., 2019) (addressing anatomical neck fractures) by more experiments and specimens, accounting for surgical neck fractures and explore CTFEA's prediction of humeri mechanical response and yield force.

METHODS

Four fresh frozen human humeri were tested in a new experimental configuration inducing surgical neck fractures. Digital image correlation (DIC) provided strains and displacements on humeri surfaces and used to validate CTFEA predictions. CTFEA were enhanced herein to improve the accuracy at the proximal neck: A cortical bone mapping (CBM) algorithm was implemented to overcome insufficient scanning resolution, and a new trabecular material mapping was investigated.

RESULTS

The new experimental setting induced impacted surgical neck fractures in all humeri. Excellent DIC to CTFEA correlation in strains was obtained at the shaft (slope 0.984, R²=0.99) and a fair agreement (slope 0.807, R²=0.73) at the neck. CBM algorithm had worsened the correlation, whereas the new material mapping had a negligible influence. Yield loads predictions improved considerably when trabecular yielding (maximum principal strain criterion) was considered instead of surface cortical yielding.

DISCUSSION

CTFEA well predicts strains on the shaft and reasonably well on the neck. This enhances former conclusions by past studies conducted using SGs, now also evident by DIC. Yield load prediction for surgical neck fractures (involving crushing of trabecular bone) is predicted better by trabecular failure laws rather than cortex ones. Further FEA studies using trabecular orthotropic constitutive models and failure laws are warrant.

Romosozumab Enhances Vertebral Bone Structure in Women With Low Bone Density

Romosozumab monoclonal antibody treatment works by binding sclerostin and causing rapid stimulation of bone formation while decreasing bone resorption. The location and local magnitude of vertebral bone accrual by romosozumab and how it compares to teriparatide remains to be investigated. Here we analyzed the data from a study collecting lumbar computed tomography (CT) spine scans at enrollment and 12 months post-treatment with romosozumab (210 mg sc monthly, n = 17), open-label daily teriparatide (20 μg sc, n = 19), or placebo (sc monthly, n = 20). For each of the 56 women, cortical thickness (Ct.Th), endocortical thickness (Ec.Th), cortical bone mineral density (Ct.bone mineral density (BMD)), cancellous BMD (Cn.BMD), and cortical mass surface density (CMSD) were measured across the first lumbar vertebral surface. In addition, color maps of the changes in the lumbar vertebrae structure were statistically analyzed and then visualized on the bone surface. At 12 months, romosozumab improved all parameters significantly over placebo and resulted in a mean vertebral Ct.Th increase of 10.3% versus 4.3% for teriparatide, an Ec.Th increase of 137.6% versus 47.5% for teriparatide, a Ct.BMD increase of 2.1% versus a -0.1% decrease for teriparatide, and a CMSD increase of 12.4% versus 3.8% for teriparatide. For all these measurements, the differences between romosozumab and teriparatide were statistically significant (p < 0.05). There was no significant difference between the romosozumab-associated Cn.BMD gains of 22.2% versus 18.1% for teriparatide, but both were significantly greater compared with the change in the placebo group (-4.6%, p < 0.05). Cortical maps showed the topographical locations of the increase in bone in fracture-prone areas of the vertebral shell, walls, and endplates. This study confirms widespread vertebral bone accrual with romosozumab or teriparatide treatment and provides new insights into how the rapid prevention of vertebral fractures is achieved in women with osteoporosis using these anabolic agents. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Protective effect of bisphosphonate on the cortical bone at key locations of the femur in aromatase inhibitor-associated bone loss: A three-dimensional cortical bone mapping study

•

Aromatase inhibitor use was associated with cortical bone loss in the hip.

•

Bisphosphonate protected hip cortical bone against aromatase inhibitor use.

•

The effect was prominent at the superior femoral neck and greater trochanter.

Aromatase inhibitor treatment in breast cancer is associated with accelerated bone loss and an increased risk of fracture. Bisphosphonates (BPs) are the mainstay treatment of aromatase inhibitor-associated bone loss (AIBL), which might improve femoral bone at key locations prone to fracture. To test this hypothesis, we performed three-dimensional cortical bone mapping based on quantitative computed tomography (QCT) scans in postmenopausal women with early breast cancer who were receiving aromatase inhibitors. Data of subjects who had both baseline and at least one follow-up QCT at Severance Hospital (South Korea) between 2005 and 2015 were analyzed (BP users, n = 93; BP non-users, n = 203). After exclusion of BP users with low medication persistence (proportion of days covered: <50%), BP users and non-users were 1:1 matched (n = 54 for each group) in terms of age, lumbar spine volumetric bone mineral density (LSvBMD), femoral neck areal BMD (FNaBMD), and total hip areal BMD (THaBMD). During a median follow-up of 2.1 years, BP use attenuated bone loss in LSvBMD (+7.2% vs. −3.8%, p < 0.001), FNaBMD (+1.3% vs. −2.7%, p < 0.001), and THaBMD (-0.3% vs. −2.5%, p = 0.024). BP had a protective effect on cortical parameters of femoral bone: estimated cortical thickness (CTh) (+3.3% vs. + 0.1%, p = 0.007) and cortical mass surface density (CMSD, cortical mass per unit surface area was calculated by multiplying cortical BMD with CTh) (+3.4% vs. −0.3%, p < 0.001). CMSD increased by up to 15% at key locations such as the superior part of the femoral neck and greater trochanter. BP prevented the thinning of average CTh of the femoral neck (-1.4% vs. −6.1%, p < 0.001), particularly at the superior anterior quadrant of femoral neck (absolute difference: +12.8% point vs. non-users). Compared to BP non-users, BP users had improved cross-sectional moment of inertia (+4.4% vs. −0.7%, p = 0.001) and less increase in buckling ratio (+1.3% vs. + 7.5%, p < 0.001). In summary, BP use prevented cortical bone deficits observed in AIBL at key locations of the proximal femur.

Spatial distribution of hip cortical thickness in postmenopausal women with different osteoporotic fractures

Few studies h ave discussed the association between cortical bone outside the fracture site and the fracture itself. Focusing on hip cortical thickness, this study revealed distinct distributions of the parameters for hip (trochanteric or femoral neck), vertebral, and peripheral osteoporotic fractures and suggested that the spatial distribution of hip cortical thickness was fracture-specific.

PURPOSE

Cortical bone is critical for bone strength. Hip cortical thickness is reported to be closely associated with the incidence of hip fractures, but its relationship with nonhip fractures is rarely studied. As the hip is a major site for fracture risk assessment, it would be of great benefit to investigate the association between hip cortical thickness and different osteoporotic fractures.

METHODS

One hundred age-matched postmenopausal women were equally assigned to 4 osteoporotic fracture groups (trochanteric, femoral neck, vertebral, and peripheral fractures) and a nonfracture group. Each subject had a clinical quantitative computed tomography scan of the bilateral hips and the lumbar spine. A cortical bone mapping algorithm was adopted to calculate hip cortical thickness. Hip and lumbar trabecular density and the hip cortical thickness distribution were compared among the groups.

RESULTS

All the fracture groups presented lower lumbar trabecular density. Compared with nonfracture controls, patients with hip or vertebral fractures but not peripheral fractures showed decreased cortical thickness and trabecular density of the hip. Fracture-specific distributions of cortical thickness were revealed, including zonal defects on the neck-intertrochanter junction, greater trochanter, and the periphery of the lesser trochanter for trochanteric fractures, a focal defect on the anterosuperior neck for femoral neck fractures, a moderate and average distribution for vertebral fractures, and focally thicker cortices on the anterosuperior greater trochanter and the periphery of the lesser trochanter for peripheral fractures.

CONCLUSION

The spatial distribution of hip cortical thickness was different for each type of osteoporotic fracture, and patients with centrally located fractures demonstrated more severe cortical deterioration. This finding needs to be validated in a larger sample.

Automated cortical thickness measurement of the mandibular condyle head on CBCT images using a deep learning method

This study proposes a deep learning model for cortical bone segmentation in the mandibular condyle head using cone-beam computed tomography (CBCT) and an automated method for measuring cortical thickness with a color display based on the segmentation results. In total, 12,800 CBCT images from 25 normal subjects, manually labeled by an oral radiologist, served as the gold-standard. The segmentation model combined a modified U-Net and a convolutional neural network for target region classification. Model performance was evaluated using intersection over union (IoU) and the Hausdorff distance in comparison with the gold standard. The second automated model measured the cortical thickness based on a three-dimensional (3D) model rendered from the segmentation results and presented a color visualization of the measurements. The IoU and Hausdorff distance showed high accuracy (0.870 and 0.928 for marrow bone and 0.734 and 1.247 for cortical bone, respectively). A visual comparison of the 3D color maps showed a similar trend to the gold standard. This algorithm for automatic segmentation of the mandibular condyle head and visualization of the measured cortical thickness as a 3D-rendered model with a color map may contribute to the automated quantification of bone thickness changes of the temporomandibular joint complex on CBCT.

Incorporating Nutrition, Vests, Education, and Strength Training (INVEST) in Bone Health: Trial Design and Methods

BACKGROUND

Achievement of 5-10% weight loss (WL) among older adults living with obesity considerably improves prognosis of health-related outcomes; however, concomitant declines in bone mineral density (BMD) limit overall benefit by increasing fracture risk. Declines in mechanical loading contribute to WL-associated BMD loss, with pilot data signaling the addition of external weight replacement (via weighted vest use) during intentional WL mitigates bone loss at weight bearing sites to a similar degree as resistance exercise training (RT). Definitive data in support of weighted vest use as a potential strategy to mitigate WL-associated bone loss in this population are needed.

METHODS

In the Incorporating Nutrition, Vests, Education, and Strength Training (INVEST) in Bone Health trial (NCT04076618), 192 older adults (60-85 years) who are overweight (BMI ≥ 27 kg/m²) with at least one obesity-related risk factor or obese (BMI = 30-40 kg/m²) will be randomly assigned to participate in one of three 12-month intervention groups: WL alone, WL + weighted vest use (WL + VEST), or WL + RT. The primary aim is to determine the effects of WL + VEST compared to WL alone and WL + RT on indicators of bone health and subsequent fracture risk.

DISCUSSION

Determining effective, translatable strategies that minimize bone loss during intentional WL among older adults holds public health potential. The INVEST in Bone Health trial offers an innovative approach for increasing mechanical stress during intentional WL in the absence of RT. If successful, findings from this study will provide evidence in support of a scalable solution to minimize bone loss during intentional WL among older adults with obesity.

Subject-specific rib finite element models with material data derived from coupon tests under bending loading

Rib fractures are common thoracic injuries in motor vehicle crashes. Several human finite element (FE) human models have been created to numerically assess thoracic injury risks. However, the accurate prediction of rib biomechanical response has shown to be challenging due to human variation and modeling approaches. The main objective of this study was to better understand the role of modeling approaches on the biomechanical response of human ribs in anterior-posterior bending. Since the development of subject specific rib models is a time-consuming process, the second objective of this study was to develop an accurate morphing approach to quickly generate high quality subject specific rib meshes. The exterior geometries and cortical-trabecular boundaries of five human 6th-level ribs were extracted from CT-images. One rib mesh was developed in a parametric fashion and the other four ribs were developed with an in-house morphing algorithm. The morphing algorithm automatically defined landmarks on both the periosteal and endosteal boundaries of the cortical layer, which were used to morph the template nodes to target geometries. Three different cortical bone material models were defined based on the stress-strain data obtained from subject-specific tensile coupon tests for each rib. Full rib anterior-posterior bending tests were simulated based on data recorded in testing. The results showed similar trends to test data with some sensitivity relative to the material modeling approach. Additionally, the FE models were substantially more resistant to failure, highlighting the need for better techniques to model rib fracture. Overall, the results of this work can be used to improve the biofidelity of human rib finite element models.

Validation of HR-pQCT against micro-CT for morphometric and biomechanical analyses: A review

High-resolution peripheral quantitative computed-tomography (HR-pQCT) has the potential to become a powerful clinical assessment and diagnostic tool. Given the recent improvements in image resolution, from 82 to 61 μm, this technology may be used to accurately quantify in vivo bone microarchitecture, a key biomarker of degenerative bone diseases. However, computational methods to assess bone microarchitecture were developed for micro computed tomography (micro-CT), a higher-resolution technology only available for ex vivo studies, and validation of these computational analysis techniques against the gold-standard micro-CT has been inconsistent and incomplete. Herein, we review methods for segmentation of bone compartments and microstructure, quantification of bone morphology, and estimation of mechanical strength using finite-element analysis, highlighting the need throughout for improved standardization across the field. Studies have relied on homogenous datasets for validation, which does not allow for robust comparisons between methods. Consequently, the adaptation and validation of novel segmentation approaches has been slow to non-existent, with most studies still using the manufacturer's segmentation for morphometric analysis despite the existence of better performing alternative approaches. The promising accuracy of HR-pQCT for capturing morphometric indices is overshadowed by considerable variability in outcomes between studies. For finite element analysis (FEA) methods, the use of disparate material models and FEA tools has led to a fragmented ability to assess mechanical bone strength with HR-pQCT. Further, the scarcity of studies comparing 62 μm HR-pQCT to the gold standard micro-CT leaves the validation of this imaging modality incomplete. This review revealed that without standardization, the capabilities of HR-pQCT cannot be adequately assessed. The need for a public, extendable, heterogeneous dataset of HR-pQCT and corresponding gold-standard micro-CT images, which would allow HR-pQCT users to benchmark existing and novel methods and select optimal methods depending on the scientific question and data at hand, is now evident. With more recent advancements in HR-pQCT, the community must learn from its past and provide properly validated technologies to ensure that HR-pQCT can truly provide value in patient diagnosis and care.

Bone and joint enhancement filtering: Application to proximal femur segmentation from uncalibrated computed tomography datasets

Methods for reliable femur segmentation enable the execution of quality retrospective studies and building of robust screening tools for bone and joint disease. An enhance-and-segment pipeline is proposed for proximal femur segmentation from computed tomography datasets. The filter is based on a scale-space model of cortical bone with properties including edge localization, invariance to density calibration, rotation invariance, and stability to noise. The filter is integrated with a graph cut segmentation technique guided through user provided sparse labels for rapid segmentation. Analysis is performed on 20 independent femurs. Rater proximal femur segmentation agreement was 0.21 mm (average surface distance), 0.98 (Dice similarity coefficient), and 2.34 mm (Hausdorff distance). Manual segmentation added considerable variability to measured failure load and volume (CV_RMS > 5%) but not density. The proposed algorithm considerably improved inter-rater reproducibility for all three outcomes (CV_RMS < 0.5%). The algorithm localized the periosteal surface accurately compared to manual segmentation but with a slight bias towards a smaller volume. Hessian-based filtering and graph cut segmentation localizes the periosteal surface of the proximal femur with comparable accuracy and improved precision compared to manual segmentation.

Can Na18F PET/CT bone scans help when deciding if early intervention is needed in patients being treated with a TSF attached to the tibia: insights from 41 patients

Purpose

To demonstrate the usefulness of positron emission tomography (PET)/computed tomography (CT) bone scans for gaining insight into healing bone status earlier than CT or X-ray alone.

Methods

Forty-one prospective patients being treated with a Taylor Spatial Frame were recruited. We registered data obtained from successive static CT scans for each patient, to align the broken bone. Radionuclide uptake was calculated over a spherical volume of interest (VOI). For all voxels in the VOI, histograms and cumulative distribution functions of the CT and PET data were used to assess the type and progress of new bone growth and radionuclide uptake. The radionuclide uptake difference per day between the PET/CT scans was displayed in a scatter plot. Superimposing CT and PET slice data and observing the spatiotemporal uptake of ¹⁸F⁻ in the region of healing bone by a time-sequenced movie allowed qualitative evaluation.

Results

Numerical evaluation, particularly the shape and distribution of Hounsfield Units and radionuclide uptake in the graphs, combined with visual evaluation and the movies enabled the identification of six patients needing intervention as well as those not requiring intervention. Every revised patient proceeded to a successful treatment conclusion.

Conclusion

Numerical and visual evaluation based on all the voxels in the VOI may aid the orthopedic surgeon to assess a patient’s progression to recovery. By identifying slow or insufficient progress at an early stage and observing the uptake of ¹⁸F⁻ in specific regions of bone, it might be possible to shorten the recovery time and avoid unnecessary late complications.

Electronic supplementary material

The online version of this article (10.1007/s00590-020-02776-2) contains supplementary material, which is available to authorized users.

Using point clouds to investigate the relationship between trabecular bone phenotype and behavior: An example utilizing the human calcaneus

OBJECTIVES

The objective of this study is to demonstrate a new method for analyzing trabecular bone volume fraction and degree of anisotropy in three dimensions.

METHODS

We use a combination of automatic mesh registration, point-cloud correspondence registration, and P-value corrected univariate statistical tests to compare bone volume fraction and degree of anisotropy on a point by point basis across the entire calcaneus of two human groups with different subsistence strategies.

RESULTS

We found that the patterns of high and low bone volume fraction and degree of anisotropy distribution between the Black Earth (hunter-gatherers) and Norris Farms (mixed-strategy agriculturalists) are very similar, but differ in magnitude. The hunter-gatherers exhibit higher levels of bone volume fraction and less anisotropic trabecular bone organization. Additionally, patterns of bone volume fraction and degree of anisotropy in the calcaneus correspond well with biomechanical expectations of relative forces experienced during walking and running.

CONCLUSIONS

We conclude that comparing site-specific, localized differences in trabecular bone variables such as bone volume fraction and degree of anisotropy in three-dimensions is a powerful analytical tool. This method makes it possible to determine where similarities and differences between groups are located within the whole skeletal element of interest. The visualization of multiple variables also provides a way for researchers to see how the trabecular bone variables interact within the morphology, and allows for a more nuanced understanding of how they relate to one another and the broader mechanical environment.

Cortical bone mapping improves finite element strain prediction accuracy at the proximal femur

Despite evidence of the biomechanical role of cortical bone, current state of the art finite element models of the proximal femur built from clinical CT data lack a subject-specific representation of the bone cortex. Our main research hypothesis is that the subject-specific modelling of cortical bone layer from CT images, through a deconvolution procedure known as Cortical Bone Mapping (CBM, validated for cortical thickness and density estimates) can improve the accuracy of CT-based FE models of the proximal femur, currently limited by partial volume artefacts. Our secondary hypothesis is that a careful choice of cortical-specific density-elasticity relationship may improve model accuracy. We therefore: (i) implemented a procedure to include subject-specific CBM estimates of both cortical thickness and density in CT-based FE models. (ii) defined alternative models that included CBM estimates and featured a cortical-specific or an independently optimised density-elasticity relationship. (iii) tested our hypotheses in terms of elastic strain estimates and failure load and location prediction, by comparing with a published cohort of 14 femurs, where strain and strength in stance and fall loading configuration were experimentally measured, and estimated through reference FE models that did not explicitly model the cortical compartment. Our findings support the main hypothesis: an explicit modelling of the proximal femur cortical bone layer including CBM estimates of cortical bone thickness and density increased the FE strains prediction, mostly by reducing peak errors (average error reduced by 30%, maximum error and 95th percentile of error distribution halved) and especially when focusing on the femoral neck locations (all error metrics at least halved). We instead rejected the secondary hypothesis: changes in cortical density-elasticity relationship could not improve validation performances. From these improved baseline strain estimates, further work is needed to achieve accurate strength predictions, as models incorporating cortical thickness and density produced worse estimates of failure load and equivalent estimates of failure location when compared to reference models. In summary, we recommend including local estimates of cortical thickness and density in FE models to estimate bone strains in physiological conditions, and especially when designing exercise studies to promote bone strength.

Reverse wedge effect following intramedullary nailing of a basicervical trochanteric fracture variant combined with a mechanically compromised greater trochanter

Background

To introduce an unreported intraoperative complication in intramedullary nailing (IN) of an anatomically reduced trochanteric fracture variant characterized by a basicervical fracture line and coronally disrupted greater trochanter (GT).

Methods

A total of 414 trochanteric fractures (TF) treated with intramedullary nails from 2013 to 2017 were included in this study. After analysis of intraoperative fluoroscopy data, 33 cases, including 21 females and 12 males, with a mean age of 72.5 years (33 to 96 years) were identified for internal rotation of the cephalocervical fragment and inferior opening at the basicervical fracture line caused by nailing a satisfactorily reduced TF. The morphological features of this group of patients were analyzed on computed tomography (CT) scan. On radiograph, the magnitude of the displacement and final femoral neck-shaft angle (NSA) were measured.

Results

CT analysis demonstrated that the basicervical fracture line and the posterolateral fragment (PLF) detached from the GT were the two dominant features of this cohort. They were classified according to the number of main fragments: a 3-fragmentary subgroup containing three consistent fragments (cephalocervical fragment, PLF and distal femoral shaft) and a 4-fragmentary subgroup embracing one additional fragment (lesser trochanter). The four subtypes were as follows: the 3-fragmentary S indicating a small PLF (6 cases), the 3-fragmentary M presenting a moderate PLF (3 cases), the 3-fragmentary L standing for the PLF involving whole lesser trochanter (LT) (4 cases) and the 4-fragmentary GL incorporating separated PLF and LT fragments (20 cases). Geological analysis demonstrated that the majority of the basicervical fracture lines (81.8%) just crossed the center of the piriformis fossa, while the others marginally involved the medial wall of the GT. Postoperatively, the mean width of the inferior opening at the basicervical region was 9.2 ± 4.6 mm. The mean NSA was 135.2 ± 7.8 degrees. The comparison between the 3- and 4-fragmentary subgroups revealed no significant differences in magnitude of displacement and NSA.

Conclusion

This unreported intraoperative complication predominantly occurred in the intramedullary nailed basicervical trochanteric fracture variant combined with a PLF from the GT. The magnitude of the secondary displacement was substantial and resulted in a relative valgus reduction. This secondary displacement was caused by an impingement of the reamer with the superolateral cortex of the cephalocervical fragment and should be addressed during the operation.

Level of evidence

Therapy IV.

An exclusion approach for addressing partial volume artifacts with quantititive computed tomography-based finite element modeling of the proximal tibia

INTRODUCTION

Quantitative computed tomography based finite element modeling (QCT-FE) has potential to clarify the role of subchondral bone stiffness in osteoarthritis. The limited spatial resolution of clinical QCT systems, however, results in partial volume (PV) artifacts and low contrast between cortical and trabecular bone, which adversely affects the accuracy of QCT-FE models. The objective of this research was to evaluate the agreement between stiffness predictions offered by QCT-FE models of proximal tibial subchondral bone (constructed with and without a new voxel-exclusion algorithm) with experimentally-derived local subchondral bone structural stiffness.

METHODS

Thirteen proximal tibial compartments were obtained and imaged using QCT. Two types of QCT-FE models were developed: (1) standard model, which employed the standard procedure for QCT-FE modeling; and (2) "voxel exclusion (VE)" model, which addressed PV artifacts by excluding low density voxels during the material mapping stage of construction. We assessed agreement between QCT-FE stiffness estimates (using standard and VE approaches) with experimental stiffness by reporting predicted variance from linear regression and mean bias with 95% Limits of Agreement (LOA).

RESULTS

The standard and VE models explained 81% and 84% of the variance in experimentally measured stiffness, respectively. The standard model showed a mean bias of -268 N/mm (LOA -1210 to 679 N/mm); the VE model showed a mean bias of +59 N/mm (LOA -762 to 910 N/mm).

INTERPRETATION

The VE model explained more variance in subchondral bone stiffness with less bias. Our findings indicate that the VE method has potential to improve QCT-FE models of bone affected by PV artifacts.

Comparison of Baumgaertner and Chang reduction quality criteria for the assessment of trochanteric fractures

Objectives

Different criteria for assessing the reduction quality of trochanteric fractures have been reported. The Baumgaertner reduction quality criteria (BRQC) are relatively common and the Chang reduction quality criteria (CRQC) are relatively new. The objectives of the current study were to compare the reliability of the BRQC and CRQC in predicting mechanical complications and to investigate the clinical implications of the CRQC.

Methods

A total of 168 patients were assessed in a retrospective observational study. Clinical information including age, sex, fracture side, American Society of Anesthesiologists (ASA) classification, tip-apex distance (TAD), fracture classification, reduction quality, blade position, BRQC, CRQC, bone quality, and the occurrence of mechanical complications were used in the statistical analysis.

Results

A total of 127 patients were included in the full analysis, and mechanical complications were observed in 26 patients. The TAD, blade position, BRQC and CRQC were significantly associated with mechanical complications in the univariate analysis. Only the TAD (p = 0.025) and the CRQC (p < 0.001) showed significant results in the multivariate analysis. In the comparison of the receiver operating characteristic curves, the CRQC also performed better than the BRQC.

Conclusion

The CRQC are reliable in predicting mechanical complications and are more reliable than the BRQC. Future studies could use the CRQC to assess fracture reduction quality. Intraoperatively, the surgeon should refer to the CRQC to achieve good reduction in trochanteric fractures.

Cite this article: Bone Joint Res 2019;8:502–508.

A Novel Blind Restoration and Reconstruction Approach for CT Images Based on Sparse Representation and Hierarchical Bayesian-MAP

Computed tomography (CT) image reconstruction and restoration are very important in medical image processing, and are associated together to be an inverse problem. Image iterative reconstruction is a key tool to increase the applicability of CT imaging and reduce radiation dose. Nevertheless, traditional image iterative reconstruction methods are limited by the sampling theorem and also the blurring of projection data will propagate unhampered artifact in the reconstructed image. To overcome these problems, image restoration techniques should be developed to accurately correct a wide variety of image degrading effects in order to effectively improve image reconstruction. In this paper, a blind image restoration technique is embedded in the compressive sensing CT image reconstruction, which can result in a high-quality reconstruction image using fewer projection data. Because a small amount of data can be obtained by radiation in a shorter time, high-quality image reconstruction with less data is equivalent to reducing radiation dose. Technically, both the blurring process and the sparse representation of the sharp CT image are first modeled as a serial of parameters. The sharp CT image will be obtained from the estimated sparse representation. Then, the model parameters are estimated by a hierarchical Bayesian maximum posteriori formulation. Finally, the estimated model parameters are optimized to obtain the final image reconstruction. We demonstrate the effectiveness of the proposed method with the simulation experiments in terms of the peak signal to noise ratio (PSNR), and structural similarity index (SSIM).

Next-generation imaging of the skeletal system and its blood supply

Bone is organized in a hierarchical 3D architecture. Traditionally, analysis of the skeletal system was based on bone mass assessment by radiographic methods or on the examination of bone structure by 2D histological sections. Advanced imaging technologies and big data analysis now enable the unprecedented examination of bone and provide new insights into its 3D macrostructure and microstructure. These technologies comprise ex vivo and in vivo methods including high-resolution computed tomography (CT), synchrotron-based imaging, X-ray microscopy, ultra-high-field magnetic resonance imaging (MRI), light-sheet fluorescence microscopy, confocal and intravital two-photon imaging. In concert, these techniques have been used to detect and quantify a novel vascular system of trans-cortical vessels in bone. Furthermore, structures such as the lacunar network, which harbours and connects osteocytes, become accessible for 3D imaging and quantification using these methods. Next-generation imaging of the skeletal system and its blood supply are anticipated to contribute to an entirely new understanding of bone tissue composition and function, from macroscale to nanoscale, in health and disease. These insights could provide the basis for early detection and precision-type intervention of bone disorders in the future.

Comparing rib cortical thickness measurements from computed tomography (CT) and Micro-CT

BACKGROUND

The objective of this study was to compare cortical thickness of rib specimens scanned with clinical computed tomography (clinical-CT) at 0.5 and 1.0 mm slice thickness versus micro-CT at 0.05 mm slice thickness. Cortical thickness variation and accuracy was explored by anatomical region (anterior vs. lateral) and cross-sectional quadrants (superior, interior, inferior, and exterior).

METHODS

A validated cortical thickness algorithm was applied to clinical-CT and micro-CT scans of 17 rib specimens from six male post mortem human subjects aged 42-81 years. Each rib specimen was segmented and the thickness measurements were partitioned into cross-sectional quadrants in the anterior and lateral regions of the rib. Within each rib quadrant, the following were calculated: average thickness ± standard deviation, mean thickness difference between clinical-CT and micro-CT, and a thickness ratio between clinical-CT and micro-CT. Correlations from linear regression and paired-t tests were determined for paired clinical-CT and micro-CT results.

RESULTS

On average, the 0.5 mm clinical-CT underestimated the micro-CT thickness by 0.005 mm, while the 1.0 mm clinical-CT overestimated the micro-CT thickness by 0.149 mm. Thickness derived from 0.5 mm clinical-CT showed greater significant linear correlations (p < 0.05) with micro-CT thickness compared to 1.0 mm clinical-CT.

CONCLUSIONS

The small mean differences and thickness ratios near 1 show validation for the cortical thickness algorithm when applied to rib clinical-CT scans. Using clinical-CT scans as way to accurately measure rib cortical thickness offers a non-invasive way to analyze millions of CT scans collected each year from males and females of all ages.

Mandibular dysmorphology due to abnormal embryonic osteogenesis in FGFR2-related craniosynostosis mice

One diagnostic feature of craniosynostosis syndromes is mandibular dysgenesis. Using three mouse models of Apert, Crouzon and Pfeiffer craniosynostosis syndromes, we investigated how embryonic development of the mandible is affected by fibroblast growth factor receptor 2 (Fgfr2) mutations. Quantitative analysis of skeletal form at birth revealed differences in mandibular morphology between mice carrying Fgfr2 mutations and their littermates that do not carry the mutations. Murine embryos with the mutations associated with Apert syndrome in humans (Fgfr2^+/S252W and Fgfr2^+/P253R) showed an increase in the size of the osteogenic anlagen and Meckel's cartilage (MC). Changes in the microarchitecture and mineralization of the developing mandible were visualized using histological staining. The mechanism for mandibular dysgenesis in the Apert Fgfr2^+/S252W mouse resulting in the most severe phenotypic effects was further analyzed in detail and found to occur to a lesser degree in the other craniosynostosis mouse models. Laser capture microdissection and RNA-seq analysis revealed transcriptomic changes in mandibular bone at embryonic day 16.5 (E16.5), highlighting increased expression of genes related to osteoclast differentiation and dysregulated genes active in bone mineralization. Increased osteoclastic activity was corroborated by TRAP assay and in situ hybridization of Csf1r and Itgb3. Upregulated expression of Enpp1 and Ank was validated in the mandible of Fgfr2^+/S252W embryos, and found to result in elevated inorganic pyrophosphate concentration. Increased proliferation of osteoblasts in the mandible and chondrocytes forming MC was identified in Fgfr2^+/S252W embryos at E12.5. These findings provide evidence that FGFR2 gain-of-function mutations differentially affect cartilage formation and intramembranous ossification of dermal bone, contributing to mandibular dysmorphogenesis in craniosynostosis syndromes.

This article has an associated First Person interview with the joint first authors of the paper.

Summary: FGFR2 gain-of-function mutations differentially affect cartilage formation and intramembranous ossification of dermal bone, resulting in abnormal embryonic osteogenesis of the mandible.

Prediction of lumbar vertebral body compressive strength of overweight and obese older adults using morphed subject-specific finite-element models to evaluate the effects of weight loss

BACKGROUND

Diet and exercise can promote weight loss in older adults; however, there is potential to increase fracture risk due to loss of bone mineral density (BMD) known to accompany weight loss. Weight loss effects on measures of bone quality and strength are currently unknown.

AIMS

The purpose of this study is to develop subject-specific finite-element (FE) models of the lumbar spine and study the effect of intentional weight loss on bone strength in a pilot data set.

METHODS

Computed tomography (CT) scans of the lumbar spine of 30 overweight and obese (mean BMI = 29.7 ± 3.9 kg/m²), older adults (mean age = 65.9 ± 4.6 years) undergoing an 18-month intentional weight loss intervention were obtained at baseline and post-intervention. Measures of volumetric BMD (vBMD) and variable cortical thickness were derived from each subject CT scan. Development of the subject-specific FE models of the lumbar spine involved model morphing techniques to accelerate the development of the models. vBMD-derived material properties and cortical thickness measures were directly mapped to baseline and post-intervention models. Bone strength was estimated through simulation of a quasi-static uniaxial compression test.

RESULTS

From baseline to 18-month post-weight loss intervention, there were statistically significant decreases in estimated bone strength (6.5% decrease; p < 0.05). Adjusting for baseline bone measures and gender revealed no statistically significant correlations between weight change and change in vBMD, cortical thickness, or bone strength.

CONCLUSION

Integration of CT-based measures and FE models with conventional areal BMD can improve the understanding of the effects of intentional weight loss on bone health.

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.

QCT of the femur: Comparison between QCTPro CTXA and MIAF Femur

QCT is commonly employed in research studies and clinical trials to measure BMD at the proximal femur. In this study we compared two analysis software options, QCTPro CTXA and MIAF-Femur, using CT scans of the semi-anthropometric European Proximal Femur Phantom (EPFP) and in vivo data from 130 Chinese elderly men and women aged 60-80 years. Integral (Int), cortical (Cort) and trabecular (Trab) vBMD, volume, and BMC of the neck (FN), trochanter (TR), inter-trochanter (IT), and total hip (TH) VOIs were compared. Accuracy was determined in the 5 mm wide central portion of the femoral neck of the EPFP. Nominal values were: cross-sectional area (CSA) 4.9 cm², cortical thickness (C.Th) 2 mm, CortBMD 723 mg/cm³ and TrabBMD 100 mg/cm³. In MIAF the so-called peeled trabecular VOI was analyzed, which excludes subcortical bone to avoid partial volume artefacts at the endocortical border that artificially increase TrabBMD. For CTXA uncorrected, so called raw cortical values were used for the analysis. QCTPro and MIAF phantom results were: CSA 5.9 cm² versus 5.1 cm²; C.Th 1.68 mm versus 1.92 mm; CortBMD 578 mg/cm³ versus 569 mg/cm³; and TrabBMD 154 mg/cm³ versus 104 mg/cm³. In vivo correlations (R²) of integral and trabecular bone parameters ranged from 0.63 to 0.96. Bland-Altman analysis for TH and FN TrabBMD showed that lower mean values were associated with higher differences, which means that TrabBMD differences between MIAF and CTXA are larger for osteoporotic than for normal patients, which can be largely explained by the inclusion of subcortical BMD in the trabecular VOI analyzed by CTXA in combination with fixed thresholds used to separate cortical from trabecular bone compartments. Correlations between CTXA corrected CortBMD and MIAF were negative, whereas raw data correlated positively with MIAF measurements for all VOIs questioning the validity of the CTXA corrections. The EPFP results demonstrated higher MIAF accuracy of cortical thickness and TrabBMD. Integral and trabecular bone parameters were highly correlated between CTXA and MIAF. Partial volume artefacts at the endocortical border artificially increased trabecular BMD by CTXA, especially for osteoporosis patients. With respect to volumetric cortical measurements with CTXA, the use raw data is recommended, because corrected data cause a negative correlation with MIAF CortBMD.

Effect of Exercise Modality During Weight Loss on Bone Health in Older Adults With Obesity and Cardiovascular Disease or Metabolic Syndrome: A Randomized Controlled Trial

The objective of this study was to determine the ability of either aerobic or resistance training to counter weight-loss-associated bone loss in older adults. There were 187 older adults (67 years, 70% women, 64% white) with obesity (BMI = 34.5 ± 3.7 kg/m² ) and cardiovascular disease and/or metabolic syndrome who were randomized to participate in an 18-month, community-based trial, with a follow-up assessment at 30 months. Intervention arms included: weight loss alone (WL; 7% to 10% baseline weight), WL plus aerobic training (WL + AT), and WL plus resistance training (WL + RT), as well as DXA-acquired total hip, femoral neck, and lumbar spine areal bone mineral density (aBMD), and trabecular bone score (TBS). Biomarkers of bone turnover (procollagen type 1 N-terminal propeptide, C-terminal telopeptide of type 1 collagen) were measured at baseline, 6, 18, and 30 (aBMD and TBS only) months. CT-acquired hip and spine volumetric BMD (vBMD), cortical thickness, and bone strength were measured in a subset at baseline (n = 55) and 18 months. Total hip aBMD was reduced by 2% in all groups at 18 months, with a primary analysis showing no significant treatment effects for any DXA, biomarker, or CT outcome. After adjustment for WL and follow-up at 30 months, secondary analyses revealed that total hip [-0.018 (-0.023 to -0.012) g/cm² versus -0.025 (-0.031 to -0.019) g/cm² ; p = 0.05] and femoral neck [-0.01 (-0.009 to 0.008) g/cm² versus -0.011 (-0.020 to -0.002) g/cm² ; p = 0.06] aBMD estimates were modestly attenuated in the WL + RT group compared with the WL group. Additionally, lumbar spine aBMD was increased in the WL [0.015 (0.007 to 0.024) g/cm² ] and the WL + RT [0.009 (0.000 to 0.017) g/cm² ] groups compared with the WL + AT [-0.003 (-0.012 to 0.005)g/cm² ] group; both p ≤ 0.01. Community-based exercise does not prevent bone loss during active WL in older adults; however, adding RT may help minimize long-term hip bone loss. © 2018 American Society for Bone and Mineral Research.

Sources of Variability in Structural Bending Response of Pediatric and Adult Human Ribs in Dynamic Frontal Impacts

Despite safety advances, thoracic injuries in motor vehicle crashes remain a significant source of morbidity and mortality, and rib fractures are the most prevalent of thoracic injuries. The objective of this study was to explore sources of variation in rib structural properties in order to identify sources of differential risk of rib fracture between vehicle occupants. A hierarchical model was employed to quantify the effects of demographic differences and rib geometry on structural properties including stiffness, force, displacement, and energy at failure and yield. Three-hundred forty-seven mid-level ribs from 182 individual anatomical donors were dynamically (~2 m/s) tested to failure in a simplified bending scenario mimicking a frontal thoracic impact. Individuals ranged in age from 4 - 108 years (mean 53 ± 23 years) and included 59 females and 123 males of diverse body sizes. Age, sex, body size, aBMD, whole rib geometry and cross-sectional geometry were explored as predictors of rib structural properties. Measures of cross-sectional rib size (Tt.Ar), bone quantity (Ct.Ar), and bone distribution (Z) generally explained more variation than any other predictors, and were further improved when normalized by rib length (e.g., robustness and WBSI). Cortical thickness (Ct.Th) was not found to be a useful predictor. Rib level predictors performed better than individual level predictors. These findings moderately explain differential risk for rib fracture and with additional exploration of the rib's role in thoracic response, may be able contribute to ATD and HBM development and alterations in addition to improvements to thoracic injury criteria and scaling methods.

Cortical Bone Mapping: Measurement and Statistical Analysis of Localised Skeletal Changes

PURPOSE OF REVIEW

Cortical bone mapping (CBM) is a technique for measuring localised skeletal changes from computed tomography (CT) images. It can provide measurements with accuracy surpassing the underlying imaging resolution. CBM can detect changes in several properties of the cortex, with no prior assumptions about the likely location of said changes. This paper summarises the theory behind CBM, discusses its strengths and limitations, and reviews some studies in which it has been applied.

RECENT FINDINGS

CBM has revealed associations between fracture risk and cortical properties in specific regions of the proximal femur which present feasible therapeutic targets. Analyses of several pharmaceutical and exercise interventions quantify effects that are distinct both in location and in the nature of the micro-architectural changes. CBM has illuminated age-related changes in the proximal femur and has recently been applied to other bones, as well as to the assessment of cartilage. The CBM processing pipeline is designed primarily for large cohort studies. Its main impact thus far has not been in the realm of clinical practice, but rather to improve our fundamental understanding of localised bone structure and changes.

Femoral neck cortical bone in female and male hip fracture cases: Differential contrasts in cortical width and sub-periosteal porosity in 112 cases and controls

OBJECTIVES

To quantitate differences between cases of hip fracture and controls in cortical width around the mid-femoral neck in men and women.

METHODS

Over 5 years, 64 (14 male) participants over age 55 (mean 79) years, who had never taken bone-active drugs and suffered intra-capsular hip fracture treated by arthroplasty, donated their routinely discarded distal intra-capsular femoral neck bone for histomorphometry. After embedding, complete femoral neck cross sections from the cut surface near the narrowest part of the neck were stained with von Kossa and cortical width was measured radially every 5 degrees of arc. Control material (n = 48, 25 male) was available through consented post mortems prior to the year 2000. Cortical widths were averaged for circumferential octants, each representing 45 degrees of arc. Divergence of individual cortical widths from their means was also examined.

RESULTS

Because sections were required to have a complete cortex, sampling was biased towards cases with sub-capital versus trans-cervical fractures. Compared to sex- and age matched controls, male cases showed larger relative differences in cortical widths than female cases. Unexpectedly, cortical widths in female but not male cases also showed marked over-representation of extremely narrow (<0.1 mm) cortical widths, located mainly posteriorly. The numbers of these very narrow cortical widths observed per subject retrospectively predicted female fracture status in logistic regression independently of mean cortical width values. Together with mean cortical width differences, the numbers of measured cortical widths <0.1 mm (out of 72 measured) raised the sensitivity of predicting fracture status in women from 48 to 80% at 80% specificity. In almost all cases, very narrow cortical widths were identified in regions enclosing a cortical pore roofed on its endosteal surface by thin structural bone defined a priori as trabecular.

CONCLUSIONS

Cortical widths <0.1 mm probably reflect zones where endosteal cortex has been trabecularised through expansion of an un-refilled sub-endosteal canal close to the periosteum. Persistent cortical defects occurring near the periosteal surface, where mechanical loading exerts its greatest stresses, are likely to result in extremes of localized concentrations of stress during a fall, unknown in young normal fallers. Such defects have the potential to help explain the excess of hip fractures among elderly women. Prevention of sub-periosteal tunnelling by osteoclasts might explain in part the additional benefits, beyond an increase in bone density, of treatments that reduce excessive bone resorption or else stimulate new bone formation on previously resorbed surfaces.

Fracture Prediction by Computed Tomography and Finite Element Analysis: Current and Future Perspectives

PURPOSE OF REVIEW

This review critiques the ability of CT-based methods to predict incident hip and vertebral fractures.

RECENT FINDINGS

CT-based techniques with concurrent calibration all show strong associations with incident hip and vertebral fracture, predicting hip and vertebral fractures as well as, and sometimes better than, dual-energy X-ray absorptiometry areal biomass density (DXA aBMD). There is growing evidence for use of routine CT scans for bone health assessment. CT-based techniques provide a robust approach for osteoporosis diagnosis and fracture prediction. It remains to be seen if further technical advances will improve fracture prediction compared to DXA aBMD. Future work should include more standardization in CT analyses, establishment of treatment intervention thresholds, and more studies to determine whether routine CT scans can be efficiently used to expand the number of individuals who undergo evaluation for fracture risk.