Individual trabecula segmentation (ITS)-based morphological analysis of microscale images of human tibial trabecular bone at limited spatial resolution

A graph model to describe the network connectivity of trabecular plates and rods

Introduction: The trabecular network is perceived as a collection of interconnected plate- (P) and rod-like (R) elements. Previous research has highlighted how these elements and their connectivity influence the mechanical properties of bone, yet further work is required to elucidate better the deeply interconnected nature of the trabecular network with distinct element formations conducting forces per their mechanical boundary conditions. Within this network, forces act through elements: a rod or plate with force applied to one end will transmit this force to a component connected to the other end, defining the boundary conditions for the loading of each element. To that end, this study has two aims: First, to investigate the connectivity of individually segmented elements of trabecular bone with respect to their local boundary conditions as defined by the surrounding trabecular network and linking them directly to the bone's overall mechanical response during loading using a mathematical graph model of the plate and rod (PR) Network. Second, we use this model to quantify side artifacts, a known artifact when testing an excised specimen of trabecular bone, where vertical trabeculae lose their load-bearing capacity due to a loss of connectivity, ultimately resulting in a change of the trabecular network topology. Resuts: Connected elements derived from our model predicted apparent elastic modulus by fitting a linear regression (R 2 = 0.81). In comparison, prediction using conventional bone volume fraction results in a lower accuracy (R 2 = 0.72), demonstrating the ability of the PR Network to estimate compressive elastic modulus independent of specimen size or loading boundary condition. Discussion: PR Network models are a novel approach to describing connectivity within the trabecular network and incorporating mechanical boundary conditions within the morphological analysis, thus enabling the study of intrinsic material properties of trabecular bone. Ultimately, PR Network models may be an early predictor or provide further insights into osteo-degenerative diseases.

Individual trabecula segmentation validation in first- and second-generation high-resolution peripheral computed tomography compared to micro-computed tomography in the distal radius and tibia

High-resolution peripheral quantitative computed tomography (HR-pQCT) has been used for in vivo 3D visualization of trabecular microstructure. Second-generation HR-pQCT (HR-pQCT II) has been shown to have good agreement with first generation HR-pQCT (HR-pQCT I). Advanced Individual Trabecula Segmentation (ITS) decomposes the trabecula network into individual plates and rods. ITS based on HR-pQCT I showed a strong correlation to ITS based on micro-computed tomography (μCT) and identified trabecular changes in metabolic bone diseases. ITS based on HR-pQCT II has new potential because of the enhanced resolution but has yet to be validated. The objective of this study was to assess the agreement between ITS based on HR-pQCT I, HR-pQCT II, and μCT to assess the capability of ITS on HR-pQCT images as a tool for studying bone structure. Freshly frozen tibia and radius bones were scanned in the distal region using HR-pQCT I at 82 μm, HR-pQCT II at 60.7 μm, and μCT at 37 μm. Images were registered, binarized, and ITS analysis was performed. Bone volume fraction (pBV/TV, rBV/TV), number density (pTb.N, rTb.N), thickness (pTb.Th, rTb.Th), and plate-to-rod (PR) ratio (pBV/rBV) of trabecular plates and rods were obtained. Paired Student's t-tests with post hoc Bonferroni analysis were used to examine the differences. Linear regression was used to determine the correlation coefficient. The HR-pQCT I parameters were different from the μCT measurements. The HR-pQCT II parameters were different from the μCT measurements except for rTb.N, and the HR-pQCT I parameters were different from the HR-pQCT II measurements except for pTb.Th. The strong correlation between HR-pQCT II and μCT microstructural analysis (R2 = 0.55-0.94) suggests that HR-pQCT II can be used to assess changes in plate and rod microstructure and that values from HR-pQCT I can be corrected.

Teriparatide Followed by Denosumab in Premenopausal Idiopathic Osteoporosis: Bone Microstructure and Strength by HR-pQCT

Premenopausal women with idiopathic osteoporosis (PreMenIOP) have marked deficits in skeletal microstructure. We have reported that sequential treatment with teriparatide and denosumab improves central skeletal bone mineral density (BMD) by dual-energy X-ray absorptiometry and central QCT in PreMenIOP. We conducted preplanned analyses of high-resolution peripheral quantitative computed tomography (HR-pQCT) scans from teriparatide and denosumab extension studies to measure effects on volumetric BMD (vBMD), microarchitecture, and estimated strength at the distal radius and tibia. Of 41 women enrolled in the parent teriparatide study (20 mcg daily), 34 enrolled in the HR-pQCT study. HR-pQCT participants initially received teriparatide (N = 24) or placebo (N = 10) for 6 months; all then received teriparatide for 24 months. After teriparatide, 26 enrolled in the phase 2B denosumab extension (60 mg q6M) for 24 months. Primary outcomes were percentage change in vBMD, microstructure, and stiffness after teriparatide and after denosumab. Changes after sequential teriparatide and denosumab were secondary outcomes. After teriparatide, significant improvements were seen in tibial trabecular number (3.3%, p = 0.01), cortical area and thickness (both 2.7%, p < 0.001), and radial trabecular microarchitecture (number: 6.8%, thickness: 2.2%, separation: -5.1%, all p < 0.02). Despite increases in cortical porosity and decreases in cortical density, whole-bone stiffness and failure load increased at both sites. After denosumab, increases in total (3.5%, p < 0.001 and 3.3%, p = 0.02) and cortical vBMD (1.7% and 3.2%; both p < 0.01), and failure load (1.1% and 3.6%; both p < 0.05) were seen at tibia and radius, respectively. Trabecular density (3.5%, p < 0.001) and number (2.4%, p = 0.03) increased at the tibia, while thickness (3.0%, p = 0.02) increased at the radius. After 48 months of sequential treatment, significant increases in total vBMD (tibia: p < 0.001; radius: p = 0.01), trabecular microstructure (p < 0.05), cortical thickness (tibia: p < 0.001; radius: p = 0.02), and whole bone strength (p < 0.02) were seen at both sites. Significant increases in total vBMD and bone strength parameters after sequential treatment with teriparatide followed by denosumab support the use of this regimen in PreMenIOP. © 2022 American Society for Bone and Mineral Research (ASBMR).

Ankle flexor torque, size and density are differential determinants of distal tibia trabecular plate-rod morphometry and bone strength: The Ankle Quality Study

HYPOTHESIS

Greater peak torque and higher myotendinous density at the ankle are associated with a more plate-like architecture at the distal tibia.

METHODS

In this cross-sectional study, women and men ≥ 50 years old with no metal implants, reconstructive surgery, muscular dystrophies, or tendinopathies in any leg were recruited by convenience. Isometric ankle dorsi-plantar flexion and inversion-eversion peak torques were measured using dynamometry. HR-pQCT distal tibia scans were completed. Both assessments were completed on the same day on the non-dominant leg. Integral and trabecular vBMD were derived from standard analyses, failure load (FL) was obtained from finite element analysis, plate-specific parameters were computed from individual trabecula segmentation (ITS) analysis, myotendinous density (MyD) and volume fraction (MyV/TV) were computed from soft tissue analysis. pQCT scans of the 66 % mid-leg were performed (500 μm at 15 mm/s) to obtain muscle density (MD) and muscle cross-sectional area (MCSA).

STATISTICAL ANALYSIS

General linear models estimated how ankle muscle group torque and muscle size and density differentially related, both separately and together, to whole-bone properties (integral vBMD, FL) and trabecular morphometry (ITS plate parameters). Models were adjusted for age, sex, BMI, use of glucocorticoids, current osteoarthritis, and participation in moderate to vigorous recreational or sport activities.

RESULTS

Among 105 participants (77 % female, mean age: 63 (10) years, BMI: 25.8 (5.4) kg/m², 25 % with OA, 17 % fracture history, 42 % falls history), all torque measures, particularly ankle dorsiflexion and eversion, were correlates of plate-plate/rod junction density and failure load. However, muscle size and density measures were further associated with vBMD. The effect of greater ankle flexor-extensor torque on more connected bone was stronger when MyD was higher (interaction p < 0.001).

CONCLUSION

Strength of muscles around the ankle are correlates of plate-like trabeculae at the distal tibia, while leaner muscle and myotendinous tissues facilitates better quality bone for stronger ankle muscle torque.

The effect of short-term high-caloric feeding and fasting on bone microarchitecture

BACKGROUND

States of chronic overnutrition and undernutrition are both associated with impaired bone health and increased fracture risk but there are no data on bone microarchitecture following short-term controlled nutritional challenges.

OBJECTIVE

The purpose of our study was to evaluate the impact of short-term high-caloric feeding and fasting on bone microarchitecture. We hypothesized that both high-caloric feeding and fasting would have negative effects on microarchitecture.

MATERIALS AND METHODS

We recruited 23 adult healthy subjects (13 males, 10 females, mean age 33.2 ± 1.4 years, mean BMI 26.0 ± 1.5 kg/m²). Subjects underwent an in-patient 10-day high-caloric visit (caloric intake with goal to achieve 7% weight gain), after which they went home to resume a normal diet for 13-18 days (stabilization period), and were then readmitted for a 10-day in-patient fasting stay (no caloric intake). All subjects underwent HRpQCT (XtremeCT, Scanco Medical AG, Brüttisellen, Switzerland) of the distal tibia and distal radius after each visit to assess volumetric bone mineral density (vBMD), trabecular and cortical microarchitecture, and strength estimates. The Wilcoxon signed rank test was used to perform within group comparisons.

RESULTS

During the high-caloric period, there was a mean increase in weight by 6.3 + 1.7% (p < 0.0001). There were no significant changes in bone parameters in the distal tibia or distal radius (p > 0.05). During the stabilization period there was a significant reduction in weight by -2.7 + 1.9% (p < 0.0001) but no change in bone parameters (p > 0.05). During the fasting period there was a further reduction in weight by -8.8 + 1.2% (p < 0.0001). In the distal tibia, there was a significant increase in total and cortical vBMD, trabecular and cortical parameters as well as strength estimates (p < 0.05). In the distal radius there was an increase in total and trabecular vBMD (p < 0.05), while there were no changes in other microarchitecture parameters or strengths estimates.

CONCLUSION

