Systematic misregistration and the statistical analysis of surface data

A Correspondence-Based Network Approach for Groupwise Analysis of Patient-Specific Spatiotemporal Data

Methods for statistically analyzing patient-specific data that vary both spatially and over time are currently either limited to summary statistics or require elaborate surface registration. We propose a new method, called correspondence-based network analysis, which leverages particle-based shape modeling to establish correspondence across a population and preserve patient-specific measurements and predictions through statistical analysis. Herein, we evaluated this method using three published datasets of the hip describing cortical bone thickness of the proximal femur, cartilage contact stress, and dynamic joint space between control and patient cohorts to evaluate activity- and group-based differences, as applicable, using traditional statistical parametric mapping (SPM) and our proposed spatially considerate correspondence-based network analysis approach. The network approach was insensitive to correspondence density, while the traditional application of SPM showed decreasing area of the region of significance with increasing correspondence density. In comparison to SPM, the network approach identified broader and more connected regions of significance for all three datasets. The correspondence-based network analysis approach identified differences between groups and activities without loss of subject and spatial specificity which could improve clinical interpretation of results.

Dynamic multi feature-class Gaussian process models

In model-based medical image analysis, three relevant features are the shape of structures of interest, their relative pose, and image intensity profiles representative of some physical properties. Often, these features are modelled separately through statistical models by decomposing the object's features into a set of basis functions through principal geodesic analysis or principal component analysis. However, analysing articulated objects in an image using independent single object models may lead to large uncertainties and impingement, especially around organ boundaries. Questions that come to mind are the feasibility of building a unique model that combines all three features of interest in the same statistical space, and what advantages can be gained for image analysis. This study presents a statistical modelling method for automatic analysis of shape, pose and intensity features in medical images which we call the Dynamic multi feature-class Gaussian process models (DMFC-GPM). The DMFC-GPM is a Gaussian process (GP)-based model with a shared latent space that encodes linear and non-linear variations. Our method is defined in a continuous domain with a principled way to represent shape, pose and intensity feature-classes in a linear space, based on deformation fields. A deformation field-based metric is adapted in the method for modelling shape and intensity variation as well as for comparing rigid transformations (pose). Moreover, DMFC-GPMs inherit properties intrinsic to GPs including marginalisation and regression. Furthermore, they allow for adding additional pose variability on top of those obtained from the image acquisition process; what we term as permutation modelling. For image analysis tasks using DMFC-GPMs, we adapt Metropolis-Hastings algorithms making the prediction of features fully probabilistic. We validate the method using controlled synthetic data and we perform experiments on bone structures from CT images of the shoulder to illustrate the efficacy of the model for pose and shape prediction. The model performance results suggest that this new modelling paradigm is robust, accurate, accessible, and has potential applications in a multitude of scenarios including the management of musculoskeletal disorders, clinical decision making and image processing.

Rib cortical bone thickness variation in adults by age and sex

Rib fractures are a common and serious outcome of blunt thoracic trauma and their likelihood is greater in older individuals. Osteoporotic bone loss is a well-documented aging phenomenon with sex-specific characteristics, but within rib bones, neither baseline maps of regional thickness nor the rates of bone thinning with age have been quantified across whole ribs. This study presents such data from 4014 ribs of 240 adult subjects aged 20-90. A validated cortical bone mapping technique was applied to clinical computed tomography scans to obtain local rib cortical bone thickness measurements over the surfaces of ribs 2 through 11. Regression models to age and sex gave rates of cortex thinning in local zones and aggregated across whole ribs. The statistical parametric mapping provided these relationships regionally as a function of rib surface location. All models showed significant reductions in bone thickness with age (p < 0.01). Average whole-rib thinning occurred at between 0.011 to 0.032 mm/decade (males) and 0.035 to 0.043 mm/decade (females), with sex and age accounting for up to 37% of population variability (R² ). Rates of thinning differed regionally and by rib, with the highest bone loss of up to 0.074 mm/decade occurring in mid-rib cutaneous and superior regions of ribs 2-6. Rates were consistently higher in females than males (significantly so across whole ribs but not all local regions) and were more pronounced in cutaneous, superior, and inferior rib aspects (average 0.025 mm/decade difference in ribs 4-8) compared to pleural aspects which had the thickest cortices but saw only minor differences in thinning rates by sex (0.045 mm/decade for females and 0.040 mm/decade for males). Regional analysis showed male and female bone thickness differences that were not statistically significant at 20 years of age (p > 0.05 across practically all regions) but subsequent cortex thinning meant that substantial pleural and cutaneous regions were thinner (p < 0.05) in females than males by 55 years of age. The techniques and results from this study can be applied to assess rib bone content loss in clinical settings across wide populations. Additionally, average cortex thickness results can be mapped directly to finite element models of the thorax, and regression results are used to modify such models to represent the ribs of men and women across their full adult lifespan.

