Independent measurement of femoral cortical thickness and cortical bone density using clinical CT

Medial cortical bone thickness of the tibial diaphysis in osteoarthritis is related to lower extremity alignment and tibial morphology

BACKGROUND

The purpose of this study was to clarify (1) the differences in cortical bone thickness (CBT) of the tibial diaphysis between healthy and osteoarthritic knees and (2) the influences of the femorotibial angle (FTA) and inclination of the medial compartment of the proximal tibia (MCT) on tibial CBT.

METHODS

The study assessed 60 subjects with varus knee osteoarthritis (OA) (22 males and 38 females; mean age, 74 ± 7 years) and 53 healthy elderly subjects (28 males and 25 females; mean age, 70 ± 6 years). Three-dimensional estimated CBT of the tibial diaphysis was automatically calculated for 2752-11,296 points using high-resolution measurements from CT. The standardized CBT was assessed in 24 regions by combining six heights and four areas. Additionally, the association between the CBT, each FTA, and MCT inclination was investigated.

RESULTS

The OA group showed a thicker CBT in the medial areas than in the lateral areas of the proximal tibia, while the healthy group had a thicker lateral CBT. The medial-to-lateral ratio of the proximal tibia was significantly higher in the OA group than in the healthy group. The proximal-medial CBT correlated with FTA and MCT inclinations in the OA group.

CONCLUSIONS

This study demonstrated that varus osteoarthritic knees showed a different trend of proximal-medial CBT with associations in FTA and MCT inclination from healthy knees, possibly due to medial load concentration.

Mapping the Spatial Evolution of Proximal Femur Osteoporosis: A Retrospective Cross-Sectional Study Based on CT Scans

Purpose

The purpose of this study was to quantify the modifications occurring in osteoporosis at the level of the human proximal femur throughout the trabecular structure, along with the identification of certain anatomic regions preferentially affected by osteoporosis. Another goal was to map the evolution of the radiodensity of the trabecular bone as osteoporosis progresses to an advanced stage.

Methods

The study included CT scans (right femur) from 51 patients, out of which 40 had various degrees of osteoporosis, but no other local pathology. Ten regions of interest in two orthogonal slices have been identified and the differences in radiodensity as well as their evolution have been statistically analyzed in terms of relative and absolute changes.

Results

A detailed spatial map showing the evolution of osteoporosis was obtained. As osteoporosis evolved, the relative decrease in radiodensity was inversely correlated to the radiodensity of the healthy bone. In particular, the region covering the Ward triangle decreased the most, by an average 61-62% in osteopenia and 101-106% in advanced osteoporosis, while the principal compressive group was affected the least, showing a decrease by an average 14-15% in osteopenia and 29-32% in advanced osteoporosis. The absolute decrease in radiodensity was not correlated to the radiodensity of the healthy bone and was shifted to the inferior-posterior edge of the femur. Inside the femoral head, the upper region was affected the most in absolute terms, while the greater trochanter was less affected than the femoral neck. The maximum metaphyseal cortical bone density was unaffected by the progression of osteoporosis.

Conclusion

Significant differences were noticed in terms of the absolute and relative osteoporotic changes in radiodensity related to different anatomical regions of the human femoral bone. These differences become more pronounced as the disease progresses.

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.

Numerical investigation of intra-abdominal pressure and spinal load-sharing upon the application of an abdominal belt

Chronic low back pain patients may experience spinal instability. Abdominal belts (ABs) have been shown to improve spine stability, trunk stiffness, and resiliency to spinal perturbations. However, research on the contributing mechanisms is inconclusive. ABs may increase intra-abdominal pressure (IAP) and reduce paraspinal soft tissue contribution to spine stability without increasing spinal compressive loads. A finite element model (FEM) of the spine inclusive of the T1-S1 vertebrae, intervertebral discs (IVDs), ribcage, pelvis, soft tissues, and abdominal cavity, without active muscle forces was developed. An identical FEM with an AB was developed. Both FEMs underwent trunk flexion. Following validation, the models' intervertebral rotation (IVR), IAP, IVD pressure, and tensile stress in the multifidus (MF), erector spinae (ES), and thoracolumbar fascia (TLF) were compared. The inclusion of an AB resulted in a 3.8 kPa IAP increase, but a decreased average soft tissue tensile stress of 0.28 kPa. The TLF withstood the majority of tension being transferred across the paraspinal soft tissues (>70 %). The average IVR in the AB model decreased by 10 %, with the lumbar spine experiencing the largest reduction. The lumbar IVDs of the AB model likewise showed a 31 % reduction in average IVD pressure. Using an AB improved trunk bending stiffness, primarily in the lumbar spine. Wearing an AB had minimal effect on reducing tensile stress in theES. The skewed stress distribution towards the TLF suggests its large contribution to spine stability and the potential advantage in unloading the structure when wearing an AB, measured herein at8 %.

Thinner femoral cortical thickness in patients with destructive rheumatoid arthritis of the knee

BACKGROUND

The examination of femoral cortical bone thickness in patients with rheumatoid arthritis (RA) has been notably limited in prior research. We aimed to compare femoral cortical thickness in patients with rheumatoid arthritis (RA) and healthy controls and to investigate the association between femoral cortical thickness and clinical parameters within the RA group.

METHODS

Forty-four patients (58 limbs) with RA who underwent total knee arthroplasty were enrolled. Preoperative computed tomography images of the lower limbs were analyzed. The femoral cortex was divided into the proximal, central, and distal diaphysis regions and further into the anterior, posterior, medial, and lateral regions. The divisions were measured using Stradwin® software and standardized by femoral length. Femoral cortical thickness was compared between RA and healthy control (n = 25) groups. Correlation analyses between standardized cortical thickness and disease parameters were performed in the RA group.

RESULTS

The RA group had significantly lower standardized femoral cortical thickness at the anterior and medial distal diaphysis than healthy controls. Standardized proximal lateral and central lateral in the RA group were significantly larger than those in the healthy control groups. Standardized femoral cortical thickness was significantly correlated with bone mineral density (BMD) in 11 areas, except the posterior central diaphysis, and with body mass index in 8 areas, except the central posterior, distal lateral, distal anterior, and distal medial diaphysis.

CONCLUSIONS

Femoral cortical thinning was noted in patients with RA complicated with destructive knee, particularly at the anterior and medial distal diaphysis. Femoral cortical thickness was significantly correlated with BMD and body mass index (BMI); thus, patients with RA and low BMD and BMI should be cared for to prevent fragility fractures.

Cortical and Trabecular Bone Deficit in Middle-Aged Men Living with HIV

