Breed and loading history influence in vivo skeletal strain patterns in pre-pubertal female chickens.

The influence of loading history on in vivo strains within a given specie remains poorly understood, and although in vivo strains have been measured at the hindlimb bones of various species, strains engendered during modes of activity other than locomotion are lacking, particularly in non-human species. For commercial egg-laying chickens specifically, there is an interest in understanding their bones' mechanical behaviour, particularly during youth, to develop early interventions to prevent the high incidence of osteoporosis in this population. We measured in vivo mechanical strains at the tibiotarsus midshaft during steady activities (ground, uphill, downhill locomotion) and non-steady activities (perching, jumping, aerial transition landing) in 48 pre-pubescent female (egg-laying) chickens from two breeds that were reared in three different housing systems, allowing varying amounts and types of physical activity. Mechanical strain patterns differed between breeds, and were dependent on the activity performed. Mechanical strains were also affected by rearing environment: chickens that were restricted from performing dynamic load bearing activity due to caged-housing generally exhibited higher mechanical strain levels during steady, but not non-steady activities, compared to chickens with prior dynamic load-bearing activity experience. Among chickens with prior experience of dynamic load bearing activity, those reared in housing systems that allowed more frequent physical activity did not exhibit lower mechanical strains. In all groups, the tibiotarsus was subjected to a loading environment consisting of a combination of axial compression, bending, and torsion, with torsion being the predominant source of strain. Aerial transition landing produced the highest strain levels with unusual strain patterns compared to other activities, suggesting it may produce the strongest anabolic response. These results exemplify how different breeds within a given specie adapt to maintain different patterns of mechanical strains, and how benefits of physical activity in terms of resistance to strain are activity-type dependent and do not necessarily increase with increased physical activity. These findings directly inform controlled loading experiments aimed at studying the bone mechanoresponse in young female chickens and can also be associated to measures of bone morphology and material properties to understand how these features influence bone mechanical properties in vivo.

3D Image Registration Marginally Improves the Precision of HR-pQCT Measurements Compared to Cross-Sectional-Area Registration in Adults With Osteogenesis Imperfecta.

Repositioning error in longitudinal high-resolution peripheral-quantitative computed tomography (HR-pQCT) imaging can lead to different bone volumes being assessed over time. To identify the same bone volumes at each time point, image registration is used. While cross-sectional area image registration corrects axial misalignment, 3D registration additionally corrects rotations. Other registration methods involving matched angle analysis (MA) or boundary transformations (3D-TB) can be used to limit interpolation error in 3D-registering micro-finite-element data. We investigated the effect of different image registration methods on short-term in vivo precision in adults with osteogenesis imperfecta, a collagen-related genetic disorder resulting in low bone mass, impaired quality, and increased fragility. The radii and tibiae of 29 participants were imaged twice on the same day with full repositioning. We compared the precision error of different image registration methods for density, microstructural, and micro-finite-element outcomes with data stratified based on anatomical site, motion status, and scanner generation. Regardless of the stratification, we found that image registration improved precision for total and trabecular bone mineral densities, trabecular and cortical bone mineral contents, area measurements, trabecular bone volume fraction, separation, and heterogeneity, as well as cortical thickness and perimeter. 3D registration marginally outperformed cross-sectional area registration for some outcomes, such as trabecular bone volume fraction and separation. Similarly, precision of micro-finite-element outcomes was improved after image registration, with 3D-TB and MA methods providing greatest improvements. Our regression model confirmed the beneficial effect of image registration on HR-pQCT precision errors, whereas motion had a detrimental effect on precision even after image registration. Collectively, our results indicate that 3D registration is recommended for longitudinal HR-pQCT imaging in adults with osteogenesis imperfecta. Since our precision errors are similar to those of healthy adults, these results can likely be extended to other populations, although future studies are needed to confirm this. © 2022 American Society for Bone and Mineral Research (ASBMR).

Fabric-elasticity relationships of tibial trabecular bone are similar in osteogenesis imperfecta and healthy individuals.