Short-term fasting after high-caloric feeding improves vBMD, bone microarchitecture and strength estimates of the distal tibia, while short-term high-caloric feeding does not change vBMD or microarchitecture. These results suggest that short-term fasting after high-caloric feeding in healthy individuals improves bone health and that these changes can be detected using HRpQCT in-vivo.

The effect of denosumab and alendronate on trabecular plate and rod microstructure at the distal tibia and radius: A post-hoc HR-pQCT study

BACKGROUND

Age-related trabecular microstructural deterioration and conversion from plate-like trabeculae to rod-like trabeculae occur because of unbalanced rapid remodeling. As denosumab achieves greater remodeling suppression and lower cortical porosity than alendronate, we hypothesized that denosumab might also preserve trabecular plate microstructure, bone stiffness and strength more effectively than alendronate.

METHODS

In this post hoc analysis of a phase 2 study, postmenopausal women randomized to placebo (P, n = 74), denosumab (D, n = 72), or alendronate (A, n = 68). HR-pQCT scans of the distal radius and tibia were performed at baseline and Month-12 (M12). Trabecular compartment was subjected to Individual Trabecula Segmentation while finite element analysis was performed to estimate stiffness and strength. Percent change from baseline at M12 of each parameter was compared between patient groups.

RESULTS

At the distal tibia, in the placebo group, plate surface area (pTb.S, -1.3%) decreased while rod bone volume fraction (rBV/TV, +4.5%) and number (rTb.N, +2.1%) increased. These changes were prevented by denosumab but persisted despite alendronate therapy (pTb.S: -1.7%; rBV/TV: +6.9%; rTb.N: +3.0%). Both treatments improved whole bone stiffness (D: +3.1%; A: +1.8%) and failure load (D: +3.0%; A: +2.2%); improvements using denosumab was significant compared to placebo (stiffness: p = 0.004; failure load: p = 0.003). At the distal radius, denosumab increased total trabecular bone volume fraction (BV/TV, +3.4%) and whole bone failure load (+4.0%), significantly different from placebo (BV/TV: p = 0.044; failure load: p = 0.046). Significantly different effects of either drug on plate and rod microstructure were not detected.

CONCLUSIONS

Denosumab preserved trabecular plate microstructure. Alendronate did not. However, estimated strength did not differ between denosumab and alendronate treated groups.

Characterization of Ultralow Density Cellular Solids: Lessons from 30 years of Bone Biomechanics Research

Advances in additive manufacturing techniques have enabled the development of micro-architectured materials displaying a combination of low-density and lightweight structures with high specific strength and toughness. The mechanical performance of micro-architectured materials can be assessed using standard techniques; however, when studying low- and ultralow density micro-architectured materials, standard characterization techniques can be subject to experimental artifacts. Additionally, quantitative assessment and comparisons of microarchitectures with distinct lattice patterns is not always straightforward. Cancellous bone is a natural, ultralow density (porosity often exceeding 90%), irregular, cellular solid that has been thoroughly characterized in terms of micro-architecture and mechanical performance over the past 30 years. However, most the literature on cancellous bone mechanical properties and micro-structure-function relationships is in the medical literature and is not immediately accessible to materials designers. Here we provide a brief review of state-of-the-art approaches for characterizing the micro-architecture and mechanical performance of ultralow density cancellous bone, including methods of addressing experimental artifacts during mechanical characterization of ultralow density cellular solids, methods of quantifying microarchitecture, and currently understood structure-function relationships.

Effects of Combination Denosumab and High-Dose Teriparatide Administration on Bone Microarchitecture and Estimated Strength: The DATA-HD HR-pQCT Study

In postmenopausal women at high risk of fracture, we previously reported that combined denosumab and high-dose (HD; 40 μg) teriparatide increased spine and hip bone mineral density (BMD) more than combination with standard-dose teriparatide (SD; 20 μg). To assess the effects of these combinations on bone microarchitecture and estimated bone strength, we performed high-resolution peripheral quantitative computed tomography (HR-pQCT) at the distal radius and distal tibia in these women, who were randomized to receive either teriparatide 20 μg (n = 39) or 40 μg (n = 37) during months 0 to 9 overlapped with denosumab 60 mg s.c. given at months 3 and 9, for a 15-month study duration. The 69 women who completed at least one study visit after baseline are included in this analysis. Over 15 months, increases in total BMD were higher in the HD-group than the SD-group at the distal tibia (5.3% versus 3.4%, p = 0.01) with a similar trend at the distal radius (2.6% versus 1.0%, p = 0.06). At 15 months, cortical porosity remained similar to baseline, with absolute differences of -0.1% and -0.7% at the distal tibia and -0.4% and -0.1% at the distal radius in the HD-group and SD-group, respectively; p = NS for all comparisons. Tibial cortical tissue mineral density increased similarly in both treatment groups (1.3% [p < 0.0001 versus baseline] and 1.5% [p < 0.0001 versus baseline] in the HD-group and SD-group, respectively; p = 0.75 for overall group difference). Improvements in trabecular microarchitecture at the distal tibia and estimated strength by micro-finite element analysis at both sites were numerically greater in the HD-group compared with SD-group but not significantly so. Together, these findings suggest that short-term treatment combining denosumab with either high- or standard-dose teriparatide improves HR-pQCT measures of bone density, microstructure, and estimated strength, with greater gains in total bone density observed in the HD-group, which may be of benefit in postmenopausal women with severe osteoporosis. © 2020 American Society for Bone and Mineral Research (ASBMR).

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

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

Differences in subchondral trabecular bone microstructure and finite element analysis-based biomechanical properties between osteoporosis and osteoarthritis

Background/Objective

The microstructure of the subchondral trabecular bone, including the composition and distribution of plates and rods, has an important influence on the disease progression and mechanical properties of osteoarthritis (OA) and osteoporosis (OP). We aimed to determine whether differences in plates and rods influence the variations in the quantities and qualities of the subchondral trabecular bone between OA and OP.

Materials and methods

Thirty-eight femoral head samples [OA, n = 13; OP, n = 17; normal control (NC), n = 8] were collected from male patients undergoing total hip arthroplasty. They were scanned using microcomputed tomography, and subchondral trabecular structures were analysed using individual trabecular segmentation. Micro-finite element analysis (μFEA) was applied to assess the mechanical property of the trabecular bone. Cartilage changes were evaluated by using histological assessment. Analysis of variance was used to compare intergroup differences in structural and mechanical properties and cartilage degradation. Pearson analysis was used to evaluate the relationship between the trabecula microstructure and biomechanical properties.

Results

Compared with the OP and NC group, there was serious cartilage damage in the OA group. With respect to the microstructure results, the OA group had the highest plate and rod trabecular microstructures including number and junction density among the three groups. For the mechanical properties detected via μFEA, the OA group had higher stiffness and failure load than did the OP group. Pearson analysis revealed that compared with OP, OA had a higher number of microstructure parameters (e.g., rod bone volume fraction and rod trabecular number) that were positively correlated with its mechanical property.

Conclusions

Compared with OP, the OA subchondral bone has both increased plate and rod microarchitecture and has more microstructures positively related with its mechanical property. These differences may help explain the variation in mechanical properties between these bone diseases.

The translational potential of this article

Our findings suggested that changes in the plates and rods of the subchondral trabecular bone play a critical role in OA and OP progression and that the improvement of the subchondral trabecular bone may be a promising treatment approach.

Dimorphism in axial and appendicular dimensions, cortical and trabecular microstructure and matrix mineral density in Chinese and Caucasian women

INTRODUCTION

Appendicular fractures are less common in Chinese than Caucasian women. Bone mineral density (BMD) is lower, not higher than in Caucasians because Chinese have smaller appendicular dimensions than Caucasians. However, smaller bones may offset the liability to fracture by being assembled with a more robust microstructure. We hypothesized that Chinese assemble an appendicular skeleton with a thicker, less porous and more mineralized cortex that is less deteriorated in advanced age than in Caucasians.

METHODS

We compared anthropometry in 477 Chinese and 278 Caucasian women and compared bone microstructure using high-resolution peripheral quantitative computed tomography in another cohort of 186 Chinese and 381 Caucasian women aged 18 to 86 years, all living in Melbourne, Australia. Trabecular plate (p) and rod (r) bone volume/total volume (BV/TV) were quantified using individual trabecula segmentation (ITS). Bone strength was estimated using micro-finite element analysis (μFEA).

RESULTS

Premenopausal Chinese were shorter than Caucasian women, mainly due to shorter leg length. Distal radial total cross sectional area (CSA) was 14.8% smaller (p < 0.001). After adjusting for age and total CSA, Chinese had similar cortical and medullary areas but 0.30 SD lower cortical porosity and 0.27 SD higher matrix mineral density (both p < 0.05). Trabecular plate-to-rod ratio was 0.55 SD higher due to a 0.41 SD higher pBV/TV and 0.36 SD lower rBV/TV (p ranging 0.001 to 0.023). Chinese also had 0.36 SD greater whole bone stiffness and 0.36 SD greater failure load than Caucasians (both p < 0.05). After adjusting for age and total CSA, postmenopausal Chinese had 3.3% smaller cortical area, medullary area was 2.1% larger, cortical porosity was no lower, matrix mineral density and pBV/TV were no higher compared with Caucasians at the distal radius. Whole bone stiffness was 0.39 SD lower and failure load was 0.40 SD lower in Chinese (both p < 0.05).

CONCLUSION

Chinese build a more robust skeleton than Caucasians during growth, an advantage not observed in advanced age due to greater bone loss or race-specific secular trends in bone morphology.

Suboptimal bone microarchitecure in adolescent girls with obesity compared to normal-weight controls and girls with anorexia nervosa

BACKGROUND

Despite their higher areal bone mineral density (aBMD), adolescents with obesity (OB) have an increase in fracture risk, particularly of the extremities, compared with normal-weight controls. Whereas bone parameters that increase fracture risk are well characterized in anorexia nervosa (AN), the other end of nutritional spectrum, these data are lacking in adolescents with obesity.

OBJECTIVE

Our objective was to compare bone parameters in adolescent girls across the nutritional spectrum, to determine whether suboptimal bone adaptation to increased body weight may explain the increased fracture risk in OB.

METHODS

