Non-invasive bone competence analysis by high-resolution pQCT: an in vitro reproducibility study on structural and mechanical properties at the human radius

What Is the Diagnostic Performance of Conventional Radiographs and Clinical Reassessment Compared With HR-pQCT Scaphoid Fracture Diagnosis?

BACKGROUND

Conventional radiographs and clinical reassessment are considered guides in managing clinically suspected scaphoid fractures. This is a unique study as it assessed the value of conventional radiographs and clinical reassessment in a cohort of patients, all of whom underwent additional imaging, regardless of the outcome of conventional radiographs and clinical reassessment.

QUESTIONS/PURPOSES

(1) What is the diagnostic performance of conventional radiographs in patients with a clinically suspected scaphoid fracture compared with high-resolution peripheral quantitative CT (HR-pQCT)? (2) What is the diagnostic performance of clinical reassessment in patients with a clinically suspected scaphoid fracture compared with HR-pQCT? (3) What is the diagnostic performance of conventional radiographs and clinical reassessment combined compared with HR-pQCT?

METHODS

Between December 2017 and October 2018, 162 patients with a clinically suspected scaphoid fracture presented to the emergency department (ED). Forty-six patients were excluded and another 25 were not willing or able to participate, which resulted in 91 included patients. All patients underwent conventional radiography in the ED and clinical reassessment 7 to 14 days later, together with CT and HR-pQCT. The diagnostic performance characteristics and accuracy of conventional radiographs and clinical reassessment were compared with those of HR-pQCT for the diagnosis of fractures since this was proven to be superior to CT scaphoid fracture detection. The cohort included 45 men and 46 women with a median (IQR) age of 52 years (29 to 67). Twenty-four patients with a median age of 44 years (35 to 65) were diagnosed with a scaphoid fracture on HR-pQCT.

RESULTS

When compared with HR-pQCT, conventional radiographs alone had a sensitivity of 67% (95% CI 45% to 84%), specificity of 85% (95% CI 74% to 93%), positive predictive value (PPV) of 62% (95% CI 46% to 75%), negative predictive value (NPV) of 88% (95% CI 80% to 93%), and a positive and negative likelihood ratio (LR) of 4.5 (95% CI 2.4 to 8.5) and 0.4 (95% CI 0.2 to 0.7), respectively. Compared with HR-pQCT, clinical reassessment alone resulted in a sensitivity of 58% (95% CI 37% to 78%), specificity of 42% (95% CI 30% to 54%), PPV of 26% (95% CI 19% to 35%), NPV of 74% (95% CI 62% to 83%), as well as a positive and negative LR of 1.0 (95% CI 0.7 to 1.5) and 1.0 (95% CI 0.6 to 1.7), respectively. Combining clinical examination with conventional radiography produced a sensitivity of 50% (95% CI 29% to 71%), specificity of 91% (95% CI 82% to 97%), PPV of 67% (95% CI 46% to 83%), NPV of 84% (95% CI 77% to 88%), as well as a positive and negative LR of 5.6 (95% CI 2.4 to 13.2) and 0.6 (95% CI 0.4 to 0.8), respectively.

CONCLUSION

The accuracy of conventional radiographs (80% compared with HR-pQCT) and clinical reassessment (46% compared with HR-pQCT) indicate that the value of clinical reassessment is limited in diagnosing scaphoid fractures and cannot be considered directive in managing scaphoid fractures. The combination of conventional radiographs and clinical reassessment does not increase the accuracy of these diagnostic tests compared with the accuracy of conventional radiographs alone and is therefore also limited in diagnosing scaphoid fractures.

LEVEL OF EVIDENCE

Level II, diagnostic study.

Hind limb muscles influence the architectural properties of long bones in frogs

The Mechanostat Theory states that osteocytes sense both the intensity and directionality of the strains induced by mechanical usage and modulate the bone design accordingly. In long bones, this process may adapt anterior-posterior and lateral-medial strength to their mechanical environment showing regional specificity. Anuran species are ideal for analyzing the muscle-bone relationships related to the different mechanical stresses induced by their many locomotor modes and habitat uses. This work aimed to explore the relationships between indicators of the force of the most relevant muscles to locomotion and the mechanical properties of femur and tibia fibula in preserved samples of three anuran species with different habitat use (aquatic, arboreal) and locomotion modes (swimmer, jumper, walker/climber). For that purpose, we measured the anatomical cross-sectional area of each dissected muscle and correlated it with the moments of inertia and bone strength indices. Significant, species-specific covariations between muscle and bone parameters were observed. Pseudis platensis, the aquatic swimmer, showed the largest muscles, followed by Boana faber, the jumper and Phyllomedusa sauvagii, the walker/climber. As we expected, bigger muscles correlate with bone parameters in all the species. Nevertheless, smaller muscles also play an important role in bone design. In aquatic species, muscle interaction enhances mostly lateral bending strength throughout the femur and lateral and antero-posterior bending strength in the tibia fibula. In the jumper species, muscles affected the femur and tibia fibula mostly in anterior-posterior bending. In the walker/climber species, responses involving both antero-posterior and lateral bending strengths were observed in the femur and tibia fibula. These results show that bones will be more or less resistant to lateral and antero-posterior bending according to the different mechanical challenges of locomotion in aquatic vs. arboreal habitats. This study provides new evidence of the muscle-bone relationships in three frog species associated with their different locomotion and habitat uses, highlighting the crucial role of muscle in determining the architectural properties of bones.

Lower bone density and microarchitecture alterations in HIV-infected Brazilian men aged 50 years and older are associated with estradiol levels

OBJECTIVE

Combination antiretroviral treatment (cART) allows for longer survival for people living with HIV and hence long-term complications of both disease and treatment are common. Our purpose was to evaluate bone alterations in men living with HIV (MLWH) and receiving cART and to identify associated factors that can be corrected or mitigated.

PATIENTS AND DESIGN

Thirty MLWH and 36 healthy controls (≥50 years) were studied for areal bone mineral density (aBMD) and body composition (dual-energy X-ray absorptiometry), volumetric bone mineral density (vBMD) and bone microstructure (high-resolution peripheral quantitative computed tomography [HR-pQCT]), serum calcium, phosphate, parathyroid hormone, 25(OH)D, testosterone (T), estradiol (E₂ ), glucose, creatinine, and albumin levels.

RESULTS

The proportion of patients classified as osteoporosis (according to the lowest aBMD T-score) was higher among MLWH as compared to controls (17.9% vs. 5.9%, p = .011). The MLWH showed significant alterations in cortical and trabecular bone on HR-pQCT, which were not associated with the duration of HIV infection or cART. These differences in vBMD and bone microstructure seen in HR-pQCT persisted in the nonosteoporotic MLWH as compared to nonosteoporotic control subjects. Body mass index (BMI) and fat mass were lower in MLWH and positively associated with total vBMD, cortical bone area, and thickness. E₂ and E₂ /T ratios were lower in MLWH than in controls and significantly correlated with several cortical and trabecular bone parameters. Multivariate regression analysis entering simultaneously age, BMI, and E₂ defined that E₂ is an independent influence on bone parameters evaluated by HR-pQCT.

CONCLUSION

MLWH have alterations in bone volumetric density and microstructure when compared with controls, irrespective of aBMD, which are associated with lower E₂ and BMI.

A computational framework for canonical holistic morphometric analysis of trabecular bone

Bone is a remarkable, living tissue that functionally adapts to external loading. Therefore, bone shape and internal structure carry information relevant to many disciplines, including medicine, forensic science, and anthropology. However, morphometric comparisons of homologous regions across different individuals or groups are still challenging. In this study, two methods were combined to quantify such differences: (1) Holistic morphometric analysis (HMA) was used to quantify morphometric values in each bone, (2) which could then be mapped to a volumetric mesh of a canonical bone created by a statistical free-form deformation model (SDM). Required parameters for this canonical holistic morphometric analysis (cHMA) method were identified and the robustness of the method was evaluated. The robustness studies showed that the SDM converged after one to two iterations, had only a marginal bias towards the chosen starting image, and could handle large shape differences seen in bones of different species. Case studies were performed on metacarpal bones and proximal femora of different primate species to confirm prior study results. The differences between species could be visualised and statistically analysed in both case studies. cHMA provides a framework for performing quantitative comparisons of different morphometric quantities across individuals or groups. These comparisons facilitate investigation of the relationship between spatial morphometric variations and function or pathology, or both.

Influence of loading conditions in finite element analysis assessed by HR-pQCT on ex vivo fracture prediction

Many fractures occur in individuals with normal areal Bone Mineral Density (aBMD) measured by Dual X-ray Absorptiometry (DXA). High Resolution peripheral Quantitative Computed Tomography (HR-pQCT) allows for non-invasive evaluation of bone stiffness and strength through micro finite element (μFE) analysis at the tibia and radius. These μFE outcomes are strongly associated with fragility fractures but do not provide clear enhancement compared with DXA measurements. The objective of this study was to establish whether a change in loading conditions in standard μFE analysis assessed by HR-pQCT enhance the discrimination of low-trauma fractured radii (n = 11) from non-fractured radii (n = 16) obtained experimentally throughout a mechanical test reproducing a forward fall. Micro finite element models were created using HR-pQCT images, and linear analyses were performed using four different types of loading conditions (axial, non-axial with two orientations and torsion). No significant differences were found between the failure load assessed with the axial and non-axial models. The different loading conditions tested presented the same area under the receiver operating characteristic (ROC) curves of 0.79 when classifying radius fractures with an accuracy of 81.5%. In comparison, the area under the curve (AUC) is 0.77 from DXA-derived ultra-distal aBMD of the forearm with an accuracy of 85.2%. These results suggest that the restricted HR-pQCT scanned region seems not sensitive to loading conditions for the prediction of radius fracture risk based on ex vivo experiments (n = 27).

