The utility of multi-stack alignment and 3D longitudinal image registration to assess bone remodeling in rheumatoid arthritis patients from second generation HR-pQCT scans

Evaluating high-resolution computed tomography derived 3-D joint space metrics of the metacarpophalangeal joints between rheumatoid arthritis and age- and sex-matched control participants

Introduction

Rheumatoid arthritis (RA) is commonly characterized by joint space narrowing. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides unparalleled in vivo visualization and quantification of joint space in extremity joints commonly affected by RA, such as the 2nd and 3rd metacarpophalangeal joints. However, age, sex, and obesity can also influence joint space narrowing. Thus, this study aimed to determine whether HR-pQCT joint space metrics could distinguish between RA patients and controls, and determine the effects of age, sex and body mass index (BMI) on these joint space metrics.

Methods

HR-pQCT joint space metrics (volume, width, standard deviation of width, maximum/minimum width, and asymmetry) were acquired from RA patients and age-and sex-matched healthy control participants 2nd and 3rd MCP joints. Joint health and functionality were assessed with ultrasound (i.e., effusion and inflammation), hand function tests, and questionnaires.

Results

HR-pQCT-derived 3D joint space metrics were not significantly different between RA and control groups (p > 0.05), despite significant differences in inflammation and joint function (p < 0.05). Joint space volume, mean joint space width (JSW), maximum JSW, minimum JSW were larger in males than females (p < 0.05), while maximum JSW decreased with age. No significant association between joint space metrics and BMI were found.

Conclusion

HR-pQCT did not detect group level differences between RA and age-and sex-matched controls. Further research is necessary to determine whether this is due to a true lack of group level differences due to well-controlled RA, or the inability of HR-pQCT to detect a difference.

A comprehensive set of ultrashort echo time magnetic resonance imaging biomarkers to assess cortical bone health: A feasibility study at clinical field strength

INTRODUCTION

Conventional bone imaging methods primarily use X-ray techniques to assess bone mineral density (BMD), focusing exclusively on the mineral phase. This approach lacks information about the organic phase and bone water content, resulting in an incomplete evaluation of bone health. Recent research highlights the potential of ultrashort echo time magnetic resonance imaging (UTE MRI) to measure cortical porosity and estimate BMD based on signal intensity. UTE MRI also provides insights into bone water distribution and matrix organization, enabling a comprehensive bone assessment with a single imaging technique. Our study aimed to establish quantifiable UTE MRI-based biomarkers at clinical field strength to estimate BMD and microarchitecture while quantifying bound water content and matrix organization.

METHODS

Femoral bones from 11 cadaveric specimens (n = 4 males 67-92 yrs of age, n = 7 females 70-95 yrs of age) underwent dual-echo UTE MRI (3.0 T, 0.45 mm resolution) with different echo times and high resolution peripheral quantitative computed tomography (HR-pQCT) imaging (60.7 μm voxel size). Following registration, a 4.5 mm HR-pQCT region of interest was divided into four quadrants and used across the multi-modal images. Statistical analysis involved Pearson correlation between UTE MRI porosity index and a signal-intensity technique used to estimate BMD with corresponding HR-pQCT measures. UTE MRI was used to calculate T1 relaxation time and a novel bound water index (BWI), compared across subregions using repeated measures ANOVA.

RESULTS

The UTE MRI-derived porosity index and signal-intensity-based estimated BMD correlated with the HR-pQCT variables (porosity: r = 0.73, p = 0.006; BMD: r = 0.79, p = 0.002). However, these correlations varied in strength when we examined each of the four quadrants (subregions, r = 0.11-0.71). T1 relaxometry and the BWI exhibited variations across the four subregions, though these differences were not statistically significant. Notably, we observed a strong negative correlation between T1 relaxation time and the BWI (r = -0.87, p = 0.0006).

CONCLUSION

UTE MRI shows promise for being an innocuous method for estimating cortical porosity and BMD parameters while also giving insight into bone hydration and matrix organization. This method offers the potential to equip clinicians with a more comprehensive array of imaging biomarkers to assess bone health without the need for invasive or ionizing procedures.

Fixed and Relative Positioning of Scans for High Resolution Peripheral Quantitative Computed Tomography

INTRODUCTION