A significant increase in the risk of hip fracture was observed in middle-aged men living with human immunodeficiency virus (MLWH), almost a decade earlier than those without infection. Data regarding cortical and trabecular bone deficit of hip, an important determinant of bone strength, in MLWH are limited. Quantitative CT was performed in consecutive MLWH aged ≥30 years between November 2017 and October 2018 at Severance Hospital, Seoul, Korea. Volumetric bone mineral density (vBMD) and cortical bone mapping parameters of hip (cortical thickness [CTh], cortical bone vBMD [CBMD], cortical mass surface density [CMSD], endocortical trabecular density [ECTD]) were compared to age-matched and body mass index (BMI)-matched controls (1:2) using a community-based healthy adults cohort. Among 83 MLWH and 166 controls (mean age: 47.2 years; BMI: 23.6 kg/m2 ), MLWH had lower total hip vBMD (280 ± 41 versus 296 ± 41 mg/cm3 ), CMSD (155 versus 160 mg/cm2 ), and ECTD (158 versus 175 mg/cm3 ) than controls that remained robust after adjustment for covariates (adjusted β: total hip vBMD, -18.8; CMSD, -7.3; ECTD, -18.0; p < 0.05 for all). Cortical bone mapping revealed localized deficit of CTh, CBMD, and CMSD in the anterolateral trochanteric region and femoral neck in MLWH compared to controls, with a more extensive ECTD deficit. In MLWH, lower CD4 T-cell count (/100 cells/mm3 decrement) and protease inhibitor (PI)-based regimen (versus non-PI regimen) at the time of antiretroviral treatment initiation were associated with lower total hip vBMD (adjusted β -7.5 for lower CD4 count; -28.3 for PI-based regimen) and CMSD (adjusted β -2.6 for lower CD4 count; -12.7 for PI-based regimen; p < 0.05 for all) after adjustment for covariates including age, BMI, smoking, alcohol use, hepatitis C virus co-infection, tenofovir exposure, and CT scanner types. MLWH had lower hip bone density with cortical and trabecular bone deficit compared to community-dwelling controls. © 2023 American Society for Bone and Mineral Research (ASBMR).

Strategies to reduce the onset of sleeve gastrectomy associated bone loss (STRONG BONES): Trial design and methods

Background

Despite recognized improvements in obesity-related comorbidities, mounting evidence implicates surgical weight loss in the onset of skeletal fragility. Sleeve gastrectomy (SG) is the most commonly performed bariatric procedure and is associated with 3-7% axial bone loss in the year following surgery. Bisphosphonates are FDA-approved medications for the prevention and treatment of age-related bone loss and may represent a strategy to reduce bone loss following SG surgery.

Methods

The Strategies to Reduce the Onset of Sleeve Gastrectomy Associated Bone Loss (STRONG BONES) trial (NCT04922333) is designed to definitively test whether monthly administration of the bisphosphonate, risedronate, for six months can effectively counter SG-associated bone loss. Approximately 120 middle-aged and older (≥40 years) SG patients will be randomized to six months of risedronate or placebo treatment, with skeletal outcomes assessed at baseline, six, and 12-months post-surgery. The primary outcome of the trial is 12-month change in total hip areal bone mineral density (aBMD), measured by dual energy x-ray absorptiometry (DXA). This will be complemented by DXA-acquired aBMD assessment at other skeletal sites and quantitative computed tomography (QCT) derived changes in bone quality. Change in muscle mass and function will also be assessed, as well as biomarkers of bone health, turnover, and crosstalk, providing mechanistic insight into intervention-related changes to the bone-muscle unit.

Discussion

Results from the STRONG BONES trial have the potential to influence current clinical practice by determining the ability of bisphosphonate use to mitigate bone loss and concomitant fracture risk in middle-aged and older SG patients.

Does Matching Femoral Size and Shape Improve Bone Fit and Patient-reported Outcomes in TKA? A Matched Controlled Study

BACKGROUND

Some researchers have suggested that achieving good component coverage over the host bone during TKA (while avoiding implant overhang) may help achieve durable implant fixation and may be associated with better outcomes scores. However, the evidence about this is limited and contradictory. Contemporary morphometric TKA includes a wide array of components with various shapes and sizes, based on large anatomic databases and specific software that simulates bone cuts. Morphometric tibial components have shown improved bone coverage and better clinical outcomes than standardized implants, but the role of morphometric femoral components in bone coverage has not been studied precisely.

QUESTIONS/PURPOSES

In a retrospective, controlled study that used patient matching, we asked: (1) Does the use of a contemporary morphometric component with more available sizes provide better femoral component fit and bone coverage than an earlier design with fewer sizes? (2) Are component fit and the presence of component overhang or underhang associated with different Knee Society Score (KSS) or Knee Injury and Osteoarthritis Outcome Score (KOOS) for Joint Replacement?

METHODS

From 2012 to 2013, we performed 403 TKAs according to the following criteria: TKA performed for primary tricompartmental arthritis of the knee; varus, valgus, and flexion deformity less than 15°; and age between 18 and 85 years on the day of surgery. Among these 403 TKAs, 237 were performed using a morphometric implant and 166 with the earlier nonmorphometric implant. At 2 years of follow-up, 3% of patients in the morphometric group and 5% in the nonmorphometric group were lost to follow-up. Based on age, BMI, gender, and preoperative KSS and KOOS, two groups of 30 patients were matched in a 1:1 ratio from this longitudinally maintained database. Clinical outcomes were measured preoperatively and at a minimum follow-up of 2 years in both groups, using the KSS and KOOS. We evaluated postoperative CT images for each patient to analyze femoral implant rotation, bone coverage, and overhang and underhang status.

RESULTS

The overhang status was similar between the two groups (23% had an overhang component in the morphometric knee group and 27% had an overhang component in the nonmorphometric knee group), and overhang was most frequently found in the lateral distal zone and medial anterior chamfer. Better cortical bone coverage was found in the morphometric knee group, with a thinner bone margin between the component edge and cortical border (morphometric group: 3 mm versus nonmorphometric knee group: 5 mm; p = 0.01). In general, there were few between-group differences in terms of patient-reported outcomes; of the seven metrics we analyzed, only the KSS favored the morphometric knee implant by a margin larger than the minimum clinically important difference (KSS mean difference: 21 points for the morphometric knee group; p < 0.05). Overhang of the femoral component of > 2 mm was associated with poorer KOOS, but not KSS, whereas a thinner bone margin had a beneficial impact on pain and global clinical scores (KOOS and KSS: p < 0.05).

CONCLUSION

The use of a morphometric femoral component design showed slightly improved bone fit and pain score according to the KSS at midterm follow-up compared with earlier implants with fewer sizes. Overhang > 2 mm was associated with worse KOOS. The tendency toward better outcomes in morphometric implants warrants longer-term evaluation before any definite conclusions about the association between bone fit and clinical results can be drawn.Level of Evidence Level III, therapeutic study.