Osteogenesis Imperfecta (OI) is an inherited form of bone fragility characterised by impaired synthesis of type I collagen, altered trabecular bone architecture and reduced bone mass. High resolution peripheral computed tomography (HR-pQCT) is a powerful method to investigate bone morphology at peripheral sites including the weight-bearing distal tibia. The resulting 3D reconstructions can be used as a basis of micro-finite element (FE) or homogenized finite element (hFE) models for bone strength estimation. The hFE scheme uses homogenized local bone volume fraction (BV/TV) and anisotropy information (fabric) to compute healthy bone strength within a reasonable computation time using fabric-elasticity relationships. However, it is unclear if these relationships quantified previously for healthy controls are valid for trabecular bone from OI patients. Thus, the aim of this study is to investigate fabric-elasticity relationships in OI trabecular bone compared to healthy controls. In the present study, the morphology of distal tibiae from 50 adults with OI were compared to 120 healthy controls using second generation HR-pQCT. Six cubic regions of interest (ROIs) were selected per individual in a common anatomical region. A first matching between OI and healthy control group was performed by selecting similar individuals to obtain identical mean and median age and sex distribution. It allowed us to perform a first morphometric analysis and compare the outcome with literature. Then, stiffness tensors of the ROIs were computed using µFE and multiple linear regressions were performed with the Zysset-Curnier orthotropic fabric-elasticity model. An initial fit was performed on both the OI group and the healthy control group using all extracted ROIs. Then, data was filtered according to a fixed threshold for a defined coefficient of variation (CV) assessing ROI heterogeneity and additional linear regressions were performed on these filtered data sets. These full and filtered data were in turn compared with previous results from µCT reconstructions obtained in other anatomical locations. Finally, the ROIs of both groups were matched according to their BV/TV and degree of anisotropy (DA). Linear regressions were performed using these matched data to detect statistical differences between the two groups. Compared to healthy controls, we found the OI samples to have significantly lower BV/TV and trabecular number (Tb.N.), significantly higher CV, trabecular separation (Tb.Sp.) and trabecular separation standard deviation (Tb.Sp.SD), but no differences in trabecular thickness (Tb.Th.). These results are in agreement with previous studies. The stiffnesses of highly heterogeneous ROIs were randomly lower with respect to the fabric-elasticity relationships, which reflects the limit of validity of the computational homogenisation methodology. This limitation does not challenge the fabric-elasticity relationship, which extrapolation to heterogeneous ROIs is probably reasonable but can simply not be evaluated with the employed homogenisation methodology. Moreover, due to their low BV/TV, the potential (unknown) errors on these heterogeneous ROIs would have negligible influence on whole bone stiffness in comparison to homogeneous ROIs which are orders of magnitude stiffer. The filtering of highly heterogeneous ROIs removed these low stiffness ROIs and led to similar correlation coefficients for both OI and healthy groups. Finally, the BV/TV and DA matched data revealed no significant differences in fabric-elasticity parameters between OI and healthy individuals. Moreover, the filtering step did not exclude a particular OI type. Compared to previous studies, the stiffness constants from the 61 µm resolution HR-pQCT ROIs were lower than for the 36 µm resolution µCT ROIs. In conclusion, OI trabecular bone of the distal tibia was shown to be significantly more heterogeneous and have a lower BV/TV than healthy controls. Despite the reduced linear regression parameters found for HR-pQCT images, the fabric-elasticity relationships between OI and healthy individuals are similar when the trabecular bone ROIs are sufficiently homogeneous to perform the computational stiffness analysis. Accordingly, the elastic properties used for FEA of healthy bones are also valid for OI bones.

Multisite longitudinal calibration of HR-pQCT scanners and precision in osteogenesis imperfecta.

BACKGROUND: For high-resolution peripheral quantitative computed tomography (HR-pQCT) to be used in longitudinal multi-center studies to assess disease and treatment effects, data must be aggregated across multiple timepoints and scanners. This requires an understanding of the factors contributing to scanner precision, and multi-scanner cross-calibration procedures, especially for clinical populations with severe phenotypes, like osteogenesis imperfecta (OI). METHODS: To address this, we first evaluated single- and multi-center short- and long-term precision errors of standard HR-pQCT parameters. Two imaging phantoms were circulated among 13 sites (7 XtremeCT and 6 XtremeCT2) and scanned in triplicate at 3 timepoints/site. Additionally, duplicate in vivo radial and tibial scans were acquired in 29 individuals with OI. Secondly, we investigated subject- and scanner-related factors that contribute to precision errors using regression analysis. Thirdly, we proposed a reference site selection criterion for multisite cross-calibration and demonstrated the external validity of phantom-based calibrations. RESULTS: Our results show excellent short-term single-site precision in both phantoms (CV % < 0.5%) and in density, microarchitecture and finite element parameters of OI participants (CV % = 0.75 to 1.2%). In vivo reproducibility significantly improved with (i) cross sectional area image registration versus no registration and (ii) scans with no motion artifacts. While reproducibility was similar across OI subtypes and anatomical sites, XtremeCT2 scanners achieved ~2.5% better precision than XtremeCT for trabecular parameters. Finally, we demonstrate that multisite longitudinal precision errors resulting from inconsistencies between scanners can be partially corrected through scanner cross-calibration. CONCLUSIONS: This study is the first to assess long-term reproducibility and cross-calibration in a study using first and second generation HR-pQCT scanners. The results presented in this context provide timely guidelines for future use of this powerful clinical imaging modality in multi-center longitudinal clinical trials.