High resolution peripheral quantitative computed tomography (HR-pQCT) imaging protocol requires defining where to position the ∼1 cm thick scan along the bone length. Discrepancies between the use of two positioning methods, the relative and fixed offset, may be problematic in the comparison between studies and participants. This study investigated how bone landmarks scale linearly with length and how this scaling affects both positioning methods aimed at providing a consistent anatomical location for scan acquisition.

METHODS

Using CT images of the radius (N = 25) and tibia (N = 42), 10 anatomical landmarks were selected along the bone length. The location of these landmarks was converted to a percent length along the bone, and the variation in their location was evaluated across the dataset. The absolute location of the HR-pQCT scan position using both offset methods was identified for all bones and converted to a percent length position relative to the HR-pQCT reference line for comparison. A secondary analysis of the location of the scan region specifically within the metaphysis was explored at the tibia.

RESULTS

The location of landmarks deviated from a linear relationship across the dataset, with a range of 3.6 % at the radius sites, and 4.5 % at the tibia sites. The consequent variation of the position of the scan at the radius was 0.6 % and 0.3 %, and at the tibia 2.4 % and 0.5 %, for the fixed and relative offset, respectively. The position of the metaphyseal junction with the epiphysis relative to the scan position was poorly correlated to bone length, with R2 = 0.06 and 0.37, for the fixed and relative offset respectively.

CONCLUSION

The variation of the scan position by either method is negated by the intrinsic variation of the bone anatomy with respect both to total bone length as well as the metaphyseal region. Therefore, there is no clear benefit of either offset method. However, the lack of difference due to the inherent variation in the underlying anatomy implies that it is reasonable to compare studies even if they are using different positioning methods.

A multi-stack registration technique to improve measurement accuracy and precision across longitudinal HR-pQCT scans

BACKGROUND

Recent applications of high-resolution peripheral quantitative computed tomography (HR-pQCT) have demonstrated that changes in local bone remodelling can be quantified in vivo using longitudinal three-dimensional image registration. However, certain emerging applications, such as fracture healing and joint analysis, require larger multi-stack scan regions that can result in stack shift image artifacts. These artifacts can be detrimental to the accurate alignment of the bone structure across multiple timepoints. The purpose of this study was to establish a multi-stack registration protocol for evaluating longitudinal HR-pQCT images and to assess the accuracy and precision error in comparison with measures obtained using previously established three-dimensional longitudinal registration.

METHODS

Three same day multi-stack HR-pQCT scans of the radius (2 stacks in length) and tibia (3 stacks in length) were obtained from 39 healthy individuals who participated in a previous reproducibility study. A fully automated multi-stack registration algorithm was developed to re-align stacks within a scan by leveraging slight offsets between longitudinal scans. Stack shift severity before and after registration was quantified using a newly proposed stack-shift severity score. The false discovery rate for bone remodelling events and precision error of bone morphology and micro-finite element analysis parameters were compared between longitudinally registered scans with and without the addition of multi-stack registration.

RESULTS

Most scans (82 %) improved in stack alignment or maintained the lowest stack shift severity score when multi-stack registration was implemented. The false discovery rate of bone remodelling events significantly decreased after multi-stack registration, resulting in median false detection of bone formation and resorption fractions between 3.2 to 7.5 % at the radius and 3.4 to 5.3 % at the tibia. Further, precision error was significantly reduced or remained unchanged in all standard bone morphology and micro-finite element analysis parameters, except for total and trabecular cross-sectional areas.

CONCLUSION

Multi-stack registration is an effective strategy for accurately aligning multi-stack HR-pQCT scans without modification of the image acquisition protocol. The algorithm presented here is a viable approach for performing accurate morphological analysis on multi-stack HR-pQCT scans, particularly for advanced application investigating local bone remodelling in vivo.

Image Registration in Longitudinal Bone Assessment Using Computed Tomography

PURPOSE OF REVIEW

Rigid image registration is an important image processing tool for the assessment of musculoskeletal chronic disease. In this paper, we critically review applications of rigid image registration in terms of similarity measurement methods over the past three years (2019-2022) in the context of monitoring longitudinal changes to bone microstructure and mechanical properties using computed tomography. This review identifies critical assumptions and trade-offs underlying different similarity measurement methods used in image registration and demonstrates the effect of using different similarity measures on registration outcomes.

RECENT FINDINGS

