ORMIR_XCT: A Python package for high resolution peripheral quantitative computed tomography image processing

High resolution peripheral quantitative computed tomography (HR-pQCT) is an imaging technique capable of imaging trabecular bone in-vivo. HR-pQCT has a wide range of applications, primarily focused on bone to improve our understanding of musculoskeletal diseases, assess epidemiological associations, and evaluate the effects of pharmaceutical interventions. Processing HR-pQCT images has largely been supported using the scanner manufacturer scripting language (Image Processing Language, IPL, Scanco Medical). However, by expanding image processing workflows outside of the scanner manufacturer software environment, users have the flexibility to apply more advanced mathematical techniques and leverage modern software packages to improve image processing. The ORMIR_XCT Python package was developed to reimplement some existing IPL workflows and provide an open and reproducible package allowing for the development of advanced HR-pQCT data processing workflows

The daily association between affect and alcohol use: a meta-analysis of individual participant data

Influential psychological theories hypothesize that people consume alcohol in response to the experience of both negative and positive emotions. Despite two decades of daily diary and ecological momentary assessment research, it remains unclear whether people consume more alcohol on days they experience higher negative and positive affect in everyday life. In this preregistered meta-analysis, we synthesized the evidence for these daily associations between affect and alcohol use. We included individual participant data from 69 studies (N = 12,394), which used daily and momentary surveys to assess affect and the number of alcoholic drinks consumed. Results indicate that people are not more likely to drink on days they experience high negative affect, but are more likely to drink and drink heavily on days high in positive affect. People self-reporting a motivational tendency to drink-to-cope and drink-to-enhance consumed more alcohol, but not on days they experienced higher negative and positive affect. Results were robust across different operationalizations of affect, study designs, study populations, and individual characteristics. These findings challenge the long-held belief that people drink more alcohol following increases in negative affect. Integrating these findings under different theoretical models and limitations of this field of research, we collectively propose an agenda for future research to explore open questions surrounding affect and alcohol use.The present study was funded by the Canadian Institutes of Health Research Grant MOP-115104 (Roisin M. O’Connor), Canadian Institutes of Health Research Grant MSH-122803 (Roisin M. O’Connor), John A. Hartford Foundation Grant (Paul Sacco), Loyola University Chicago Research Support Grant (Tracy De Hart), National Institute for Occupational Safety and Health Grant T03OH008435 (Cynthia Mohr), National Institutes of Health (NIH) Grant F31AA023447 (Ryan W. Carpenter), NIH Grant R01AA025936 (Kasey G. Creswell), NIH Grant R01AA025969 (Catharine E. Fairbairn), NIH Grant R21AA024156 (Anne M. Fairlie), NIH Grant F31AA024372 (Fallon Goodman), NIH Grant R01DA047247 (Kevin M. King), NIH Grant K01AA026854 (Ashley N. Linden-Carmichael), NIH Grant K01AA022938 (Jennifer E. Merrill), NIH Grant K23AA024808 (Hayley Treloar Padovano), NIH Grant P60AA11998 (Timothy Trull), NIH Grant MH69472 (Timothy Trull), NIH Grant K01DA035153 (Nisha Gottfredson), NIH Grant P50DA039838 (Ashley N. Linden-Carmichael), NIH Grant K01DA047417 (David M. Lydon-Staley), NIH Grant T32DA037183 (M. Kushner), NIH Grant R21DA038163 (A. Moore), NIH Grant K12DA000167 (M. Potenza, Stephanie S. O’Malley), NIH Grant R01AA025451 (Bruce Bartholow, Thomas M. Piasecki), NIH Grant P50AA03510 (V. Hesselbrock), NIH Grant K01AA13938 (Kristina M. Jackson), NIH Grant K02AA028832 (Kevin M. King), NIH Grant T32AA007455 (M. Larimer), NIH Grant R01AA025037 (Christine M. Lee, M. Patrick), NIH Grant R01AA025611 (Melissa Lewis), NIH Grant R01AA007850 (Robert Miranda), NIH Grant R21AA017273 (Robert Miranda), NIH Grant R03AA014598 (Cynthia Mohr), NIH Grant R29AA09917 (Cynthia Mohr), NIH Grant T32AA07290 (Cynthia Mohr), NIH Grant P01AA019072 (P. Monti), NIH Grant R01AA015553 (J. Morgenstern), NIH Grant R01AA020077 (J. Morgenstern), NIH Grant R21AA017135 (J. Morgenstern), NIH Grant R01AA016621 (Stephanie S. O’Malley), NIH Grant K99AA029459 (Marilyn Piccirillo), NIH Grant F31AA022227 (Nichole Scaglione), NIH Grant R21AA018336 (Katie Witkiewitz), Portuguese State Budget Foundation for Science and Technology Grant UIDB/PSI/01662/2020 (Teresa Freire), University of Washington Population Health COVID-19 Rapid Response Grant (J. Kanter, Adam M. Kuczynski), U.S. Department of Defense Grant W81XWH-13-2-0020 (Cynthia Mohr), SANPSY Laboratory Core Support Grant CNRS USR 3413 (Marc Auriacombe), Social Sciences and Humanities Research Council of Canada Grant (N. Galambos), and Social Sciences and Humanities Research Council of Canada Grant (Andrea L. Howard)