Multiparametric 3-D analysis of bone and joint space width at the knee from weight bearing computed tomography

OBJECTIVE

Computed tomography (CT) can deliver multiple parameters relevant to osteoarthritis. In this study we demonstrate that a 3-D multiparametric approach at the weight bearing knee with cone beam CT is feasible, can include multiple parameters from across the joint space, and can reveal stronger relationships with disease status in combination.

DESIGN

33 participants with knee weight bearing CT (WBCT) were analysed with joint space mapping and cortical bone mapping to deliver joint space width (JSW), subchondral bone plate thickness, endocortical thickness, and trabecular attenuation at both sides of the joint. All data were co-localised to the same canonical surface. Statistical parametric mapping (SPM) was applied in uni- and multivariate models to demonstrate significant dependence of parameters on Kellgren & Lawrence grade (KLG). Correlation between JSW and bony parameters and 2-week test-retest repeatability were also calculated.

RESULTS

SPM revealed that the central-to-posterior medial tibiofemoral joint space was significantly narrowed by up to 0.5 mm with significantly higher tibial trabecular attenuation up to 50 units for each increment in KLG as single features, and in a wider distribution when combined (p<0.05). These were also more strongly correlated with worsening KLG grade category. Test-retest repeatability was subvoxel (0.37 mm) for nearly all thickness parameters.

CONCLUSIONS

3-D JSW and tibial trabecular attenuation are repeatable and significantly dependent on radiographic disease severity at the weight bearing knee joint not just alone, but more strongly in combination. A quantitative multiparametric approach with WBCT may have potential for more sensitive investigation of disease progression in osteoarthritis.

Quantifying Differences Between Affine and Nonlinear Spatial Normalization of FP-CIT Spect Images

Spatial normalization helps us to compare quantitatively two or more input brain scans. Although using an affine normalization approach preserves the anatomical structures, the neuroimaging field is more common to find works that make use of nonlinear transformations. The main reason is that they facilitate a voxel-wise comparison, not only when studying functional images but also when comparing MRI scans given that they fit better to a reference template. However, the amount of bias introduced by the nonlinear transformations can potentially alter the final outcome of a diagnosis especially when studying functional scans for neurological disorders like Parkinson's Disease. In this context, we have tried to quantify the bias introduced by the affine and the nonlinear spatial registration of FP-CIT SPECT volumes of healthy control subjects and patients with PD. For that purpose, we calculated the deformation fields of each participant and applied these deformation fields to a 3D-grid. As the space between the edges of small cubes comprising the grid change, we can quantify which parts from the brain have been enlarged, compressed or just remain the same. When the nonlinear approach is applied, scans from PD patients show a region near their striatum very similar in shape to that of healthy subjects. This artificially increases the interclass separation between patients with PD and healthy subjects as the local intensity is decreased in the latter region, and leads machine learning systems to biased results due to the artificial information introduced by these deformations.

Physical activity induced adaptation can increase proximal femur strength under loading from a fall onto the greater trochanter

Physical activity enhances proximal femur bone mass, but it remains unclear whether the benefits translate into an enhanced ability to resist fracture related loading. We recently used baseball pitchers as a within-subject controlled model to demonstrate physical activity induced proximal femur adaptation in regions associated with weight bearing and increased strength under single-leg stance loading. However, there was no measurable benefit to resisting common injurious loading (e.g. a fall onto the greater trochanter). A lack of power and a small physical activity effect size may have contributed to the latter null finding. Softball pitchers represent an alternative within-subject controlled model to explore adaptation of the proximal femur to physical activity, exhibiting greater dominant-to-nondominant (D-to-ND) leg differences than baseball pitchers. The current study used quantitative computed tomography, statistical parametric mapping, and subject-specific finite element (FE) modeling to explore adaptation of the proximal femur to physical activity in female softball pitchers (n = 25). Female cross-country runners (n = 15) were included as symmetrically loaded controls, showing very limited D-to-ND leg differences. Softball pitchers had D-to-ND leg differences in proximal femur, femoral neck, and trochanteric volumetric bone mineral density and content, and femoral neck volume. Voxel-based morphometry analyses and cortical bone mapping showed D-to-ND leg differences within a large region connecting the superior femoral head, inferior femoral neck and medial intertrochanteric region, and within the greater trochanter. FE modeling revealed pitchers had 19.4% (95%CI, 15.0 to 23.9%) and 4.9% (95%CI, 1.7 to 8.2%) D-to-ND leg differences in predicted ultimate strength under single-leg stance loading and a fall onto the greater trochanter, respectively. These data affirm the spatial and strength adaptation of the proximal femur to weight bearing directed loading and demonstrate that the changes can also have benefits, albeit smaller, on resisting loads associated with a sideways fall onto the greater trochanter.

