Structural analysis of cortical porosity applied to HR‐pQCT data

PurposeThe investigation of cortical porosity is an important aspect of understanding biological, pathoetiological, and biomechanical processes occurring within the skeleton. With the emergence of HR-pQCT as a noninvasive tool suitable for clinical use, cortical porosity at appendicular sites can be directly visualized in vivo. The aim of this study was to introduce a novel topological analysis of the cortical pore network for HR-pQCT data and determine the influence of resolution on measures of cortical pore network microstructure and topology.MethodsCadaveric radii were scanned using HR-pQCT at two different voxel sizes (41 and 82 μm) and also using μCT at a voxel size of 18 μm. HR-pQCT and μCT image sets were spatially coregistered. Segmentation and quantification of cortical porosity (Ct.Po) and mean pore diameter (Ct.Po.Dm) were achieved using an established extended cortical analysis technique. Topological classification of individual pores was performed using topology-preserving skeletonization and multicolor dilation algorithms. Based on the pore skeleton topological classification, the following parameters were quantified: total number of planar surface-skeleton canals (N.Slabs), tubular curve-skeleton canals (N.Tubes), and junction elements (N.Junctions), mean slab volume (Slab.Vol), mean tube volume (Tube.Vol), mean slab orientation (Slab.θ), mean tube orientation (Tube.θ), N.Slabs/N.Tubes, and integral (total) slab volume/integral tube volume (iSlab.Vol/iTube.Vol). An in vivo reproducibility study was also conducted to assess short-term precision of the topology parameters. Precision error was characterized using root mean square coefficient of variation (RMSCV%).ResultsCorrelations to μCT values for Ct.Po were significant for both the 41 and 82 μm HR-pQCT data (41: r(2) = 0.82, p < 0.001, 82: r(2) = 0.75, p < 0.001). For Ct.Po.Dm, only the 41 μm data were significantly predictive of μCT values (r(2) = 0.72, p < 0.01) Data at both HR-pQCT voxel sizes were strongly predictive of the μCT values for N.Slabs (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.84, p < 0.001), N.Tubes (41: r(2) = 0.94, p < 0.001; 82: r(2) = 0.84, p < 0.001), and N.Junctions (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.78, p < 0.001), though proportional bias was evident in these correlations. Weak correlations were seen for iSlab.Vol/iTube.Vol at both voxel sizes (41: r(2) = 0.52, p < 0.01; 82: r(2) = 0.39, p < 0.05). Slab.Vol was significantly correlated to μCT data at 41 μm (r(2) = 0.60, p < 0.01) but not at 82 μm, while Tube.Vol was significantly correlated at both voxel sizes (41: r(2) = 0.79, p < 0.001; 82: r(2) = 0.68, p < 0.01). In vivo precision error for these parameters ranged from 2.31 to 9.68 RMSCV%.ConclusionsStrong correlations between μCT- and HR-pQCT-derived measurements were found, particularly in HR-pQCT images obtained at 41 μm. These data are in agreement with our previous study investigating the effect of voxel size on standard HR-pQCT metrics of trabecular and cortical microstructure, and extend our previous findings to include topological descriptors of the cortical pore network

The influence of disuse on bone microstructure and mechanics assessed by HR-pQCT.

Numerous clinical cohorts are exposed to reduced skeletal loading and associated bone loss, including surgical patients, stroke and spinal cord injury victims, and women on bed rest during pregnancy. In this context, understanding disuse-related bone loss is critical to developing interventions to prevent fractures and the associated morbidity, mortality, and cost to the health care system. The aim of this pilot study was to use high-resolution peripheral QCT (HR-pQCT) to examine changes in trabecular and cortical microstructure and biomechanics during a period of non weight bearing (WB) and during recovery following return to normal WB. Surgical patients requiring a 6-week non WB period (n=12, 34.8±7.7 yrs) were scanned at the affected and contralateral tibia prior to surgery, after the 6-week non WB period, and 6 and 13 weeks after returning to full WB. At the affected ultradistal tibia, integral vBMD (including both trabecular and cortical compartments) decreased with respect to baseline (-1.2%), trabecular number increased (+5.6%), while trabecular thickness (-5.4%), separation (-4.6%), and heterogeneity (-7.2%) decreased (all p<0.05). Six weeks after return to full WB, trabecular structure measures reverted to baseline levels. In contrast, integral vBMD continued to decrease after 6 (-2.0%, p<0.05) and 13 weeks (-2.5%, p=0.07) of full WB. At the affected distal site, the disuse period resulted in increased porosity (+16.1%, p<0.005), which remained elevated after 6 weeks (+16.8%, p<0.01) and after 13 weeks (+16.2%, p<0.05). A novel topological analysis applied to the distal tibia cortex demonstrated increased number of canals with surface topology ("slabs" +21.7%, p<0.01) and curve topology ("tubes" +15.0%, p<0.05) as well as increased number of canal junctions (+21.4%, p<0.05) following the disuse period. Porosity increased uniformly through increases in both pore size and number. Finite element analysis at the ultradistal tibia showed decreased stiffness and failure load (-2.8% and -2.4%, p<0.01) following non WB. These biomechanical predictions remained depressed following 6 and 13 weeks of full WB. Finite element analysis at the distal site followed similar trends. Our results suggest that detectable microstructural and biomechanical degradation occurs--particularly within the cortical compartment--as a result of non WB and persists following return to normal loading. A better understanding of these microstructural changes and their short- and long-term influence on biomechanics may have clinical relevance in the context of disuse-related fracture prevention