Reproducibility and repeatability of a semi-automated pipeline to quantify trapeziometacarpal joint angles using dynamic computed tomography

Abstract Background The trapeziometacarpal (TMC) joint is a mechanically complex joint and is commonly affected by musculoskeletal diseases such as osteoarthritis. Quantifying in vivo TMC joint biomechanics, such as joint angles, with traditional reflective marker-based methods can be difficult due to the joint’s location in the hand. Dynamic computed tomography (CT) can facilitate the quantification of TMC joint motion by continuously capturing three-dimensional volumes over time. However, post-processing of dynamic CT datasets can be time intensive and automated methods are needed to reduce processing times to allow for application to larger clinical studies. The purpose of this work is to introduce a fast, semi-automated pipeline to quantify joint angles from dynamic CT scans of the TMC joint and evaluate the associated error in joint angle and translation computation by means of a reproducibility and repeatability study. Methods Ten cadaveric hands were scanned with dynamic CT using a passive motion device to move thumbs in a radial abduction–adduction motion. Static CT scans and high-resolution peripheral quantitative CT scans were also acquired to generate high-resolution bone meshes. Abduction–adduction, flexion–extension, and axial rotation angles were computed using a joint coordinate system. Reproducibility and repeatability were assessed using intraclass correlation coefficients, Bland–Altman analysis, and root mean square errors. Target registration errors were computed to evaluate errors associated with image registration. Results We found good repeatability for flexion–extension, abduction–adduction, and axial rotation angles. Reproducibility was moderate for all three angles. Joint translations exhibited greater repeatability than reproducibility. Specimens with greater joint degeneration had lower repeatability and reproducibility. We found that the difference in resulting joint angles and translations were likely due to differences in segment coordinate system definition between multiple raters, rather than due to registration errors. Conclusions The proposed semi-automatic processing pipeline was fast, repeatable, and moderately reproducible when quantifying TMC joint angles and translations. This work provides a range of errors for TMC joint angles from dynamic CT scans using manually selected anatomical landmarks

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

Abstract Background Medical imaging plays an important role in determining the progression of joint damage in rheumatoid arthritis (RA). High resolution peripheral quantitative computed tomography (HR-pQCT) is a sensitive tool capable of evaluating bone microarchitecture and erosions, and 3D rigid image registration can be used to visualize and quantify bone remodeling over time. However, patient motion during image acquisition can cause a “stack shift” artifact resulting in loss of information and reducing the number of erosions that can be analyzed using HR-pQCT. The purpose of this study was to use image registration to improve the number of useable HR-pQCT scans and to apply image-based bone remodeling assessment to the metacarpophalangeal (MCP) joints of RA patients. Methods Ten participants with RA completed HR-pQCT scans of the 2nd and 3rd MCP joints at enrolment to the study and at a 6-month follow-up interval. At 6-months, an additional repeat scan was acquired to evaluate reliability. HR-pQCT images were acquired in three individual 1 cm acquisitions (stacks) with a 25% overlap. We completed analysis first using standard evaluation methods, and second with multi-stack registration. We assessed whether additional erosions could be evaluated after multi-stack registration. Bone remodeling analysis was completed using registration and transformation of baseline and follow-up images. We calculated the bone formation and resorption volume fractions with 6-month follow-up, and same-day repositioning as a negative control. Results 13/57 (23%) of erosions could not be analyzed from raw images due to a stack shift artifact. All erosions could be volumetrically assessed after multi-stack registration. We observed that there was a median bone formation fraction of 2.1% and resorption fraction of 3.8% in RA patients over the course of 6 months. In contrast to the same-day rescan negative control, we observed median bone formation and resorption fractions of 0%. Conclusions Multi-stack image registration is a useful tool to improve the number of useable scans when analyzing erosions using HR-pQCT. Further, image registration can be used to longitudinally assess bone remodeling. These methods could be implemented in future studies to provide important pathophysiological information on the progression of bone damage

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 2 = 0.85) and XtremeCTII (r 2 = 0.87). Erosion volume assessment by BAM-Erosion agreed strongly (r 2 = 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

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 (r2 = 0.85) and XtremeCTII (r2 = 0.87). Erosion volume assessment by BAM-Erosion agreed strongly (r2 = 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