We assessed bone endpoints in 153 adolescent girls 14-21 years old: 50 OB, 48 controls and 55 AN. We used (i) DXA to assess aBMD at the lumbar spine, proximal femur and whole body, and body composition, (ii) high resolution peripheral quantitative CT (HRpQCT) to assess bone geometry, microarchitecture and volumetric BMD (vBMD), and (iii) finite element analysis to assess failure load (a strength estimate) at the distal radius and tibia. All aBMD, microarchitecture and FEA analyses were controlled for age and race.

RESULTS

Groups did not differ for age or height. Areal BMD Z-scores at all sites were highest in OB, intermediate in controls and lowest in AN (p < 0.0001). At the radius, cortical area and thickness were higher in OB compared to AN and control groups (p = 0.001) while trabecular area did not differ across groups. Compared to controls, OB had higher cortical porosity (p = 0.003), higher trabecular thickness (p = 0.024), and higher total, cortical and trabecular vBMD and rod BV/TV (p < 0.04). Plate BV/TV did not differ in OB vs. controls, but was higher than in AN (p = 0.001). At the tibia, total, cortical, and trabecular area and cortical thickness were higher in OB vs. controls and AN (p < 0.005). OB also had higher cortical porosity (p < 0.007) and lower trabecular thickness (p < 0.02) than the other two groups. Trabecular number, total and trabecular vBMD, and rod BV/TV were higher in OB vs. controls and AN (p < 0.02), while cortical vBMD and plate BV/TV did not differ in OB vs. the other two groups. Finally, failure load (a strength estimate) was higher in OB at the radius and tibia compared to controls and AN (p < 0.004 for all). However, after adjusting for body weight, failure load was lower in OB vs. controls at both sites (p < 0.05), and lower than in AN at the distal tibia.

CONCLUSION

Not all bone parameters demonstrate appropriate adaptation to higher body weight. Cortical porosity and plate BV/TV at the radius and tibia, and cortical vBMD and trabecular thickness at the tibia are particularly at risk. These effects may contribute to the higher risk for fracture reported in OB vs. controls.

Trabecular microstructure is influenced by race and sex in Black and White young adults

Lower fracture rates in Black men and women compared to their White counterparts are incompletely understood. High-resolution imaging specific to trabecular bone may provide insight. Black participants have enhanced trabecular morphology. These differences may contribute to the lower fracture risk in Black versus White individuals.

INTRODUCTION

Lower fracture rates in Black men and women compared to their White counterparts may be explained by favorable bone microstructure in Black individuals. Individual trabecular segmentation (ITS) analysis, which characterizes the alignment and plate- and rod-like nature of trabecular bone using high-resolution peripheral quantitative computed tomography (HR-pQCT), may provide insight into trabecular differences by race/ethnic origin.

PURPOSE

We determined differences in trabecular bone microarchitecture, connectivity, and alignment according to race/ethnic origin and sex in young adults.

METHODS

We analyzed HR-pQCT scans of 184 adult (24.2 ± 3.4 years) women (n = 51 Black, n = 50 White) and men (n = 34 Black, n = 49 White). We used ANCOVA to compare bone outcomes, and adjusted for age, height, and weight.

RESULTS

Overall, the effect of race on bone outcomes did not differ by sex, and the effect of sex on bone outcomes did not differ by race. After adjusting for covariates, Black participants and men of both races had greater trabecular plate volume fraction, plate thickness, plate number density, plate surface area, and greater axial alignment of trabeculae, leading to higher trabecular bone stiffness compared to White participants and women, respectively (p < 0.05 for all).

CONCLUSION

These findings demonstrate that more favorable bone microarchitecture in Black individuals compared to White individuals and in men compared to women is not unique to the cortical bone compartment. Enhanced plate-like morphology and greater trabecular axial alignment, established in young adulthood, may contribute to the improved bone strength and lower fracture risk in Black versus White individuals and in men compared to women.

Sexual Dimorphism in Cortical and Trabecular Bone Microstructure Appears During Puberty in Chinese Children

Distal forearm fractures during growth are more common in males than females. Because metaphyseal cortical bone is formed by coalescence of trabeculae emerging from the periphery of the growth plate, we hypothesized that the later onset of puberty in males produces a longer delay in trabecular bone formation and coalescence, which leaves a transient phase of high cortical porosity, low matrix mineral density, and high trabecular density relative to females. We quantified the nondominant distal radial microstructure using high-resolution peripheral quantitative computed tomography in 214 healthy Chinese boys and 219 Chinese girls aged between 7 and 17 years living in Hong Kong. Measurements of 110 slices (9.02 mm) were acquired 5 mm proximal to the growth plate of the nondominant distal radius. Porosity was measured using StrAx1.0 (Straxcorp, Melbourne, VIC, Australia) and trabecular plate and rod structure were measured using individual trabecula segmentation (ITS). Mechanical properties were estimated using finite element analysis (FEA). Results were adjusted for age, total bone cross-sectional area (CSA), dietary calcium intake, and physical activity. In boys, total bone CSA was 17.2% to 22.9% larger throughout puberty, cortical/total bone CSA was 5.1% smaller in Tanner stage 2 only, cortical porosity was 9.4% to 17.5% higher, and matrix mineral density was 1.0% to 2.5% lower in Tanner stage 2 to 5, than girls. Boys had higher trabecular rod BV/TV in Tanner stage 3 and 4, but higher trabecular plate BV/TV and plate to rod ratio in Tanner stage 5, than girls. Boys had 17.0% lower apparent modulus than girls in Tanner stage 2. A transient phase of higher porosity due to dissociation between bone mineral accrual and linear growth may contribute to higher distal radial bone fragility in Chinese boys compared to girls. © 2018 American Society for Bone and Mineral Research.

Robust Trabecular Microstructure in Type 2 Diabetes Revealed by Individual Trabecula Segmentation Analysis of HR-pQCT Images

Type 2 diabetes (T2D) patients have an increased fracture risk, which may be partly explained by compromised bone microarchitecture within the cortical bone compartment. Data on trabecular bone parameters in T2D are contradictory. By high-resolution peripheral quantitative computed tomography (HR-pQCT), trabecular microarchitecture is preserved, yet larger trabecular holes are detected in T2D by MRI and DXA-based trabecular bone scores are abnormal. To determine if there are differences in trabecular microstructure, connectivity, and alignment in postmenopausal women with T2D as compared with controls, we performed an individual trabecula segmentation (ITS) analysis on HR-pQCT scans of the distal radius and tibia in 92 women with (n = 42) and without (n = 50) T2D. Unadjusted analyses showed that T2D subjects had greater total trabecular bone volume, trabecular plate volume fraction, plate number density, plate junction density, and axial alignment at the radius and tibia, and increased plate tissue fraction, but decreased rod tissue fraction and rod length at the radius (p < 0.05 for all). After adjustments for clinical covariates, plate number density and plate junction density remained higher at the radius and tibia, whereas total trabecular bone volume was increased and trabecular rod length was decreased at the radius. These differences remained significant after adjustment for hip BMD and trabecular volumetric bone density. Notably, the increased plate-like ITS qualities were seen in those with T2D duration of <10 years, whereas ITS parameters in subjects with T2D duration ≥10 years did not differ from those of control subjects. In conclusion, postmenopausal women with early T2D had a greater plate-like and less rod-like trabecular network. This early advantage in trabecular plate quality does not explain the well-established increased fracture risk in these patients and does not persist in the later stage of T2D. © 2018 American Society for Bone and Mineral Research.

Commonality in the microarchitecture of trabecular bone: A preliminary study

Understanding the relationship between the microstructure and mechanical function of trabecular bone is critical for prediction and prevention of bone fragility fractures. However, a detailed understanding of the structural design of trabecular microarchitecture is still missing. This study hypothesized that there exists a commonality in the underlying probabilistic distributions of microstructural features of trabecular bones, whereas the microstructural differences among individuals are primarily describe by a set of scalar parameters. To test the hypothesis, twenty-three trabecular bone specimens were obtained from two anatomic locations (i.e., femoral neck and vertebral body) and a diverse group of seventeen donors of different age and sex. The number, size, spatial location, and orientation of individual plates and rods in the trabecular bone specimens were determined via volumetric decomposition of 3D μCT images using the Individual Trabecula Segmentation (ITS) technique. Then, m/n bootstrap Kolmogorov-Smirnov tests were performed to compare the normalized distributions of size, orientation, and spatial arrangement of trabecular plates and rods in the specimens. The results showed that 100% of the twenty-three normalized distributions of each microstructural feature were statistically equivalent irrespective of individual differences among the bone specimens, except the distributions of rod spatial arrangement (<100%). On the other hand, nonparametric Mann-Whitney U tests showed that a set of scalar parameters (i.e., the number, average size, and average nearest neighbor distance of trabecular plates and rods) were statistically different among the individual specimens (p<0.05). Due to the commonality of the underlying distributions, the individual differences in the trabecular microstructure among the specimens seemed to be reflected primarily by changes in the scalar parameters. The above results strongly support the hypothesis of this study and may shed more light on understanding the natural design of trabecular bone microstructures.

Adults with cystic fibrosis have deficits in bone structure and strength at the distal tibia despite similar size and measuring standard and relative sites

Individuals with cystic fibrosis (CF) have lower bone mineral density (BMD) by DXA and are at higher risk of fracture than healthy controls. However, the 2-dimensional measurement of areal BMD (aBMD) provided by DXA is influenced by bone size and the true extent of the bone deficit is unclear. Our objective was to use high-resolution peripheral quantitative computed tomography (HR-pQCT) and individual trabecula segmentation (ITS) analysis to compare volumetric BMD (vBMD), microarchitecture and estimated strength at the distal radius and tibia in 26 young adults with CF and 26 controls matched for age, gender, and race. To assess the effect of limb length and minimize the confounding effects of size on HR-pQCT outcomes, we scanned participants at both the standard fixed HR-pQCT measurement sites and at a subject-specific relative site that varied according to limb length. CF participants did not differ significantly in age, height, weight, or BMI from controls. Ulnar and tibial lengths were 9mm shorter in CF patients, though differences were not significant. CF patients had significantly lower BMI-adjusted aBMD by DXA at the lumbar spine (8.9%, p<0.01), total hip (11.5%, p<0.01) and femoral neck (14.5%, p<0.01), but not at the forearm. At the fixed radius site, thickness of trabecular plates and torsional stiffness were significantly lower in CF participants than controls. At the relative radius site, only torsional stiffness was significantly lower in CF participants. At the tibia, total, trabecular and cortical vBMD were significantly lower at both fixed and relative sites in CF participants, with fewer, more widely-spaced trabecular plates, lower trabecular connectivity, and lower axial and torsional stiffness. Our results confirm that aBMD is lower at the spine and hip in young adults with CF, independent of BMI and body size. We also conclude that vBMD and stiffness are lower at the weight-bearing tibia. The pathogenesis of these differences in bone density and strength at the tibia appear to be related to trabecular drop-out and reduced trabecular connectivity and to be independent of differences in limb length, as assessed by scanning participants at both standard and relative sites. We concluded that significant deficits in bone structure and strength persist in young adults with CF, despite advances in care that permit them to attain relatively normal height and weight.