Distal radius sections offer accurate and precise estimates of forearm fracture load.

BACKGROUND: Forearm fracture risk can be estimated via factor-of-risk: the ratio of applied impact force to forearm fracture load. Simple techniques are available for estimating impact force associated with a fall; estimating forearm fracture load is more challenging. Our aim was to assess whether failure load estimates of sections of the distal radius (acquired using High-Resolution peripheral Quantitative Computed Tomography and finite element modeling) offer accurate and precise estimates of forearm fracture load. METHODS: We scanned a section of the distal radius of 19 cadaveric forearms (female, mean age 83.7, SD 8.3), and 34 women (75.0, 7.7). Sections were converted to finite element models and failure loads were acquired for different failure criteria. We assessed forearm fracture load using experimental testing simulating a fall on the outstretched hand. We used linear regression to derive relationships between ex vivo forearm fracture load and finite element derived distal radius failure load. We used derived regression coefficients to estimate forearm fracture load, and assessed explained variance and prediction error. We used root-mean-squared coefficients of variation to assess in vivo precision errors of estimated forearm fracture load. FINDINGS: Failure load estimates of sections of the distal radius, used in conjunction with derived regression coefficients, explained 89-90% of the variance in experimentally-measured forearm fracture load with prediction errors <6.8% and precision errors <5.0%. INTERPRETATION: Failure load estimates of distal radius sections can reliably estimate forearm fracture load experienced during a fall. Forearm fracture load estimates can be used to improve factor-of-risk predictions for forearm fracture.

Reliability of Annual Changes and Monitoring Time Intervals for Bone Strength, Size, Density, and Microarchitectural Development at the Distal Radius and Tibia in Children: A 1-Year HR-pQCT Follow-Up.

High-resolution peripheral quantitative computed tomography (HR-pQCT) imaging, together with computational finite element analysis (FEA), offers an attractive, noninvasive tool to quantify bone strength development in pediatric studies. Evidence of annual changes and errors in repeated HR-pQCT measures is limited, and time intervals required to reliably capture changes in children's bone strength or microarchitecture have not yet been defined. Our objectives were: (1) to quantify annual changes in bone strength and microarchitectural properties; (2) to define precision errors for pediatric bone strength outcomes; (3) to characterize annual changes in contrast to pediatric precision errors; and (4) to estimate monitoring time intervals (MTIs) required to reliably characterize bone development at the distal radius and tibia. We obtained distal radius (7% of ulnar length) and tibia (8%) bone properties using HR-pQCT and FEA from 38 follow-up study participants (21 girls) at baseline (mean age 10.6 years, SD 1.7 years) and after 1 year; and from 32 precision study participants (16 girls) at baseline (mean age 11.3 years, SD 1.6 years) and after 1 week. We characterized mean annual changes (paired t tests) contrasted to pediatric precision errors (CV%RMS ) and estimated MTIs. Annual increases in bone strength, total area, cortical thickness, and density ranged between 3.0% and 25.3% and 2.4% and 15.6% at the distal radius and tibia, respectively. Precision errors for all bone strength outcomes were ≤6.8% and ≤5.1% at the distal radius and tibia, respectively, and appeared lower than annual gains in bone strength at both sites. Cortical porosity decreased 19.6% at the distal radius and 6.6% at the distal tibia; these changes exceeded respective precision errors, indicating cortical bone consolidation. MTIs ranged between 0.5 years and infinity at the distal radius and 0.5 and 5.9 years at the distal tibia. Estimated MTIs suggest that pediatric bone strength, cortical bone density, and porosity development can be reliably monitored with annual measurements. © 2019 American Society for Bone and Mineral Research.