Using 3D image registration to maximize the reproducibility of longitudinal bone strength assessment by HR-pQCT and finite element analysis

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D image registration to account for misalignment between images. This reduced variability in longitudinal FE estimates and improved our ability to measure in vivo changes in HR-pQCT studies of bone strength.

INTRODUCTION

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D rigid-body registration (3DR) to maximize reproducibility by accounting for misalignment between images.

METHODS

In our proposed approach, we used the full common bone volume defined by 3DR to estimate standard FE parameters. Using standard HR-pQCT imaging protocols, we validated the 3DR approach with ex vivo samples of the distal radius (n = 10, four repeat scans) by assessing whether 3DR can reduce measurement variability from repositioning error. We used in vivo data (n = 40, five longitudinal scans) to assess the sensitivity of 3DR to detect changes in bone strength at the distal radius by the standard deviation of the rate of change (σ), where the ideal value of σ is minimized to define true change. FE estimates by 3DR were compared to estimates by no registration (NR) and slice-matching (SM).

RESULTS

Group-wise comparisons of ex vivo variation (CV_RMS, %) found that FE measurement precision was improved by SM (CV_RMS < 0.80%) and 3DR (CV_RMS < 0.62%) compared to NR (CV_RMS~2%), and 3DR was advantageous as repositioning error increased. Longitudinal in vivo reproducibility was minimized by 3DR for failure load estimates (σ = 0.008 kN/month).

CONCLUSION

Although 3D registration cannot negate motion artifacts, it plays an important role in detecting and reducing variability in FE estimates for longitudinal HR-pQCT data and is well suited for estimating effects of interventions in in vivo longitudinal studies of bone strength.

Opportunistic CT screening predicts individuals at risk of major osteoporotic fracture

Millions of CT scans are performed annually and could be also used to opportunistically assess musculoskeletal health; however, it is unknown how well this secondary assessment relates to osteoporotic fracture. This study demonstrates that opportunistic CT screening is a promising tool to predict individuals with previous osteoporotic fracture.

INTRODUCTION

Opportunistic computed tomography (oCT) screening for osteoporosis and fracture risk determination complements current dual X-ray absorptiometry (DXA) diagnosis. This study determined major osteoporotic fracture prediction by oCT at the spine and hip from abdominal CT scans.

METHODS

Initial 1158 clinical abdominal CT scans were identified from administrative databases and were the basis to generate a cohort of 490 men and women with suitable abdominal CT scans. Participant CT scans met the following criteria: over 50 years of age, the scan had no image artifacts, and the field-of-view included the L4 vertebra and proximal femur. A total of 123 participants were identified as having previously suffered a fracture within 5 years of CT scan date. Fracture cause was identified from clinical data and used to create a low-energy fracture sub-cohort. At each skeletal site, bone mineral density (BMD) and finite element (FE)-estimated bone strength were determined. Logistic regression predicted fracture and receiver-operator characteristic curves analyzed prediction capabilities.

RESULTS

In participants with a fracture, low-energy fractures occurred in 88% of women and 79% of men. Fracture prediction by combining both BMD and FE-estimated bone strength was not statistically different than using either BMD or FE-estimated bone strength alone. Predicting low-energy fractures in women determined the greatest AUC of 0.710 by using both BMD and FE-estimated bone strength.

CONCLUSIONS

oCT screening using abdominal CT scans is effective at predicting individuals with previous fracture at major osteoporotic sites and offers a promising screening tool for skeletal health assessment.

Effect of high-dose vitamin D supplementation on peripheral arterial calcification: secondary analysis of a randomized controlled trial

Although high-dose vitamin D supplementation is common, effects on arterial calcification remain unexplored. Tibial artery calcification was identified and quantified over 3 years in participants randomized to 400, 4000, or 10,000 IU vitamin D₃ daily. High-dose vitamin D supplementation did not affect the development or progression of arterial calcification.

INTRODUCTION

To determine whether vitamin D supplementation has a dose-dependent effect on development and progression of arterial calcification.

METHODS

This was a secondary analysis of the Calgary Vitamin D Study, a 3-year, double-blind, randomized controlled trial conducted at a single-center in Calgary, Canada. Participants were community-dwelling adults aged 55-70 years with serum 25-hydroxyvitamin D 30-125 nmol/L. Participants were randomized 1:1:1 to receive vitamin D₃ 400, 4000, or 10,000 IU/day for 3 years. Tibial artery calcification was identified and quantified (in milligrams of hydroxyapatite, mgHA) using high-resolution peripheral quantitative computed tomography (HR-pQCT) at baseline and 6, 12, 24, and 36 months. Changes in calcification over time and treatment group interaction were evaluated using a constrained linear mixed effects model.

RESULTS