Validation of calcaneus trabecular microstructure measurements by HR-pQCT

OBJECTIVE

Assessment of calcaneus microstructure using high-resolution peripheral quantitative computed tomography (HR-pQCT) might be used to improve fracture risk predictions or to assess responses to pharmacological and physical interventions. To develop a standard clinical protocol for the calcaneus, we validated calcaneus trabecular microstructure measured by HR-pQCT against 'gold-standard' micro-CT measurements.

METHODS

Ten human cadaveric feet were scanned in situ using HR-pQCT (isotropic 82μm voxel size) at 100, 150 and 200ms integration times, and at 100ms integration time following removal of the calcaneus from the foot (ex vivo). Dissected portions of these bones were scanned using micro-computed tomography (micro-CT) at an isotropic 17.4μm voxel size. HR-pQCT images were rigidly registered to those obtained with micro-CT and divided into multiple 5mm sided cubes to evaluate and compare morphometric parameters between the modalities. Standard HR-pQCT measurements (derived bone volume fraction (BV/TVd); trabecular number, Tb.N; derived trabecular thickness, Tb.Thd; derived trabecular spacing, Tb.Spd) and corresponding micro-CT voxel-based measurements (BV/TV, Tb.N, Tb.Th, Tb.Sp) were compared.

RESULTS

A total of 108 regions of interest were analysed across the 10 specimens. At all integration times HR-pQCT BV/TVd was strongly correlated with micro-CT BV/TV (r2=0.95-0.98, RMSE=1%), but BV/TVd was systematically lower than that measured by micro-CT (mean bias=5%). In contrast, HR-pQCT systematically overestimated Tb.N at all integration times; of the in situ scans, 200ms yielded the lowest mean bias and the strongest correlation with micro-CT (r2=0.61, RMSE=0.15mm-1). Regional analysis revealed greater accuracy for Tb.N in the superior regions of the calcaneus at all integration times in situ (mean bias=0.44-0.85mm-1; r2=0.70-0.88, p<0.001 versus mean bias=0.63-1.46mm-1; r2≤0.08, p≥0.21 for inferior regions). Tb.Spd was underestimated by HR-pQCT compared to micro-CT, but showed similar trends with integration time and the region evaluated as Tb.N. HR-pQCT Tb.Thd was also underestimated and moderately correlated (r2=0.53-0.59) with micro-CT Tb.Th, independently from the integration time. Stronger correlations, smaller biases and error were found in the scans of the calcaneus ex vivo compared to in situ.

CONCLUSION

Calcaneus trabecular BV/TVd and trabecular microstructure, particularly in the superior region of the calcaneus, can be assessed by HR-pQCT. The highest integration time examined, 200ms, compared best with micro-CT. Weaker correlations for microstructure at inferior regions, and also with lower integration times, might limit the use of the proposed protocol, which warrants further investigation in vivo.

Trabecular Bone Morphology Correlates With Skeletal Maturity and Body Composition in Healthy Adolescent Girls

Context

Growth in healthy children is associated with changes in bone density and microarchitecture. Trabecular morphology is an additional important determinant of bone strength, but little is currently known about trabecular morphology in healthy young people.

Objective

To investigate associations of trabecular morphology with increasing maturity and with body composition in healthy girls.

Design

Cross-sectional study.

Setting

Academic research center.

Participants

Eighty-six healthy girls aged 9 to 18 years.

Main Outcome Measures

High-resolution peripheral quantitative computed tomography and individual trabecula segmentation were used to assess volumetric bone density, microarchitecture, and trabecular morphology (plate-like vs rod-like) at the distal radius and tibia.

Results

Plate-like bone volume divided by total volume (pBV/TV) increased statistically significantly at the tibia (R = 0.41, P < 0.001), whereas rod-like BV/TV (rBV/TV) decreased statistically significantly at both the radius and tibia (R = -0.34, P = 0.003 and R = -0.28, P = 0.008, respectively) with increasing bone age. In multivariable models, lean mass positively correlated with pBV/TV and plate number at the radius and with plate thickness at both sites. In contrast, fat mass negatively correlated with plate thickness at the tibia and plate surface at both sites. In addition, fat mass positively correlated with rBV/TV and number at the tibia. pBV/TV at both the distal radius and tibia was positively correlated with spine bone mineral density.

Conclusions

Increasing maturity across late childhood and adolescence is associated with changes in trabecular morphology anticipated to contribute to bone strength. Body composition correlates with trabecular morphology, suggesting that muscle mass and adiposity in youth may contribute to long-term skeletal health.

Cortical microstructure compensates for smaller bone size in young Caribbean Hispanic versus non-Hispanic white men

Hispanic men have smaller bone size but thicker and denser cortices compared to white men, leading to similar mechanical competence.

INTRODUCTION

The purpose of this study was to assess differences in vBMD and microarchitecture in young Caribbean Hispanic (n = 30) and non-Hispanic Caucasian (n = 30) men.

METHODS

We measured areal bone mineral density (aBMD) at the spine, total hip (TH), femoral neck (FN), and forearm by dual-energy X-ray absorptiometry (DXA) and bone geometry, mass, microarchitecture, and mechanical competence by high-resolution peripheral quantitative computed tomography (HRpQCT), individual trabecula segmentation (ITS), and finite element analysis (FEA).

RESULTS

Hispanic men were slightly older, shorter, and heavier and had higher BMI compared with white men. aBMD, measured by DXA, did not differ at the spine, TH, or forearm before or after adjustment for age, height, weight, and the interaction of height and weight. At the FN, marginally significant higher BMD in Hispanics prior to adjustment was attenuated and no longer differed after adjustment for covariates. Adjusted HRpQCT indices indicated smaller total and trabecular area at the radius but greater total volumetric density and cortical thickness in Hispanic versus white men. The adjusted difference in cortical density at the radius was of borderline significance. Trabecular and ITS microstructure tended not to differ at the radius. At the tibia, results were similar. Bone size tended to be smaller and covariate-adjusted cortical density and cortical thickness were greater in Hispanic versus white men. Additionally, cortical porosity was lower at the tibia in Hispanic compared to white men. Stiffness and failure load did not differ at either skeletal site by ethnicity.

CONCLUSION

In conclusion, greater cortical thickness and density as well as lower cortical porosity tend to compensate for smaller bone size in Hispanic men, leading to similar mechanical competence compared with white men.

Differences in Trabecular Microstructure Between Black and White Women Assessed by Individual Trabecular Segmentation Analysis of HR-pQCT Images

Black women have lower fracture risk compared with white women, which may be partly explained by improved volumetric bone mineral density (vBMD) and bone microarchitecture primarily within the cortical bone compartment. To determine if there are differences in trabecular microstructure, connectivity, and alignment according to race/ethnicity, we performed individual trabecular segmentation (ITS) analyses on high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of the distal radius and tibia in 273 peri- and postmenopausal black (n = 100) and white (n = 173) women participating in the Study of Women's Health Across the Nation in Boston. Unadjusted analyses showed that black women had greater trabecular plate volume fraction, plate thickness, plate number density, and plate surface area along with greater axial alignment of trabeculae, whereas white women had greater trabecular rod tissue fraction (p < 0.05 for all). Adjustment for clinical covariates augmented these race/ethnicity-related differences in plates and rods, such that white women had greater trabecular rod number density and rod-rod connectivity, whereas black women continued to have superior plate structural characteristics and axial alignment (p < 0.05 for all). These differences remained significant after adjustment for hip BMD and trabecular vBMD. In conclusion, black women had more plate-like trabecular morphology and higher axial alignment of trabeculae, whereas white women had more rod-like trabeculae. These differences may contribute to the improved bone strength and lower fracture risk observed in black women. © 2016 American Society for Bone and Mineral Research.

Differences in bone quality and strength between Asian and Caucasian young men

This is a cross-sectional study to assess differences in bone quality in young Asian and Caucasian (n = 30/group) men between 25 and 35 years. We found that Asians had smaller bones, thicker and denser cortices, and more plate-like trabeculae, but stiffness did not differ between groups.

INTRODUCTION

We conducted a cross-sectional study to assess differences in bone quality in young Asian and Caucasian (n = 30/group) men between 25 and 35 years.

METHODS

We measured bone mineral density (BMD) at the spine, total hip (TH), femoral neck (FN), and forearm by dual energy X-ray absorptiometry (DXA), and bone geometry, density, microarchitecture, and mechanical competence at the radius and tibia by high-resolution peripheral quantitative computed tomography (HR-pQCT) with application of individual trabecula segmentation (ITS) and trabecular and whole bone finite element analysis (FEA). We measured load-to-strength ratio to account for differences in bone size and height, respectively. We used Wilcoxon rank sum and generalized linear models adjusted for height, weight, and their interaction for comparisons.

RESULTS

Asians were 3.9 % shorter and weighed 6.5 % less than Caucasians. In adjusted models: by DXA, there were no significant race-based differences in areal BMD; by HR-pQCT, at the radius, Asians had smaller total and trabecular area (p = 0.003 for both), and denser (p = 0.01) and thicker (p = 0.04) cortices at the radius; by ITS, at the radius Asians, had more plate-like than rod-like trabeculae (PR ratio p = 0.01), greater plate trabecular surface (p = 0.009) and longer rod length (p = 0.002). There were no significant race-based differences in FEA or the load-to-strength ratio.

CONCLUSIONS

Asians had smaller bones, thicker and denser cortices, and more plate-like trabeculae, but biomechanical estimates of bone strength did not differ between groups. Studies are needed to determine whether these differences persist later in life.

Young adults with cystic fibrosis have altered trabecular microstructure by ITS-based morphological analysis