Indications of musculoskeletal health in deceased male individuals with lower-limb amputations: comparison to non-amputee and diabetic controls

Individuals with lower-limb amputations, many of whom have type 2 diabetes, experience impaired musculoskeletal health. This study: (1) compared residual and intact limbs of diabetic and non-diabetic post-mortem individuals with amputation to identify structures vulnerable to injury, and (2) compared findings to diabetic and healthy control groups to differentiate influences of amputation and diabetes on musculoskeletal health. Postmortem CT scans of three groups, ten individuals each, were included: (1) individuals with transtibial or transfemoral amputations, half with diabetes (2) diabetic controls, and (3) healthy controls. Hip and knee joint spaces, cross-sectional thigh muscle and fat areas, and cross-sectional bone properties (e.g. area, thickness, geometry) were measured. Wilcoxon Signed-Rank and Kruskal-Wallis tests assessed statistical significance. Asymmetry percentages between limbs assessed clinical significance. Residual limbs of individuals with amputation, particularly those with diabetes, had significantly less thigh muscle area and thinner distal femoral cortical bone compared to intact limbs. Compared to control groups, individuals with amputation had significantly narrower joint spaces, less thigh muscle area bilaterally, and thinner proximal femoral cortical bone in the residual limb. Diabetic individuals with amputation had the most clinically significant asymmetry. Findings tended to align with those of living individuals. However, lack of available medical information and small sample sizes reduced the anticipated clinical utility. Larger sample sizes of living individuals are needed to assess generalizability of findings. Quantifying musculoskeletal properties and differentiating influences of amputation and diabetes could eventually help direct rehabilitation techniques.

Cross-sectional properties of rib geometry from an adult population

Introduction: Human body models (HBMs) play a key role in improving modern vehicle safety systems to protect broad populations. However, their geometry is commonly derived from single individuals chosen to meet global anthropometric targets, thus their internal anatomy may not fully represent the HBM's target demographic. Past studies show sixth rib cross-sectional geometry differences between HBM ribs and population-derived ribs, and corrections to HBM ribs based on these data have improved HBM's abilities to predict rib fracture locations. Methods: We measure and report average and standard deviations (SDs) in rib cross-sectional geometric properties derived from live subject CT scans of 240 adults aged 18-90. Male and female results are given as functions of rib number and rib lengthwise position for ribs 2 through 11. Population means/SDs are reported for measures of rib total area, rib cortical bone area, and rib endosteal area, as well as inertial moment properties of these rib sections. These population corridors are compared between males and females, and against the baseline rib geometries defined in six current HBMs. Results: Total cross-sectional area results found average males ribs to be larger than those of females by between approximately 1-2 SDs depending on rib number and position, and larger in cortical bone cross-sectional area by between 0-1 SDs. Inertial moment ratios showed female ribs being between approximately 0-1 SDs more elongated than male ribs, dependent again on rib number and position. Rib cross-sectional areas from 5 of the 6 HBMs were found to be overly large along substantial portions of most ribs when compared to average population corridors. Similarly, rib aspect ratios in HBMs deviated from average population data by up to 3 SDs in regions towards sternal rib ends. Discussion: Overall, while most HBMs capture overall trends such as reductions in cross-section along shaft lengths, many also exhibit local variation that deviates from population trends. This study's results provide the first reference values for assessing the cross-sectional geometry of human ribs across a wide range of rib levels. Results also further provide clear guidelines to improve rib geometry definitions present in current HBMs in order to better represent their target demographic.

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

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

Prediction of femoral strength of elderly men based on quantitative computed tomography images using machine learning

Hip fracture is the most common complication of osteoporosis, and its major contributor is compromised femoral strength. This study aimed to develop practical machine learning models based on clinical quantitative computed tomography (QCT) images for predicting proximal femoral strength. Eighty subjects with entire QCT data of the right hip region were randomly selected from the full MrOS cohorts, and their proximal femoral strengths were calculated by QCT-based finite element analysis (QCT/FEA). A total of 50 parameters of each femur were extracted from QCT images as the candidate predictors of femoral strength, including grayscale distribution, regional cortical bone mapping (CBM) measurements, and geometric parameters. These parameters were simplified by using feature selection and dimensionality reduction. Support vector regression (SVR) was used as the machine learning algorithm to develop the prediction models, and the performance of each SVR model was quantified by the mean squared error (MSE), the coefficient of determination (R² ), the mean bias, and the SD of bias. For feature selection, the best prediction performance of SVR models was achieved by integrating the grayscale value of 30% percentile and specific regional CBM measurements (MSE ≤ 0.016, R² ≥ 0.93); and for dimensionality reduction, the best prediction performance of SVR models was achieved by extracting principal components with eigenvalues greater than 1.0 (MSE ≤ 0.014, R² ≥ 0.93). The femoral strengths predicted from the well-trained SVR models were in good agreement with those derived from QCT/FEA. This study provided effective machine learning models for femoral strength prediction, and they may have great potential in clinical bone health assessments.

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.

Spinal stability analysis of lumbar interbody fusion according to pelvic type and cage angle based on simplified spinal model with various pelvic indices

Recently, the objectives of lumbar interbody fusion (LIF) have been extended to include the correction of broader/relative indications in addition to spinal fixation. Accordingly, LIF must be optimized for sagittal alignment while simultaneously achieving decompression. Therefore, a representative model classified into three pelvic types, i.e., neutral pelvis (NP), anterior pelvis (AP), and retroverted pelvis (RP), was selected according to the pelvic index, and LIF was performed on each representative model to analyze Lumbar lordosis (LL) and the corresponding equivalent stress. The finite element (FE) model was based on a sagittal 2D X-ray image. The calculation efficiency and convergence were improved by simplifying the modeling of the vertebral body in general and its posterior portion in particular. Based on the position of the pelvis, according to the pelvic shape, images of patients were classified into three types: AP, RP, and NP. Subsequently, representative images were selected for each type. The fixation device used in the fusion model was a pedicle screw and a spinal rod of a general type. PEEK was used as the cage material, and the cage shape was varied by using three different cage angles: 0°, 4°, and 8°. Spinal mobility: The pelvic type with the highest range of motion (ROM) for the spine was the NP type; the AP type had the highest LL. Under a combination load, the NP type exhibited the highest lumbar flexibility (LF), which was 2.46° lower on average compared to the case where a pure moment was applied. Equivalent stress on the spinal fixation device: The equivalent stress acting on the vertebrae was lowest when cage 0 was used for the NP and AP type. For the RP type, the lowest equivalent stress on the vertebrae was observed when cage 4 was used. Finally, for the L5 upper endplate, the stress did not vary significantly for a given type of cage. In conclusion, there was no significant difference in ROM according to cage angle, and the highest ROM, LL and LF were shown in the pelvic shape of NP type. However, when comparing the results with other pelvic types, it was not possible to confirm that LF is completely dependent on LL and ROM.