Automated segmentation of fractured distal radii by 3D geodesic active contouring of in vivo HR-pQCT images

Radius fractures are among the most common fracture types; however, there is limited consensus on the standard of care. A better understanding of the fracture healing process could help to shape future treatment protocols and thus improve functional outcomes of patients. High-resolution peripheral quantitative computed tomography (HR-pQCT) allows monitoring and evaluation of the radius on the micro-structural level, which is crucial to our understanding of fracture healing. However, current radius fracture studies using HR-pQCT are limited by the lack of automated contouring routines, hence only including small number of patients due to the prohibitively time-consuming task of manually contouring HR-pQCT images. In the present study, a new method to automatically contour images of distal radius fractures based on 3D morphological geodesic active contours (3D-GAC) is presented. Contours of 60 HR-pQCT images of fractured and conservatively treated radii spanning the healing process up to one year post-fracture are compared to the current gold standard, hand-drawn 2D contours, to assess the accuracy of the algorithm. Furthermore, robustness was established by applying the algorithm to HR-pQCT images of intact radii of 73 patients and comparing the resulting morphometric indices to the gold standard patient evaluation including a threshold- and dilation-based contouring approach. Reproducibility was evaluated using repeat scans of intact radii of 19 patients. The new 3D-GAC approach offers contours within inter-operator variability for images of fractured distal radii (mean Dice score of 0.992 ± 0.005 versus median operator Dice score of 0.992 ± 0.006). The generated contours for images of intact radii yielded morphometric indices within the in vivo reproducibility limits compared to the current gold standard. Additionally, the 3D-GAC approach shows an improved robustness against failure (n = 5) when dealing with cortical interruptions, fracture fragments, etc. compared with the automatic, default manufacturer pipeline (n = 40). Using the 3D-GAC approach assures consistent results, while reducing the need for time-consuming hand-contouring.

Quantitative measures of bone shape, cartilage morphometry and joint alignment are associated with disease in an ACLT and MMx rat model of osteoarthritis

Multi-scale, subject-specific quantitative methods to characterize and monitor osteoarthritis in animal models and therapeutic treatments could help reveal causal relationships in disease development and distinguish treatment strategies. In this work, we demonstrate a reproducible and sensitive quantitative image analysis to characterize bone, cartilage and joint measures describing a rat model of post-traumatic osteoarthritis. Eleven 3-month-old male Wistar rats underwent medial anterior cruciate ligament (ACL) transection and medial meniscectomy on the right knee to destabilise the right tibiofemoral joint. They were sacrificed 6 weeks post-surgery and a silicon-based micro-bead contrast agent was injected in the joint space, before scanning with micro-computed tomography (microCT). Subsequently, 3D quantitative morphometric analysis (QMA), previously developed for rabbit joints, was performed. This included cartilage, subchondral cortical and epiphyseal bone measures, as well as novel tibiofemoral joint metrics. Semi-quantitative evaluation was performed on matching two-dimensional (2D) histology and microCT images. Reproducibility of the QMA was tested on eleven age-matched additional joints. The results indicate the QMA method is accurate and reproducible and that microCT-derived cartilage measurements are valid for the analysis of rat joints. The pathologic changes caused by transection of the ACL and medial meniscectomy were reflected in measurements of bone shape, cartilage morphology, and joint alignment. Furthermore, we were able to identify model-specific predictive parameters based on morphometric parameters measured with the QMA.

3D computational anatomy of the scaphoid and its waist for use in fracture treatment

BACKGROUND

A detailed understanding of scaphoid anatomy helps anatomic fracture reduction, and optimal screw position. Therefore, we analysed (1) the size and shape variations of the cartilage and osseous surface, (2) the distribution of volumetric bone mineral density (vBMD) and (3) if the vBMD values differ between a peripheral and a central screw pathway?

METHODS

Forty-three fresh frozen hand specimens (17 females, 26 males) were analysed with high-resolution peripheral quantitative computed tomography (HR-pQCT) and dissected to compute a 3D-statistical osseous and cartilage surface model and a 3D-averaged vBMD model of the scaphoid. 3D patterns were analysed using principal component analysis (PCA). vBMD was analysed via averaging HR-pQCT grey values and virtual bone probing along a central and peripheral pathway.

RESULTS

(1) PCA displayed most notable variation in length ranging from 1.7 cm (- 2SD) to 2.6 cm (mean) and 3.7 cm (+ 2SD) associated with differences of the width and configuration of the dorsal surface (curved and narrow (4 mm) to a wider width (9 mm)). (2) High vBMD was located in the peripheral zone. Lowest vBMD was observed in the centre and waist. (3) Virtual probing along a peripheral pathway near to the cartilage surfaces for the capitate and lunate allowed the center region to be bypassed, resulting in increased vBMD compared to a central pathway.

CONCLUSION

High anatomical variations regarding the osseous and cartilage surfaces were associated with three distinct concentrically arranged zones with notable different vBMD. The complex scaphoid anatomy with its waist might alter the strategy of fracture fixation, education and research.

Bone health assessment via digital wrist tomosynthesis in the mammography setting

Bone fractures attributable to osteoporosis are a significant problem. Though preventative treatment options are available for individuals who are at risk of a fracture, a substantial number of these individuals are not identified due to lack of adherence to bone screening recommendations. The issue is further complicated as standard diagnosis of osteoporosis is based on bone mineral density (BMD) derived from dual energy x-ray absorptiometry (DXA), which, while helpful in identifying many at risk, is limited in fully predicting risk of fracture. It is reasonable to expect that bone screening would become more prevalent and efficacious if offered in coordination with digital breast tomosynthesis (DBT) exams, provided that osteoporosis can be assessed using a DBT modality. Therefore, the objective of the current study was to explore the feasibility of using digital tomosynthesis imaging in a mammography setting. To this end, we measured density, cortical thickness and microstructural properties of the wrist bone, correlated these to reference measurements from microcomputed tomography and DXA, demonstrated the application in vivo in a small group of participants, and determined the repeatability of the measurements. We found that measurements from digital wrist tomosynthesis (DWT) imaging with a DBT scanner were highly repeatable ex vivo (error = 0.05%-9.62%) and in vivo (error = 0.06%-10.2%). In ex vivo trials, DWT derived BMDs were strongly correlated with reference measurements (R = 0.841-0.980), as were cortical thickness measured at lateral and medial cortices (R = 0.991 and R = 0.959, respectively) and the majority of microstructural measures (R = 0.736-0.991). The measurements were quick and tolerated by human patients with no discomfort, and appeared to be different between young and old participants in a preliminary comparison. In conclusion, DWT is feasible in a mammography setting, and informative on bone mass, cortical thickness, and microstructural qualities that are known to deteriorate in osteoporosis. To our knowledge, this study represents the first application of DBT for imaging bone. Future clinical studies are needed to further establish the efficacy for diagnosing osteoporosis and predicting risk of fragility fracture using DWT.

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.

In vivo repeatability of homogenized finite element analysis based on multiple HR-pQCT sections for assessment of distal radius and tibia strength

INTRODUCTION

Micro finite element analysis (μFE) is a widely applied tool in biomedical research for assessing in vivo mechanical properties of bone at measurement sites, including the ultra-distal radius and tibia. A finite element approach (hFE) based on homogenized constitutive models for trabecular bone offers an attractive alternative for clinical use, as it is computationally less expensive than traditional μFE. The respective patient-specific models for in vivo bone strength estimation are usually based on standard clinical high-resolution peripheral quantitative CT (HR-pQCT) measurements. They include a scan region of roughly 10 mm in height and are referred to as single-sections. It has been shown, that these small peripheral bone sections don't reliably cover the fracture line in Colles' fractures and therefore the weakest region at the radius. Recently introduced multiple section (multiple adjacent single-sections) measurements might improve the evaluation of bone strength, but little is known about the repeatability of hFE estimations in general, and especially for multiple section measurement protocols. Accordingly, the aim of the present work is to quantify repeatability of clinical in vivo bone strength measurement by hFE on multiple section HR-pQCT reconstructions at the distal radius and tibia.

METHODS

Nineteen healthy Swiss women (43.6y ± 17.8y) and twenty men (48.2y ± 19.4y) were examined with HR-pQCT at 61 μm isotropic voxel resolution. Each subject was first scanned three times using a double-section (336 slices) at the distal radius and then three times using a triple-section (504 slices) at the distal tibia. The multiple section HR-pQCT reconstructions were graded for motion artefacts and non-linear hFE models (radius and tibia) and linear μFE models (only radius) were generated for estimation of stiffness and ultimate load. Then in vivo repeatability errors were computed in terms of root mean square coefficients of variation (CV).

RESULTS

In vivo repeatability errors of non-linear hFE stiffness (S) and ultimate load (F) were significantly higher at the radius (S: 2.71% and F: 2.97%) compared to the tibia (S: 1.21%, F: 1.45%). Multiple section linear μFE at the radius resulted in substantially higher repeatability errors (S: 5.38% and F: 10.80%) compared to hFE.

DISCUSSION/CONCLUSION