Structural analysis of cortical porosity applied to HR-pQCT data

Purpose: The investigation of cortical porosity is an important aspect of understanding biological, pathoetiological, and biomechanical processes occurring within the skeleton. With the emergence of HR-pQCT as a noninvasive tool suitable for clinical use, cortical porosity at appendicular sites can be directly visualized in vivo. The aim of this study was to introduce a novel topological analysis of the cortical pore network for HR-pQCT data and determine the influence of resolution on measures of cortical pore network microstructure and topology. Methods: Cadaveric radii were scanned using HR-pQCT at two different voxel sizes (41 and 82 μm) and also using μCT at a voxel size of 18 μm. HR-pQCT and μCT image sets were spatially coregistered. Segmentation and quantification of cortical porosity (Ct.Po) and mean pore diameter (Ct.Po.Dm) were achieved using an established extended cortical analysis technique. Topological classification of individual pores was performed using topology-preserving skeletonization and multicolor dilation algorithms. Based on the pore skeleton topological classification, the following parameters were quantified: total number of planar surface-skeleton canals (N.Slabs), tubular curve-skeleton canals (N.Tubes), and junction elements (N.Junctions), mean slab volume (Slab.Vol), mean tube volume (Tube.Vol), mean slab orientation (Slab.θ), mean tube orientation (Tube.θ), N.Slabs/N.Tubes, and integral (total) slab volume/integral tube volume (iSlab.Vol/iTube.Vol). An in vivo reproducibility study was also conducted to assess short-term precision of the topology parameters. Precision error was characterized using root mean square coefficient of variation (RMSCV%). Results: Correlations to μCT values for Ct.Po were significant for both the 41 and 82 μm HR-pQCT data (41: r(2) = 0.82, p < 0.001, 82: r(2) = 0.75, p < 0.001). For Ct.Po.Dm, only the 41 μm data were significantly predictive of μCT values (r(2) = 0.72, p < 0.01) Data at both HR-pQCT voxel sizes were strongly predictive of the μCT values for N.Slabs (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.84, p < 0.001), N.Tubes (41: r(2) = 0.94, p < 0.001; 82: r(2) = 0.84, p < 0.001), and N.Junctions (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.78, p < 0.001), though proportional bias was evident in these correlations. Weak correlations were seen for iSlab.Vol/iTube.Vol at both voxel sizes (41: r(2) = 0.52, p < 0.01; 82: r(2) = 0.39, p < 0.05). Slab.Vol was significantly correlated to μCT data at 41 μm (r(2) = 0.60, p < 0.01) but not at 82 μm, while Tube.Vol was significantly correlated at both voxel sizes (41: r(2) = 0.79, p < 0.001; 82: r(2) = 0.68, p < 0.01). In vivo precision error for these parameters ranged from 2.31 to 9.68 RMSCV%. Conclusions: Strong correlations between μCT- and HR-pQCT-derived measurements were found, particularly in HR-pQCT images obtained at 41 μm. These data are in agreement with our previous study investigating the effect of voxel size on standard HR-pQCT metrics of trabecular and cortical microstructure, and extend our previous findings to include topological descriptors of the cortical pore network

The influence of disuse on bone microstructure and mechanics assessed by HR-pQCT

Numerous clinical cohorts are exposed to reduced skeletal loading and associated bone loss, including surgical patients, stroke and spinal cord injury victims, and women on bed rest during pregnancy. In this context, understanding disuse-related bone loss is critical to developing interventions to prevent fractures and the associated morbidity, mortality, and cost to the health care system. The aim of this pilot study was to use high-resolution peripheral QCT (HR-pQCT) to examine changes in trabecular and cortical microstructure and biomechanics during a period of non weight bearing (WB) and during recovery following return to normal WB. Surgical patients requiring a 6-week non-WB period (n = 12, 34.8 ± 7.7 yrs) were scanned at the affected and contralateral tibia prior to surgery, after the 6-week non-WB period, and 6 and 13 weeks after returning to full-WB. At the affected ultradistal tibia, integral vBMD (including both trabecular and cortical compartments) decreased with respect to baseline (−1.2%), trabecular number increased (+5.6%), while trabecular thickness (−5.4%), separation (−4.6%), and heterogeneity (−7.2%) decreased (all p<0.05). Six weeks after return to full-WB, trabecular structure measures reverted to baseline levels. In contrast, integral vBMD continued to decrease after 6 (−2.0%, p < 0.05) and 13 weeks (−2.5%, p = 0.07) of full-WB. At the affected distal site, the disuse period resulted in increased porosity (+16.1%, p < 0.005), which remained elevated after 6 weeks (+16.8%, p < 0.01) and after 13 weeks (+16.2%, p < 0.05). A novel topological analysis applied to the distal tibia cortex demonstrated increased number of canals with surface topology (“slabs” +21.7%, p < 0.01) and curve topology (“tubes” +15.0%, p < 0.05) as well as increased number of canal junctions (+21.4%, p < 0.05) following the disuse period. Porosity increased uniformly through increases in both pore size and number. Finite element analysis at the ultradistal tibia showed decreased stiffness and failure load (−2.8% and −2.4%, p < 0.01) following non-WB. These biomechanical predictions remained depressed following 6 and 13 weeks of full-WB. Finite element analysis at the distal site followed similar trends. Our results suggest that detectable microstructural and biomechanical degradation occurs – particularly within the cortical compartment – as a result of non-WB and persists following return to normal loading. A better understanding of these microstructural changes and their short- and long-term influence on biomechanics may have clinical relevance in the context of disuse-related fracture prevention