Of 311 randomized participants, 302 (400: 105, 4000: 96, 10,000: 101) were eligible for analysis of arterial calcification (54% male, mean (SD) age 62 (4) years, mean (SD) 25-hydroxyvitamin D 78.9 (19.9) nmol/L). At baseline, 85 (28%) had tibial artery calcification, and mean (95% CI) calcification quantity was 2.8 mgHA (95% CI 1.7-3.9). In these 85 participants, calcification quantity increased linearly by 0.020 mgHA/month (95% CI 0.012-0.029) throughout the study, with no evidence of a treatment-group effect (p = 0.645 for interaction). No participants developed new arterial calcifications during the study.

CONCLUSIONS

In this population of community-dwelling adults who were vitamin D replete at baseline, supplementation with vitamin D 400, 4000, or 10,000 IU/day did not have differential effects on the development or progression of arterial calcification over 3 years.

TRIAL REGISTRATION

clinicaltrials.gov (NCT01900860).

Longitudinal bone microarchitectural changes are best detected using image registration

Longitudinal studies of bone using high-resolution medical imaging may result in non-physiological measurements of longitudinal changes. In this study, we determined that three-dimensional image processing techniques best capture realistic longitudinal changes in bone density and should therefore be used with high-resolution imaging when studying bone changes over time.

INTRODUCTION

The purpose of this study was to determine which longitudinal analysis technique (no registration (NR), slice-match (SM) registration, or three-dimensional registration (3DR)) produced the most realistic longitudinal changes in a 3-year study of bone density and structure using high-resolution peripheral quantitative computed tomography (HR-pQCT).

METHODS

We assessed HR-pQCT scans of the distal radius and tibia for men and women (N = 40) aged 55-70 years at baseline and 6, 12, 24, and 36 months. To evaluate which longitudinal analysis technique (NR, SM, or 3DR) best captured physiologically reasonable 3-year changes, we calculated the standard deviation of the absolute rate of change in each bone parameter. The data were compared between longitudinal analysis techniques using repeated measures ANOVA and post hoc analysis.

RESULTS

As expected, both SM and 3DR better captured physiological longitudinal changes than NR. At the tibia, there were no differences between SM and 3DR; however, at the radius where precision was lower, 3DR produced better results for total bone mineral density.

CONCLUSIONS

At least SM or 3DR should be implemented in longitudinal studies using HR-pQCT. 3DR is preferable, particularly at the radius, to ensure that physiological changes in bone density are observed.

Optimizing HR-pQCT workflow: a comparison of bias and precision error for quantitative bone analysis

Manual correction of automatically generated contours for high-resolution peripheral quantitative computed tomography can be time consuming and introduces precision error. However, bias related to the automated protocol is unknown. This study provides insight into error bias that is present when using uncorrected contours and inter-operator precision error based on operator training.

INTRODUCTION

High-resolution peripheral quantitative computed tomography workflow includes manually correcting contours generated by the manufacturer's automated protocol. There is interest in minimizing corrections to save time and reduce precision error; however, bias related to the automated protocol is unknown. This study quantifies error bias when contours are uncorrected and identifies the impact of operator training on bias and precision error.

METHODS

Forty-five radii and tibiae scans across a representative range of bone density were analyzed using the automated and manually corrected contours of three operators, with training ranging from beginner to expert, and compared with a "ground truth" to estimate bias. Inter-operator precision was measured across operators.

RESULTS

The tibia had greater error bias than the radius when contours were uncorrected, with compartmental bone mineral densities and cortical microarchitecture having greatest biases, which could have significant implications for interpretation of studies using this skeletal site. Bias and precision error were greatest when contours were corrected by the beginner operator; however, when this operator was removed, bias was no longer present and inter-operator precision was between 0.01 and 3.74% for all parameters except cortical porosity.

CONCLUSION

These findings establish the need for manual correction and provide guidance on operator training needed to maximize workflow efficiency.

Highly stable fullerene-based porous molecular crystals with open metal sites

The synthesis of conventional porous crystals involves building a framework using reversible chemical bond formation, which can result in hydrolytic instability. In contrast, porous molecular crystals assemble using only weak intermolecular interactions, which generally do not provide the same environmental stability. Here, we report that the simple co-crystallization of a phthalocyanine derivative and a fullerene (C₆₀ or C₇₀) forms porous molecular crystals with environmental stability towards high temperature and hot aqueous base or acid. Moreover, by using diamond anvil cells and synchrotron single-crystal measurements, stability towards extreme pressure (>4 GPa) is demonstrated, with the stabilizing fullerene held between two phthalocyanines and the hold tightening at high pressure. Access to open metal centres within the porous molecular co-crystal is demonstrated by in situ crystallographic analysis of the chemisorption of pyridine, oxygen and carbon monoxide. This suggests strategies for the formation of highly stable and potentially functional porous materials using only weak van der Waals intermolecular interactions.