Quantitative Three-dimensional Assessment of Knee Joint Space Width from Weight-bearing CT

Background Imaging of structural disease in osteoarthritis has traditionally relied on MRI and radiography. Joint space mapping (JSM) can be used to quantitatively map joint space width (JSW) in three dimensions from CT images. Purpose To demonstrate the reproducibility, repeatability, and feasibility of JSM of the knee using weight-bearing CT images. Materials and Methods Two convenience samples of weight-bearing CT images of left and right knees with radiographic Kellgren-Lawrence grades (KLGs) less than or equal to 2 were acquired from 2014 to 2018 and were analyzed retrospectively with JSM to deliver three-dimensional JSW maps. For reproducibility, images of three sets of knees were used for novice training, and then the JSM output was compared against an expert's assessment. JSM was also performed on 2-week follow-up images in the second cohort, yielding three-dimensional JSW difference maps for repeatability. Statistical parametric mapping was performed on all knee imaging data (KLG, 0-4) to show the feasibility of a surface-based analysis in three dimensions. Results Reproducibility (in 20 individuals; mean age, 58 years ± 7 [standard deviation]; mean body mass index, 28 kg/m² ± 6; 14 women) and repeatability (in nine individuals; mean age, 53 years ± 6; mean body mass index, 26 kg/m² ± 4; seven women) reached their lowest performance at a smallest detectable difference less than ±0.1 mm in the central medial tibiofemoral joint space for individuals without radiographically demonstrated disease. The average root mean square coefficient of variation was less than 5% across all groups. Statistical parametric mapping (33 individuals; mean age, 57 years ± 7; mean body mass index, 27 kg/m² ± 6; 23 women) showed that the central-to-posterior medial joint space was significantly narrower by 0.5 mm for each incremental increase in the KLG (threshold P < .05). One knee (KLG, 2) demonstrated a baseline versus 24-month change in its three-dimensional JSW distribution that was beyond the smallest detectable difference across the lateral joint space. Conclusion Joint space mapping of the knee using weight-bearing CT images is feasible, demonstrating a relationship between the three-dimensional joint space width distribution and structural joint disease. It is reliably learned by novice users, can be personalized for disease phenotypes, and can be used to achieve a smallest detectable difference that is at least 50% smaller than that reported to be achieved at the highest performance level in radiography. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Roemer in this issue.

Using point clouds to investigate the relationship between trabecular bone phenotype and behavior: An example utilizing the human calcaneus

OBJECTIVES

The objective of this study is to demonstrate a new method for analyzing trabecular bone volume fraction and degree of anisotropy in three dimensions.

METHODS

We use a combination of automatic mesh registration, point-cloud correspondence registration, and P-value corrected univariate statistical tests to compare bone volume fraction and degree of anisotropy on a point by point basis across the entire calcaneus of two human groups with different subsistence strategies.

RESULTS

We found that the patterns of high and low bone volume fraction and degree of anisotropy distribution between the Black Earth (hunter-gatherers) and Norris Farms (mixed-strategy agriculturalists) are very similar, but differ in magnitude. The hunter-gatherers exhibit higher levels of bone volume fraction and less anisotropic trabecular bone organization. Additionally, patterns of bone volume fraction and degree of anisotropy in the calcaneus correspond well with biomechanical expectations of relative forces experienced during walking and running.

CONCLUSIONS

We conclude that comparing site-specific, localized differences in trabecular bone variables such as bone volume fraction and degree of anisotropy in three-dimensions is a powerful analytical tool. This method makes it possible to determine where similarities and differences between groups are located within the whole skeletal element of interest. The visualization of multiple variables also provides a way for researchers to see how the trabecular bone variables interact within the morphology, and allows for a more nuanced understanding of how they relate to one another and the broader mechanical environment.

Heterogeneous Spatial and Strength Adaptation of the Proximal Femur to Physical Activity: A Within-Subject Controlled Cross-Sectional Study