A Review of CT-Based Fracture Risk Assessment with Finite Element Modeling and Machine Learning

PURPOSE OF REVIEW

We reviewed advances over the past 3 years in assessment of fracture risk based on CT scans, considering methods that use finite element models, machine learning, or a combination of both.

RECENT FINDINGS

Several studies have demonstrated that CT-based assessment of fracture risk, using finite element modeling or biomarkers derived from machine learning, is equivalent to currently used clinical tools. Phantomless calibration of CT scans for bone mineral density enables accurate measurements from routinely taken scans. This opportunistic use of CT scans for fracture risk assessment is facilitated by high-quality automated segmentation with deep learning, enabling workflows that do not require user intervention. Modeling of more realistic and diverse loading conditions, as well as improved modeling of fracture mechanisms, has shown promise to enhance our understanding of fracture processes and improve the assessment of fracture risk beyond the performance of current clinical tools. CT-based screening for fracture risk is effective and, by analyzing scans that were taken for other indications, could be used to expand the pool of people screened, therefore improving fracture prevention. Finite element modeling and machine learning both provide valuable tools for fracture risk assessment. Future approaches should focus on including more loading-related aspects of fracture risk.

Rib cortical bone thickness variation in adults by age and sex

Rib fractures are a common and serious outcome of blunt thoracic trauma and their likelihood is greater in older individuals. Osteoporotic bone loss is a well-documented aging phenomenon with sex-specific characteristics, but within rib bones, neither baseline maps of regional thickness nor the rates of bone thinning with age have been quantified across whole ribs. This study presents such data from 4014 ribs of 240 adult subjects aged 20-90. A validated cortical bone mapping technique was applied to clinical computed tomography scans to obtain local rib cortical bone thickness measurements over the surfaces of ribs 2 through 11. Regression models to age and sex gave rates of cortex thinning in local zones and aggregated across whole ribs. The statistical parametric mapping provided these relationships regionally as a function of rib surface location. All models showed significant reductions in bone thickness with age (p < 0.01). Average whole-rib thinning occurred at between 0.011 to 0.032 mm/decade (males) and 0.035 to 0.043 mm/decade (females), with sex and age accounting for up to 37% of population variability (R² ). Rates of thinning differed regionally and by rib, with the highest bone loss of up to 0.074 mm/decade occurring in mid-rib cutaneous and superior regions of ribs 2-6. Rates were consistently higher in females than males (significantly so across whole ribs but not all local regions) and were more pronounced in cutaneous, superior, and inferior rib aspects (average 0.025 mm/decade difference in ribs 4-8) compared to pleural aspects which had the thickest cortices but saw only minor differences in thinning rates by sex (0.045 mm/decade for females and 0.040 mm/decade for males). Regional analysis showed male and female bone thickness differences that were not statistically significant at 20 years of age (p > 0.05 across practically all regions) but subsequent cortex thinning meant that substantial pleural and cutaneous regions were thinner (p < 0.05) in females than males by 55 years of age. The techniques and results from this study can be applied to assess rib bone content loss in clinical settings across wide populations. Additionally, average cortex thickness results can be mapped directly to finite element models of the thorax, and regression results are used to modify such models to represent the ribs of men and women across their full adult lifespan.

Average Lumbar Hounsfield Units Predicts Osteoporosis-Related Complications Following Lumbar Spine Fusion

STUDY DESIGN

Retrospective Study.

OBJECTIVE

To compare methods of assessing pre-operative bone density to predict risk for osteoporosis related complications (ORC), defined as proximal junctional kyphosis, pseudarthrosis, accelerated adjacent segment disease, reoperation, compression fracture, and instrument failure following spine fusions.

METHODS

Chart review of primary posterior thoracolumbar or lumbar fusion patients during a 7 year period. Inclusion criteria: preoperative dual-energy x-ray absorptiometry (DXA) test within 1 year and lumbar CT scan within 6 months prior to surgery with minimum of 1 year follow-up. Exclusion criteria: <18 years at time of index procedure, infection, trauma, malignancy, skeletal dysplasia, neuromuscular disorders, or anterior-posterior procedures.

RESULTS

140 patients were included. The average age was 67.9 years, 83 (59.3%) were female, and 45 (32%) had an ORC. There were no significant differences in patient characteristics between those with and without an ORC. Multilevel fusions were associated with ORCs (46.7% vs 26.3%, p = 0.02). Patients with ORCs had lower DXA t-scores (-1.62 vs -1.10, p = 0.003) and average Hounsfield units (HU) (112.1 vs 148.1, p ≤ 0.001). Multivariable binary logistic regression analysis showed lower average HU (Adj. OR 0.00 595% CI 0.0001-0.1713, p = 0.001) was an independent predictor of an ORC. The odds of an ORC increased by 1.7-fold for every 25 point decrease in average HU.

CONCLUSIONS

The gold standard for assessing bone mineral density has been DXA t-scores, but the best predictor of ORC remains unclear. While both lower t-scores and average HU were associated with ORC, only HU was an independent predictor of ORC.

Spatial assessment of femoral neck bone density and microstructure in hip osteoarthritis

Osteoarthritis (OA) is known to involve profound changes in bone density and microstructure near to, and even distal to, the joint. Critically, however, a full, spatial picture of these abnormalities has not been well documented in a quantitative fashion in hip OA. Here, micro-computed tomography (44.8 μm/voxel) and data-driven computational anatomy were used to generate 3-D maps of the distribution of bone density and microstructure in human femoral neck samples with early (6F/4M, mean age = 51.3 years), moderate (14F/8M, mean age = 60 years), and severe (16F/6M, mean age = 63.3 years) radiographic OA. With increasing severity of radiographic OA, there was decreased cortical bone mineral density (BMD) (p=0.003), increased cortical thickness (p=0.001), increased cortical porosity (p=0.0028), and increased cortical cross-sectional area (p=0.0012, due to an increase in periosteal radius (p=0.018)), with no differences detected in the total femoral neck or trabecular compartment measures. No OA-related region-specific differences were detected through Statistical Parametric Mapping, but there were trends towards decreased tissue mineral density (TMD) in the inferior femoral neck with increasing OA severity (0.050 < p ≤ 0.091), possibly due to osteophytes. Overall, the lack of differences in cortical TMD among radiographic OA groups indicated that the decrease in cortical BMD with increasing OA severity was largely due to the increased cortical porosity rather than decreased tissue mineralization. As porosity is inversely associated with stiffness and strength in cortical bone, increased porosity may offset the effect that increased cortical cross-sectional area would be expected to have on reducing stresses within the femoral neck. The use of high-resolution imaging and quantitative spatial assessment in this study provide insight into the heterogeneous and multi-faceted changes in density and microstructure in hip OA, which have implications for OA progression and fracture risk.