Young adults with cystic fibrosis have compromised plate-like trabecular microstructure, altered axial alignment of trabeculae, and reduced connectivity between trabeculae that may contribute to the reduced bone strength and increased fracture risk observed in this patient population.

INTRODUCTION

The risk of fracture is increased in patients with cystic fibrosis (CF). Individual trabecular segmentation (ITS)-based morphological analysis of high-resolution peripheral quantitative computed tomography (HR-pQCT) images segments trabecular bone into individual plates and rods of different alignment and connectivity, which are important determinants of trabecular bone strength. We sought to determine whether alterations in ITS variables are present in patients with CF and may help explain their increased fracture risk.

METHODS

Thirty patients with CF ages 18-40 years underwent DXA scans of the hip and spine and HR-pQCT scans of the radius and tibia with further assessment of trabecular microstructure by ITS. These CF patients were compared with 60 healthy controls matched for age (±2 years), race, and gender.

RESULTS

Plate volume fraction, thickness, and density as well as plate-plate and plate-rod connectivity were reduced, and axial alignment of trabeculae was lower in subjects with CF at both the radius and the tibia (p < 0.05 for all). At the radius, adjustment for BMI eliminated most of these differences. At the tibia, however, reductions in plate volume fraction and number, axially aligned trabeculae, and plate-plate connectivity remained significant after adjustment for BMI alone and for BMI and aBMD (p < 0.05 for all).

CONCLUSIONS

Young adults with CF have compromised plate-like and axially aligned trabecular morphology and reduced connectivity between trabeculae. ITS analysis provides unique information about bone integrity, and these trabecular deficits may help explain the increased fracture risk in adults with CF not accounted for by BMD and/or traditional bone microarchitecture measurements.

Deterioration of trabecular plate-rod and cortical microarchitecture and reduced bone stiffness at distal radius and tibia in postmenopausal women with vertebral fractures

Postmenopausal women with vertebral fractures have abnormal bone microarchitecture at the distal radius and tibia by HR-pQCT, independent of areal BMD. However, whether trabecular plate and rod microarchitecture is altered in women with vertebral fractures is unknown. This study aims to characterize the abnormalities of trabecular plate and rod microarchitecture, cortex, and bone stiffness in postmenopausal women with vertebral fractures. HR-pQCT images of distal radius and tibia were acquired from 45 women with vertebral fractures and 45 control subjects without fractures. Trabecular and cortical compartments were separated by an automatic segmentation algorithm and subjected to individual trabecula segmentation (ITS) analysis for measuring trabecular plate and rod morphology and cortical bone evaluation for measuring cortical thickness and porosity, respectively. Whole bone and trabecular bone stiffness were estimated by finite element analysis. Fracture and control subjects did not differ according to age, race, body mass index, osteoporosis risk factors, or medication use. Women with vertebral fractures had thinner cortices, and larger trabecular area compared to the control group. By ITS analysis, fracture subjects had fewer trabecular plates, less axially aligned trabeculae and less trabecular connectivity at both the radius and the tibia. Fewer trabecular rods were observed at the radius. Whole bone stiffness and trabecular bone stiffness were 18% and 22% lower in women with vertebral fractures at the radius, and 19% and 16% lower at the tibia, compared with controls. The estimated failure load of the radius and tibia were also reduced in the fracture subjects by 13% and 14%, respectively. In summary, postmenopausal women with vertebral fractures had both trabecular and cortical microstructural deterioration at the peripheral skeleton, with a preferential loss of trabecular plates and cortical thinning. These microstructural deficits translated into lower whole bone and trabecular bone stiffness at the radius and tibia. Our results suggest that abnormalities in trabecular plate and rod microstructure may be important mechanisms of vertebral fracture in postmenopausal women.

An Advanced Quantitative Echosound Methodology for Femoral Neck Densitometry

The aim of this paper was to investigate the clinical feasibility and the accuracy in femoral neck densitometry of the Osteoporosis Score (O.S.), an ultrasound (US) parameter for osteoporosis diagnosis that has been recently introduced for lumbar spine applications. A total of 377 female patients (aged 61-70 y) underwent both a femoral dual X-ray absorptiometry (DXA) and an echographic scan of the proximal femur. Recruited patients were sub-divided into a reference database used for ultrasound spectral model construction and a study population for repeatability assessments and accuracy evaluations. Echographic images and radiofrequency signals were analyzed through a fully automatic algorithm that performed a series of combined spectral and statistical analyses, providing as a final output the O.S. value of the femoral neck. Assuming DXA as a gold standard reference, the accuracy of O.S.-based diagnoses resulted 94.7%, with k = 0.898 (p < 0.0001). Significant correlations were also found between O.S.-estimated bone mineral density and corresponding DXA values, with r(2) up to 0.79 and root mean square error = 5.9-7.4%. The reported accuracy levels, combined with the proven ease of use and very good measurement repeatability, provide the adopted method with a potential for clinical routine application in osteoporosis diagnosis.

High-resolution peripheral quantitative computed tomography (HR-pQCT) can assess microstructural and biomechanical properties of both human distal radius and tibia: Ex vivo computational and experimental validations

High-resolution peripheral quantitative computed tomography (HR-pQCT) provides in vivo three-dimensional (3D) imaging at the distal radius and tibia and has been increasingly used to characterize cortical and trabecular bone morphology in clinical studies. In this study, we comprehensively examined the accuracy of HR-pQCT and HR-pQCT based micro finite element (μFE) analysis predicted bone elastic stiffness and strength through comparisons with gold-standard micro computed tomography (μCT) based morphological/μFE measures and direct mechanical testing results. Twenty-six sets of human cadaveric distal radius and tibia segments were imaged by HR-pQCT and μCT. Microstructural analyses were performed for the registered HR-pQCT and μCT images. Bone stiffness and yield strength were determined by both HR-pQCT and μCT based linear and nonlinear μFE predictions and mechanical testing. Our results suggested that strong and significant correlations existed between the HR-pQCT standard, model-independent and corresponding μCT measurements. HR-pQCT based nonlinear μFE overestimated stiffness and yield strength while the linear μFE underestimated yield strength, but both were strongly correlated with those predicted by μCT μFE and measured by mechanical testing at both radius and tibia (R(2)>0.9). The microstructural differences between HR-pQCT and μCT were also examined by the Bland-Altman plots. Our results showed HR-pQCT morphological measurements of BV/TV(d), Tb.Th, and Tb.Sp, can be adjusted by correction values to approach true values measured by gold-standard μCT. In addition, we observed moderate correlations of HR-pQCT biomechanical and microstructural parameters between the distal radius and tibia. We concluded that morphological and mechanical properties of human radius and tibia bone can be assessed by HR-pQCT based measures.

Altered trabecular bone morphology in adolescent and young adult athletes with menstrual dysfunction

CONTEXT

Young amenorrheic athletes (AA) have lower bone mineral density (BMD) and an increased prevalence of fracture compared with eumenorrheic athletes (EA) and non-athletes. Trabecular morphology is a determinant of skeletal strength and may contribute to fracture risk.

OBJECTIVES

To determine the variation in trabecular morphology among AA, EA, and non-athletes and to determine the association of trabecular morphology with fracture among AA.

DESIGN AND SETTING

A cross-sectional study performed at an academic clinical research center.

PARTICIPANTS

161 girls and young women aged 14-26 years (97 AA, 32 EA, and 32 non-athletes).

MAIN OUTCOME MEASURE

We measured volumetric BMD (vBMD) and skeletal microarchitecture using high-resolution peripheral quantitative computed tomography. We evaluated trabecular morphology (plate-like vs. rod-like), orientation, and connectivity by individual trabecula segmentation.

RESULTS

At the non-weight-bearing distal radius, the groups did not differ for trabecular vBMD. However, plate-like trabecular bone volume fraction (pBV/TV) was lower in AA vs. EA (p=0.03), as were plate number (p=0.03) and connectivity (p=0.03). At the weight-bearing distal tibia, trabecular vBMD was higher in athletes vs. non-athletes (p=0.05 for AA and p=0.009 for EA vs. non-athletes, respectively). pBV/TV was higher in athletes vs. non-athletes (p=0.04 AA and p=0.005 EA vs. non-athletes), as were axially-aligned trabeculae, plate number, and connectivity. Among AA, those with a history of recurrent stress fracture had lower pBV/TV, axially-aligned trabeculae, plate number, plate thickness, and connectivity at the distal radius.

CONCLUSIONS

Trabecular morphology and alignment differ among AA, EA, and non-athletes. These differences may be associated with increased fracture risk.

Three-dimensional evaluation of human jaw bone microarchitecture: correlation between the microarchitectural parameters of cone beam computed tomography and micro-computer tomography

OBJECTIVE

To evaluate the potential feasibility of cone beam computed tomography (CBCT) in the assessment of trabecular bone microarchitecture.

STUDY DESIGN

Sixty-eight specimens from four pairs of human jaw were scanned using both micro-computed tomography (micro-CT) of 19.37-μm voxel size and CBCT of 100-μm voxel size. The correlation of 3-dimensional parameters between CBCT and micro-CT was evaluated.

RESULTS

All parameters, except bone-specific surface and trabecular thickness, showed linear correlations between the 2 imaging modalities (P < .05). Among the parameters, bone volume, percent bone volume, trabecular separation, and degree of anisotropy (DA) of CBCT images showed strong correlations with those of micro-CT images. DA showed the strongest correlation (r = 0.693).

CONCLUSIONS

Most microarchitectural parameters from CBCT were correlated with those from micro-CT. Some microarchitectural parameters, especially DA, could be used as strong predictors of bone quality in the human jaw.

In Vivo Precision of Digital Topological Skeletonization Based Individual Trabecula Segmentation (ITS) Analysis of Trabecular Microstructure at the Distal Radius and Tibia by HR-pQCT