Cortical porosity exhibits accelerated rate of change in peri- compared with post-menopausal women

The rate of change in bone density was not different between peri- and post-menopausal women. Differences in rate of change were observed in bone microarchitecture, specifically cortical porosity (Ct.Po), where peri-menopausal women increased +9% per year compared with the +6% per year for post-menopausal women.

INTRODUCTION

The purpose of this study was to compare changes in bone density and microarchitecture in peri- and post-menopausal women over 6 years.

METHODS

Peri- (n = 26) and post- (n = 65) menopausal women were selected from the Canadian Multicenter Osteoporosis Study. Caucasian women were scanned on dual x-ray absorptiometry (DXA) and high-resolution peripheral quantitative computed tomography (HR-pQCT) at baseline and follow-up, an average 6 years later. To compare repeat scans, automated 3D image registration was conducted. At the radius and tibia, total volumetric BMD (Tt.BMD), total bone area (Tt.Ar) and cortical porosity (Ct.Po) were assessed, and finite element analysis estimated apparent bone strength.

RESULTS

At the tibia, the rate of change for Ct.Po and Tt.Ar was different between groups. Peri-menopausal women had a + 9% per year increase in Ct.Po, but this increase was slower for post-menopausal women at +6% per year (p = 0.049). In addition, post-menopausal women had an increase in Tt.Ar of +0.13% per year compared with a slower increase of +0.06% per year for peri-menopausal women (p = 0.017). The rate of change of density between groups was not significantly different and was approximately -1% per year at the hip by DXA, and -1% per year at the radius and -0.5% per year tibia by HR-pQCT.

CONCLUSION

This is a 6-year prospective HR-pQCT study exploring rate of change in Caucasian peri- and post-menopausal women. The microarchitectural features represented by Ct.Po and Tt.Ar changed at a significantly different rate between groups, but group differences were not detected by density measures.

A hindered subphthalocyanine that forms crystals with included aromatic solvent but will not play ball with C60

A novel subphthalocyanine containing bulky substituents placed at its peripheral sites (i.e. 2,3,9,10,16,17-hexa(2′,6′-di-iso-propylphenoxy)boron subphthalocyanine) was prepared and assessed for supramolecular binding with C[Formula: see text], through crystallisation and fluorescence studies. Three different crystal polymorphs of the subphthalocyanine were obtained that showed inclusion of a single aromatic solvent molecule within the well-defined cavity of the molecule but complete exclusion of C[Formula: see text]. Analysis of the crystal structures indicated that the bowl-shaped cavity of the subphthalocyanine molecule was only accessible to small molecules due to the steric congestion surrounding the macrocycle, which results in hindered rotation of the substituents on the NMR timescale. Enhanced solubility, up to 42.0 g/L in common organic solvents, was demonstrated consistent with the crystal structures which are dominated by relatively weak intermolecular interactions, which allow solvent molecules to play a role in crystallisation.

The relationship between serum 25(OH)D and bone density and microarchitecture as measured by HR-pQCT

The relation between serum 25-hydroxy vitamin D [25(OH)D] and bone quality is not well understood, particularly for high levels. We measured bone microarchitecture in three groups of people stratified by their serum 25(OH)D. There was a weak association of serum 25(OH)D and microarchitecture for this cross-sectional population, suggesting possible benefits to bone quality.

INTRODUCTION

Vitamin D plays an important role in bone and mineral metabolism, but the relation between serum 25(OH)D and bone quality is not well understood. Here, we present a cross-sectional study that investigated a convenience group of participants from an ongoing health initiative in Alberta, Canada, who have been receiving daily vitamin D supplementation.

METHODS

A total of 105 participants were organized into three groups based on their serum 25(OH)D levels: low (<75 nmol/L), medium (75-175 nmol/L), and high (>175 nmol/L). They were also assessed with 25(OH)D as a continuous variable. Average daily supplementation was 7670 ± 438 IU, and the change in 25(OH)D ranged from 22 to 33 % during the period of receiving supplements. We used high-resolution peripheral quantitative computed tomography measurements at the radius and tibia to assess bone microarchitecture.

RESULTS

Microarchitectural parameters were not strongly associated with serum 25(OH)D. In the tibia, there were fewer trabeculae (TbN; p = 0.015) and a non-significant trend toward thicker trabeculae (p = 0.067) of the high group. Body mass index (BMI) was negatively associated with serum 25(OH)D levels (p < 0.001) and PTH levels (p < 0.001). There was no clinically significant relationship detected between high serum 25(OH)D and high serum calcium.

CONCLUSION

These data suggest a weak relationship between serum 25(OH)D and bone microarchitecture in this population of mostly vitamin-D-sufficient participants, and there were no indications of negative effects related to the high supplementation levels. These data provided a basis to design and implement our 3-year dose-dependent randomized controlled trial investigating the effects of vitamin D supplementation on bone health outcomes.