Image registration has been used in recent studies for: correcting positional shifts between longitudinal scans to quantify changes to bone microstructural and mechanical properties over time, developing registration-based workflows for longitudinal assessment of bone properties in pre-clinical and clinical studies, and developing and validating registration techniques for longitudinal studies. In evaluating the recent literature, it was found that the assumptions at the root of different similarity measures used in rigid image registration are not always confirmed and reported. Each similarity measurement has its advantages and disadvantages, as well as underlying assumptions. Breaking these assumptions can lead to poor and inaccurate registration results. Thus, care must be taken with regards to the choice of similarity measurement and interpretation of results. We propose that understanding and verifying the assumptions of similarity measurements will enable more accurate and efficient quantitative assessments of structural changes over time.

Open-source image analysis tool for the identification and quantification of cortical interruptions and bone erosions in high-resolution peripheral quantitative computed tomography images of patients with rheumatoid arthritis

Identification of bone erosions and quantification of erosion volume is important for rheumatoid arthritis diagnosis, and can add important information to evaluate disease progression and treatment effects. High-resolution peripheral quantitative computed tomography (HR-pQCT) is well suited for this purpose, however analysis methods are not widely available. The purpose of this study was to develop an open-source software tool for the identification and quantification of bone erosions using images acquired by HR-pQCT. The collection of modules, Bone Analysis Modules (BAM) - Erosion, implements previously published erosion analysis techniques as modules in 3D Slicer, an open-source image processing and visualization tool. BAM includes a module to automatically identify cortical interruptions, from which erosions are manually selected, and a hybrid module that combines morphological and level set operations to quantify the volume of bone erosions. HR-pQCT images of the second and third metacarpophalangeal (MCP) joints were acquired in patients with RA (XtremeCT, n = 14, XtremeCTII, n = 22). The number of cortical interruptions detected by BAM-Erosion agreed strongly with the previously published cortical interruption detection algorithm for both XtremeCT (r² = 0.85) and XtremeCTII (r² = 0.87). Erosion volume assessment by BAM-Erosion agreed strongly (r² = 0.95) with the Medical Image Analysis Framework. BAM-Erosion provides an open-source erosion analysis tool that produces comparable results to previously published algorithms, with improved options for visualization. The strength of the tool is that it implements multiple image processing algorithms for erosion analysis on a single, widely available, open-source platform that can accommodate future updates.

Clinical observation of diminished bone quality and quantity through longitudinal HR-pQCT-derived remodeling and mechanoregulation

High resolution peripheral quantitative computed tomography (HR-pQCT) provides methods for quantifying volumetric bone mineral density and microarchitecture necessary for early diagnosis of bone disease. When combined with a longitudinal imaging protocol and finite element analysis, HR-pQCT can be used to assess bone formation and resorption (i.e., remodeling) and the relationship between this remodeling and mechanical loading (i.e., mechanoregulation) at the tissue level. Herein, 25 patients with a contralateral distal radius fracture were imaged with HR-pQCT at baseline and 9-12 months follow-up: 16 patients were prescribed vitamin D3 with/without calcium supplement based on a blood biomarker measures of bone metabolism and dual-energy X-ray absorptiometry image-based measures of normative bone quantity which indicated diminishing (n = 9) or poor (n = 7) bone quantity and 9 were not. To evaluate the sensitivity of this imaging protocol to microstructural changes, HR-pQCT images were registered for quantification of bone remodeling and image-based micro-finite element analysis was then used to predict local bone strains and derive rules for mechanoregulation. Remodeling volume fractions were predicted by both average values of trabecular and cortical thickness and bone mineral density (R² > 0.8), whereas mechanoregulation was affected by dominance of the arm and group classification (p < 0.05). Overall, longitudinal, extended HR-pQCT analysis enabled the identification of changes in bone quantity and quality too subtle for traditional measures.

Advancements in Osteoporosis Imaging, Screening, and Study of Disease Etiology

PURPOSE OF REVIEW

The purpose of this review is to inform researchers and clinicians with the most recent imaging techniques that are employed (1) to opportunistically screen for osteoporosis and (2) to provide a better understanding into the disease etiology of osteoporosis.

RECENT FINDINGS

Phantomless calibration techniques for computed tomography (CT) may pave the way for better opportunistic osteoporosis screening and the retroactive analysis of imaging data. Additionally, hardware advances are enabling new applications of dual-energy CT and cone-beam CT to the study of bone. Advances in MRI sequences are also improving imaging evaluation of bone properties. Finally, the application of image registration techniques is enabling new uses of imaging to investigate soft tissue-bone interactions as well as bone turnover. While DXA remains the most prominent imaging tool for osteoporosis diagnosis, new imaging techniques are becoming more widely available and providing additional information to inform clinical decision-making.