Physical activity (PA) enhances proximal femur bone mass, as assessed using projectional imaging techniques. However, these techniques average data over large volumes, obscuring spatially heterogeneous adaptations. The current study used quantitative computed tomography, statistical parameter mapping, and subject-specific finite element (FE) modeling to explore spatial adaptation of the proximal femur to PA. In particular, we were interested in adaptation occurring at the superior femoral neck and improving strength under loading from a fall onto the greater trochanter. High/long jump athletes (n = 16) and baseball pitchers (n = 16) were utilized as within-subject controlled models as they preferentially load their take-off leg and leg contralateral to their throwing arm, respectively. Controls (n = 15) were included but did not show any dominant-to-nondominant (D-to-ND) leg differences. Jumping athletes showed some D-to-ND leg differences but less than pitchers. Pitchers had 5.8% (95% confidence interval [CI] 3.9%-7.6%) D-to-ND leg differences in total hip volumetric bone mineral density (vBMD), with increased vBMD in the cortical compartment of the femoral neck and trochanteric cortical and trabecular compartments. Voxel-based morphometry analyses and cortical bone mapping showed pitchers had D-to-ND leg differences within the regions of the primary compressive trabeculae, inferior femoral neck, and greater trochanter but not the superior femoral neck. FE modeling revealed pitchers had 4.1% (95% CI 1.4%-6.7%) D-to-ND leg differences in ultimate strength under single-leg stance loading but no differences in ultimate strength to a fall onto the greater trochanter. These data indicate the asymmetrical loading associated with baseball pitching induces proximal femur adaptation in regions associated with weight bearing and muscle contractile forces and increases strength under single-leg stance loading. However, there were no benefits evident at the superior femoral neck and no measurable improvement in ultimate strength to common injurious loading during aging (ie, fall onto the greater trochanter), raising questions as to how to better target these variables with PA. © 2019 American Society for Bone and Mineral Research.

Quantitative 3D imaging parameters improve prediction of hip osteoarthritis outcome

Osteoarthritis is an increasingly important health problem for which the main treatment remains joint replacement. Therapy developments have been hampered by a lack of biomarkers that can reliably predict disease, while 2D radiographs interpreted by human observers are still the gold standard for clinical trial imaging assessment. We propose a 3D approach using computed tomography—a fast, readily available clinical technique—that can be applied in the assessment of osteoarthritis using a new quantitative 3D analysis technique called joint space mapping (JSM). We demonstrate the application of JSM at the hip in 263 healthy older adults from the AGES-Reykjavík cohort, examining relationships between 3D joint space width, 3D joint shape, and future joint replacement. Using JSM, statistical shape modelling, and statistical parametric mapping, we show an 18% improvement in prediction of joint replacement using 3D metrics combined with radiographic Kellgren & Lawrence grade (AUC 0.86) over the existing 2D FDA-approved gold standard of minimum 2D joint space width (AUC 0.73). We also show that assessment of joint asymmetry can reveal significant differences between individuals destined for joint replacement versus controls at regions of the joint that are not captured by radiographs. This technique is immediately implementable with standard imaging technologies.

Adaptation of the proximal humerus to physical activity: A within-subject controlled study in baseball players

The proximal humerus is a common, yet understudied site for osteoporotic fracture. The current study explored the impact of prolonged physical activity on proximal humerus bone health by comparing bone properties between the throwing and nonthrowing arms within professional baseball players. The proximal humerus in throwing arms had 28.1% (95% CI, 17.8 to 38.3%) greater bone mass compared to nonthrowing arms, as assessed using dual-energy x-ray absorptiometry. At the level of the surgical neck, computed tomography revealed 12.0% (95% CI, 8.2 to 15.8%) greater total cross-sectional area and 31.0% (95% CI, 17.8 to 44.2%) greater cortical thickness within throwing arms, which contributed to 56.8% (95% CI, 44.9 to 68.8%) greater polar moment of inertia (i.e., estimated ability to resist torsional forces) compared to nonthrowing arms. Within the humeral head and greater tubercle regions, throwing arms had 3.1% (95% CI, 1.1 to 5.1%) more trabecular bone, as assessed using high-resolution magnetic resonance imaging. Three-dimensional mapping of voxel- and vertex-wise differences between arms using statistical parametric mapping techniques revealed throwing arms had adaptation within much of the proximal diaphysis, especially the posterolateral cortex. The pattern of proximal diaphysis adaptation approximated the pattern of strain energy distribution within the proximal humerus during a fastball pitch derived from a musculoskeletal and finite element model in a representative player. These data demonstrate the adaptive ability of the proximal humerus to physical activity-related mechanical loads. It remains to be established how they translate to exercise prescription to improve bone health within the proximal humerus; however, they provide unique insight into the relationship between prolonged loading and skeletal adaptation at a clinically relevant osteoporotic site.

