ABSTRACT
High-resolution peripheral quantitative computed tomography (HR-pQCT) provides a new opportunity for the noninvasive measurement of bone microarchitecture in patient research of rheumatic disease in joints of the extremities. Rheumatoid arthritis of the hand and osteoarthritis of the knee are highly prevalent and there are no cures, so an improved understanding of the etiology of these diseases, particularly when combining HR-pQCT with other imaging modalities, is important to develop new disease management strategies. Analysis of HR-pQCT data is challenging due to the large size of the datasets, but the field continues to rapidly develop, and there is excellent opportunity to integrate sophisticated analyses for this novel research tool.
Subject(s)
Osteoarthritis, Knee , Tomography, X-Ray Computed , Humans , Tomography, X-Ray Computed/methods , Osteoarthritis, Knee/diagnostic imaging , Rheumatic Diseases/diagnostic imaging , Knee Joint/diagnostic imagingABSTRACT
INTRODUCTION: Traumatic bone marrow lesions (BML) are frequently identified on knee MRI scans in patients following an acute full-thickness, complete ACL tear. BMLs coincide with regions of elevated localized bone loss, and studies suggest these may act as a precursor to the development of post-traumatic osteoarthritis. This study addresses the labour-intensive manual assessment of BMLs by using a 3D U-Net for automated identification and segmentation from MRI scans. METHODS: A multi-task learning approach was used to segment both bone and BML from T2 fat-suppressed (FS) fast spin echo (FSE) MRI sequences for BML assessment. Training and testing utilized datasets from individuals with complete ACL tears, employing a five-fold cross-validation approach and pre-processing involved image intensity normalization and data augmentation. A post-processing algorithm was developed to improve segmentation and remove outliers. Training and testing datasets were acquired from different studies with similar imaging protocol to assess the model's performance robustness across different populations and acquisition conditions. RESULTS: The 3D U-Net model exhibited effectiveness in semantic segmentation, while post-processing enhanced segmentation accuracy and precision through morphological operations. The trained model with post-processing achieved a Dice similarity coefficient (DSC) of 0.75 ± 0.08 (mean ± std) and a precision of 0.87 ± 0.07 for BML segmentation on testing data. Additionally, the trained model with post-processing achieved a DSC of 0.93 ± 0.02 and a precision of 0.92 ± 0.02 for bone segmentation on testing data. This demonstrates the approach's high accuracy for capturing true positives and effectively minimizing false positives in the identification and segmentation of bone structures. CONCLUSION: Automated segmentation methods are a valuable tool for clinicians and researchers, streamlining the assessment of BMLs and allowing for longitudinal assessments. This study presents a model with promising clinical efficacy and provides a quantitative approach for bone-related pathology research and diagnostics.
Subject(s)
Anterior Cruciate Ligament Injuries , Bone Marrow , Deep Learning , Magnetic Resonance Imaging , Humans , Magnetic Resonance Imaging/methods , Anterior Cruciate Ligament Injuries/diagnostic imaging , Bone Marrow/diagnostic imaging , Male , Female , Adult , Image Interpretation, Computer-Assisted/methodsABSTRACT
Standard microarchitectural analysis of bone using micro-computed tomography produces a large number of parameters that quantify the structure of the trabecular network. Analyses that perform statistical tests on many parameters are at elevated risk of making Type I errors. However, when multiple testing correction procedures are applied, the risk of Type II errors is elevated if the parameters being tested are strongly correlated. In this article, we argue that four commonly used trabecular microarchitectural parameters (thickness, separation, number, and bone volume fraction) are interdependent and describe only two independent properties of the trabecular network. We first derive theoretical relationships between the parameters based on their geometric definitions. Then, we analyze these relationships with an aggregated in vivo dataset with 2987 images from 1434 participants and a synthetically generated dataset with 144 images using principal component analysis (PCA) and linear regression analysis. With PCA, when trabecular thickness, separation, number, and bone volume fraction are combined, we find that 92 % to 97 % of the total variance in the data is explained by the first two principal components. With linear regressions, we find high coefficients of determination (0.827-0.994) and fitted coefficients within expected ranges. These findings suggest that to maximize statistical power in future studies, only two of trabecular thickness, separation, number and bone volume fraction should be used for statistical testing.
ABSTRACT
BACKGROUND: Weight bearing computed tomography (WBCT) utilizes cone beam CT technology to provide assessments of lower limb joint structures while they are functionally loaded. Grey-scale values indicative of X-ray attenuation that are output from cone beam CT are challenging to calibrate, and their use for bone mineral density (BMD) measurement remains debatable. To determine whether WBCT can be reliably used for cortical and trabecular BMD assessment, we sought to establish the accuracy of BMD measurements at the knee using modern WBCT by comparing them to measurements from conventional CT. METHODS: A hydroxyapatite phantom with three inserts of varying densities was used to systematically quantify signal uniformity and BMD accuracy across the acquisition volume. We evaluated BMD in vivo (nâ¯=â¯5, female) using synchronous and asynchronous calibration techniques in WBCT and CT. To account for variation in attenuation along the height (z-axis) of acquisition volumes, we tested a height-dependent calibration approach for both WBCT and CT images. RESULTS: Phantom BMD measurement error in WBCT was as high as 15.3% and consistently larger than CT (up to 5.6%). Phantom BMD measures made under synchronous conditions in WBCT improved measurement accuracy by up to 3% but introduced more variability in measured BMD. We found strong correlations (Râ¯=â¯0.96) as well as wide limits of agreement (-324 mgHA/cm3 to 183 mgHA/cm3) from Bland-Altman analysis between WBCT and CT measures in vivo that were not improved by height-dependent calibration. CONCLUSION: Whilst BMD accuracy from WBCT was found to be dependent on apparent density, accuracy was independent of the calibration technique (synchronous or asynchronous) and the location of the measurement site within the field of view. Overall, we found strong correlations between BMD measures from WBCT and CT and in vivo measures to be more accurate in trabecular bone regions. Importantly, WBCT can be used to distinguish between anatomically relevant differences in BMD, however future work is necessary to determine the repeatability and sensitivity of BMD measures in WBCT.
Subject(s)
Bone Density , Cone-Beam Computed Tomography , Phantoms, Imaging , Weight-Bearing , Cone-Beam Computed Tomography/methods , Humans , Female , Calibration , Middle Aged , Knee Joint/diagnostic imaging , Knee Joint/physiology , Reproducibility of Results , Cancellous Bone/diagnostic imaging , Aged , Absorptiometry, Photon/methodsABSTRACT
Older men with high bone turnover have faster bone loss. We assessed the link between the baseline levels of bone turnover markers (BTMs) and the prospectively assessed bone microarchitecture decline in men. In 825 men aged 60-87 yr, we measured the serum osteocalcin (OC), bone alkaline phosphatase (BAP), N-terminal propeptide of type I procollagen (PINP), and C-terminal telopeptide of type I collagen (CTX-I), and urinary total deoxypyridinoline (tDPD). Bone microarchitecture and strength (distal radius and distal tibia) were estimated by high-resolution pQCT (XtremeCT, Scanco Medical) at baseline and then after 4 and 8 yr. Thirty-seven men took medications affecting bone metabolism. Statistical models were adjusted for age and BMI. At the distal radius, the decrease in the total bone mineral density (Tt.BMD), cortical BMD (Ct.BMD), cortical thickness (Ct.Thd), and cortical area (Ct.Ar) and failure load was faster in the highest vs the lowest CTX-I quartile (failure load: -0.94 vs -0.31% yr-1, P < .001). Patterns were similar for distal tibia. At the distal tibia, bone decline (Tt.BMD, Ct.Thd, Ct.Ar, Ct.BMD, and failure load) was faster in the highest vs the lowest tDPD quartile. At each skeletal site, the rate of decrease in Tb.BMD differed between the extreme OC quartiles (P < .001). Men in the highest BAP quartile had a faster loss of Tt.BMD, Tb.BMD, reaction force, and failure load vs the lowest quartile. The link between PINP and bone decline was poor. The BTM score is the sum of the nos. of the quartiles for each BTM. Men in the highest quartile of the score had a faster loss of cortical bone and bone strength vs the lowest quartile. Thus, in the older men followed prospectively for 8 yr, the rate of decline in bone microarchitecture and estimated bone strength was 50%-215% greater in men with high bone turnover (highest quartile, CTX-I above the median) compared to the men with low bone turnover (lowest quartile, CTX-I below the median).
Subject(s)
Bone Density , Bone and Bones , Male , Humans , Aged , Female , Prospective Studies , Bone Remodeling , RadiusSubject(s)
Bone Density , Osteoporosis , Humans , Calibration , Osteoporosis/diagnosis , Tomography, X-Ray Computed/methodsABSTRACT
Fracture risk increases with lower areal bone mineral density (aBMD); however, aBMD-related estimate of risk may decrease with age. This may depend on technical limitations of 2-dimensional (2D) dual energy X-ray absorptiometry (DXA) which are reduced with 3D high-resolution peripheral quantitative computed tomography (HR-pQCT). Our aim was to examine whether the predictive utility of HR-pQCT measures with fracture varies with age. We analyzed associations of HR-pQCT measures at the distal radius and distal tibia with two outcomes: incident fractures and major osteoporotic fractures. We censored follow-up time at first fracture, death, last contact or 8 years after baseline. We estimated hazard ratios (HR) and 95%CI for the association between bone traits and fracture incidence across age quintiles. Among 6835 men and women (ages 40-96) with at least one valid baseline HR-pQCT scan who were followed prospectively for a median of 48.3 months, 681 sustained fractures. After adjustment for confounders, bone parameters at both the radius and tibia were associated with higher fracture risk. The estimated HRs for fracture did not vary significantly across age quintiles for any HR-pQCT parameter measured at either the radius or tibia. In this large cohort, the homogeneity of the associations between the HR-pQCT measures and fracture risk across age groups persisted for all fractures and for major osteoporotic fractures. The patterns were similar regardless of the HR-pQCT measure, the type of fracture, or the statistical models. The stability of the associations between HR-pQCT measures and fracture over a broad age range shows that bone deficits or low volumetric density remain major determinants of fracture risk regardless of age group. The lower risk for fractures across measures of aBMD in older adults in other studies may be related to factors which interfere with DXA but not with HR-pQCT measures.
Subject(s)
Tomography, X-Ray Computed , Humans , Aged , Male , Female , Prospective Studies , Middle Aged , Aged, 80 and over , Risk Factors , Bone Density , Adult , Fractures, Bone/diagnostic imaging , Fractures, Bone/epidemiology , Aging , Radius/diagnostic imaging , Tibia/diagnostic imaging , Tibia/pathologyABSTRACT
INTRODUCTION: The continued development of high-resolution peripheral quantitative computed tomography (HR-pQCT) has led to a second-generation scanner with higher resolution and longer scan region. However, large multicenter prospective cohorts were collected with first-generation HR-pQCT and have been used to develop bone phenotyping and fracture risk prediction (µFRAC) models. This study establishes whether there is sufficient universality of these first-generation trained models for use with second-generation scan data. METHODS: HR-pQCT data were collected for a cohort of 60 individuals, who had been scanned on both first- and second-generation scanners on the same day to establish the universality of the HR-pQCT models. These data were each used as input to first-generation trained bone microarchitecture models for bone phenotyping and fracture risk prediction, and their outputs were compared for each study participant. Reproducibility of the models were assessed using same-day repeat scans obtained from first-generation (n = 37) and second-generation (n = 74) scanners. RESULTS: Across scanner generations, the bone phenotyping model performed with an accuracy of 93.1%. Similarly, the 5-year fracture risk assessment by µFRAC was well correlated with a Pearson's (r) correlation coefficient of r > 0.83 for the three variations of µFRAC (varying inclusion of clinical risk factors, finite element analysis, and dual X-ray absorptiometry). The first-generation reproducibility cohort performed with an accuracy for categorical assignment of 100% (bone phenotyping) and a correlation coefficient of 0.99 (µFRAC), whereas the second-generation reproducibility cohort performed with an accuracy of 96.4% (bone phenotyping) and a correlation coefficient of 0.99 (µFRAC). CONCLUSION: We demonstrated that bone microarchitecture models trained using first-generation scan data generalize well to second-generation scans, performing with a high level of accuracy and reproducibility. Less than 4% of individuals' estimated fracture risk led to a change in treatment threshold, and in general, these dissimilar outcomes using second-generation data tended to be more conservative.
Establishing the universality of first-generation-trained HR-pQCT prediction models on second-generation scan data is important to move the bone microarchitecture field forward. We found that despite the difference in resolutions between the two HR-pQCT generations, models developed with first-generation data generalized well to second-generation systems. This avoids unnecessarily repeating complex studies.
Subject(s)
Fractures, Bone , Phenotype , Tomography, X-Ray Computed , Humans , Female , Male , Risk Assessment , Fractures, Bone/diagnostic imaging , Aged , Middle Aged , Bone and Bones/diagnostic imaging , Adult , Bone DensityABSTRACT
Opportunistic computed tomography (CT) scans, which can assess relevant osteoporotic bones of interest, offer a potential solution for identifying osteoporotic individuals. CT scans usually do not contain calibration phantoms, so internal calibration methods have been developed to create a voxel-specific density calibration that can be used in opportunistic CT. It remains a challenge, however, to account for potential sources of error in internal calibration, such as beam hardening or heterogeneous internal reference tissues. The purpose of this work was to introduce our internal calibration method that accounts for these variations and to estimate error bounds for the bone mineral density (BMD) measurements taken from internally calibrated scans. The error bounds are derived by incorporating a combination of a Monte Carlo simulation and standard error propagation into our previously established internal calibration method. A cohort of 138 clinical abdominal CT scans were calibrated for BMD assessment with a phantom placed in the field of view and used as the ground truth. Our modified internal calibration method provided error bounds on the same images and was tested to contain the ground truth phantom-calibrated BMD. This was repeated using 10 different internal reference tissue combinations to explore how error bounds are affected by the choice of internal tissue referents. We found that the tissue combination of air, skeletal muscle, and cortical bone yielded the most accurate BMD estimates while maintaining error bounds that were sufficiently conservative to account for sources of error such as beam hardening and heterogeneous tissue samples. The mean difference between the phantom BMD and the BMD resulting from the tissue combination of air, skeletal muscle and cortical bone was 2.12 mg/cc (0.06% BMD error) and 1.13 mg/cc (0.02 % BMD error) for the left and right femur, respectively. Providing error bounds for internal calibration provides a method to explore the influence of internal reference tissues and confidence for BMD estimates.
Subject(s)
Bone Density , Tomography, X-Ray Computed , Humans , Calibration , Tomography, X-Ray Computed/methods , Bone Density/physiology , Bone and Bones/diagnostic imaging , Femur , Phantoms, ImagingABSTRACT
INTRODUCTION: High resolution peripheral quantitative computed tomography (HR-pQCT) imaging protocol requires defining where to position the â¼1 cm thick scan along the bone length. Discrepancies between the use of two positioning methods, the relative and fixed offset, may be problematic in the comparison between studies and participants. This study investigated how bone landmarks scale linearly with length and how this scaling affects both positioning methods aimed at providing a consistent anatomical location for scan acquisition. METHODS: Using CT images of the radius (Nâ¯=â¯25) and tibia (Nâ¯=â¯42), 10 anatomical landmarks were selected along the bone length. The location of these landmarks was converted to a percent length along the bone, and the variation in their location was evaluated across the dataset. The absolute location of the HR-pQCT scan position using both offset methods was identified for all bones and converted to a percent length position relative to the HR-pQCT reference line for comparison. A secondary analysis of the location of the scan region specifically within the metaphysis was explored at the tibia. RESULTS: The location of landmarks deviated from a linear relationship across the dataset, with a range of 3.6â¯% at the radius sites, and 4.5â¯% at the tibia sites. The consequent variation of the position of the scan at the radius was 0.6â¯% and 0.3â¯%, and at the tibia 2.4â¯% and 0.5â¯%, for the fixed and relative offset, respectively. The position of the metaphyseal junction with the epiphysis relative to the scan position was poorly correlated to bone length, with R2â¯=â¯0.06 and 0.37, for the fixed and relative offset respectively. CONCLUSION: The variation of the scan position by either method is negated by the intrinsic variation of the bone anatomy with respect both to total bone length as well as the metaphyseal region. Therefore, there is no clear benefit of either offset method. However, the lack of difference due to the inherent variation in the underlying anatomy implies that it is reasonable to compare studies even if they are using different positioning methods.
Subject(s)
Radius , Tomography, X-Ray Computed , Humans , Tomography, X-Ray Computed/methods , Radius/diagnostic imaging , Tibia/diagnostic imaging , Upper Extremity , Epiphyses , Bone DensityABSTRACT
The study objective was to assess bone quality measured by high resolution peripheral quantitative computed tomography (HR-pQCT) in competitive athletes. Medline, EMBASE and Sport Discus were searched through May 2022. Prior to submission, a follow-up database search was performed (January 2023). Studies of competitive athletes using HR-pQCT to assess bone quality were included. Athletes were aged between 14 and 45 years. Data extraction included study design and location (country), skeletal imaging modality and site, bone variables and any additional musculoskeletal-related outcome. Information identifying sports and athletes were also extracted. This review included 14 manuscripts and a total of 928 individuals (male: n=75; female: n=853). Athletes comprised 78% (n=722) of the included individuals and 93% of athletes were female. Assessment scores indicate the studies were good to fair quality. The athletes included in this review can be categorized into three groups: 1) healthy athletes, 2) athletes with compromised menstrual function (e.g., amenorrhoea), and 3) athletes with compromised bone health (e.g., bone stress injuries). When assessing bone quality using HR-pQCT, healthy competitive athletes had denser, stronger and larger bones with better microarchitecture, compared with controls. However, the same cannot be said for athletes with amenorrhoea or bone stress injuries.
Subject(s)
Amenorrhea , Bone Density , Humans , Male , Female , Adolescent , Young Adult , Adult , Middle Aged , Bone and Bones/diagnostic imaging , Tomography, X-Ray Computed/methods , Athletes , RadiusABSTRACT
A method was previously developed to identify participant-specific parameters in a model of trabecular bone adaptation from longitudinal computed tomography (CT) imaging. In this study, we use these numerical methods to estimate changes in astronaut bone health during the distinct phases of spaceflight and recovery on Earth. Astronauts (N = 16) received high-resolution peripheral quantitative CT (HR-pQCT) scans of their distal tibia prior to launch (L), upon their return from an approximately six-month stay on the international space station (R+0), and after six (R+6) and 12 (R+12) months of recovery. To model trabecular bone adaptation, we determined participant-specific parameters at each time interval and estimated their bone structure at R+0, R+6, and R+12. To assess the fit of our model to this population, we compared static and dynamic bone morphometry as well as the Dice coefficient and symmetric distance at each measurement. In general, modeled and observed static morphometry were highly correlated (R2> 0.94) and statistically different (p < 0.0001) but with errors close to HR-pQCT precision limits. Dynamic morphometry, which captures rates of bone adaptation, was poorly estimated by our model (p < 0.0001). The Dice coefficient and symmetric distance indicated a reasonable local fit between observed and predicted bone volumes. This work applies a general and versatile computational framework to test bone adaptation models. Future work can explore and test increasingly sophisticated models (e.g., those including load or physiological factors) on a participant-specific basis.
ABSTRACT
Relative Energy Deficiency in Sport (REDs) is a syndrome describing the relationship between prolonged and/or severe low energy availability and negative health and performance outcomes. The high energy expenditures incurred during training and competition put endurance athletes at risk of REDs. The objective of this study was to investigate differences in bone quality in winter endurance athletes classified as either low-risk versus at-risk for REDs. Forty-four participants were recruited (M = 18; F = 26). Bone quality was assessed at the distal radius and tibia using high resolution peripheral quantitative computed tomography (HR-pQCT), and at the hip and spine using dual X-ray absorptiometry (DXA). Finite element analysis was used to estimate bone strength. Participants were grouped using modified criteria from the REDs Clinical Assessment Tool Version 1. Fourteen participants (M = 3; F = 11), were classified as at-risk of REDs (≥ 3 risk factors). Measured with HR-pQCT, cortical bone area (radius) and bone strength (radius and tibia) were 6.8%, 13.1% and 10.3% lower (p = 0.025, p = 0.033, p = 0.027) respectively, in at-risk compared with low-risk participants. Using DXA, femoral neck areal bone density was 9.4% lower in at-risk compared with low-risk participants (p = 0.005). At-risk male participants had 21.9% lower femoral neck areal bone density (via DXA) than low-risk males (p = 0.020) with no significant differences in females. Overall, 33.3% of athletes were at-risk for REDs and had lower bone quality than those at low-risk.
Subject(s)
Bone Density , Bone and Bones , Female , Humans , Male , Cross-Sectional Studies , Absorptiometry, Photon , Tibia , Risk Factors , Radius/diagnostic imaging , AthletesABSTRACT
Computed tomography (CT) offers advanced biomedical imaging of the body and is broadly utilized for clinical diagnosis. Traditionally, clinical CT scans have not been used for volumetric bone mineral density (vBMD) assessment; however, computational advances can now leverage clinically obtained CT data for the secondary analysis of bone, known as opportunistic CT analysis. Initial applications focused on using clinically acquired CT scans for secondary osteoporosis screening, but opportunistic CT analysis can also be applied to answer research questions related to vBMD changes in response to various disease states. There are several considerations for opportunistic CT analysis, including scan acquisition, contrast enhancement, the internal calibration technique, and bone segmentation, but there remains no consensus on applying these methods. These factors may influence vBMD measures and therefore the robustness of the opportunistic CT analysis. Further research and standardization efforts are needed to establish a consensus and optimize the application of opportunistic CT analysis for accurate and reliable assessment of vBMD in clinical and research settings. This review summarizes the current state of opportunistic CT analysis, highlighting its potential and addressing the associated challenges.
ABSTRACT
In a 36-month randomized controlled trial examining the effect of high-dose vitamin D3 on radial and tibial total bone mineral density (TtBMD), measured by high-resolution peripheral quantitative tomography (HR-pQCT), participants (311 healthy males and females aged 55-70 years with dual-energy X-ray absorptiometry T-scores > -2.5 without vitamin D deficiency) were randomized to receive 400 IU (N = 109), 4000 IU (N = 100), or 10,000 IU (N = 102) daily. Participants had HR-pQCT radius and tibia scans and blood sampling at baseline, 6, 12, 24, and 36 months. This secondary analysis examined the effect of vitamin D dose on plasma measurements of the vitamin D metabolome by liquid chromatography-tandem mass spectrometry (LC-MS/MS), exploring whether the observed decline in TtBMD was associated with changes in four key metabolites [25-(OH)D3 ; 24,25-(OH)2 D3 ; 1,25-(OH)2 D3 ; and 1,24,25-(OH)3 D3 ]. The relationship between peak values in vitamin D metabolites and changes in TtBMD over 36 months was assessed using linear regression, controlling for sex. Increasing vitamin D dose was associated with a marked increase in 25-(OH)D3 , 24,25-(OH)2 D3 and 1,24,25-(OH)3 D3 , but no dose-related change in plasma 1,25-(OH)2 D3 was observed. There was a significant negative slope for radius TtBMD and 1,24,25-(OH)3 D3 (-0.05, 95% confidence interval [CI] -0.08, -0.03, p < 0.001) after controlling for sex. A significant interaction between TtBMD and sex was seen for 25-(OH)D3 (female: -0.01, 95% CI -0.12, -0.07; male: -0.04, 95% CI -0.06, -0.01, p = 0.001) and 24,25-(OH)2 D3 (female: -0.75, 95% CI -0.98, -0.52; male: -0.35, 95% CI -0.59, -0.11, p < 0.001). For the tibia there was a significant negative slope for 25-(OH)D3 (-0.03, 95% CI -0.05, -0.01, p < 0.001), 24,25-(OH)2 D3 (-0.30, 95% CI -0.44, -0.16, p < 0.001), and 1,24,25-(OH)3 D3 (-0.03, 95% CI -0.05, -0.01, p = 0.01) after controlling for sex. These results suggest vitamin D metabolites other than 1,25-(OH)2 D3 may be responsible for the bone loss seen in the Calgary Vitamin D Study. Although plasma 1,25-(OH)2 D3 did not change with vitamin D dose, it is possible rapid catabolism to 1,24,25-(OH)3 D3 prevented the detection of a dose-related rise in plasma 1,25-(OH)2 D3 . © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).
Subject(s)
Bone Diseases, Metabolic , Vitamin D , Male , Female , Humans , Chromatography, Liquid , Tandem Mass Spectrometry , Cholecalciferol/pharmacology , Ergocalciferols , Metabolome , Dietary SupplementsABSTRACT
BACKGROUND: Data-driven development of medical biomarkers of bone requires a large amount of image data but physical measurements are generally too restricted in size and quality to perform a robust training. PURPOSE: This study aims to provide a reliable in silico method for the generation of realistic bone microstructure with defined microarchitectural properties. Synthetic bone samples may improve training of neural networks and serve for the development of new diagnostic parameters of bone architecture and mineralization. METHODS: One hundred-fifty cadaveric lumbar vertebrae from 48 different male human spines were scanned with a high resolution peripheral quantitative CT. After prepocessing the scans, we extracted 10,795 purely spongeous bone patches, each with a side length of 32 voxels (5 mm) and isotropic voxel size of 164 µm. We trained a volumetric generative adversarial network (GAN) in a progressive manner to create synthetic microstructural bone samples. We then added a style transfer technique to allow the generation of synthetic samples with defined microstructure and gestalt by simultaneously optimizing two entangled loss functions. Reliability testing was performed by comparing real and synthetic bone samples on 10 well-understood microstructural parameters. RESULTS: The method was able to create synthetic bone samples with visual and quantitative properties that effectively matched with the real samples. The GAN contained a well-formed latent space allowing to smoothly morph bone samples by their microstructural parameters, visual appearance or both. Optimum performance has been obtained for bone samples with voxel size 32 × 32 × 32, but also samples of size 64 × 64 × 64 could be synthesized. CONCLUSIONS: Our two-step-approach combines a parameter-agnostic GAN with a parameter-specific style transfer technique. It allows to generate an unlimited anonymous database of microstructural bone samples with sufficient realism to be used for the development of new data-driven methods of bone-biomarkers. Particularly, the style transfer technique can generate datasets of bone samples with specific conditions to simulate certain bone pathologies.
Subject(s)
Bone and Bones , Tomography, X-Ray Computed , Humans , Male , Reproducibility of Results , Tomography, X-Ray Computed/methods , Bone and Bones/diagnostic imaging , Neural Networks, Computer , Biomarkers , Image Processing, Computer-Assisted/methodsABSTRACT
Most fracture risk assessment tools use clinical risk factors combined with bone mineral density (BMD) to improve assessment of osteoporosis; however, stratifying fracture risk remains challenging. This study developed a fracture risk assessment tool that uses information about volumetric bone density and three-dimensional structure, obtained using high-resolution peripheral quantitative compute tomography (HR-pQCT), to provide an alternative approach for patient-specific assessment of fracture risk. Using an international prospective cohort of older adults (n = 6802) we developed a tool to predict osteoporotic fracture risk, called µFRAC. The model was constructed using random survival forests, and input predictors included HR-pQCT parameters summarizing BMD and microarchitecture alongside clinical risk factors (sex, age, height, weight, and prior adulthood fracture) and femoral neck areal BMD (FN aBMD). The performance of µFRAC was compared to the Fracture Risk Assessment Tool (FRAX) and a reference model built using FN aBMD and clinical covariates. µFRAC was predictive of osteoporotic fracture (c-index = 0.673, p < 0.001), modestly outperforming FRAX and FN aBMD models (c-index = 0.617 and 0.636, respectively). Removal of FN aBMD and all clinical risk factors, except age, from µFRAC did not significantly impact its performance when estimating 5-year and 10-year fracture risk. The performance of µFRAC improved when only major osteoporotic fractures were considered (c-index = 0.733, p < 0.001). We developed a personalized fracture risk assessment tool based on HR-pQCT that may provide an alternative approach to current clinical methods by leveraging direct measures of bone density and structure. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).
Subject(s)
Osteoporotic Fractures , Humans , Aged , Adult , Osteoporotic Fractures/diagnostic imaging , Prospective Studies , Tomography, X-Ray Computed , Bone Density , Risk AssessmentABSTRACT
Rapid loss of areal bone mineral density (aBMD) is associated with higher fracture risk after adjustment for confounders including initial aBMD. However, the link between bone microarchitecture decline and fracture is not clear. We studied the association between bone microarchitecture deterioration assessed prospectively over 4 years and the subsequent fracture risk in older men. Bone microarchitecture at the distal radius and tibia was assessed by high-resolution peripheral QCT (HR-pQCT; XtremeCT, Scanco Medical) (baseline, 4 years) in 732 men aged 60-87 years. During the 8-year follow-up, 109 men had fragility fractures. Areal BMD was assessed by dual-energy X-ray absorptiometry. After adjustment for age, weight, prior falls and fractures, distal radius aBMD (baseline, slope), and baseline distal radius total volumetric BMD (Tt.BMD), a faster decrease in distal radius Tt.BMD was associated with higher fracture risk (hazard ratio [HR] = 1.54/SD, 95% confidence interval: 1.20-1.95, p < .005). Rapid cortical bone loss was associated with higher fracture risk (cortical thickness: HR = 1.48; 1.15-1.90, p < .01; cortical BMD: HR = 1.38; 1.11-1.72, p < .01). The rate of trabecular bone loss at the distal radius and the rate of bone microarchitecture decline at the distal tibia were not associated with fracture risk. After adjustment for aBMD and distal radius HR-pQCT measures assessed after 4 years, changes in Tt.BMD were associated with higher fracture risk (e.g., Tt.BMD: HR = 1.37; 1.11-1.69, p < .005). Compared with the reference model (age, weight, prior fractures and falls, baseline and slope of aBMD, baseline HR-pQCT value), further addition of the slope of the HR-pQCT measure did not improve the fracture prediction. Thus, rapid cortical bone loss at the distal radius is associated with higher fracture risk in the multivariable models including baseline values of the HR-pQCT measure. However, repeated HR-pQCT measurements did not improve the assessment of the fracture risk in older men (compared with the reference model defined earlier). © 2023 American Society for Bone and Mineral Research (ASBMR).
Subject(s)
Bone Diseases, Metabolic , Fractures, Bone , Male , Humans , Aged , Radius/diagnostic imaging , Tomography, X-Ray Computed , Fractures, Bone/diagnostic imaging , Fractures, Bone/epidemiology , Bone and Bones , Bone Density , Cortical Bone/diagnostic imaging , Absorptiometry, Photon , Tibia/diagnostic imagingABSTRACT
High-resolution peripheral quantitative computed tomography (HR-pQCT) is an emerging in vivo imaging modality for quantification of bone microarchitecture. However, extraction of quantitative microarchitectural parameters from HR-pQCT images requires an accurate segmentation of the image. The current standard protocol using semi-automated contouring for HR-pQCT image segmentation is laborious, introduces inter-operator biases into research data, and poses a barrier to streamlined clinical implementation. In this work, we propose and validate a fully automated algorithm for segmentation of HR-pQCT radius and tibia images. A multi-slice 2D U-Net produces initial segmentation predictions, which are post-processed via a sequence of traditional morphological image filters. The U-Net was trained on a large dataset containing 1822 images from 896 unique participants. Predicted segmentations were compared to reference segmentations on a disjoint dataset containing 386 images from 190 unique participants, and 156 pairs of repeated images were used to compare the precision of the novel and current protocols. The agreement of morphological parameters obtained using the predicted segmentation relative to the reference standard was excellent (R2 between 0.938 and > 0.999). Precision was significantly improved for several outputs, most notably cortical porosity. This novel and robust algorithm for automated segmentation will increase the feasibility of using HR-pQCT in research and clinical settings.
Subject(s)
Bone and Bones , Tomography, X-Ray Computed , Humans , Tomography, X-Ray Computed/methods , Tibia/diagnostic imaging , Radius , Upper Extremity , Bone DensityABSTRACT
Prolonged use of antiresorptives such as the bisphosphonate alendronate (ALN) and the RANKL inhibitor denosumab (DMAb) are associated with rare cases of atypical femoral fracture (AFF). The etiology of AFF is unclear, but it has been hypothesized that potent osteoclast inhibitors may reduce bone fatigue resistance. The purpose of this study was to quantify the relationship between antiresorptive treatment and fatigue life (cycles to failure) in bone from ovariectomized cynomolgus monkeys. We analyzed humeral bone from 30 animals across five treatment groups. Animals were treated for 12 months with subcutaneous (sc) vehicle (VEH), sc DMAb (25 mg/kg/month), or intravenous (iv) ALN (50 µg/kg/month). Another group received 6 months VEH followed by 6 months DMAb (VEH-DMAb), and the final group received 6 months ALN followed by 6 months DMAb (ALN-DMAb). A total of 240 cortical beam samples were cyclically tested in four-point bending at 80, 100, 120, or 140 MPa peak stress. High-resolution imaging and density measurements were performed to evaluate bone microstructure and composition. Samples from the ALN (p = 0.014), ALN-DMAb (p = 0.008), and DMAb (p < 0.001) groups illustrated higher fatigue-life measurements than VEH. For example, at 140 MPa the VEH group demonstrated a median ± interquartile range (IQR) fatigue life of 1987 ± 10593 cycles, while animals in the ALN, ALN-DMAb, and DMAb groups survived 9850 ± 13648 (+395% versus VEH), 10493 ± 16796 (+428%), and 14495 ± 49299 (+629%) cycles, respectively. All antiresorptive treatment groups demonstrated lower porosity, smaller pore size, greater pore spacing, and lower number of canals versus VEH (p < 0.001). Antiresorptive treatment was also associated with greater apparent density, dry density, and ash density (p ≤ 0.03). We did not detect detrimental changes following antiresorptive treatments that would explain their association with AFF. In contrast, 12 months of treatment may have a protective effect against fatigue fractures. © 2022 American Society for Bone and Mineral Research (ASBMR).