Formation Dominates Resorption With Increasing Mineralized Density and Time Postfracture in Cortical but Not Trabecular Bone: A Longitudinal HRpQCT Imaging Study in the Distal Radius

Clinical evaluation of fracture healing is often limited to an assessment of fracture bridging from radiographic images, without consideration for other aspects of bone quality. However, recent advances in HRpQCT offer methods to accurately monitor microstructural bone remodeling throughout the healing process. In this study, local bone formation and resorption were investigated during the first year post fracture in both the fractured (n = 22) and contralateral (n = 19) radii of 34 conservatively treated patients (24 female, 10 male) who presented with a unilateral radius fracture at the Innsbruck University Hospital, Austria. HRpQCT images and clinical metrics were acquired at six time points for each patient. The standard HRpQCT image acquisition was captured for all radii, with additional distal and proximal image acquisitions for the fractured radii. Measured radial bone densities were isolated with a voxel‐based mask and images were rigidly registered to images from the previous imaging session using a pyramid‐based approach. From the registered images, bone formation and resorption volume fractions were quantified for multiple density‐based thresholds and compared between the fractured and contralateral radius and relative to demographics, bone morphometrics, and fracture metrics using regression. Compared with the contralateral radius, both bone formation and resorption were significantly increased in the fractured radius throughout the study for nearly all evaluated thresholds. Higher density cortical bone formation continually increased throughout the duration of the study and was significantly greater than resorption during late‐stage healing in both the fractured and intact regions of the radius. With the small and diverse study population, only weak relationships between fracture remodeling and patient‐specific parameters were unveiled. However this study provides methods for the analysis of local bone remodeling during fracture healing and highlights relevant considerations for future studies, specifically that remodeling postfracture is likely to continue beyond 12‐months postfracture. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Heterogenous bone response to biologic DMARD therapies in rheumatoid arthritis patients and their relationship to functional indices

Objectives: Previous studies of high-resolution peripheral quantitative computed tomography (HR-pQCT) imaging of hand joints in patients with rheumatoid arthritis (RA) have suggested that erosion healing may occur. Our objective was to examine changes in erosion volume, joint space width (JSW), bone mineral density (BMD), and bone remodelling, and their association with clinical outcomes and measures of patient hand function.Method: We examined 48 patients who achieved a good response to a newly initiated biologic therapy. HR-pQCT images of the dominant hands' second and third metacarpophalangeal joints were obtained 3 and 12 months after therapy initiation. Bone erosion volume, JSW, BMD, and bone remodelling were quantified from HR-pQCT images, with improvement, no change (unchanged), or progression in these measures determined by least significant change. Disease activity and hand function measures were collected.Results: There were no significant group changes in HR-pQCT outcomes over the 9 month period. Twenty-two patients had total erosion volumes that remained unchanged, nine showed improvement, and two progressed. The majority of JSW and BMD measures remained unchanged. There was a significant association between the baseline Health Assessment Questionnaire score and the change in minimum JSW, but no other significant associations between HR-pQCT outcomes and function were observed.Conclusions: The vast majority of patients maintained unchanged JSW and BMD over the course of follow-up. Significant improvements in total erosion volume occurred in 27% of patients, suggesting that biologic therapies may lead to erosion healing in some patients, although this did not have an impact on self-reported and demonstrated hand function.

Longitudinal time-lapse in vivo micro-CT reveals differential patterns of peri-implant bone changes after subclinical bacterial infection in a rat model