Cortical Bone Mapping: Measurement and Statistical Analysis of Localised Skeletal Changes

PURPOSE OF REVIEW

Cortical bone mapping (CBM) is a technique for measuring localised skeletal changes from computed tomography (CT) images. It can provide measurements with accuracy surpassing the underlying imaging resolution. CBM can detect changes in several properties of the cortex, with no prior assumptions about the likely location of said changes. This paper summarises the theory behind CBM, discusses its strengths and limitations, and reviews some studies in which it has been applied.

RECENT FINDINGS

CBM has revealed associations between fracture risk and cortical properties in specific regions of the proximal femur which present feasible therapeutic targets. Analyses of several pharmaceutical and exercise interventions quantify effects that are distinct both in location and in the nature of the micro-architectural changes. CBM has illuminated age-related changes in the proximal femur and has recently been applied to other bones, as well as to the assessment of cartilage. The CBM processing pipeline is designed primarily for large cohort studies. Its main impact thus far has not been in the realm of clinical practice, but rather to improve our fundamental understanding of localised bone structure and changes.

Sex differences in the spatial distribution of bone in relation to incident hip fracture: Findings from the AGES-Reykjavik study

In this case-cohort study, we used data-driven computational anatomy approaches to assess within and between sex spatial differences in proximal femoral bone characteristics in relation to incident hip fracture. One hundred male and 234 female incident hip fracture cases, and 1047 randomly selected noncase subcohort participants (562 female) were chosen from the population-based AGES-Reykjavik study (mean age of 77 years). The baseline -i.e. before hip fracture- hip quantitative computed tomography scans of these subjects were analyzed using voxel-based morphometry, tensor-based morphometry, and surface-based statistical parametric mapping to assess the spatial distribution of volumetric bone mineral density (vBMD), internal structure, and cortical bone properties (thickness, vBMD and trabecular vBMD adjacent to the endosteal surface) of the proximal femur, respectively, in relation to incident hip fracture. Results showed that in both men and women: 1) the superior aspect of the femoral neck and the trochanteric region (except for cortical bone thickness) were consistently identified as being associated with incident hip fracture, and 2) differences in bone properties between noncases and incident hip fracture cases followed similar trends, were located at compatible regions, and manifested heterogeneity in the spatial distribution of their magnitude with focal regions showing larger differences. With respect to sex differences, most of the regions with a significant interaction between fracture group and sex showed: 1) differences of greater magnitude in men between noncases and incident hip fracture cases with different spatial distributions for all bone properties with the exception of cortical bone thickness, and 2) that while most of these regions showed better bone quality in male cases than in female cases, female cases showed higher vBMD in the principal compressive group and higher endotrabecular vBMD at several regions including the anterior, posterior, and lateral aspects of the proximal femur. These findings indicate the value of these image analysis techniques by providing unique information about the specific patterns of bone deterioration associated with incident hip fracture and their sex differences, highlighting the importance of looking to men and women separately in the assessment of hip fracture risk.

Trabecular bone patterning across the human hand

Hand bone morphology is regularly used to link particular hominin species with behaviors relevant to cognitive/technological progress. Debates about the functional significance of differing hominin hand bone morphologies tend to rely on establishing phylogenetic relationships and/or inferring behavior from epigenetic variation arising from mechanical loading and adaptive bone modeling. Most research focuses on variation in cortical bone structure, but additional information about hand function may be provided through the analysis of internal trabecular structure. While primate hand bone trabecular structure is known to vary in ways that are consistent with expected joint loading differences during manipulation and locomotion, no study exists that has documented this variation across the numerous bones of the hand. We quantify the trabecular structure in 22 bones of the human hand (early/extant modern Homo sapiens) and compare structural variation between two groups associated with post-agricultural/industrial (post-Neolithic) and foraging/hunter-gatherer (forager) subsistence strategies. We (1) establish trabecular bone volume fraction (BV/TV), modulus (E), degree of anisotropy (DA), mean trabecular thickness (Tb.Th) and spacing (Tb.Sp); (2) visualize the average distribution of site-specific BV/TV for each bone; and (3) examine if the variation in trabecular structure is consistent with expected joint loading differences among the regions of the hand and between the groups. Results indicate similar distributions of trabecular bone in both groups, with those of the forager sample presenting higher BV/TV, E, and lower DA, suggesting greater and more variable loading during manipulation. We find indications of higher loading along the ulnar side of the forager sample hand, with high site-specific BV/TV distributions among the carpals that are suggestive of high loading while the wrist moves through the 'dart-thrower's' motion. These results support the use of trabecular structure to infer behavior and have direct implications for refining our understanding of human hand evolution and fossil hominin hand use.