Treatment Failure in Femoral Neck Fractures in Adults Less Than 50 Years of Age: Analysis of 492 Patients Repaired at 26 North American Trauma Centers

OBJECTIVES

To assess the operative results of femoral neck fractures (FNFs) in young adults in a large multicenter series, specifically focusing on risk factors for treatment failure.

DESIGN

Large multicenter retrospective cohort series.

SETTING

Twenty-six North American Level 1 trauma centers.

PATIENTS

Skeletally mature patients younger than 50 years with displaced and nondisplaced FNFs treated between 2005 and 2017.

INTERVENTION

Operative repair of FNF.

MAIN OUTCOME MEASUREMENTS

The main outcome measure is treatment failure: nonunion and/or failed fixation, osteonecrosis, malunion, and need for subsequent major reconstructive surgery (arthroplasty or proximal femoral osteotomy). Logistic regression models were conducted to examine factors associated with treatment failure.

RESULTS

Of 492 patients with FNFs studied, a major complication and/or subsequent major reconstructive surgery occurred in 45% (52% of 377 displaced fractures and 21% of 115 nondisplaced fractures). Overall, 23% of patients had nonunion/failure of fixation, 12% osteonecrosis type 2b or worse, 15% malunion (>10 mm), and 32% required major reconstructive surgery. Odds of failure were increased with fair-to-poor reduction [odds ratio (OR) = 5.29, 95% confidence interval (CI) = 2.41-13.31], chronic alcohol misuse (OR = 3.08, 95% CI = 1.59-6.38), comminution (OR = 2.63, 95% CI = 1.69-4.13), multiple screw constructs (vs. fixed-angle devices, OR = 1.95, 95% CI = 1.30-2.95), metabolic bone disease (OR = 1.77, 95% CI = 1.17-2.67), and increasing age (OR = 1.03, 95% CI = 1.01-1.06). Women (OR = 0.57, 95% CI = 0.37-0.88), Pauwels angle ≤50 degrees (type 1 or 2; OR = 0.64, 95% CI = 0.41-0.98), or associated femoral shaft fracture (OR = 0.19, 95% CI = 0.10-0.33) had lower odds of failure.

CONCLUSIONS

FNFs in adults <50 years old remain a difficult clinical and surgical problem, with 45% of patients experiencing major complications and 32% undergoing subsequent major reconstructive surgery. Risk factors for complications after treatment of displaced FNFs were numerous.

LEVEL OF EVIDENCE

Therapeutic Level III. See Instructions for Authors for a complete description of levels of evidence.

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).

Development of an open-source measurement system to assess the areal bone mineral density of the proximal femur from clinical CT images

Commercial software is generally needed to measure the areal bone mineral density (aBMD) of the proximal femur from clinical computed tomography (CT) images. This study developed and verified an open-source reproducible system to quantify CT-aBMD to screen osteoporosis using clinical CT images.

PURPOSE

For existing CT images acquired for various reasons other than osteoporosis, it might be beneficial to estimate areal BMD as assessed by dual-energy X-ray absorptiometry (DXA-based BMD) to ascertain the bone status based on DXA. In this study, we aimed to (1) develop an open-source reproducible measurement system to quantify DXA-based BMD from CT images and (2) validate its accuracy.

METHODS

This study analyzed 75 pairs of hip CT and DXA images of women that were acquired for the preoperative assessment of total hip arthroplasty. From the CT images, the femur and a calibration phantom were automatically segmented using pre-trained codes/models available at https://github.com/keisuke-uemura . The proximal femoral region was isolated by manually selected landmarks and was projected onto the coronal plane to measure the areal density (CT-aHU). The calibration phantom was employed to convert the CT-aHU into CT-aBMD. Each parameter was correlated with DXA-based BMD, and the residual errors of CT images to estimate the T-scores in DXA were calculated using the standard error of estimate (SEE).

RESULTS

The correlation coefficients of DXA-based BMD with CT-aHU and CT-aBMD were 0.947 and 0.950, respectively (both p < 0.001). The SEE for quantifying the T-scores in DXA were 0.51 and 0.50 for CT-aHU and CT-aBMD, respectively.

CONCLUSION

With the method developed herein, CT permits estimation of the DXA-based BMD of the proximal femur within the standard DXA total hip region of interest with an SEE of 0.5 in T-scores. The radiation dose for CT acquisition needs consideration; therefore, our data do not provide a rationale for performing CT for screening osteoporosis. However, on CT images already acquired for clinical indications other than osteoporosis, researchers may use this open-source system to investigate osteoporosis status through the estimated DXA-based BMD of the proximal femur.

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.

Finite Element Assessment of Bone Fragility from Clinical Images

PURPOSE OF REVIEW

We re-evaluated clinical applications of image-to-FE models to understand if clinical advantages are already evident, which proposals are promising, and which questions are still open.

RECENT FINDINGS

CT-to-FE is useful in longitudinal treatment evaluation and groups discrimination. In metastatic lesions, CT-to-FE strength alone accurately predicts impending femoral fractures. In osteoporosis, strength from CT-to-FE or DXA-to-FE predicts incident fractures similarly to DXA-aBMD. Coupling loads and strength (possibly in dynamic models) may improve prediction. One promising MRI-to-FE workflow may now be tested on clinical data. Evidence of artificial intelligence usefulness is appearing. CT-to-FE is already clinical in opportunistic CT screening for osteoporosis, and risk of metastasis-related impending fractures. Short-term keys to improve image-to-FE in osteoporosis may be coupling FE with fall risk estimates, pool FE results with other parameters through robust artificial intelligence approaches, and increase reproducibility and cross-validation of models. Modeling bone modifications over time and bone fracture mechanics are still open issues.

Statistical Shape and Appearance Models: Development Towards Improved Osteoporosis Care

PURPOSE OF REVIEW

Statistical models of shape and appearance have increased their popularity since the 1990s and are today highly prevalent in the field of medical image analysis. In this article, we review the recent literature about how statistical models have been applied in the context of osteoporosis and fracture risk estimation.

RECENT FINDINGS

Recent developments have increased their ability to accurately segment bones, as well as to perform 3D reconstruction and classify bone anatomies, all features of high interest in the field of osteoporosis and fragility fractures diagnosis, prevention, and treatment. An increasing number of studies used statistical models to estimate fracture risk in retrospective case-control cohorts, which is a promising step towards future clinical application. All the reviewed application areas made considerable steps forward in the past 5-6 years. Heterogeneities in validation hinder a thorough comparison between the different methods and represent one of the future challenges to be addressed to reach clinical implementation.

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.