Subclinical infection associated with orthopedic devices can be challenging to diagnose. The goal of this study was to evaluate longitudinal, microcomputed tomography (microCT) imaging in a rat model of subclinical orthopedic device-related infection caused by Staphylococcus epidermidis and four different Cutibacterium (previously Propionibacterium) acnes strains, and compare outcomes with non-inoculated and historical S. aureus-inoculated controls. Sterile screws or screws colonized with bacteria were placed in the tibia of 38 adult Wistar rats [n = 6 sterile screws; n = 6 S. epidermidis-colonized screws; n = 26 C. acnes-colonized screws (covering all three main subspecies)]. Regular microCT scans were taken over 28 days and processed for quantitative time-lapse imaging with dynamic histomorphometry. At euthanasia, tissues were processed for semiquantitative histopathology or quantitative bacteriology. All rats receiving sterile screws were culture-negative at euthanasia and displayed progressive bony encapsulation of the screw. All rats inoculated with S. epidermidis-colonized screws were culture-positive and displayed minor changes in peri-implant bone, characteristic of subclinical infection. Five of the 17 rats in the C. acnes inoculated group were culture positive at euthanasia and displayed bone changes at the interface of the screw and bone, but not deeper in the peri-implant bone. Dynamic histomorphometry revealed significant differences in osseointegration, bone remodeling and periosteal reactions between groups that were not measurable by visual observation of still microCT images. Our study illustrates the added value of merging 3D microCT data from subsequent timepoints and producing inherently richer 4D data for the detection and characterization of subclinical orthopedic infections, whilst also reducing animal use.

Multi-Modal Imaging to Assess the Interaction Between Inflammation and Bone Damage Progression in Inflammatory Arthritis

Combining results from multiple imaging techniques (i.e., multi-modal imaging) through image registration can result in the better characterization of joint tissue characteristics. In the context of inflammatory arthritis conditions, high-resolution peripheral quantitative computed tomography (HR-pQCT) provides excellent bone contrast while magnetic resonance imaging (MRI) provides superior contrast and resolution of soft tissue and inflammatory characteristics. Superimposing these imaging results upon each other provides a robust characterization of the joint. In a preliminary study of nine rheumatoid arthritis (RA) participants in clinical remission, we acquired HR-pQCT and MR images of their 2nd and 3rd metacarpophalangeal (MCP) joints at two timepoints 6 months apart. We present the benefits of a multi-modal imaging approach, in which we demonstrate the ability to localize regions of inflammation with subtle changes in bone erosion volume. Using HR-pQCT and MRI to visualize bone damage and inflammation, respectively, will improve our understanding of the impact that subclinical inflammation has on bone damage progression, and demonstrating if bone repair occurs where inflammation is resolved. The presented multi-modal imaging technique has the potential to study the progression of bone damage in relation to inflammation that otherwise would not be possible with either imaging technique alone. The multi-modal image registration technique will be helpful to understanding the development and pathogenesis of RA-associated bone erosions. Additionally, multi-modal imaging may provide a technique to probe the tissue-level changes that occur as a result of treatment regimes.

High-Resolution Peripheral Quantitative Computed Tomography for Bone Evaluation in Inflammatory Rheumatic Disease

High resolution peripheral quantitative computed tomography (HR-pQCT) is a 3-dimensional imaging modality with superior sensitivity for bone changes and abnormalities. Recent advances have led to increased use of HR-pQCT in inflammatory arthritis to report quantitative volumetric measures of bone density, microstructure, local anabolic (e.g., osteophytes, enthesiophytes) and catabolic (e.g., erosions) bone changes and joint space width. These features may be useful for monitoring disease progression, response to therapy, and are responsive to differentiating between those with inflammatory arthritis conditions and healthy controls. We reviewed 69 publications utilizing HR-pQCT imaging of the metacarpophalangeal (MCP) and/or wrist joints to investigate arthritis conditions. Erosions are a marker of early inflammatory arthritis progression, and recent work has focused on improvement and application of techniques to sensitively identify erosions, as well as quantifying erosion volume changes longitudinally using manual, semi-automated and automated methods. As a research tool, HR-pQCT may be used to detect treatment effects through changes in erosion volume in as little as 3 months. Studies with 1-year follow-up have demonstrated progression or repair of erosions depending on the treatment strategy applied. HR-pQCT presents several advantages. Combined with advances in image processing and image registration, individual changes can be monitored with high sensitivity and reliability. Thus, a major strength of HR-pQCT is its applicability in instances where subtle changes are anticipated, such as early erosive progression in the presence of subclinical inflammation. HR-pQCT imaging results could ultimately impact decision making to uptake aggressive treatment strategies and prevent progression of joint damage. There are several potential areas where HR-pQCT evaluation of inflammatory arthritis still requires development. As a highly sensitive imaging technique, one of the major challenges has been motion artifacts; motion compensation algorithms should be implemented for HR-pQCT. New research developments will improve the current disadvantages including, wider availability of scanners, the field of view, as well as the versatility for measuring tissues other than only bone. The challenge remains to disseminate these analysis approaches for broader clinical use and in research.