How does the femoral cortex depend on bone shape? A methodology for the joint analysis of surface texture and shape

In humans, there is clear evidence of an association between hip fracture risk and femoral neck bone mineral density, and some evidence of an association between fracture risk and the shape of the proximal femur. Here, we investigate whether the femoral cortex plays a role in these associations: do particular morphologies predispose to weaker cortices? To answer this question, we used cortical bone mapping to measure the distribution of cortical mass surface density (CMSD, mg/cm2) in a cohort of 125 females. Principal component analysis of the femoral surfaces identified three modes of shape variation accounting for 65% of the population variance. We then used statistical parametric mapping (SPM) to locate regions of the cortex where CMSD depends on shape, allowing for age. Our principal findings were increased CMSD with increased gracility over much of the proximal femur; and decreased CMSD at the superior femoral neck, coupled with increased CMSD at the calcar femorale, with increasing neck-shaft angle. In obtaining these results, we studied the role of spatial normalization in SPM, identifying systematic misregistration as a major impediment to the joint analysis of CMSD and shape. Through a series of experiments on synthetic data, we evaluated a number of registration methods for spatial normalization, concluding that only those predicated on an explicit set of homologous landmarks are suitable for this kind of analysis. The emergent methodology amounts to an extension of Geometric Morphometric Image Analysis to the domain of textured surfaces, alongside a protocol for labelling homologous landmarks in clinical CT scans of the human proximal femur.

Ricci-flow based conformal mapping of the proximal femur to identify exercise loading effects

The causal relationship between habitual loading and adaptive response in bone morphology is commonly explored by analysing the spatial distribution of mechanically relevant features. In this study, 3D distribution of features in the proximal femur of 91 female athletes (5 exercise loading groups representing habitual loading) is contrasted with 20 controls. A femur specific Ricci-flow based conformal mapping procedure was developed for establishing correspondence among the periosteal surfaces. The procedure leverages the invariance of the conformal mapping method to isometric shape differences to align surfaces in the 2D parametric domain, to produce dense correspondences across an isotopological set of surfaces. This is implemented through a multi-parametrisation approach to detect surface features and to overcome the issue of inconsistency in the anatomical extent present in the data. Subsequently, the group-wise distribution of two mechanically relevant features was studied – cortical thickness and surface principal strains (simulation results of a sideways fall). Statistical inferences over the surfaces were made by contrasting the athlete groups with the controls through statistical parametric mapping. With the aid of group-wise and composite-group maps, proximal femur regions affected by specific loading groups were identified with a high degree of spatial localisation.

Spatial Differences in the Distribution of Bone Between Femoral Neck and Trochanteric Fractures

There is little knowledge about the spatial distribution differences in volumetric bone mineral density and cortical bone structure at the proximal femur between femoral neck fractures and trochanteric fractures. In this case-control study, a total of 93 women with fragility hip fractures, 72 with femoral neck fractures (mean ± SD age: 70.6 ± 12.7 years) and 21 with trochanteric fractures (75.6 ± 9.3 years), and 50 control subjects (63.7 ± 7.0 years) were included for the comparisons. Differences in the spatial distributions of volumetric bone mineral density, cortical bone thickness, cortical volumetric bone mineral density, and volumetric bone mineral density in a layer adjacent to the endosteal surface were investigated using voxel-based morphometry (VBM) and surface-based statistical parametric mapping (SPM). We compared these spatial distributions between controls and both types of fracture, and between the two types of fracture. Using VBM, we found spatially heterogeneous volumetric bone mineral density differences between control subjects and subjects with hip fracture that varied by fracture type. Interestingly, femoral neck fracture subjects, but not subjects with trochanteric fracture, showed significantly lower volumetric bone mineral density in the superior aspect of the femoral neck compared with controls. Using surface-based SPM, we found that compared with controls, both fracture types showed thinner cortices in regions in agreement with the type of fracture. Most outcomes of cortical and endocortical volumetric bone mineral density comparisons were consistent with VBM results. Our results suggest: 1) that the spatial distribution of trabecular volumetric bone mineral density might play a significant role in hip fracture; 2) that focal cortical bone thinning might be more relevant in femoral neck fractures; and 3) that areas of reduced cortical and endocortical volumetric bone mineral density might be more relevant for trochanteric fractures in Chinese women. © 2017 American Society for Bone and Mineral Research.