Bone mineral density modeling via random field: Normality, stationarity, sex and age dependence

BACKGROUND AND OBJECTIVE

Capturing the population variability of bone properties is of paramount importance to biomedical engineering. The aim of the present paper is to describe variability and correlations in bone mineral density with a spatial random field inferred from routine computed tomography data.

METHODS

Random fields were simulated by transforming pairwise uncorrelated Gaussian random variables into correlated variables through the spectral decomposition of an age-detrended correlation matrix. The validity of the random field model was demonstrated in the spatiotemporal analysis of bone mineral density. The similarity between the computed tomography samples and those generated via random fields was analyzed with the energy distance metric.

RESULTS

The random field of bone mineral density was found to be approximately Gaussian/slightly left-skewed/strongly right-skewed at various locations. However, average bone density could be simulated well with the proposed Gaussian random field for which the energy distance, i.e., a measure that quantifies discrepancies between two distribution functions, is convergent with respect to the number of correlation eigenpairs.

CONCLUSIONS

The proposed random field model allows the enhancement of computational biomechanical models with variability in bone mineral density, which could increase the usability of the model and provides a step forward in in-silico medicine.

Bone Analysis Using Trabecular Bone Score and Dual-Energy X-Ray Absorptiometry-Based 3-Dimensional Modeling in Postmenopausal Women With Primary Hyperparathyroidism

OBJECTIVE

Predominance of bone loss in cortical sites with relative preservation of trabecular bone, even in postmenopausal women, has been described in primary hyperparathyroidism (PHPT). The aim of this study was to evaluate bone microarchitectural differences using dual-energy x-ray absorptiometry (DXA), trabecular bone score (TBS), and DXA-based 3-dimensional (3D) modeling (3D-DXA) between postmenopausal women diagnosed with PHPT (PM-PHPT) and healthy postmenopausal controls.

METHODS

This retrospective study included 44 women with PM-PHPT (9 of whom had fractures) and 48 healthy women matched by age, body mass index, and years since menopause treated at Hospital Universitario Fundación Jiménez Díaz between 2008 and 2017. The bone mineral density (BMD) of the lumbar spine (LS), femoral neck, total hip (TH), and 1/3 radius was assessed using DXA, and trabecular volumetric BMD (vBMD), cortical vBMD, integral vBMD, cortical thickness, and cortical surface BMD at TH were assessed using a 3D-DXA software and TBS at LS.

RESULTS

The mean adjusted BMD values at LS, the femoral neck, and TH; TBS at LS; and TH 3D-DXA parameters (trabecular vBMD, integral vBMD, cortical thickness, and cortical surface BMD) were significantly reduced in women with PM-PHPT compared with those in the controls. However, differences in mean cortical vBMD were not statistically significant (P = .078). There were no significant differences in mean BMD, TBS, or the 3D-DXA parameters between patients with fractures and those without fractures. The 25-hydroxyvitamin D level appeared to be associated with TBS but not with DXA and 3D-DXA measurements.

CONCLUSION

PM-PHPT has significant involvement of the trabecular and cortical compartments of the bone, as determined by DXA, TBS, and 3D-DXA.

A new approach to analyzing regenerated bone quality in the mouse digit amputation model using semi-automatic processing of microCT data

Bone regeneration is a critical area of research impacting treatment of diseases such as osteoporosis, age-related decline, and orthopaedic implants. A crucial question in bone regeneration is that of bone architectural quality, or how "good" is the regenerated bone tissue structurally? Current methods address typical long bone architecture, however there exists a need for improved ability to quantify structurally relevant parameters of bone in non-standard bone shapes. Here we present a new analysis approach based on open-source semi-automatic methods combining image processing, solid modeling, and numerical calculations to analyze bone tissue at a more granular level using μCT image data from a mouse digit model of bone regeneration. Examining interior architecture, growth patterning, spatial mineral content, and mineral density distribution, these methods are then applied to two types of 6-month old mouse digits - 1) those prior to amputation injury (unamputated) and 2) those 42 days after amputation when bone has regenerated. Results show regenerated digits exhibit increased inner void fraction, decreased patterning, different patterns of spatial mineral distribution, and increased mineral density values when compared to unamputated bone. Our approach demonstrates the utility of this new analysis technique in assessment of non-standard bone models, such as the regenerated bone of the digit, and aims to bring a deeper level of analysis with an open-source, integrative platform to the greater bone community.

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.

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.

Cortical thickness of the tibial diaphysis reveals age- and sex-related characteristics between non-obese healthy young and elderly subjects depending on the tibial regions

PURPOSE

This study aimed to evaluate the age- and sex-related characteristics in cortical thickness of the tibial diaphysis between non-obese healthy young and elderly subjects as reference data.

METHODS

The study investigated 31 young subjects (12 men and 19 women; mean age, 25 ± 8 years) and 54 elderly subjects (29 men and 25 women; mean age, 70 ± 6 years). Three-dimensional estimated cortical thickness of the tibial diaphysis was automatically calculated for 5000-9000 measurement points using the high-resolution cortical thickness measurement from clinical computed tomography data. In 12 assessment regions created by combining three heights (proximal, central, and distal diaphysis) and four areas of the axial plane at 90° (medial, anterior, lateral, and posterior areas) in the tibial coordinate system, the standardized thickness was assessed using the tibial length.

RESULTS

As structural characteristics, there were no differences in the medial and lateral thicknesses, while the anterior thickness was greater than the posterior thickness in all groups. The sex-related difference was not shown. As an age-related difference, elderly subjects showed greater or lesser cortical thickness than the young subjects, depending on the regions of the tibia.

CONCLUSIONS

Cortical thickness was different depending on sex, age, and regions in the tibia. The results of this study are of clinical relevance as reference points to clarify the causes of various pathological conditions for diseases.

LEVEL OF EVIDENCE

Level 3.

Three-Dimensional Surface-Based Analysis of Cartilage MRI Data in Knee Osteoarthritis: Validation and Initial Clinical Application

BACKGROUND

Traditional quantitative analysis of cartilage with MRI averages measurements (eg, thickness) across regions-of-interest (ROIs) which may reduce responsiveness.

PURPOSE

To validate and describe clinical application of a semiautomated surface-based method for analyzing cartilage relaxation times ("composition") and morphology on MRI, 3D cartilage surface mapping (3D-CaSM).

STUDY TYPE

Validation study in cadaveric knees and prospective observational (cohort) study in human participants.

POPULATION

Four cadaveric knees and 14 participants aged 40-60 with mild-moderate knee osteoarthritis (OA) and 6 age-matched healthy volunteers, imaged at baseline, 1, and 6 months.

FIELD STRENGTH/SEQUENCE