Gymnastics participation is associated with skeletal benefits in the distal forearm: a 6-month study using peripheral Quantitative Computed Tomography

OBJECTIVES

Musculoskeletal development of the upper limbs during exposure to weight-bearing loading is under-researched during early pubescent growth. The purpose was to assess the changes in upper body musculoskeletal strength in young girls following 6 months of non-elite gymnastics participation.

METHODS

Eighty-four girls, 6-12 years were divided into groups based on gymnastics participation: high-training (HGYM, 6-16 hr/wk), low-training (LGYM, 1-5 hr/wk), and non-gymnasts (NONGYM). Volumetric BMD, bone geometry, estimated bone strength and muscle size were assessed at the non-dominant forearm (4% and 66% radius and ulna) with pQCT. DXA assessed aBMD and body composition. Tests for explosive power, muscle strength, and endurance were also performed.

RESULTS

Interaction effects were observed in all variables at the 4% radius. At the 66% ulna, HGYM and LGYM had greater bone mass, size and bone strength than NONGYM, furthermore a dose-response relationship was observed at this location. Body composition was better for HGYM than LGYM and NONGYM, however muscle function was better for HGYM and LGYM than NONGYM.

CONCLUSION

The greatest changes were obtained with more than one gymnastics class per week. Separating gymnastics participation-related changes from those associated with normal growth and development remains difficult, particularly at the 4% radius.

Non-elite gymnastics participation is associated with greater bone strength, muscle size, and function in pre- and early pubertal girls

Recent reports indicate an increase in forearm fractures in children. Bone geometric properties are an important determinant of bone strength and therefore fracture risk. Participation in non-elite gymnastics appears to contribute to improving young girls' musculoskeletal health, more specifically in the upper body.

INTRODUCTION

The primary aim of this study was to determine the association between non-elite gymnastics participation and upper limb bone mass, geometry, and strength in addition to muscle size and function in young girls.

METHODS

Eighty-eight pre- and early pubertal girls (30 high-training gymnasts [HGYM, 6-16 hr/ wk], 29 low-training gymnasts [LGYM, 1-5 h r/wk] and 29 non-gymnasts [NONGYM]), aged 6-11 years were recruited. Upper limb lean mass, BMD and BMC were derived from a whole body DXA scan. Forearm volumetric BMD, bone geometry, estimated strength, and muscle CSA were determined using peripheral QCT. Upper body muscle function was investigated with muscle strength, explosive power, and muscle endurance tasks.

RESULTS

HGYM showed greater forearm bone strength compared with NGYM, as well as greater arm lean mass, BMC, and muscle function (+5% to +103%, p < 0.05). LGYM displayed greater arm lean mass, BMC, muscle power, and endurance than NGYM (+4% to +46%, p < 0.05); however, the difference in bone strength did not reach significance. Estimated fracture risk at the distal radius, which accounted for body weight, was lower in both groups of gymnasts. Compared with NONGYM, HGYM tended to show larger skeletal differences than LGYM; yet, the two groups of gymnasts only differed for arm lean mass and muscle CSA.

CONCLUSION

Non-elite gymnastics participation was associated with musculoskeletal benefits in upper limb bone geometry, strength and muscle function. Differences between the two gymnastic groups emerged for arm lean mass and muscle CSA, but not for bone strength.

Skeletal differences at the ulna and radius between pre-pubertal non-elite female gymnasts and non-gymnasts

OBJECTIVE

To compare skeletal parameters between the ulna and radius in pre-pubertal non-elite gymnasts and non-gymnasts.

METHODS

Fifty-eight non-elite artistic gymnasts, aged 6-11 years, were compared with 28 non-gymnasts for bone mineral content (BMC), total and cortical bone area (ToA, CoA), trabecular and cortical volumetric density (TrD, CoD) and estimated bone strength (BSI and SSIp), obtained by pQCT at the distal and proximal forearm.

RESULTS

Gymnasts had greater estimated bone strength than non-gymnasts at both sites of the forearm. At the distal forearm, the gymnastics-induced skeletal benefits were greater at the radius than ulna (Z-scores for BMC, TrD and BSI +0.40 to +0.61 SD, p<0.05 vs. +0.15 to +0.48 SD, NS). At the proximal forearm, the skeletal benefits were greater at the ulna than the radius (Z-scores for BMC, ToA, CoA and SSIp +0.59 to +0.82 SD, p<0.01 vs. +0.35 (ToA) and +0.43 SD (SSIp), p<0.01).

CONCLUSION

Skeletal benefits at the distal and proximal forearm emerged in young non-elite gymnasts. Benefits were larger when considering skeletal parameters at both the ulna and radius, than the radius alone as traditionally performed with pQCT. These findings suggest the ulna is worth investigating in future studies aiming to accurately quantify exercise-induced skeletal adaptations.