Trabecular and cortical bone structure of the talus and distal tibia in Pan and Homo

OBJECTIVES

Internal bone structure, both cortical and trabecular bone, remodels in response to loading and may provide important information regarding behavior. The foot is well suited to analysis of internal bone structure because it experiences the initial substrate reaction forces, due to its proximity to the substrate. Moreover, as humans and apes differ in loading of the foot, this region is relevant to questions concerning arboreal locomotion and bipedality in the hominoid fossil record.

MATERIALS AND METHODS

We apply a whole-bone/epiphysis approach to analyze trabecular and cortical bone in the distal tibia and talus of Pan troglodytes and Homo sapiens. We quantify bone volume fraction (BV/TV), degree of anisotropy (DA), trabecular thickness (Tb.Th), bone surface to volume ratio (BS/BV), and cortical thickness and investigate the distribution of BV/TV and cortical thickness throughout the bone/epiphysis.

RESULTS

We find that Pan has a greater BV/TV, a lower BS/BV and thicker cortices than Homo in both the talus and distal tibia. The trabecular structure of the talus is more divergent than the tibia, having thicker, less uniformly aligned trabeculae in Pan compared to Homo. Differences in dorsiflexion at the talocrural joint and in degree of mobility at the talonavicular joint are reflected in the distribution of cortical and trabecular bone.

DISCUSSION

Overall, quantified trabecular parameters represent overall differences in bone strength between the two species, however, DA may be directly related to joint loading. Cortical and trabecular bone distributions correlate with habitual joint positions adopted by each species, and thus have potential for interpreting joint position in fossil hominoids.

Focal osteoporosis defects play a key role in hip fracture

BACKGROUND

Hip fractures are mainly caused by accidental falls and trips, which magnify forces in well-defined areas of the proximal femur. Unfortunately, the same areas are at risk of rapid bone loss with ageing, since they are relatively stress-shielded during walking and sitting. Focal osteoporosis in those areas may contribute to fracture, and targeted 3D measurements might enhance hip fracture prediction. In the FEMCO case-control clinical study, Cortical Bone Mapping (CBM) was applied to clinical computed tomography (CT) scans to define 3D cortical and trabecular bone defects in patients with acute hip fracture compared to controls. Direct measurements of trabecular bone volume were then made in biopsies of target regions removed at operation.

METHODS

The sample consisted of CT scans from 313 female and 40 male volunteers (158 with proximal femoral fracture, 145 age-matched controls and 50 fallers without hip fracture). Detailed Cortical Bone Maps (c.5580 measurement points on the unfractured hip) were created before registering each hip to an average femur shape to facilitate statistical parametric mapping (SPM). Areas where cortical and trabecular bone differed from controls were visualised in 3D for location, magnitude and statistical significance. Measures from the novel regions created by the SPM process were then tested for their ability to classify fracture versus control by comparison with traditional CT measures of areal Bone Mineral Density (aBMD). In women we used the surgical classification of fracture location ('femoral neck' or 'trochanteric') to discover whether focal osteoporosis was specific to fracture type. To explore whether the focal areas were osteoporotic by histological criteria, we used micro CT to measure trabecular bone parameters in targeted biopsies taken from the femoral heads of 14 cases.

RESULTS

Hip fracture patients had distinct patterns of focal osteoporosis that determined fracture type, and CBM measures classified fracture type better than aBMD parameters. CBM measures however improved only minimally on aBMD for predicting any hip fracture and depended on the inclusion of trabecular bone measures alongside cortical regions. Focal osteoporosis was confirmed on biopsy as reduced sub-cortical trabecular bone volume.

CONCLUSION

Using 3D imaging methods and targeted bone biopsy, we discovered focal osteoporosis affecting trabecular and cortical bone of the proximal femur, among men and women with hip fracture.

Association between femur size and a focal defect of the superior femoral neck