3D spoiled gradient echo, T1 rho/T2 magnetization-prepared 3D fast spin echo for mapping of T1 rho/T2 relaxation times and delayed gadolinium enhanced MRI of cartilage (dGEMRIC) using variable flip angle T1 relaxation time mapping at 3T.

ASSESSMENT

3D-CaSM was validated against high-resolution peripheral quantitative computed tomography (HRpQCT) in cadaveric knees, with comparison to expert manual segmentation. The clinical study assessed test-retest repeatability and sensitivity to change over 6 months for cartilage thickness and relaxation times.

STATISTICAL TESTS

Bland-Altman analysis was performed for the validation study and evaluation of test-retest repeatability. Six-month changes were assessed via calculation of the percentage of each cartilage surface affected by areas of significant change (%SC), defined using thresholds based on area and smallest detectable difference (SDD).

RESULTS

Bias and precision (0.06 ± 0.25 mm) of 3D-CaSM against reference HRpQCT data were comparable to expert manual segmentation (-0.13 ± 0.26 mm). 3D-CaSM demonstrated significant (>SDD) 6-month changes in cartilage thickness and relaxation times in both OA participants and healthy controls. The parameter demonstrating the greatest 6-month change was T2 relaxation time (OA median %SC [IQR] = 8.8% [5.5 to 12.6]).

DATA CONCLUSION

This study demonstrates the construct validity and potential clinical utility of 3D-CaSM, which may offer advantages to conventional ROI-based methods.

LEVEL OF EVIDENCE

2.

TECHNICAL EFFICACY STAGE

2. J. Magn. Reson. Imaging 2020;52:1139-1151.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Age- and sex-related characteristics in cortical thickness of femoral diaphysis for young and elderly subjects

INTRODUCTION

Cortical thickness of the femoral diaphysis is assumed to be a preferred parameter in the assessment of the structural adaptation by mechanical use and biological factors. This study aimed to investigate the age- and sex-specific characteristics in cortical thickness of the femoral diaphysis between young and elderly non-obese people.

MATERIALS AND METHODS

This study investigated 34 young subjects (21 men and 13 women; mean age: 27 ± 8 years) and 52 elderly subjects (29 men and 23 women; mean age: 70 ± 6 years). Three-dimensional (3D) cortical thickness of the femoral diaphysis was automatically calculated for 5000-8000 measurement points using the high-resolution cortical thickness measurement from clinical CT data. In 12 assessment regions created by combining three heights (proximal, central, and distal diaphysis) and four areas of the axial plane at 90° (medial, anterior, lateral, and posterior areas) in the femoral coordinate system, the standardized thickness was assessed using the femoral length.

RESULTS

As per the trends, (1) there were no differences in medial and lateral thicknesses, while the posterior thickness was greater than the anterior thickness, (2) the thickness in men was higher than that in women, and (3) the thickness in young subjects was higher than that in elderly subjects.

CONCLUSIONS

The results of this study are of clinical relevance as reference points to clarify the causes of various pathological conditions for diseases of the lower extremities.

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.

Clinical Contrast-Enhanced Computed Tomography With Semi-Automatic Segmentation Provides Feasible Input for Computational Models of the Knee Joint

Computational models can provide information on joint function and risk of tissue failure related to progression of osteoarthritis (OA). Currently, the joint geometries utilized in modeling are primarily obtained via manual segmentation, which is time-consuming and hence impractical for direct clinical application. The aim of this study was to evaluate the applicability of a previously developed semi-automatic method for segmenting tibial and femoral cartilage to serve as input geometry for finite element (FE) models. Knee joints from seven volunteers were first imaged using a clinical computed tomography (CT) with contrast enhancement and then segmented with semi-automatic and manual methods. In both segmentations, knee joint models with fibril-reinforced poroviscoelastic (FRPVE) properties were generated and the mechanical responses of articular cartilage were computed during physiologically relevant loading. The mean differences in the absolute values of maximum principal stress, maximum principal strain, and fibril strain between the models generated from semi-automatic and manual segmentations were <1 MPa, <0.72% and <0.40%, respectively. Furthermore, contact areas, contact forces, average pore pressures, and average maximum principal strains were not statistically different between the models (p >0.05). This semi-automatic method speeded up the segmentation process by over 90% and there were only negligible differences in the results provided by the models utilizing either manual or semi-automatic segmentations. Thus, the presented CT imaging-based segmentation method represents a novel tool for application in FE modeling in the clinic when a physician needs to evaluate knee joint function.

Quantitative 3D imaging parameters improve prediction of hip osteoarthritis outcome

Osteoarthritis is an increasingly important health problem for which the main treatment remains joint replacement. Therapy developments have been hampered by a lack of biomarkers that can reliably predict disease, while 2D radiographs interpreted by human observers are still the gold standard for clinical trial imaging assessment. We propose a 3D approach using computed tomography—a fast, readily available clinical technique—that can be applied in the assessment of osteoarthritis using a new quantitative 3D analysis technique called joint space mapping (JSM). We demonstrate the application of JSM at the hip in 263 healthy older adults from the AGES-Reykjavík cohort, examining relationships between 3D joint space width, 3D joint shape, and future joint replacement. Using JSM, statistical shape modelling, and statistical parametric mapping, we show an 18% improvement in prediction of joint replacement using 3D metrics combined with radiographic Kellgren & Lawrence grade (AUC 0.86) over the existing 2D FDA-approved gold standard of minimum 2D joint space width (AUC 0.73). We also show that assessment of joint asymmetry can reveal significant differences between individuals destined for joint replacement versus controls at regions of the joint that are not captured by radiographs. This technique is immediately implementable with standard imaging technologies.

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.

Automated segmentation of cortical and trabecular bone to generate finite element models for femoral bone mechanics

Finite element (FE) models based on quantitative computed tomography (CT) images are better predictors of bone strength than conventional areal bone mineral density measurements. However, FE models require manual segmentation of the femur, which is not clinically applicable. This study developed a method for automated FE analyses from clinical CT images. Clinical in-vivo CT images of 13 elderly female subjects were collected to evaluate the method. Secondly, proximal cadaver femurs were harvested and imaged with clinical CT (N = 17). Of these femurs, 14 were imaged with µCT and three had earlier been tested experimentally in stance-loading, while collecting surface deformations with digital image correlation. Femurs were segmented from clinical CT images using an automated method, based on the segmentation tool Stradwin. The method automatically distinguishes trabecular and cortical bone, corrects partial volume effect and generates input for FE analysis. The manual and automatic segmentations agreed within about one voxel for in-vivo subjects (0.99 ± 0.23 mm) and cadaver femurs (0.21 ± 0.07 mm). The strains from the FE predictions closely matched with the experimentally measured strains (R² = 0.89). The method can automatically generate meshes suitable for FE analysis. The method may bring us one step closer to enable clinical usage of patient-specific FE analyses.