Repeatability errors of hFE outcomes based on multiple section measurements at the distal radius and tibia were generally lower compared to respective reported single-section μFE repeatability errors. Therefore, hFE is an attractive alternative to today's gold standard of μFE models and should especially be encouraged when analyzing multiple section measurements.

Best Practices for Conducting Observational Research to Assess the Relation between Nutrition and Bone: An International Working Group Summary

Diet is a modifiable factor that can affect bone strength and integrity, and the risk of fractures. Currently, a hierarchy of scientific evidence contributes to our understanding of the role of diet on bone health and fracture risk. The strength of evidence is generally based on the type of study conducted, the quality of the methodology employed, the rigor and integrity of the data collected and analysis plan, and the transparency and completeness of the results. Randomized controlled trials (RCTs) are considered to be the gold standard from a clinical research paradigm, but there is a dearth of high-quality diet-related intervention trials with bone as the primary outcome, forcing the use of observational research to inform research and clinical practices. However, for observational research to be of the most utility, standardization and optimization of the study design, accurate and reliable measurement of key variables, and appropriate data analysis and data reporting are paramount. Although there have been recommendations made in relation to RCTs in the field of nutrition, no clear rubric exists for best practices in conducting observational research with regard to nutrition and bone health. Therefore, the purpose of this paper is to describe the best practices and considerations for designing, conducting, analyzing, interpreting, and reporting observational research specifically for understanding the role of nutrition in bone health, amassed by a global panel of scientific experts with strengths in bone, nutrition epidemiology, physical activity, public health, clinical and translational trials, and observational study methods. The global panel of scientific experts represents the leadership and selected participants from the 10th annual International Symposium for the Nutritional Aspects of Osteoporosis. The topics selected and best practices presented reflect expert opinion and areas of scientific expertise of the authors rather than a systematic or comprehensive literature review or professional reporting guidelines.

Virtual supersampling as post-processing step preserves the trabecular bone morphometry in human peripheral quantitative computed tomography scans

In the clinical field of diagnosis and monitoring of bone diseases, high-resolution peripheral quantitative computed tomography (HR-pQCT) is an important imaging modality. It provides a resolution where quantitative bone morphometry can be extracted in vivo on patients. It is known that HR-pQCT provides slight differences in morphometric indices compared to the current standard approach micro-computed tomography (micro-CT). The most obvious reason for this is the restriction of the radiation dose and with this a lower image resolution. With advances in micro-CT evaluation techniques such as patient-specific remodeling simulations or dynamic bone morphometry, a higher image resolution would potentially also allow the application of such novel evaluation techniques to clinical HR-pQCT measurements. Virtual supersampling as post-processing step was considered to increase the image resolution of HR-pQCT scans. The hypothesis was that this technique preserves the structural bone morphometry. Supersampling from 82 μm to virtual 41 μm by trilinear interpolation of the grayscale values of 42 human cadaveric forearms resulted in strong correlations of structural parameters (R²: 0.96–1.00). BV/TV was slightly overestimated (4.3%, R²: 1.00) compared to the HR-pQCT resolution. Tb.N was overestimated (7.47%; R²: 0.99) and Tb.Th was slightly underestimated (-4.20%; R²: 0.98). The technique was reproducible with PE_%CV between 1.96% (SMI) and 7.88% (Conn.D). In a clinical setting with 205 human forearms with or without fracture measured at 82 μm resolution HR-pQCT, the technique was sensitive to changes between groups in all parameters (p < 0.05) except trabecular thickness. In conclusion, we demonstrated that supersampling preserves the bone morphometry from HR-pQCT scans and is reproducible and sensitive to changes between groups. Supersampling can be used to investigate on the resolution dependency of HR-pQCT images and gain more insight into this imaging modality.

Computerized anatomy of the distal radius and its relevance to volar plating, research, and teaching

Distal radius fractures are common and fracture patterns and fixation can be complex. Computerized anatomy evaluation (CAE) might offer non-invasive and enhanced anatomy assessment that might help with implant selection and placement and screw length determination. Our goal was to test the accuracy of two CAE methods for anatomical volar plate positioning and screw lengths measurement of the distal radius. We included 56 high-resolution peripheral quantitative computed tomography scans of intact, human distal radii. Plates were placed manually onto 3D printed models (method 1), which was compared with automated computerized plate placement onto the 3D computer models (method 2). Subsequently, screw lengths were determined digitally for both methods. Screw lengths evaluations were compared via Bland-Altman plots. Both CAE methods resulted in identical volar plate selection and in anatomical plate positioning. For screw length the concordance correlation coefficient was ≥0.91, the location shift ≤0.22 mm, and the scale shift ≤0.16. The differences were smaller than ±1 mm in all samples. Both CAE methods allow for comparable plate positioning and subsequent screw length measurement in distal radius volar plating. Both can be used as a non-invasive teaching environment for volar plate fixation. Method 2 even offers fully computerized assessments. Future studies could compare our models to other anatomical areas, post-operative volar plate positioning, and model performance in actual distal radius fracture instead of intact radii. Clin. Anat. 32:361-368, 2019. © 2018 The Authors. Clinical Anatomy published by Wiley Periodicals, Inc. on behalf of American Association of Clinical Anatomists.

Mechanical Metrics of the Proximal Tibia are Precise and Differentiate Osteoarthritic and Normal Knees: A Finite Element Study

Our objective was to identify precise mechanical metrics of the proximal tibia which differentiated OA and normal knees. We developed subject-specific FE models for 14 participants (7 OA, 7 normal) who were imaged three times each for assessing precision (repeatability). We assessed various mechanical metrics (minimum principal and von Mises stress and strain as well as structural stiffness) across the proximal tibia for each subject. In vivo precision of these mechanical metrics was assessed using CV%RMS. We performed parametric and non-parametric statistical analyses and determined Cohen's d effect sizes to explore differences between OA and normal knees. For all FE-based mechanical metrics, average CV%RMS was less than 6%. Minimum principal stress was, on average, 75% higher in OA versus normal knees while minimum principal strain values did not differ. No difference was observed in structural stiffness. FE modeling could precisely quantify and differentiate mechanical metrics variations in normal and OA knees, in vivo. This study suggests that bone stress patterns may be important for understanding OA pathogenesis at the knee.

The role of imaging in early diagnosis and prevention of joint damage in inflammatory arthritis

INTRODUCTION

Inflammatory arthritis is characterized by chronic inflammation in the synovium, associated with degradation of cartilage and erosion of juxta-articular bone. The bone loss and joint destruction mediated by aberrant immunological responses resulting in proinflammatory cytokine release and various immune cell activation are known as osteoimmunology. Areas covered: A structured literature search including Medline and PubMed, Cochrane meta-analyses and abstracts of international congresses was performed to review joint damage in inflammatory arthritis in terms of pathogenesis, novel imaging assessment, and prevention. Expert commentary: Deeper understanding of the integration of the skeletal and immune as well as inflammatory system is paving the way to prevent bone loss and bone destruction in inflammatory arthritis. With the availability of various imaging modalities such as ultrasound, magnetic resonance imaging (MRI) and high-resolution peripheral quantitative computed tomography (HR-pQCT), we are now able to detect early joint damage, early diagnosis of inflammatory arthritis, monitor the progression or even ascertain whether the inflammatory process is effectively suppressed to allow repair of joint damage by novel therapeutic agents.

Imaging of metabolic bone disease

Osteoporosis is the most important metabolic bone disease, with a wide distribution among the elderly. It is characterized by low bone mass and micro architectural deterioration of bone tissue, leading to enhanced bone fragility and a consequent increase in fracture risk. Identify bone weakening with an appropriate and accurate use of diagnostic imaging is of critical importance in the diagnosis and follow-up of osteoporotic patients. The aim of this review is to evaluate the detection rates of the different imaging modalities in the evaluation of bone strength, in the assessment of fracture risk and in the management of fragility fractures. (www.actabiomedica.it)

The SPECTRA Collaboration OMERACT Special Interest Group: Current Research and Future Directions

OBJECTIVE

High-resolution peripheral quantitative computed tomography (HR-pQCT) has the potential to improve radiographic progression determination in clinical trials and longitudinal observational studies. The goal of this work was to describe the current state of research presented at Outcome Measures in Rheumatology (OMERACT) 2016 and ensuing future directions outlined during discussion among attendees.

METHODS

At OMERACT 2016, SPECTRA (Study grouP for xtrEme-Computed Tomography in Rheumatoid Arthritis) introduced efforts to (1) validate the HR-pQCT according to OMERACT guidelines, focusing on rheumatoid arthritis (RA), and (2) find alternatives for automated joint space width (JSW) analysis. The Special Interest Group (SIG) was presented to patient research partners, physicians/researchers, and SIG leaders followed by a 40-min discussion on future directions.

RESULTS

A consensus definition for RA erosion using HR-pQCT was demonstrated through a systematic literature review and a Delphi exercise. Histopathology and perfusion studies were presented that analyzed the true characteristics of cortical breaks in HR-pQCT images, and to provide criterion validity. Results indicate that readers were able to discriminate between erosion and small vascular channels. Moderate reliability (ICC 0.206-0.871) of direct erosion size measures was shown, which improved (> 0.9) only when experienced readers were considered. Quantification of erosion size was presented for scoring, direct measurement, and volumetric approaches, as well as a reliability exercise for direct measurement. Three methods for JSW measurement were compared, all indicating excellent reproducibility with differences at the extremes (i.e., near-zero and joint edge thickness).

CONCLUSION