Trabecular plate and rod microstructure plays a dominant role in the apparent mechanical properties of trabecular bone. With high-resolution computed tomography (CT) images, digital topological analysis (DTA) including skeletonization and topological classification was applied to transform the trabecular three-dimensional (3D) network into surface and curve skeletons. Using the DTA-based topological analysis and a new reconstruction/recovery scheme, individual trabecula segmentation (ITS) was developed to segment individual trabecular plates and rods and quantify the trabecular plate- and rod-related morphological parameters. High-resolution peripheral quantitative computed tomography (HR-pQCT) is an emerging in vivo imaging technique to visualize 3D bone microstructure. Based on HR-pQCT images, ITS was applied to various HR-pQCT datasets to examine trabecular plate- and rod-related microstructure and has demonstrated great potential in cross-sectional and longitudinal clinical applications. However, the reproducibility of ITS has not been fully determined. The aim of the current study is to quantify the precision errors of ITS plate-rod microstructural parameters. In addition, we utilized three different frequently used contour techniques to separate trabecular and cortical bone and to evaluate their effect on ITS measurements. Overall, good reproducibility was found for the standard HR-pQCT parameters with precision errors for volumetric BMD and bone size between 0.2%-2.0%, and trabecular bone microstructure between 4.9%-6.7% at the radius and tibia. High reproducibility was also achieved for ITS measurements using all three different contour techniques. For example, using automatic contour technology, low precision errors were found for plate and rod trabecular number (pTb.N, rTb.N, 0.9% and 3.6%), plate and rod trabecular thickness (pTb.Th, rTb.Th, 0.6% and 1.7%), plate trabecular surface (pTb.S, 3.4%), rod trabecular length (rTb.ℓ, 0.8%), and plate-plate junction density (P-P Junc.D, 2.3%) at the tibia. The precision errors at the radius were similar to those at the tibia. In addition, precision errors were affected by the contour technique. At the tibia, precision error by the manual contour method was significantly different from automatic and standard contour methods for pTb.N, rTb.N and rTb.Th. Precision error using the manual contour method was also significantly different from the standard contour method for rod trabecular number (rTb.N), rod trabecular thickness (rTb.Th), rod-rod and plate-rod junction densities (R-R Junc.D and P-R Junc.D) at the tibia. At the radius, the precision error was similar between the three different contour methods. Image quality was also found to significantly affect the ITS reproducibility. We concluded that ITS parameters are highly reproducible, giving assurance that future cross-sectional and longitudinal clinical HR-pQCT studies are feasible in the context of limited sample sizes.

Trabecular plates and rods determine elastic modulus and yield strength of human trabecular bone

The microstructure of trabecular bone is usually perceived as a collection of plate-like and rod-like trabeculae, which can be determined from the emerging high-resolution skeletal imaging modalities such as micro-computed tomography (μCT) or clinical high-resolution peripheral quantitative CT (HR-pQCT) using the individual trabecula segmentation (ITS) technique. It has been shown that the ITS-based plate and rod parameters are highly correlated with elastic modulus and yield strength of human trabecular bone. In the current study, plate-rod (PR) finite element (FE) models were constructed completely based on ITS-identified individual trabecular plates and rods. We hypothesized that PR FE can accurately and efficiently predict elastic modulus and yield strength of human trabecular bone. Human trabecular bone cores from proximal tibia (PT), femoral neck (FN) and greater trochanter (GT) were scanned by μCT. Specimen-specific ITS-based PR FE models were generated for each μCT image and corresponding voxel-based FE models were also generated in comparison. Both types of specimen-specific models were subjected to nonlinear FE analysis to predict the apparent elastic modulus and yield strength using the same trabecular bone tissue properties. Then, mechanical tests were performed to experimentally measure the apparent modulus and yield strength. Strong linear correlations for both elastic modulus (r(2) = 0.97) and yield strength (r(2) = 0.96) were found between the PR FE model predictions and experimental measures, suggesting that trabecular plate and rod morphology adequately captures three-dimensional (3D) microarchitecture of human trabecular bone. In addition, the PR FE model predictions in both elastic modulus and yield strength were highly correlated with the voxel-based FE models (r(2) = 0.99, r(2) = 0.98, respectively), resulted from the original 3D images without the PR segmentation. In conclusion, the ITS-based PR models predicted accurately both elastic modulus and yield strength determined experimentally across three distinct anatomic sites. Trabecular plates and rods accurately determine elastic modulus and yield strength of human trabecular bone.

A novel ultrasound methodology for estimating spine mineral density

We investigated the possible clinical feasibility and accuracy of an innovative ultrasound (US) method for diagnosis of osteoporosis of the spine. A total of 342 female patients (aged 51-60 y) underwent spinal dual X-ray absorptiometry and abdominal echographic scanning of the lumbar spine. Recruited patients were subdivided into a reference database used for US spectral model construction and a study population for repeatability and accuracy evaluation. US images and radiofrequency signals were analyzed via a new fully automatic algorithm that performed a series of spectral and statistical analyses, providing a novel diagnostic parameter called the osteoporosis score (O.S.). If dual X-ray absorptiometry is assumed to be the gold standard reference, the accuracy of O.S.-based diagnoses was 91.1%, with k = 0.859 (p < 0.0001). Significant correlations were also found between O.S.-estimated bone mineral densities and corresponding dual X-ray absorptiometry values, with r(2) values up to 0.73 and a root mean square error of 6.3%-9.3%. The results obtained suggest that the proposed method has the potential for future routine application in US-based diagnosis of osteoporosis.

The Ellipsoid Factor for Quantification of Rods, Plates, and Intermediate Forms in 3D Geometries

The ellipsoid factor (EF) is a method for the local determination of the rod- or plate-like nature of porous or spongy continua. EF at a point within a 3D structure is defined as the difference in axis ratios of the greatest ellipsoid that fits inside the structure and that contains the point of interest, and ranges from -1 for strongly oblate (discus-shaped) ellipsoids, to +1 for strongly prolate (javelin-shaped) ellipsoids. For an ellipsoid with axes a ≤ b ≤ c, EF = a/b - b/c. Here, EF is demonstrated in a Java plugin, "Ellipsoid Factor" for ImageJ, distributed in the BoneJ plugin collection. Ellipsoid Factor utilizes an ellipsoid optimization algorithm, which assumes that maximal ellipsoids are centered on the medial axis, then dilates, rotates, and translates slightly each ellipsoid until it cannot increase in volume any further. EF successfully identifies rods, plates, and intermediate structures within trabecular bone, and summarizes the distribution of geometries with an overall EF mean and SD, EF histogram, and Flinn diagram displaying a/b versus b/c. EF is released to the community for testing, use, and improvement.

Bone health in HIV-infected children, adolescents and young adults: a systematic review

BACKGROUND

Children and adolescents, who either acquire HIV infection perinatally, from contaminated blood products or via sexual transmission early in life, have the greatest cumulative exposure to the negative direct and indirect effects of HIV infection and ART on bone, which may lead to increased lifetime risk for osteoporosis and fracture. We conducted a systematic review to evaluate the literature on bone health in children and adolescents with HIV.

METHODS

We performed a comprehensive search of the Medline, Scopus, and Cochrane Library databases (up to April 1, 2014) for studies that reported on bone imaging or bone fractures in HIV-infected children, adolescents, or young adults.

RESULTS

A total of 32 publications met our inclusion criteria. Seventeen studies were cross-sectional and 15 were longitudinal. The majority of studies were conducted in high-income countries, three in middle-income countries and none in low-income countries. Overall, the studies we reviewed indicate that measures of bone mass are reduced, with increased prevalence of low BMD in children and adolescents with HIV. However, the studies are highly variable with respect to comparison sources, measurement methods, adjustment techniques for body size or growth retardation, and highlighted risk factors, including aspects related to medication exposures as well as the effects of HIV infection per se.

CONCLUSION

HIV infection appears to be associated with decreased bone accrual throughout childhood and adolescence. Initial studies indicate that sub-optimal bone accrual may be persistent and result in reduced peak bone mass, an important determinant of future risk of osteoporosis and fracture. Important areas for future research include evaluation of bone mass, bone quality and fracture risk across the life course among those with early-life infection with HIV, particularly in resource-limited settings where the majority of children with HIV live.

Bone density, microarchitecture and stiffness in Caucasian and Caribbean Hispanic postmenopausal American women

Hispanic Americans of Caribbean origin are a fast-growing subset of the US population, but there are no studies on bone density, microstructure and biomechanical integrity in this minority group. In this study, we aimed to compare Caucasian and Caribbean Hispanic postmenopausal American women with respect to these characteristics. Thirty-three Caribbean Hispanics were age-matched to thirty-three Caucasian postmenopausal women. At the lumbar spine, the Hispanic women had significantly lower areal bone mineral density (aBMD). At the radius by high-resolution peripheral quantitative computed tomography (HR-pQCT), there were minimal differences between Hispanic and Caucasian women. At the tibia, Hispanic women had lower trabecular volumetric bone density and trabecular number, and higher trabecular separation. Individual trabecula segmentation (ITS) analyses indicated that at the tibia, Hispanic women not only had significantly lower bone volume fraction, but also had significantly lower rod bone volume fraction, plate trabecular number, rod trabecular number and lower plate–plate, plate–rod and rod–rod junction densities compared to Caucasian women. The differences in bone quantity and quality contributed to lower whole bone stiffness at the radius, and both whole bone and trabecular bone stiffness at the tibia in Hispanic women. In conclusion, Hispanic women had poorer bone mechanical and microarchitectural properties than Caucasian women, especially at the load-bearing distal tibia.

The developing juvenile ischium: macro-radiographic insights

Despite the importance of the human pelvis as a weight-bearing structure, there is a paucity of literature that discusses the development of the juvenile innominate from a biomechanical perspective. This study aims to add to the limited body of literature pertaining to this topic through the qualitative analysis of the gross architecture of the human ischium during the juvenile period. Macro-radiographs of 55 human ischia ranging from 28 intra-uterine weeks to 14 years of age were examined using intensity-gradient color mapping to highlight changes in gross structural morphology with increasing age. A clear pattern of maturation was observed in the juvenile ischium with increasing age. The acetabular component and ramus of the ischium consistently displayed low bone intensity in the postnatal skeletal material. Conversely the posterior body of the ischium, and in particular the ischial spine and lesser sciatic notch, exhibited increasing bone intensity which first arose at 1-2 years of age and became more expansive in older cohorts. The intensity patterns observed within the developing juvenile ischium are indicative of the potential factors influencing the maturation of this skeletal element. While the low intensity acetabular fossa indicates a lack of significant biomechanical interactions, the posterior increase in bone intensity may be related to the load-bearing nature of the posterior ischium.

Teriparatide increases strength of the peripheral skeleton in premenopausal women with idiopathic osteoporosis: a pilot HR-pQCT study

CONTEXT