Regional maps of rib cortical bone thickness and cross-sectional geometry

Here we present detailed regional bone thickness and cross-sectional measurements from full adult ribs using high resolution CT scans processed with a cortical bone mapping technique. Sixth ribs from 33 subjects ranging from 24 to 99 years of age were used to produce average cortical bone thickness maps and to provide average ± 1SD corridors for expected cross-section properties (cross-sectional areas and inertial moments) as a function of rib length. Results obtained from CT data were validated at specific rib locations using direct measurements from cut sections. Individual thickness measurements from CT had an accuracy (mean error) and precision (SD error) of -0.013 ± 0.167 mm (R² coefficient of determination of 0.84). CT-based measurement errors for rib cross-sectional geometry were -0.1 ± 13.1% (cortical bone cross-sectional area) and 4.7 ± 1.8% (total cross-sectional area). Rib cortical bone thickness maps show the expected regional variation across a typical rib's surface. The local mid-rib maxima in cortical thickness along the pleural rib aspect ranged from range 0.9 to 2.6 mm across the study population with an average map maximum of 1.4 mm. Along the cutaneous aspect, rib cortical bone thickness ranged from 0.7 to 1.9 mm with an average map thickness of 0.9 mm. Average cross-sectional properties show a steady reduction in total cortical bone area from 10% along the rib's length through to the sternal end, whereas overall cross-sectional area remains relatively constant along the majority of the rib's length before rising steeply towards the sternal end. On average, male ribs contained more cortical bone within a given cross-section than was seen for female ribs. Importantly, however, this difference was driven by male ribs having larger overall cross-sectional areas, rather than by sex differences in the bone thickness observed at specific local cortex sites. The cortical bone thickness results here can be used directly to improve the accuracy of current human body and rib models. Furthermore, the measurement corridors obtained from adult subjects across a wide age range can be used to validate future measurements from more widely available image sources such as clinical CT where gold standard reference measures (e.g. such as direct measurements obtained from cut sections) are otherwise unobtainable.

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.

Trabecular architecture of the great ape and human femoral head

Studies of femoral trabecular structure have shown that the orientation and volume of bone are associated with variation in loading and could be informative about individual joint positioning during locomotion. In this study, we analyse for the first time trabecular bone patterns throughout the femoral head using a whole-epiphysis approach to investigate how potential trabecular variation in humans and great apes relates to differences in locomotor modes. Trabecular architecture was analysed using microCT scans of Pan troglodytes (n = 20), Gorilla gorilla (n = 14), Pongo sp. (n = 5) and Homo sapiens (n = 12) in medtool 4.1. Our results revealed differences in bone volume fraction (BV/TV) distribution patterns, as well as overall trabecular parameters of the femoral head between great apes and humans. Pan and Gorilla showed two regions of high BV/TV in the femoral head, consistent with hip posture and loading during two discrete locomotor modes: knuckle-walking and climbing. Most Pongo specimens also displayed two regions of high BV/TV, but these regions were less discrete and there was more variability across the sample. In contrast, Homo showed only one main region of high BV/TV in the femoral head and had the lowest BV/TV, as well as the most anisotropic trabeculae. The Homo trabecular structure is consistent with stereotypical loading with a more extended hip compared with great apes, which is characteristic of modern human bipedalism. Our results suggest that holistic evaluations of femoral head trabecular architecture can reveal previously undetected patterns linked to locomotor behaviour in extant apes and can provide further insight into hip joint loading in fossil hominins and other primates.

Discrimination of osteoporosis-related vertebral fractures by DXA-derived 3D measurements: a retrospective case-control study

A retrospective case-control study assessing the association of DXA-derived 3D measurements with osteoporosis-related vertebral fractures was performed. Trabecular volumetric bone mineral density was the measurement that best discriminates between fracture and control groups.

INTRODUCTION

The aim of the present study was to evaluate the association of DXA-derived 3D measurements at the lumbar spine with osteoporosis-related vertebral fractures.

METHODS

We retrospectively analyzed a database of 74 postmenopausal women: 37 subjects with incident vertebral fractures and 37 age-matched controls without any type of fracture. DXA scans at the lumbar spine were acquired at baseline (i.e., before the fracture event for subjects in the fracture group), and areal bone mineral density (aBMD) was measured. DXA-derived 3D measurements, such as volumetric BMD (vBMD), were assessed using a DXA-based 3D modeling software (3D-SHAPER). vBMD was computed at the trabecular, cortical, and integral bone. Cortical thickness and cortical surface BMD were also measured. Differences in DXA-derived measurements between fracture and control groups were evaluated using unpaired t test. Odds ratio (OR) and area under the receiver operating curve (AUC) were also computed. Subgroup analyses according to fractured vertebra were performed.

RESULTS

aBMD of fracture group was 9.3% lower compared with control group (p < 0.01); a higher difference was found for trabecular vBMD in the vertebral body (- 16.1%, p < 0.001). Trabecular vBMD was the measurement that best discriminates between fracture and control groups, with an AUC of 0.733, against 0.682 for aBMD. Overall, similar findings were observed within the subgroup analyses. The L1 vertebral fractures subgroup had the highest AUC at trabecular vBMD (0.827), against aBMD (0.758).

CONCLUSION

This study showed the ability of cortical and trabecular measurements from DXA-derived 3D models to discriminate between fracture and control groups. Large cohorts need to be analyzed to determine if these measurements could improve fracture risk prediction in clinical practice.

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.

A cross-sectional study on the age-related cortical and trabecular bone changes at the femoral head in elderly female hip fracture patients

Bone is in a continuous state of remodeling whereby old bone is absorbed and new bone is formed in its place. During this process, new formations reinforce the bone in the direction of the dominant stress trajectories through a functional adaptation. In normal aging, the balance between bone resorption and formation can be shifted. How this affects the functional adaptation remains to be investigated. Furthermore, how or whether the bone continues to change beyond the age of 85 is not yet studied in detail. In this study we examined the age-related changes in the cortical and trabecular bone in old age, and assessed whether we can find evidence of the presence of functional adaptation. We measured cortical and trabecular parameters from micro-computed tomography scans of the femoral head extracted from hip fracture patients between the age of 70 and 93 years. A significant decrease in global trabecular bone mineral density (38.1%) and cortical thickness (13.0%) was seen from the 9th to the 10th decade of life. The degree of anisotropy was maintained globally as well as locally in both high and low stress regions. The local trabecular bone mineral density decreased in both high stress and low stress regions between the 9th and 10th decade of life with similar trends. This suggests that the role of functional adaptation in maintaining the bone structural integrity in old age may be limited. This study highlights the need for a controlled clinical trial examining the cause of the continued bone degradation throughout old age.