Initial reports on HR-pQCT are promising; however, to consider its use in clinical trials and longitudinal observational studies, it is imperative to assess the responsiveness of erosion measurement quantification.

Fast estimation of Colles' fracture load of the distal section of the radius by homogenized finite element analysis based on HR-pQCT

Fractures of the distal section of the radius (Colles' fractures) occur earlier in life than other osteoporotic fractures. Therefore, they can be interpreted as a warning signal for later, more deleterious fractures of vertebral bodies or the femoral neck. In the past decade, the advent of HR-pQCT allowed a detailed architectural analysis of the distal radius and an automated but time-consuming estimation of its strength with linear micro-finite element (μFE) analysis. Recently, a second generation of HR-pQCT scanner (XtremeCT II, SCANCO Medical, Switzerland) with a resolution beyond 61 μm became available for even more refined biomechanical investigations in vivo. This raises the question how biomechanical outcome variables compare between the original (LR) and the new (HR) scanner resolution. Accordingly, the aim of this work was to validate experimentally a patient-specific homogenized finite element (hFE) analysis of the distal section of the human radius for the fast prediction of Colles' fracture load based on the last generation HR-pQCT. Fourteen pairs of fresh frozen forearms (mean age = 77.5±9) were scanned intact using the high (61 μm) and the low (82 μm) resolution protocols that correspond to the new and original HR-pQCT systems. From each forearm, the 20mm most distal section of the radius were dissected out, scanned with μCT at 16.4 μm and tested experimentally under compression up to failure for assessment of stiffness and ultimate load. Linear and nonlinear hFE models together with linear micro finite element (μFE) models were then generated based on the μCT and HR-pQCT reconstructions to predict the aforementioned mechanical properties of 24 sections. Precision errors of the short term reproducibility of the FE analyses were measured based on the repeated scans of 12 sections. The calculated failure loads correlated strongly with those measured in the experiments: accounting for donor as a random factor, the nonlinear hFE provided a marginal coefficient of determination (R_m²) of 0.957 for the high resolution (HR) and 0.948 for the low resolution (LR) protocols, the linear hFE with R_m² of 0.957 for the HR and 0.947 for the LR protocols. Linear μFE predictions of the ultimate load were similar with an R_m² of 0.950 for the HR and 0.954 for the LR protocols, respectively. Nonlinear hFE strength computation led to precision errors of 2.2 and 2.3% which were higher than the ones calculated based on the linear hFE (1.6 and 1.9%) and linear μFE (1.2 and 1.6%) for the HR and LR protocols respectively. Computation of the fracture load with nonlinear hFE demanded in average 6h of CPU time which was 3 times faster than with linear μFE, while computation with linear hFE took only a few minutes. This study delivers an extensive experimental and numerical validation for the application of an accurate and fast hFE diagnostic tool to help in identifying individuals who may be at risk of an osteoporotic wrist fracture and to follow up pharmacological and other treatments in such patients.

Validation of a new multiscale finite element analysis approach at the distal radius

High-resolution peripheral computed tomography is commonly used to evaluate mechanical behavior of the distal radius microstructure using micro-finite element (FE) modeling. However, only a 9mm section is considered and boundary conditions (BCs) are usually simplified (platen-compression), and may not represent physiologic loading. Regardless, these methods are increasingly being used for clinical evaluations. Our goal was to develop and validate a novel multiscale solution that allows for physiologically relevant loading simulations (such as bracing during a fall), and show that mechanical behavior in the distal radius is different under platen BCs. Our approach incorporated bone microstructure together with organ-level radius geometry, by replacing matching continuum regions with micro-FE sections in user-defined regions of interest. Multiscale model predicted strains showed a strong correlation and a significant relationship with measured strains (r=0.836, p<0.001; slope=0.881, intercept=-12.17 µε, p<0.001). Interestingly, platen BC simulated strains were almost 50% lower than measured strains (r=0.835, p<0.001), and strain distributions were clearly different. Our multiscale method demonstrated excellent potential as a computationally efficient alternative for observing true mechanical environment within distal radius microstructure under physiologically accurate loading.

Bone Mass Distribution of the Distal Tibia in Normal, Osteopenic, and Osteoporotic Conditions: An Ex Vivo Assessment Using HR-pQCT, DXA, and Computational Modelling

Osteoporosis leads to bone loss and structural deterioration, which increase the risk of fractures. The aim of this study was to characterize the three-dimensional (3D) bone mass distributions of the distal tibia in normal, osteopenic, and osteoporotic conditions. High-resolution peripheral quantitative computed tomography (HR-pQCT) of the 33 % of the distal tibia and local dual-energy X-ray absorptiometry were applied to 53 intact, fresh-frozen tibiae. The HR-pQCTs were graded to assign local T-scores and merged into three equally sized average normal, osteopenic, and osteoporotic surface models. Volumetric bone mineral density (vBMD) was determined using categorized T-scores, volumetric visualization, and virtual bore probes at the dia-, meta-, and epiphyseal sites (T-DIA, T-META, and T-EPI). We observed a distinct 3D bone mass distribution that was gradually uninfluenced by T-score categories. T-DIA was characterized by the lowest bone mass located in the medullary cavity and a wide homogenous cortex containing the maximum vBMD. The T-META showed decreased cortical thickness and maximal vBMD. At the T-EPI, the relatively low vBMD of the mostly trabecular bone was similar to the maximal cortical vBMD in this sub-region. Four trabecular regions of low bone mass were identified in the recesses. The bone content gradually decreased at all sites, whereas the pattern of bone mass distribution remained essentially unchanged, with the exception of disproportionate losses at T-DIA, T-META, and T-EPI that consistently showed increased endocortical, intracortical, and trabecular bone loss. Extra information can be obtained from the specific pattern of bone mass distribution, potential disproportionate bone losses, and method used.

A STUDY ON THE EVALUATION OF THE PLGA-CURCUMIN TREATED RAT MODELS FOR OSTEOPOROSIS IMPROVEMENT USING FINITE ELEMENT ANALYSIS

Osteoporosis induced bone fracture is not limited to simple damage according to its fractured region and the age of the damaged person. Therefore, there has been more necessity to understand and study the biomechanical aspects of osteoporosis, and relevant research has actively been conducted. With the development of medical imaging technology, the reconstruction of CT images into a 3D model has been made possible, and it has been conducted to analyze 3-dimensionally and inhomogeneously distributed bone mineral density according to CT values. The 3D bone strength computational analysis based on 2D bone images utilized for evaluating the new treatment using the mixture of PLGA and curcumin for rat osteoporosis model. After the treatment, bone mineral density recovered by 65% compared to the osteoporosis model, moreover, bone strength improved by 113%. This study suggested that the limitations of conventional bone monitoring method using 2D bone images can be overcome by considering 3D bone computational analysis. It is considered that introducing the newly proposed Deformation Per Volume (DPW) and Stress Per Weight (SPW) can be also used as effective parameters to more accurately predict the change in bone stiffness.

Voxel size dependency, reproducibility and sensitivity of an in vivo bone loading estimation algorithm

A bone loading estimation algorithm was previously developed that provides in vivo loading conditions required for in vivo bone remodelling simulations. The algorithm derives a bone's loading history from its microstructure as assessed by high-resolution (HR) computed tomography (CT). This reverse engineering approach showed accurate and realistic results based on micro-CT and HR-peripheral quantitative CT images. However, its voxel size dependency, reproducibility and sensitivity still need to be investigated, which is the purpose of this study. Voxel size dependency was tested on cadaveric distal radii with micro-CT images scanned at 25 µm and downscaled to 50, 61, 75, 82, 100, 125 and 150 µm. Reproducibility was calculated with repeated in vitro as well as in vivo HR-pQCT measurements at 82 µm. Sensitivity was defined using HR-pQCT images from women with fracture versus non-fracture, and low versus high bone volume fraction, expecting similar and different loading histories, respectively. Our results indicate that the algorithm is voxel size independent within an average (maximum) error of 8.2% (32.9%) at 61 µm, but that the dependency increases considerably at voxel sizes bigger than 82 µm. In vitro and in vivo reproducibility are up to 4.5% and 10.2%, respectively, which is comparable to other in vitro studies and slightly higher than in other in vivo studies. Subjects with different bone volume fraction were clearly distinguished but not subjects with and without fracture. This is in agreement with bone adapting to customary loading but not to fall loads. We conclude that the in vivo bone loading estimation algorithm provides reproducible, sensitive and fairly voxel size independent results at up to 82 µm, but that smaller voxel sizes would be advantageous.

In vivo precision of three HR-pQCT-derived finite element models of the distal radius and tibia in postmenopausal women

Background

The distal radius is the most common osteoporotic fracture site occurring in postmenopausal women. Finite element (FE) modeling is a non-invasive mathematical technique that can estimate bone strength using inputted geometry/micro-architecture and tissue material properties from computed tomographic images. Our first objective was to define and compare in vivo precision errors for three high-resolution peripheral quantitative computed tomography (HR-pQCT, XtremeCT; Scanco) based FE models of the distal radius and tibia in postmenopausal women. Our second objective was to assess the role of scan interval, scan quality, and common region on precision errors of outcomes for each FE model.

Methods

Models included: single-tissue model (STM), cortical-trabecular dual-tissue model (DTM), and one scaled model using imaged bone mineral density (E-BMD). Using HR-pQCT, we scanned the distal radius and tibia of 34 postmenopausal women (74 ± 7 years), at two time points. Primary outcomes included: tissue stiffness, apparent modulus, average von Mises stress, and failure load. Precision errors (root-mean-squared coefficient of variation, CV%_RMS) were calculated. Multivariate ANOVA was used to compare the mean of individual CV% among the 3 HR-pQCT-based FE models. Spearman correlations were used to characterize the associations between precision errors of all FE model outcomes and scan/time interval, scan quality, and common region. Significance was accepted at P < 0.05.

Results

At the distal radius, CV%_RMS precision errors were <9 % (Range STM: 2.8–5.3 %; DTM: 2.9–5.4 %; E-BMD: 4.4–8.7 %). At the distal tibia, CV%_RMS precision errors were <6 % (Range STM: 2.7–4.8 %; DTM: 2.9–3.8 %; E-BMD: 1.8–2.5 %). At the radius, Spearman correlations indicated associations between the common region and associated precision errors of the E-BMD-derived apparent modulus (ρ = −0.392; P < 0.001) and von Mises stress (ρ = −0.297; P = 0.007).

Conclusion

Results suggest that the STM and DTM are more precise for modeling apparent modulus, average von Mises stress, and failure load at the distal radius. Precision errors were comparable for all three models at the distal tibia. Results indicate that the noted differences in precision error at the distal radius were associated with the common scan region, illustrating the importance of participant repositioning within the cast and reference line placement in the scout view during the scanning process.

Microarchitecture and Peripheral BMD are Impaired in Postmenopausal White Women With Fracture Independently of Total Hip T-Score: An International Multicenter Study

Because single-center studies have reported conflicting associations between microarchitecture and fracture prevalence, we included high-resolution peripheral quantitative computed tomography (HR-pQCT) data from five centers worldwide into a large multicenter analysis of postmenopausal women with and without fracture. Volumetric BMD (vBMD) and microarchitecture were assessed at the distal radius and tibia in 1379 white postmenopausal women (age 67 ± 8 years); 470 (34%) had at least one fracture including 349 with a major fragility fracture. Age, height, weight, and total hip T-score differed across centers and were employed as covariates in analyses. Women with fracture had higher BMI, were older, and had lower total hip T-score, but lumbar spine T-score was similar between groups. At the radius, total and trabecular vBMD and cortical thickness were significantly lower in fractured women in three out of five centers, and trabecular number in two centers. Similar results were found at the tibia. When data from five centers were combined, however, women with fracture had significantly lower total, trabecular, and cortical vBMD (2% to 7%), lower trabecular number (4% to 5%), and thinner cortices (5% to 6%) than women without fracture after adjustment for covariates. Results were similar at the radius and tibia. Similar results were observed with analysis restricted to major fragility fracture, vertebral and hip fractures, and peripheral fracture (at the radius). When focusing on osteopenic women, each SD decrease of total and trabecular vBMD was associated with a significantly increased risk of major fragility fracture (OR = 1.55 to 1.88, p < 0.01) after adjustment for covariates. Moreover, trabecular architecture modestly improved fracture discrimination beyond peripheral total vBMD. In conclusion, we observed differences by center in the magnitude of fracture/nonfracture differences at both the distal radius and tibia. However, when data were pooled across centers and the sample size increased, we observed significant and consistent deficits in vBMD and microarchitecture independent of total hip T-score in all postmenopausal white women with fracture and in the subgroup of osteopenic women, compared to women who never had a fracture. © 2016 American Society for Bone and Mineral Research.

Magnetic resonance imaging of the ear for patient-specific reconstructive surgery

Introduction

Like a fingerprint, ear shape is a unique personal feature that should be reconstructed with a high fidelity during reconstructive surgery. Ear cartilage tissue engineering (TE) advantageously offers the possibility to use novel 3D manufacturing techniques to reconstruct the ear, thus allowing for a detailed auricular shape. However it also requires detailed patient-specific images of the 3D cartilage structures of the patient’s intact contralateral ear (if available). Therefore the aim of this study was to develop and evaluate an imaging strategy for acquiring patient-specific ear cartilage shape, with sufficient precision and accuracy for use in a clinical setting.

Methods and Materials

Magnetic resonance imaging (MRI) was performed on 14 volunteer and six cadaveric auricles and manually segmented. Reproducibility of cartilage volume (Cg.V), surface (Cg.S) and thickness (Cg.Th) was assessed, to determine whether raters could repeatedly define the same volume of interest. Additionally, six cadaveric auricles were harvested, scanned and segmented using the same procedure, then dissected and scanned using high resolution micro-CT. Correlation between MR and micro-CT measurements was assessed to determine accuracy.

Results

Good inter- and intra-rater reproducibility was observed (precision errors <4% for Cg.S and <9% for Cg.V and Cg.Th). Intraclass correlations were good for Cg.V and Cg.S (>0.82), but low for Cg.Th (<0.23) due to similar average Cg.Th between patients. However Pearson’s coefficients showed that the ability to detect local cartilage shape variations is unaffected. Good correlation between clinical MRI and micro-CT (r>0.95) demonstrated high accuracy.

Discussion and Conclusion

This study demonstrated that precision and accuracy of the proposed method was high enough to detect patient-specific variation in ear cartilage geometry. The present study provides a clinical strategy to access the necessary information required for the production of 3D ear scaffolds for TE purposes, including detailed patient-specific shape. Furthermore, the protocol is applicable in daily clinical practice with existing infrastructure.

Assessment of the healing process in distal radius fractures by high resolution peripheral quantitative computed tomography

In clinical practice, fracture healing is evaluated by clinical judgment in combination with conventional radiography. Due to limited resolution, radiographs don't provide detailed information regarding the bone micro-architecture and bone strength. Recently, assessment of in vivo bone density, architectural and mechanical properties at the microscale became possible using high resolution peripheral quantitative computed tomography (HR-pQCT) in combination with micro finite element analysis (μFEA). So far, such techniques have been used mainly to study intact bone. The aim of this study was to explore whether these techniques can also be used to assess changes in bone density, micro-architecture and bone stiffness during fracture healing. Therefore, the fracture region in eighteen women, aged 50 years or older with a stable distal radius fracture, was scanned using HR-pQCT at 1-2 (baseline), 3-4, 6-8 and 12weeks post-fracture. At 1-2 and 12 weeks post-fracture the distal radius at the contra-lateral side was also scanned as control. Standard bone density, micro-architectural and geometric parameters were calculated and bone stiffness in compression, torsion and bending was assessed using μFEA. A linear mixed effect model with time post-fracture as fixed effect was used to detect significant (p-value ≤0.05) changes from baseline. Wrist pain and function were scored using the patient-rated wrist evaluation (PRWE) questionnaire. Correlations between the bone parameters and the PRWE score were calculated by Spearman's correlation coefficient. At the fracture site, total and trabecular bone density increased by 11% and 20%, respectively, at 6-8 weeks, whereas cortical density was decreased by 4%. Trabecular thickness increased by 23-31% at 6-8 and 12 weeks and the intertrabecular area became blurred, indicating intertrabecular bone formation. Compared to baseline, calculated bone stiffness in compression, torsion and bending was increased by 31% after 12 weeks. A moderate negative correlation was found between the stiffness and the PRWE score. No changes were observed at the contra-lateral side. The results demonstrate that it is feasible to assess clinically relevant and significant longitudinal changes in bone density, micro-architecture and mechanical properties at the fracture region during the healing process of stable distal radius fractures using HR-pQCT.

Challenges in longitudinal measurements with HR-pQCT: evaluation of a 3D registration method to improve bone microarchitecture and strength measurement reproducibility

Definition of identical regions between repeated computed tomography (CT) scans is a key factor to monitor changes in bone microarchitecture. In longitudinal studies, accurate determination of the volume of interest (VOI), using three dimensional (3D) registration may improve precision. Therefore, the aim of our study was to investigate the short-term reproducibility of bone geometry, density, microstructure and biomechanical parameters assessed by HR-pQCT and micro-finite element (μFE) derived analyses, using the cross-sectional area (CSA) registration method in comparison with the use of 3D registration, to find overlapping regions between scans. Fifteen healthy individuals (aged 21-47 years) underwent 3 separate scans at the distal radius and tibia, within a one-month interval. Reproducibility was assessed after double contouring the cortical compartment and after applying three different methods to determine the common region between repeated scans: (i) the VOI was determined with no registration, i.e., on 110 slices, (ii) the VOI was determined after CSA-based registration, and (iii) the VOI was determined after 3D registration. Both pre- and post-registration short-term reproducibility for each subject was determined. With no registration, CVrms of geometry parameters ranged from 0.5 to 3.7%, showing a slight variation in the CSA between scans. When the CSA registration method was employed, the variability of geometry (CVrms<1.8%) and density parameters (CVrms<1.8%), was better than that obtained without registration. By removing the effect of repositioning, the 3D registration further improved the reproducibility of cortical bone measurements compared to other methods. Indeed, significant improvements were found for cortical geometry and microstructure measurements (CVrms ranged from 0.4% to 10.7% at both sites; p<0.05), whereas the impact on trabecular bone measurements was restricted to its geometry parameter. The repositioning error was significantly reduced, most markedly at the radius compared to the tibia. For μFE measures, the impact of 3D registration on whole bone stiffness was negligible, indicating adequate assessment of longitudinal changes in estimated biomechanical properties, even without registration. In conclusion, we have shown that the 3D registration improved the identification of the common region retained for longitudinal analysis, contributing to improve the reproducibility of cortical bone parameter measurements. We also quantified the minimally detectable bone changes to help designing future studies with HR-pQCT.

Change of mechanical vertebrae properties due to progressive osteoporosis: combined biomechanical and finite-element analysis within a rat model

For assessing mechanical properties of osteoporotic bone, biomechanical testing combined with in silico modeling plays a key role. The present study focuses on microscopic mechanical bone properties in a rat model of postmenopausal osteoporosis. Female Sprague-Dawley rats were (1) euthanized without prior interventions, (2) sham-operated, and (3) subjected to ovariectomy combined with a multi-deficiencies diet. Rat vertebrae (corpora vertebrae) were imaged by micro-CT, their stiffness was determined by compression tests, and load-induced stress states as well as property changes due to the treatment were analyzed by finite-element modeling. By comparing vertebra stiffness measurements with finite-element calculations of stiffness, an overall microscopic Young's modulus of the bone was determined. Macroscopic vertebra stiffness as well as the microscopic modulus diminish with progression of osteoporosis by about 70 %. After strong initial changes of bone morphology, further decrease in macroscopic stiffness is largely due to decreasing microscopic Young's modulus. The micromechanical stress calculations reveal particularly loaded vertebra regions prone to failure. Osteoporosis-induced changes of the microscopic Young's modulus alter the fracture behavior of bone, may influence bone remodeling, and should be considered in the design of implant materials.

A customized protocol to assess bone quality in the metacarpal head, metacarpal shaft and distal radius: a high resolution peripheral quantitative computed tomography precision study

Background

High Resolution-Peripheral Quantitative Computed Tomography (HR-pQCT) is an emerging technology for evaluation of bone quality in Rheumatoid Arthritis (RA). However, there are limitations with standard HR-pQCT imaging protocols for examination of regions of bone commonly affected in RA. We developed a customized protocol for evaluation of volumetric bone mineral density (vBMD) and microstructure at the metacarpal head (MH), metacarpal shaft (MS) and ultra-ultra-distal (UUD) radius; three sites commonly affected in RA. The purpose was to evaluate short-term measurement precision for bone density and microstructure at these sites.

Methods

12 non-RA participants, individuals likely to have no pre-existing bone damage, consented to participate [8 females, aged 23 to 71 y [median (IQR): 44 (28) y]. The custom protocol includes more comfortable/stable positioning and adapted cortical segmentation and direct transformation analysis methods. Dominant arm MH, MS and UUD radius scans were completed on day one; repeated twice (with repositioning) three to seven days later. Short-term precision for repeated measures was explored using intraclass correlational coefficient (ICC), mean coefficient of variation (CV%), root mean square coefficient of variation (RMSCV%) and least significant change (LSC%₉₅).

Results

Bone density and microstructure precision was excellent: ICCs varied from 0.88 (MH₂ trabecular number) to .99 (MS₃ polar moment of inertia); CV% varied from < 1 (MS₂ vBMD) to 6 (MS₃ marrow space diameter); RMSCV% varied from < 1 (MH₂ full bone vBMD) to 7 (MS₃ marrow space diameter); and LSC% ₉₅varied from 2 (MS₂ full bone vBMD to 21 (MS₃ marrow space diameter). Cortical porosity measures were the exception; RMSCV% varying from 19 (MS₃) to 42 (UUD). No scans were stopped for discomfort. 5% (5/104) were repeated due to motion during imaging. 8% (8/104) of final images had motion artifact graded > 3 on 5 point scale.

Conclusion

In our facility, this custom protocol extends the potential for in vivo HR-pQCT imaging to assess, with high precision, regional differences in bone quality at three sites commonly affected in RA. Our methods are easy to adopt and we recommend other users of HR-pQCT consider this protocol for further evaluations of its precision and feasibility in their imaging facilities.

Improved fracture risk assessment based on nonlinear micro-finite element simulations from HRpQCT images at the distal radius

More accurate techniques to estimate fracture risk could help reduce the burden of fractures in postmenopausal women. Although micro-finite element (µFE) simulations allow a direct assessment of bone mechanical performance, in this first clinical study we investigated whether the additional information obtained using geometrically and materially nonlinear µFE simulations allows a better discrimination between fracture cases and controls. We used patient data and high-resolution peripheral quantitative computed tomography (HRpQCT) measurements from our previous clinical study on fracture risk, which compared 100 postmenopausal women with a distal forearm fracture to 105 controls. Analyzing these data with the nonlinear µFE simulations, the odds ratio (OR) for the factor-of-risk (yield load divided by the expected fall load) was marginally higher (1.99; 95% confidence interval [CI], 1.41-2.77) than for the factor-of-risk computed from linear µFE (1.89; 95% CI, 1.37-2.69). The yield load and the energy absorbed up to the yield point as computed from nonlinear µFE were highly correlated with the initial stiffness (R(2)  = 0.97 and 0.94, respectively) and could therefore be derived from linear simulations with little loss in precision. However, yield deformation was not related to any other measurement performed and was itself a good predictor of fracture risk (OR, 1.89; 95% CI, 1.39-2.63). Moreover, a combined risk score integrating information on relative bone strength (yield load-based factor-of-risk), bone ductility (yield deformation), and the structural integrity of the bone under critical loads (cortical plastic volume) improved the separation of cases and controls by one-third (OR, 2.66; 95% CI, 1.84-4.02). We therefore conclude that nonlinear µFE simulations provide important additional information on the risk of distal forearm fractures not accessible from linear µFE nor from other techniques assessing bone microstructure, density, or mass.

Clinical imaging of bone microarchitecture with HR-pQCT

Osteoporosis, a disease characterized by loss of bone mass and structural deterioration, is currently diagnosed by dual-energy x-ray absorptiometry (DXA). However, DXA does not provide information about bone microstructure, which is a key determinant of bone strength. Recent advances in imaging permit the assessment of bone microstructure in vivo using high-resolution peripheral quantitative computed tomography (HR-pQCT). From these data, novel image processing techniques can be applied to characterize bone quality and strength. To date, most HR-pQCT studies are cross-sectional comparing subjects with and without fracture. These studies have shown that HR-pQCT is capable of discriminating fracture status independent of DXA. Recent longitudinal studies present new challenges in terms of analyzing the same region of interest and multisite calibrations. Careful application of analysis techniques and educated clinical interpretation of HR-pQCT results have improved our understanding of various bone-related diseases and will no doubt continue to do so in the future.

Computational identification and quantification of trabecular microarchitecture classes by 3-D texture analysis-based clustering

High resolution peripheral quantitative computed tomography (HR-pQCT) permits the non-invasive assessment of cortical and trabecular bone density, geometry, and microarchitecture. Although researchers have developed various post-processing algorithms to quantify HR-pQCT image properties, few of these techniques capture image features beyond global structure-based metrics. While 3D-texture analysis is a key approach in computer vision, it has been utilized only infrequently in HR-pQCT research. Motivated by high isotropic spatial resolution and the information density provided by HR-pQCT scans, we have developed and evaluated a post-processing algorithm that quantifies microarchitecture characteristics via texture features in HR-pQCT scans. During a training phase in which clustering was applied to texture features extracted from each voxel of trabecular bone, three distinct clusters, or trabecular microarchitecture classes (TMACs) were identified. These TMACs represent trabecular bone regions with common texture characteristics. The TMACs were then used to automatically segment the voxels of new data into three regions corresponding to the trained cluster features. Regional trabecular bone texture was described by the histogram of relative trabecular bone volume covered by each cluster. We evaluated the intra-scanner and inter-scanner reproducibility by assessing the precision errors (PE), intra class correlation coefficients (ICC) and Dice coefficients (DC) of the method on 14 ultradistal radius samples scanned on two HR-pQCT systems. DC showed good reproducibility in intra-scanner set-up with a mean of 0.870±0.027 (no unit). Even in the inter-scanner set-up the ICC showed high reproducibility, ranging from 0.814 to 0.964. In a preliminary clinical test application, the TMAC histograms appear to be a good indicator, when differentiating between postmenopausal women with (n=18) and without (n=18) prevalent fragility fractures. In conclusion, we could demonstrate that 3D-texture analysis and feature clustering seems to be a promising new HR-pQCT post-processing tool with good reproducibility, even between two different scanners.

A new fracture assessment approach coupling HR-pQCT imaging and fracture mechanics-based finite element modeling

A new fracture assessment approach that combines HR-pQCT imaging with fracture mechanics-based finite element modeling was developed to evaluate distal radius fracture load. Twenty distal radius images obtained from postmenopausal women (fracture, n=10; nonfracture, n=10) were processed to obtain a cortical and a whole bone model for each subject. The geometrical properties of each model were evaluated and the corresponding fracture load was determined under realistic fall conditions using cohesive finite element modeling. The results showed that the whole bone fracture load can be estimated based on the cortical fracture load for nonfracture (R(2)=0.58, p=0.01) and pooled data (R(2)=0.48, p<0.001) but not for the fracture group. The portion of the whole bone fracture load carried by the cortical bone increased with increasing cortical fracture load (R(2)≥0.5, p<0.05) indicating that a more robust cortical bone carries a larger percentage of whole bone fracture load. Cortical thickness was found to be the best predictor of both cortical and whole bone fracture load for all groups (R(2) range: 0.49-0.96, p<0.02) with the exception of fracture group whole bone fracture load showing the predictive capability of cortical geometrical properties in determining whole bone fracture load. Fracture group whole bone fracture load was correlated with trabecular thickness (R(2)=0.4, p<0.05) whereas the nonfracture and the pooled group did not show any correlation with the trabecular parameters. In summary, this study introduced a new modeling approach that coupled HR-pQCT imaging with fracture mechanics-based finite element simulations, incorporated fracture toughness and realistic fall loading conditions in the models, and showed the significant contribution of the cortical compartment to the overall fracture load of bone. Our results provide more insight into the fracture process in bone and may lead to improved fracture load predictions.

Multicenter precision of cortical and trabecular bone quality measures assessed by high-resolution peripheral quantitative computed tomography

High-resolution peripheral quantitative computed tomography (HR-pQCT) has recently been introduced as a clinical research tool for in vivo assessment of bone quality. The utility of this technology to address important skeletal health questions requires translation to standardized multicenter data pools. Our goal was to evaluate the feasibility of pooling data in multicenter HR-pQCT imaging trials. Reproducibility imaging experiments were performed using structure and composition-realistic phantoms constructed from cadaveric radii. Single-center precision was determined by repeat scanning over short-term (<72 hours), intermediate-term (3-5 months), and long-term intervals (28 months). Multicenter precision was determined by imaging the phantoms at nine different HR-pQCT centers. Least significant change (LSC) and root mean squared coefficient of variation (RMSCV) for each interval and across centers was calculated for bone density, geometry, microstructure, and biomechanical parameters. Single-center short-term RMSCVs were <1% for all parameters except cortical thickness (Ct.Th) (1.1%), spatial variability in cortical thickness (Ct.Th.SD) (2.6%), standard deviation of trabecular separation (Tb.Sp.SD) (1.8%), and porosity measures (6% to 8%). Intermediate-term RMSCVs were generally not statistically different from short-term values. Long-term variability was significantly greater for all density measures (0.7% to 2.0%; p < 0.05 versus short-term) and several structure measures: cortical thickness (Ct.Th) (3.4%; p < 0.01 versus short-term), cortical porosity (Ct.Po) (15.4%; p < 0.01 versus short-term), and trabecular thickness (Tb.Th) (2.2%; p < 0.01 versus short-term). Multicenter RMSCVs were also significantly higher than short-term values: 2% to 4% for density and micro-finite element analysis (µFE) measures (p < 0.0001), 2.6% to 5.3% for morphometric measures (p < 0.001), whereas Ct.Po was 16.2% (p < 0.001). In the absence of subject motion, multicenter precision errors for HR-pQCT parameters were generally less than 5%. Phantom-based multicenter precision was comparable to previously reported in in vivo single-center precision errors, although this was approximately two to five times worse than ex vivo short-term precision. The data generated from this study will contribute to the future design and validation of standardized procedures that are broadly translatable to multicenter study designs.

Longitudinal in vivo imaging of bone formation and resorption using fluorescence molecular tomography

Bone research often focuses on anatomical imaging of the bone microstructure, but in order to gain better understanding in how bone remodeling is modulated through interventions also bone formation and resorption processes should be investigated. With this in mind, the purpose of this study was to establish a longitudinal in vivo imaging approach of bone formation and resorption using fluorescence molecular tomography (FMT). In this study the reproducibility, accuracy and sensitivity of FMT for bone imaging were assessed by performing longitudinal measurements with FMT and comparing it to in vivo micro-computed tomography on a set of control mice, and mice in which load-adaptation was induced in the sixth caudal vertebra. The precision error for FMT measurements, expressed as coefficient of variation, was smaller than 16%, indicating acceptable reproducibility. A correlation was found between bone resorption measured with FMT and bone resorption rate measured with in vivo micro-computed tomography only over the first 14days (R=0.81, p<0.01), but not between bone formation measured with FMT and bone formation rate measured with in vivo micro-CT. Bone formation measured by FMT was 89-109% greater (p<0.05) for mice subjected to mechanical loading than control mice. Bone resorption was 5-8% lower, but did not reach a significant difference between groups, indicating moderate sensitivity for FMT. In conclusion, in vivo FMT in mouse tail bones is feasible but needs to be optimized for monitoring load adaptation in living mice.

Age- and gender-related differences in cortical geometry and microstructure: Improved sensitivity by regional analysis

OBJECTIVE

While the importance of cortical structure quantification is increasingly underscored by recent literature, conventional analysis techniques obscure potentially important regional variations in cortical structure. The objective of this study was to characterize the spatial variability in cortical geometry and microstructure at the distal radius and tibia using high resolution peripheral quantitative computed tomography (HR-pQCT). We show that spatially-resolved analysis is able to identify cortical sub-regions with increased sensitivity to the effects of gender and aging.

METHODS

HR-pQCT scans of 146 volunteers (92 female/54 male) spanning a wide range of ages (20-78years) were analyzed. For each subject, radius and tibia scans were obtained using a clinical HR-pQCT system. Measures describing geometry (cortical bone thickness (Ct.Th)), microstructure (porosity (Ct.Po), pore diameter (Ct.Po.Dm), and pore size heterogeneity (Ct.Po.Dm SD)), and cortical bone density were calculated from the image data. Biomechanical parameters describing load and stress distribution were calculated using linear finite element analysis. Cortical quadrants were defined based on anatomic axes to quantify regional parameter variation. Subjects were categorized by gender, and age, and menopausal status for analysis.

RESULTS

Significant regional variation was found in all geometric and microstructural parameters in both the radius and tibia. In general, the radius showed more pronounced and significant variations in all parameters as compared with the tibia. At both sites, Ct.Po displayed the greatest regional variations. Correlation coefficients for Ct.Po and Ct.Th with respect to load and stress distribution provided evidence of an association between regional cortical structure and biomechanics in the tibia. Comparing women to men, differences in Ct.Po were most pronounced in the anterior quadrant of the radius (36% lower in women (p<0.01)) and the posterior quadrant of the tibia (27% lower in women (p<0.01)). Comparing elderly to young women, differences in Ct.Po were most pronounced in the lateral quadrant of the radius (328% higher in elderly women (p<0.001)) and the anterior quadrant of the tibia (433% higher in elderly women (p<0.001)). Comparing elderly to young men, the most pronounced age differences were found in the anterior radius (205% higher in elderly men, (p<0.001)) and the anterior tibia (190% higher in elderly men (p<0.01)). All subregional Ct.Po differences provided greater sensitivity to gender and age effects than those based on the global means.

CONCLUSION

These results show significant regional variation in all geometric and microarchitectural parameters studied in both the radius and tibia. Quantification of region-specific parameters provided increased sensitivity in the analysis of age- and gender-related differences, in many cases providing statistically significant differentiation of groups where conventional global analysis failed to detect differences. These results suggest that regional analysis may be important in studies of disease and therapeutic effects, particularly where microstructural parameters based on global analyses have thus far failed to identify a response in bone quality.

Patient-specific bone modelling and remodelling simulation of hypoparathyroidism based on human iliac crest biopsies

We previously developed a load-adaptive bone modelling and remodelling simulation model that can predict changes in the bone micro-architecture as a result of changes in mechanical loading or cell activity. In combination with a novel algorithm to estimate loading conditions, this offers the possibility for patient-specific predictions of bone modelling and remodelling. Based on such models, the underlying mechanisms of bone diseases and/or the effects of certain drugs and their influence on the bone micro-architecture can be investigated. In the present study we test the ability of this approach to predict changes in bone micro-architecture during hypoparathyroidism (HypoPT), as an illustrative example. We hypothesize that, apart from reducing bone turnover, HypoPT must also lead to increased osteocyte mechanosensitivity in order to explain the changes in bone mass seen in patients. Healthy human iliac crest biopsies were used as the starting point for the simulations that mimic HypoPT conditions and the resultant micro-architectures were compared to age-matched clinical HypoPT biopsies. Simulation results were in good agreement with the clinical data when osteocyte mechanosensitivity was increased by 40%. In conclusion, the results confirm our hypothesis, and also demonstrate that patient-specific bone modelling and remodelling simulations are feasible.

Short-term in vivo precision of BMD and parameters of trabecular architecture at the distal forearm and tibia

In vivo hr-pQCT precision was determined in 42 postmenopausal women using double baseline measurements from a multicenter trial of odanacatib. Errors, e.g., at the radius below 1.3% for BMD and below 6.3% for trabecular structure, were comparable to single-center results. Motion artifacts remain a challenge, particularly at the forearm.

INTRODUCTION

The short-term in vivo precision of BMD, trabecular bone structure, cortical thickness and porosity of the forearm and tibia was measured by hr-pQCT. Also the effect of image quality on precision was evaluated.

METHODS

In 42 postmenopausal women (age 64.4 ± 6.8 years) out of 214 subjects enrolled in a multi center advanced imaging phase III study of odanacatib (DXA spine or hip T-scores between -1.5 and -3.5), double baseline hr-pQCT (XtremeCT) measurements with repositioning were performed. The standard ultradistal location and a second, more proximally located VOI were measured at the radius and tibia to better assess cortical thickness and porosity. Image analysis and quality grading (grades: perfect, slight artifacts, pronounced artifacts, unacceptable) were performed centrally.

RESULTS

At the radius RMS%CV values varied from 0.7% to 1.3% for BMD and BV/TV and from 5.6% to 6.3% for Tb.Sp, Tb.Th, Tb.N, and cortical porosity. Numerically at the tibia, precision errors were approx. 0.5% lower for BMD and 1% to 2% lower for structural parameters although most differences were insignificant. In the radius but not in the tibia, precision errors for cortical thickness were smaller at the distal compared to the ultradistal location (1% versus 2%).

CONCLUSIONS

BMD precision errors were lower than those for trabecular architecture and cortical porosity. Motion artifacts remain a challenge, particularly at the forearm. Quality grading remains subjective, and more objective evaluation methods are needed. Precision in the context of a multicenter clinical trial, with centralized training and scan analysis, was comparable to single-center results previously reported.

Image interpolation allows accurate quantitative bone morphometry in registered micro-computed tomography scans

Time-lapsed in vivo micro-computed tomography is a powerful tool to analyse longitudinal changes in the bone micro-architecture. Registration can overcome problems associated with spatial misalignment between scans; however, it requires image interpolation which might affect the outcome of a subsequent bone morphometric analysis. The impact of the interpolation error itself, though, has not been quantified to date. Therefore, the purpose of this ex vivo study was to elaborate the effect of different interpolator schemes [nearest neighbour, tri-linear and B-spline (BSP)] on bone morphometric indices. None of the interpolator schemes led to significant differences between interpolated and non-interpolated images, with the lowest interpolation error found for BSPs (1.4%). Furthermore, depending on the interpolator, the processing order of registration, Gaussian filtration and binarisation played a role. Independent from the interpolator, the present findings suggest that the evaluation of bone morphometry should be done with images registered using greyscale information.

Quality control for bone quality parameters affected by subject motion in high-resolution peripheral quantitative computed tomography

Subject motion during high-resolution peripheral quantitative computed tomography (HR-pQCT) causes image artifacts that affect morphological analysis of bone quality. The aim of our study was to determine effectiveness of techniques for quality control in the presence of motion in vivo including automated and manual approaches. First, repeatability of manual grading was determined within and between laboratories. Given proper training using a standardized scale and training images (provided by the manufacturer), we found that manual grading is suitable for repeatable image quality grading within and across sites (ICC>0.7). Both a new automated technique providing motion measures based on projection moments, and traditional manual grading (1=best, 5=worst) were subsequently used to assess subject data for motion in N=137 image pairs (scan/re-scan) from the Canadian Multicentre Osteoporosis Study (CaMos) Calgary cohort. High quality image pairs were selected and measurement precision was estimated by calculating the coefficient of variation (CV). Consistent with previous data, density parameters (e.g. total bone mineral density) are more robust than structural (e.g. trabecular number) or finite element parameters (e.g. failure load). To obtain acceptable measurement precision, images should not exceed a manual grading of 3 (on a scale from 1 to 5) or an automatic (ε(T)) grading of 1.2. Automatic and manual grading provide comparable quality control, but the advantage of the automated technique is its ability to provide a motion value at scan time (providing a basis for real time decision regarding re-scan requirements), and the assessment is objective. Notably, automatic motion measurement can be performed retrospectively based on original scan data, and is therefore well suited for multi-center studies as well as any research where objective quality control is paramount.

QCT-based finite element models predict human vertebral strength in vitro significantly better than simulated DEXA

SUMMARY

While dual energy X-ray absorptiometry (DXA) is considered the gold standard to evaluate fracture risk in vivo, in the present study, the quantitative computed tomography (QCT)-based finite element modeling has been found to provide a quantitative and significantly improved prediction of vertebral strength in vitro. This technique might be used in vivo considering however the much larger doses of radiation needed for QCT.

INTRODUCTION

Vertebral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. QCT-based finite element (FE) modeling is an engineering method to predict vertebral strength. The aim of this study was to compare the ability of FE and clinical diagnostic tools to predict vertebral strength in vitro using an improved testing protocol.

METHODS

Thirty-seven vertebral sections were scanned with QCT and high resolution peripheral QCT (HR-pQCT). Bone mineral content (BMC), total BMD (tBMD), areal BMD from lateral (aBMD-lat), and anterior-posterior (aBMD-ap) projections were evaluated for both resolutions. Wedge-shaped fractures were then induced in each specimen with a novel testing setup. Nonlinear homogenized FE models (hFE) and linear micro-FE (μFE) were generated from QCT and HR-pQCT images, respectively. For experiments and models, both structural properties (stiffness, ultimate load) and material properties (apparent modulus and strength) were computed and compared.

RESULTS

Both hFE and μFE models predicted material properties better than structural ones and predicted strength significantly better than aBMD computed from QCT and HR-pQCT (hFE: R² = 0.79, μFE: R² = 0.88, aBMD-ap: R² = 0.48-0.47, aBMD-lat: R² = 0.41-0.43). Moreover, the hFE provided reasonable quantitative estimations of the experimental mechanical properties without fitting the model parameters.

CONCLUSIONS

The QCT-based hFE method provides a quantitative and significantly improved prediction of vertebral strength in vitro when compared to simulated DXA. This superior predictive power needs to be verified for loading conditions that simulate even more the in vivo case for human vertebrae.

Integrated imaging approach to osteoporosis: state-of-the-art review and update

Osteoporosis is the most common of all metabolic bone disorders. It is characterized by low bone mass and microarchitectural deterioration of bone tissue, with a consequent increase in bone fragility and susceptibility to fractures. Because of the increasing aging of the world population, the number of persons affected by osteoporosis is also increasing. Complications related to osteoporosis can create social and economic burdens. For these reasons, the early diagnosis of osteoporosis is crucial. Conventional radiography allows qualitative and semiquantitative evaluation of osteoporosis, whereas other imaging techniques allow quantification of bone loss (eg, dual-energy x-ray absorptiometry and quantitative computed tomography [CT]), assessment for the presence of fractures (morphometry), and the study of bone properties (ultrasonography). In recent years, new imaging modalities such as micro-CT and high-resolution magnetic resonance imaging have been developed in an attempt to help diagnose osteoporosis in its early stages, thereby reducing social and economic costs and preventing patient suffering. The correct diagnosis of osteoporosis results in better management in terms of prevention and adequate pharmacologic or surgical treatment.

Determination of vertebral and femoral trabecular morphology and stiffness using a flat-panel C-arm-based CT approach

The importance of assessing trabecular architecture together with bone mineral density to determine bone stiffness and fracture risk in osteoporosis has been well established. However, no imaging modalities are available to assess trabecular architecture at clinically relevant sites in the axial skeleton. Recently developed flat-panel CT devices, however, offer resolutions that are potentially good enough to resolve bone architecture at these sites. The goal of the present study was to investigate how accurate trabecular architecture and stiffness can be determined based on images from such a device (XperCT, Philips Healthcare). Ten cadaver human C3 vertebrae, twelve T12 vertebrae and 12 proximal femora were scanned with XperCT while mimicking in-vivo scanning conditions and compared to scans of the same bones with microCT. Standard segmentation and morphology quantification algorithms were applied as well as finite element (FE) simulation based on segmented and gray value images. Results showed that mean trabecular separation (Tb.Sp) and number (Tb.N) can be accurately determined at all sites. The accuracy of other parameters, however, depended on the site. For T12 no other structural parameters could be accurately quantified and no FE-results could be obtained from segmented images. When using gray-level images, however, accurate determination of cancellous bone stiffness was possible. For the C3 vertebrae and proximal femora, mean bone volume fraction (BV/TV), Tb.Sp, Tb.N, and anisotropy (C3 only) could be determined accurately. For Tb.Th, structure model index (SMI, femur only), and anisotropy good correlations were obtained but the values were not determined accurately. FE simulations based on segmented images were accurate for the C3 vertebrae, but severely underestimated bone stiffness for the femur. Here also, this was improved by using the gray value models. In conclusion, XperCT does provide a resolution that is good enough to determine trabecular architecture, but the signal to noise ratio is key to the accuracy of the morphology measurement. When the trabeculae are thick e.g. in the femur or the noise is low, e.g. cervical spine, architecture and stiffness could be determined accurately, but when the trabeculae are thin and the noise is high, e.g. thoracic spine, architecture could not be determined accurately and the connectivity was lost and hence no mechanical properties could be calculated directly.

In vivo precision of a depth-specific topographic mapping technique in the CT analysis of osteoarthritic and normal proximal tibial subchondral bone density

OBJECTIVE

To test the in vivo precision of a depth-specific topographic mapping technique (CT-TOMASD, computed tomography topographic mapping of subchondral density) for the 3D assessment of subchondral cortical and trabecular bone density in normal and osteoarthritic (OA) human tibiae.

METHODS

Fourteen participants were recruited (3 men:11 women; mean age: 49.9, SD: 11.9 years) and categorized as normal (n = 7) or OA (n = 7). Each participant was scanned using clinical quantitative CT (QCT) three times over 2 days. We assessed average subchondral bone mineral density (BMD) across three layers (0-2.5 mm, 2.5-5 mm and 5-10 mm) measured in relation to depth from the subchondral surface. Regional analyses included: medial plateau BMD; lateral plateau BMD; anterior/central/posterior compartment BMD; medial:lateral (M:L) BMD ratio; and average BMD of a 10-mm diameter "focal spot," which searched each medial and lateral plateau for the highest focal densities present within each plateau. Precision was assessed using root mean square coefficients of variation (CV%(RMS)) and intraclass correlation coefficients (ICC).

RESULTS

Average CV%(RMS) precision errors for BMD measures were 2.3%, reaching a maximum CV%(RMS) of 3.9%. ICC showed high repeatability above 0.98.

CONCLUSIONS

CT-TOMASD offered precise 3D measures of subchondral BMD. This method has the potential to identify and quantify changes in subchondral BMD associated with OA in vivo.