Within each sex, there is an association between hip fracture risk and the size of the proximal femur, with larger femurs apparently more susceptible to fracture. Here, we investigate whether the thickness and density of the femoral cortex play a role in this association: might larger femurs harbour focal, cortical defects? To answer this question, we used cortical bone mapping to measure the distribution of cortical mass surface density (CMSD, mg/cm(2)) in cohorts of 308 males and 125 females. Principal component analysis of the various femoral surfaces led to a measure of size that is linearly independent from shape. After mapping the data onto a canonical femur surface, we used statistical parametric mapping to identify any regions where CMSD depends on size, allowing for other confounding covariates including shape. Our principal finding was a focal patch on the superior femoral neck, where CMSD is reduced by around 1% for each 1% increase in proximal-distal size (p<0.000005 in the males, p<0.001 in the females). This finding appears to be consistent with models of functional adaptation, and may help with the design of interventional strategies for reducing fracture risk.

Predicting Hip Fracture Type With Cortical Bone Mapping (CBM) in the Osteoporotic Fractures in Men (MrOS) Study

Hip fracture risk is known to be related to material properties of the proximal femur, but fracture prediction studies adding richer quantitative computed tomography (QCT) measures to dual-energy X-ray (DXA)-based methods have shown limited improvement. Fracture types have distinct relationships to predictors, but few studies have subdivided fracture into types, because this necessitates regional measurements and more fracture cases. This work makes use of cortical bone mapping (CBM) to accurately assess, with no prior anatomical presumptions, the distribution of properties related to fracture type. CBM uses QCT data to measure the cortical and trabecular properties, accurate even for thin cortices below the imaging resolution. The Osteoporotic Fractures in Men (MrOS) study is a predictive case-cohort study of men over 65 years old: we analyze 99 fracture cases (44 trochanteric and 55 femoral neck) compared to a cohort of 308, randomly selected from 5994. To our knowledge, this is the largest QCT-based predictive hip fracture study to date, and the first to incorporate CBM analysis into fracture prediction. We show that both cortical mass surface density and endocortical trabecular BMD are significantly different in fracture cases versus cohort, in regions appropriate to fracture type. We incorporate these regions into predictive models using Cox proportional hazards regression to estimate hazard ratios, and logistic regression to estimate area under the receiver operating characteristic curve (AUC). Adding CBM to DXA-based BMD leads to a small but significant (p < 0.005) improvement in model prediction for any fracture, with AUC increasing from 0.78 to 0.79, assessed using leave-one-out cross-validation. For specific fracture types, the improvement is more significant (p < 0.0001), with AUC increasing from 0.71 to 0.77 for trochanteric fractures and 0.76 to 0.82 for femoral neck fractures. In contrast, adding DXA-based BMD to a CBM-based predictive model does not result in any significant improvement.

Independent measurement of femoral cortical thickness and cortical bone density using clinical CT

The local structure of the proximal femoral cortex is of interest since both fracture risk, and the effects of various interventions aimed at reducing that risk, are associated with cortical properties focused in particular regions rather than dispersed over the whole bone. Much of the femoral cortex is less than 3mm thick, appearing so blurred in clinical CT that its actual density is not apparent in the data, and neither thresholding nor full-width half-maximum techniques are capable of determining its width. Our previous work on cortical bone mapping showed how to produce more accurate estimates of cortical thickness by assuming a fixed value of the cortical density for each hip. However, although cortical density varies much less over the proximal femur than thickness, what little variation there is leads to errors in thickness measurement. In this paper, we develop the cortical bone mapping technique by exploiting local estimates of imaging blur to correct the global density estimate, thus providing a local density estimate as well as more accurate estimates of thickness. We also consider measurement of cortical mass surface density and the density of trabecular bone immediately adjacent to the cortex. Performance is assessed with ex vivo clinical QCT scans of proximal femurs, with true values derived from high resolution HRpQCT scans of the same bones. We demonstrate superior estimation of thickness than is possible with alternative techniques (accuracy 0.12 ± 0.39 mm for cortices in the range 1-3mm), and that local cortical density estimation is feasible for densities >800 mg/cm(3).

Spatially varying registration using Gaussian processes

In this paper we propose a new approach for spatially-varying registration using Gaussian process priors. The method is based on the idea of spectral tempering, i.e. the spectrum of the Gaussian process is modified depending on a user defined tempering function. The result is a non-stationary Gaussian process, which induces different amount of smoothness in different areas. In contrast to most other schemes for spatially-varying registration, our approach does not require any change in the registration algorithm itself, but only affects the prior model. Thus we can obtain spatially-varying versions of any registration method whose deformation prior can be formulated in terms of a Gaussian process. This includes for example most spline-based models, but also statistical shape or deformation models. We present results for the problem of atlas based skull-registration of cone beam CT images. These datasets are difficult to register as they contain a large amount of noise around the teeth. We show that with our method we can become robust against noise, but still obtain accurate correspondence where the data is clean.