In premenopausal women with idiopathic osteoporosis (IOP), treatment with teriparatide leads to substantial improvement in bone density and quality at central skeletal sites. The effects of teriparatide may differ on cortical and trabecular bone and also at the central and the peripheral skeleton.

OBJECTIVE

The objective of the study was to determine whether teriparatide was associated with improvements in compartmental volumetric bone mineral density (BMD), bone microarchitecture, and estimated bone strength of the distal radius and tibia as assessed by high-resolution peripheral quantitative computed tomography.

DESIGN, SETTING, AND PARTICIPANTS

Premenopausal women (n = 20, age 41 ± 5 y) with IOP (low trauma fractures and/or Z-scores ≤ -2.0) were scanned with high-resolution peripheral quantitative computed tomography at baseline and after 18 months of teriparatide treatment. Cortical and trabecular volumetric BMD and microarchitecture were measured by both standard and advanced techniques, including individual trabecula segmentation, and bone strength was estimated by finite element analysis.

MAIN OUTCOME MEASURES

The total volumetric BMD and homogeneous bone stiffness were measured.

RESULTS

Trabecular volumetric BMD increased significantly by 2.6% (1.8, 6.2) [median (interquartile range)] at the radius and 2.5% (1.1, 3.6) at the tibia. In addition, trabecular plate bone volume fraction increased by 9.1% (2.1, 17.1) at the radius and 7.6% (1.0, 9.7) at the tibia. Cortical thickness and volumetric density did not change; however, cortical porosity increased at the radius but not at the tibia. Despite these changes, whole-bone stiffness and failure load estimated by finite element analysis increased at both the radius and tibia.

CONCLUSIONS

In premenopausal women with IOP, 18 months of teriparatide was associated with increases in trabecular volumetric BMD, improved trabecular microarchitecture, and estimated bone strength at both the radius and tibia.

Lower peak bone mass and abnormal trabecular and cortical microarchitecture in young men infected with HIV early in life

INTRODUCTION

HIV infection and antiretroviral therapy (ART) early in life may interfere with acquisition of peak bone mass, thereby increasing fracture risk in adulthood.

METHODS

We conducted a cross-sectional study of dual-energy X-ray absorptiometry (DXA) and high-resolution peripheral quantitative computed tomography (HR-pQCT) in 30 HIV-infected African-American or Hispanic Tanner stage 5 men aged 20-25 on ART (15 perinatally infected and 15 infected during adolescence) and 15 HIV-uninfected controls.

RESULTS

HIV-infected men were similar in age and BMI, but were more likely to be African-American (P = 0.01) than uninfected men. DXA-derived areal bone mineral density (aBMD) Z-scores were 0.4-1.2 lower in HIV-infected men at the spine, hip, and radius (all P < 0.05). At the radius and tibia, total and trabecular volumetric BMD (vBMD), and cortical and trabecular thickness were between 6 and 19% lower in HIV-infected than uninfected men (P <0.05). HIV-infected men had dramatic deficiencies in plate-related parameters by individual trabeculae segmentation (ITS) analyses and 14-17% lower bone stiffness by finite element analysis. Differences in most HR-pQCT parameters remained significant after adjustment for race/ethnicity. No DXA or HR-pQCT parameters differed between men infected perinatally or during adolescence.

CONCLUSION

At an age by which young men have typically acquired peak bone mass, HIV-infected men on ART have lower BMD, markedly abnormal trabecular plate and cortical microarchitecture, and decreased whole bone stiffness, whether infected perinatally or during adolescence. Reduced bone strength in young adults infected with HIV early in life may place them at higher risk for fractures as they age.

Are patterns of bone loss in anorexic and postmenopausal women similar? Preliminary results using high resolution peripheral computed tomography

This study intended to compare bone density and architecture in three groups of women: young women with anorexia nervosa (AN), an age-matched control group of young women, and healthy late postmenopausal women. Three-dimensional peripheral quantitative high resolution computed-tomography (HR-pQCT) at the ultradistal radius, a technology providing measures of cortical and trabecular bone density and microarchitecture, was performed in the three cohorts. Thirty-six women with AN aged 18-30 years (mean duration of AN: 5.8 years), 83 healthy late postmenopausal women aged 70-81 as well as 30 age-matched healthy young women were assessed. The overall cortical and trabecular bone density (D100), the absolute thickness of the cortical bone (CTh), and the absolute number of trabecules per area (TbN) were significantly lower in AN patients compared with healthy young women. The absolute number of trabecules per area (TbN) in AN and postmenopausal women was similar, but significantly lower than in healthy young women. The comparison between AN patients and post-menopausal women is of interest because the latter reach bone peak mass around the middle of the fertile age span whereas the former usually lose bone before reaching optimal bone density and structure. This study shows that bone mineral density and bone compacta thickness in AN are lower than those in controls but still higher than those in postmenopause. Bone compacta density in AN is similar as in controls. However, bone inner structure in AN is degraded to a similar extent as in postmenopause. This last finding is particularly troubling.

Dependence of mechanical properties of trabecular bone on plate-rod microstructure determined by individual trabecula segmentation (ITS)

Individual trabecula segmentation (ITS) technique can decompose the trabecular bone network into individual trabecular plates and rods and is capable of quantifying the plate/rod-related microstructural characteristics of trabecular bone. This novel technique has been shown to be able to provide in-depth insights into micromechanics and failure mechanisms of human trabecular bone, as well as to distinguish the fracture status independent of area bone mineral density in clinical applications. However, the plate/rod microstructural parameters from ITS have never been correlated to experimentally determined mechanical properties of human trabecular bone. In this study, on-axis cylindrical trabecular bone samples from human proximal tibia (n=22), vertebral body (n=10), and proximal femur (n=21) were harvested, prepared, scanned using micro computed-tomography (µCT), analyzed with ITS and mechanically tested. Regression analyses showed that the plate bone volume fraction (pBV/TV) and axial bone volume fraction (aBV/TV) calculated by ITS analysis correlated the best with elastic modulus (R(2)=0.96-0.97) and yield strength (R(2)=0.95-0.96). Trabecular plate-related microstructural parameters correlated highly with elastic modulus and yield strength, while most rod-related parameters were found inversely and only moderately correlated with the mechanical properties. In addition, ITS analysis also identified that trabecular bone at human femoral neck had the highest trabecular plate-related parameters while the other sites were similar with each other in terms of plate-rod microstructure.

Fast trabecular bone strength predictions of HR-pQCT and individual trabeculae segmentation-based plate and rod finite element model discriminate postmenopausal vertebral fractures

Although high-resolution peripheral quantitative computed tomography (HR-pQCT) has advanced clinical assessment of trabecular bone microstructure, nonlinear microstructural finite element (µFE) prediction of yield strength using a HR-pQCT voxel model is impractical for clinical use due to its prohibitively high computational costs. The goal of this study was to develop an efficient HR-pQCT-based plate and rod (PR) modeling technique to fill the unmet clinical need for fast bone strength estimation. By using an individual trabecula segmentation (ITS) technique to segment the trabecular structure into individual plates and rods, a patient-specific PR model was implemented by modeling each trabecular plate with multiple shell elements and each rod with a beam element. To validate this modeling technique, predictions by HR-pQCT PR model were compared with those of the registered high-resolution micro-computed tomography (HR-µCT) voxel model of 19 trabecular subvolumes from human cadaveric tibia samples. Both the Young's modulus and yield strength of HR-pQCT PR models strongly correlated with those of µCT voxel models (r²  = 0.91 and 0.86). Notably, the HR-pQCT PR models achieved major reductions in element number (>40-fold) and computer central processing unit (CPU) time (>1200-fold). Then, we applied PR model µFE analysis to HR-pQCT images of 60 postmenopausal women with (n = 30) and without (n = 30) a history of vertebral fracture. HR-pQCT PR model revealed significantly lower Young's modulus and yield strength at the radius and tibia in fracture subjects compared to controls. Moreover, these mechanical measurements remained significantly lower in fracture subjects at both sites after adjustment for areal bone mineral density (aBMD) T-score at the ultradistal radius or total hip. In conclusion, we validated a novel HR-pQCT PR model of human trabecular bone against µCT voxel models and demonstrated its ability to discriminate vertebral fracture status in postmenopausal women. This accurate nonlinear µFE prediction of the HR-pQCT PR model, which requires only seconds of desktop computer time, has tremendous promise for clinical assessment of bone strength.

High-resolution peripheral quantitative computed tomography for the assessment of bone strength and structure: a review by the Canadian Bone Strength Working Group

Bone structure is an integral determinant of bone strength. The availability of high resolution peripheral quantitative computed tomography (HR-pQCT) has made it possible to measure three-dimensional bone microarchitecture and volumetric bone mineral density in vivo, with accuracy previously unachievable and with relatively low-dose radiation. Recent studies using this novel imaging tool have increased our understanding of age-related changes and sex differences in bone microarchitecture, as well as the effect of different pharmacological therapies. One advantage of this novel tool is the use of finite element analysis modelling to non-invasively estimate bone strength and predict fractures using reconstructed three-dimensional images. In this paper, we describe the strengths and limitations of HR-pQCT and review the clinical studies using this tool.

Primary hyperparathyroidism is associated with abnormal cortical and trabecular microstructure and reduced bone stiffness in postmenopausal women

Typically, in the milder form of primary hyperparathyroidism (PHPT), now seen in most countries, bone density by dual-energy X-ray absorptiometry (DXA) and detailed analyses of iliac crest bone biopsies by histomorphometry and micro-computed tomography (µCT) show detrimental effects in cortical bone, whereas the trabecular site (lumbar spine by DXA) and the trabecular compartment (by bone biopsy) appear to be relatively well preserved. Despite these findings, fracture risk at both vertebral and nonvertebral sites is increased in PHPT. Emerging technologies, such as high-resolution peripheral quantitative computed tomography (HRpQCT), may provide additional insight into microstructural features at sites such as the forearm and tibia that have heretofore not been easily accessible. Using HRpQCT, we determined cortical and trabecular microstructure at the radius and tibia in 51 postmenopausal women with PHPT and 120 controls. Individual trabecula segmentation (ITS) and micro-finite element (µFE) analyses of the HRpQCT images were also performed to further understand how the abnormalities seen by HRpQCT might translate into effects on bone strength. Women with PHPT showed, at both sites, decreased volumetric densities at trabecular and cortical compartments, thinner cortices, and more widely spaced and heterogeneously distributed trabeculae. At the radius, trabeculae were thinner and fewer in PHPT. The radius was affected to a greater extent in the trabecular compartment than the tibia. ITS analyses revealed, at both sites, that plate-like trabeculae were depleted, with a resultant reduction in the plate/rod ratio. Microarchitectural abnormalities were evident by decreased plate-rod and plate-plate junctions at the radius and tibia, and rod-rod junctions at the radius. These trabecular and cortical abnormalities resulted in decreased whole-bone stiffness and trabecular stiffness. These results provide evidence that in PHPT, microstructural abnormalities are pervasive and not limited to the cortical compartment, which may help to account for increased global fracture risk in PHPT.

Accuracy of individual trabecula segmentation based plate and rod finite element models in idealized trabecular bone microstructure

Currently, specimen-specific micro finite element (μFE) analysis based micro computed tomography (μCT) images have become a major computational tool for the assessment of the mechanical properties of human trabecular bone. Despite the fine characterization of the three-dimensional (3D) trabecular microstructure based on high-resolution μCT images, conventional μFE models with each voxel converted to an element are not efficient in predicting the nonlinear failure behavior of bone due to a prohibitive computational cost. Recently, a highly efficient individual trabecula segmentation (ITS)-based plate and rod (PR) modeling technique has been developed by substituting individual plates and rods with shell and beam elements, respectively. In this technical brief, the accuracy of novel PR μFE models was examined in idealized microstructure models over a broad range of trabecular thicknesses. The Young's modulus and yield strength predicted by simplified PR models strongly correlated with those of voxel models at various voxel sizes. The conversion from voxel models to PR models resulted in an ∼762-fold reduction in the largest model size and significantly accelerated the nonlinear FE analysis. The excellent predictive power of the PR μFE models, demonstrated in an idealized trabecular microstructure, provided a quantitative mechanical basis for this promising tool for an accurate and efficient assessment of trabecular bone mechanics and fracture risk.

Local topological analysis at the distal radius by HR-pQCT: Application to in vivo bone microarchitecture and fracture assessment in the OFELY study

High-resolution peripheral quantitative computed tomography (HR-pQCT) is an in-vivo technique used to analyze the distal radius and tibia. It provides a voxel size of 82μm. In addition to providing the usual microarchitecture parameters, local topological analysis (LTA) depicting rod- and plate-like trabeculae may improve prediction of bone fragility. Thirty-three women with prevalent wrist fractures from the OFELY cohort were compared with age-matched controls. Bone microarchitecture, including the structural model index (SMI), was assessed by HR-pQCT, and micro-finite element analysis (μFE) was computed on trabecular bone images of the distal radius (XtremeCT, Scanco Medical AG). A new LTA method was applied to label each bone voxel as a rod, plate or node. Then the bone volume fraction (BV/TV*), the rod, plate and node ratios over bone volume (RV/BV*, PV/BV*, NV/BV*) or total volume (RV/TV*, PV/TV*, NV/TV*) and the rod to plate ratio (RV/PV*) were calculated. Associations between LTA parameters and wrist fractures were computed in a conditional logistic regression model. Multivariate models were tested to predict the μFE-derived trabecular bone stiffness. RV/TV* (OR=4.41 [1.05-18.62]) and BV/TV* (OR=6.45 [1.06-39.3]), were significantly associated with prevalent wrist fracture, after adjustment for ultra distal radius aBMD. Multivariate linear models including PV/TV* or BV/TV*+RV/PV* predicted trabecular stiffness with the same magnitude as those including SMI. Conversion from plates into rods was significantly associated with bone fragility, with a negative correlation between RV/PV* and trabecular bone stiffness (r=-0.63, p<0.0001). We conclude that our local topological analysis is feasible for a voxel size of 82μm. After further validation, it may improve bone fragility description.

Noninvasive imaging of bone microarchitecture

The noninvasive quantification of peripheral compartment-specific bone microarchitecture is feasible with high-resolution peripheral quantitative computed tomography (HR-pQCT) and high-resolution magnetic resonance imaging (HR-MRI). In addition to classic morphometric indices, both techniques provide a suitable basis for virtual biomechanical testing using finite element (FE) analyses. Methodical limitations, morphometric parameter definition, and motion artifacts have to be considered to achieve optimal data interpretation from imaging studies. With increasing availability of in vivo high-resolution bone imaging techniques, special emphasis should be put on quality control including multicenter, cross-site validations. Importantly, conclusions from interventional studies investigating the effects of antiosteoporotic drugs on bone microarchitecture should be drawn with care, ideally involving imaging scientists, translational researchers, and clinicians.

Individual trabecula segmentation (ITS)-based morphological analyses and microfinite element analysis of HR-pQCT images discriminate postmenopausal fragility fractures independent of DXA measurements

Osteoporosis is typically diagnosed by dual-energy X-ray absorptiometry (DXA) measurements of areal bone mineral density (aBMD). Emerging technologies, such as high-resolution peripheral quantitative computed tomography (HR-pQCT), may increase the diagnostic accuracy of DXA and enhance our mechanistic understanding of decreased bone strength in osteoporosis. Women with (n = 68) and without (n = 101) a history of postmenopausal fragility fracture had aBMD measured by DXA, trabecular plate and rod microarchitecture measured by HR-pQCT image-based individual trabecula segmentation (ITS) analysis, and whole bone and trabecular bone stiffness by microfinite element analysis (µFEA) of HR-pQCT images at the radius and tibia. DXA T-scores were similar in women with and without fractures at the spine, hip, and 1/3 radius, but lower in fracture subjects at the ultradistal radius. Trabecular microarchitecture of fracture subjects was characterized by preferential reductions in trabecular plate bone volume, number, and connectivity over rod trabecular parameters, loss of axially aligned trabeculae, and a more rod-like trabecular network. In addition, decreased thickness and size of trabecular plates were observed at the tibia. The differences between groups were greater at the radius than the tibia for plate number, rod bone volume fraction and number, and plate-rod and rod-rod junction densities. Most differences between groups remained after adjustment for T-score by DXA. At a fixed bone volume fraction, trabecular plate volume, number, and connectivity were directly associated with bone stiffness. In contrast, rod volume, number, and connectivity were inversely associated with bone stiffness. In summary, HR-pQCT-based ITS and µFEA measurements discriminate fracture status in postmenopausal women independent of DXA measurements. Moreover, these results suggest that preferential loss of plate-like trabeculae contribute to lower trabecular bone and whole bone stiffness in women with fractures. We conclude that HR-pQCT-based ITS and µFEA measurements increase our understanding of the microstructural pathogenesis of fragility fracture in postmenopausal women.

Quantitative characterization of trabecular bone micro-architecture using tensor scale and multi-detector CT imaging

Osteoporosis, characterized by low bone mineral density (BMD) and micro-architectural deterioration of trabecular bone (TB), increases risk of fractures associated with substantial morbidity, mortality, and financial costs. A quantitative measure of TB micro-architecture with high reproducibility, large between-subjects variability and strong association with bone strength that may be computed via in vivo imaging would be an important indicator of bone quality for clinical trials evaluating fracture risks under different clinical conditions. Previously, the notion of tensor scale (t-scale) was introduced using an ellipsoidal model that yields a unified representation of structure size, orientation and anisotropy. Here, we develop a new 3-D t-scale algorithm for fuzzy objects and investigate its application to compute quantitative measures characterizing TB micro-architecture acquired by in vivo multi-row detector CT (MD-CT) imaging. Specifically, new measures characterizing individual trabeculae on the continuum of a perfect plate and a perfect rod and their orientation are directly computed in a volumetric BMD representation of a TB network. Reproducibility of these measures is evaluated using repeat MD-CT scans and also by comparing their correlation between MD-CT and micro-CT imaging. Experimental results have demonstrated that the t-scale-based TB micro-architectural measures are highly reproducible with strong association of their values at MD-CT and micro-CT resolutions. Results of an experimental mechanical study have proved these measures' ability to predict TB's bone strength.

3D shape-dependent thinning method for trabecular bone characterization

PURPOSE

Curve and surface thinning are widely-used skeletonization techniques for modeling objects in three dimensions. In the case of trabecular bone analysis, however, neither curve nor surface thinning is really efficient since the internal geometry of the object is usually composed of both rod and plate shapes. The purpose of this paper is to propose an original method called hybrid skeleton which better matches the geometry of the data compared to curve and surface skeletons. In the hybrid skeleton algorithm, 1D curves represent rod-shaped zones whereas 2D surfaces represent plate-shaped elements.

METHODS

The proposed hybrid skeleton algorithm is based on a combination of three methods. (1) A new variant of the method proposed by Bonnassie et al. for the classification of voxels as belonging to plate-like or rod-like structures, where the medial axis (MA) algorithm is replaced by a fast and connected skeletonization algorithm. In addition, the reversibility of the MA algorithm is replaced by an isotropic region-growth method to spread the rod and plate labels back to the original object. (2) A well chosen surface thinning method applied on the plate voxels set. (3) A well chosen curve skeleton thinning method applied on the rod voxels set. The efficiency and the robustness of the proposed algorithm were evaluated using synthesis test vectors. A clinical study was led on micro-CT (computed tomography) images of two different populations of osteoarthritic and osteoporotic trabecular bone samples. The morphological and topological characteristics of the two populations were evaluated using the proposed hybrid skeleton as well as the classification algorithm.

RESULTS

When evaluated on test vectors and compared to Bonnassie's algorithm, the proposed classification algorithm gives a slightly better rate of classification. The hybrid skeleton preserves the shape information of the processed objects. Interesting morphological and topological features as well as volumetric ones were extracted from the skeleton and from the classified volumes, respectively. The extracted features enable the two populations of osteoarthritic and osteoporotic trabecular bone samples to be distinguished.

CONCLUSIONS

Compared to curve-based or surface-based skeletons, the hybrid skeleton better matches the geometry of the data. Each rod is represented by a one-voxel-thick arc and each plate is represented by a one-voxel-thick surface. The hybrid skeleton as well as the proposed classification algorithm introduce relevant parameters linked to the presence of plates in the trabecular bone data, showing that rods and plates contain independent information about trabeculae. The hybrid skeleton offers a new opportunity for precise studies of porous media such as trabecular bone.