Publications by Year: 2021

2021

Nitsch J, Sack J, Halle MW, Moltz JH, Wall A, Rutherford AE, Kikinis R, Meine H. MRI-Based Radiomic Feature Analysis of End-Stage Liver Disease for Severity Stratification. Int J Comput Assist Radiol Surg. 2021;16(3):457–66.
PURPOSE: We aimed to develop a predictive model of disease severity for cirrhosis using MRI-derived radiomic features of the liver and spleen and compared it to the existing disease severity metrics of MELD score and clinical decompensation. The MELD score is compiled solely by blood parameters, and so far, it was not investigated if extracted image-based features have the potential to reflect severity to potentially complement the calculated score. METHODS: This was a retrospective study of eligible patients with cirrhosis ([Formula: see text]) who underwent a contrast-enhanced MR screening protocol for hepatocellular carcinoma (HCC) screening at a tertiary academic center from 2015 to 2018. Radiomic feature analyses were used to train four prediction models for assessing the patient’s condition at time of scan: MELD score, MELD score [Formula: see text] 9 (median score of the cohort), MELD score [Formula: see text] 15 (the inflection between the risk and benefit of transplant), and clinical decompensation. Liver and spleen segmentations were used for feature extraction, followed by cross-validated random forest classification. RESULTS: Radiomic features of the liver and spleen were most predictive of clinical decompensation (AUC 0.84), which the MELD score could predict with an AUC of 0.78. Using liver or spleen features alone had slightly lower discrimination ability (AUC of 0.82 for liver and AUC of 0.78 for spleen features only), although this was not statistically significant on our cohort. When radiomic prediction models were trained to predict continuous MELD scores, there was poor correlation. When stratifying risk by splitting our cohort at the median MELD 9 or at MELD 15, our models achieved AUCs of 0.78 or 0.66, respectively. CONCLUSIONS: We demonstrated that MRI-based radiomic features of the liver and spleen have the potential to predict the severity of liver cirrhosis, using decompensation or MELD status as imperfect surrogate measures for disease severity.
Wu C, Diao B, Ungi T, Sedghi A, Kikinis R, Mousavi P, Fichtinger G. Development of an Open-Source System for Prostate Biopsy Training in Senegal. SPIE Medical Imaging. 2021;.
PURPOSE: Prostate cancer is the second most common cancer diagnosed in men. The rate is disproportionately high among men in sub-Saharan Africa where, unlike in North America and Western Europe, the screening process for prostate cancer has historically not been routine. Currently, as awareness regarding prostate health increases, more patients in this region are being referred to trans-rectal ultrasound guided prostate biopsy, a diagnosis procedure which requires a strong understanding of prostate zonal anatomy. To aid in the instruction of this procedure, prostate biopsy training programs need to be implemented. Unfortunately, current TRUS-guided training tools are not ideal for reproducibility in these Western African countries. To answer this challenge, we are developing an affordable and open-source training simulator for TRUS-guided prostate biopsy, for use in Senegal.  In this paper, we present the implementation of the training simulator’s virtual interface, highlighting the generation and evaluation of the critical training component of zonal anatomy overlaid on TRUS. METHODS: For the simulator’s dataset, we registered TRUS and MRI volumes together to obtain the zonal segmentation from the MRI volumes. After generating ten pairings of TRUS overlaid with zonal segmentation, we designed and implemented a virtual TRUS training system, developed in open-source software. The objective of our simulator is to teach trainees to accurately identify the prostate’s anatomical zones in TRUS. To confirm the system’s usability for training zonal identification, we conducted a two-part survey on the quality of the zonal overlays with 7 urology experts. In the first part, they assessed the zonal overlay for visual correctness by rating 10 images from one patient’s TRUS with registered overlay on a 5-point Likert scale. For the second part, they labelled 10 plain TRUS volumes with zonal anatomy and the labels were compared to the labels of our overlay. RESULTS: On average, experts rated the zonal overlay’s visual accuracy at 4 out of 5. Furthermore, 7 out of 7 experts labelled the peripheral, anterior, and transitional zones in the same regions we overlaid them, and 5 out of 7 labelled the central zone in the same region we overlaid it. CONCLUSION: We created the prototype of a TRUS imaging simulator in open-source software. A vital training component, zonal overlay, was generated using publicly accessible data and validated by expert urologists for prostate zone identification, confirming the concept.
Steinmann S, Lyall AE, Langhein M, Nägele FL, Rauh J, Cetin-Karayumak S, Zhang F, Mussmann M, Billah T, Makris N, Pasternak O, Donnell LJO, Rathi Y, Kubicki M, Leicht G, Shenton ME, Mulert C. Sex-Related Differences in White Matter Asymmetry and Its Implications for Verbal Working Memory in Psychosis High-Risk State. Front Psychiatry. 2021;12:686967.
Objective: Sexual dimorphism has been investigated in schizophrenia, although sex-specific differences among individuals who are at clinical high-risk (CHR) for developing psychosis have been inconclusive. This study aims to characterize sexual dimorphism of language areas in the brain by investigating the asymmetry of four white matter tracts relevant to verbal working memory in CHR patients compared to healthy controls (HC). HC typically show a leftward asymmetry of these tracts. Moreover, structural abnormalities in asymmetry and verbal working memory dysfunctions have been associated with neurodevelopmental abnormalities and are considered core features of schizophrenia. Methods: Twenty-nine subjects with CHR (17 female/12 male) for developing psychosis and twenty-one HC (11 female/10 male) matched for age, sex, and education were included in the study. Two-tensor unscented Kalman filter tractography, followed by an automated, atlas-guided fiber clustering approach, were used to identify four fiber tracts related to verbal working memory: the superior longitudinal fasciculi (SLF) I, II and III, and the superior occipitofrontal fasciculus (SOFF). Using fractional anisotropy (FA) of tissue as the primary measure, we calculated the laterality index for each tract. Results: There was a significantly greater right>left asymmetry of the SLF-III in CHR females compared to HC females, but no hemispheric difference between CHR vs. HC males. Moreover, the laterality index of SLF-III for CHR females correlated negatively with Backward Digit Span performance, suggesting a greater rightward asymmetry was associated with poorer working memory functioning. Conclusion: This study suggests increased rightward asymmetry of the SLF-III in CHR females. This finding of sexual dimorphism in white matter asymmetry in a language-related area of the brain in CHR highlights the need for a deeper understanding of the role of sex in the high-risk state. Future work investigating early sex-specific pathophysiological mechanisms, may lead to the development of novel personalized treatment strategies aimed at preventing transition to a more chronic and difficult-to-treat disorder.
Fichtinger G, Mousavi P, Ungi T, Fenster A, Abolmaesumi P, Kronreif G, Ruiz-Alzola J, Ndoye A, Diao B, Kikinis R. Design of an Ultrasound-Navigated Prostate Cancer Biopsy System for Nationwide Implementation in Senegal. J Imaging. 2021;7(8):154.
This paper presents the design of NaviPBx, an ultrasound-navigated prostate cancer biopsy system. NaviPBx is designed to support an affordable and sustainable national healthcare program in Senegal. It uses spatiotemporal navigation and multiparametric transrectal ultrasound to guide biopsies. NaviPBx integrates concepts and methods that have been independently validated previously in clinical feasibility studies and deploys them together in a practical prostate cancer biopsy system. NaviPBx is based entirely on free open-source software and will be shared as a free open-source program with no restriction on its use. NaviPBx is set to be deployed and sustained nationwide through the Senegalese Military Health Service. This paper reports on the results of the design process of NaviPBx. Our approach concentrates on "frugal technology", intended to be affordable for low-middle income (LMIC) countries. Our project promises the wide-scale application of prostate biopsy and will foster time-efficient development and programmatic implementation of ultrasound-guided diagnostic and therapeutic interventions in Senegal and beyond.
Ye Y, Barapatre S, Davis MK, Elliston KO, Davatzikos C, Fedorov A, Fillion-Robin JC, Foster I, Gilbertson JR, Lasso A, Miller J V, Morgan M, Pieper S, Raumann BE, Sarachan BD, Savova G, Silverstein JC, Taylor DP, Zelnis JB, Zhang GQ, Cuticchia J, Becich MJ. Open-source Software Sustainability Models: Initial White Paper From the Informatics Technology for Cancer Research Sustainability and Industry Partnership Working Group. J Med Internet Res. 2021;23(12):e20028.
BACKGROUND: The National Cancer Institute Informatics Technology for Cancer Research (ITCR) program provides a series of funding mechanisms to create an ecosystem of open-source software (OSS) that serves the needs of cancer research. As the ITCR ecosystem substantially grows, it faces the challenge of the long-term sustainability of the software being developed by ITCR grantees. To address this challenge, the ITCR sustainability and industry partnership working group (SIP-WG) was convened in 2019. OBJECTIVE: The charter of the SIP-WG is to investigate options to enhance the long-term sustainability of the OSS being developed by ITCR, in part by developing a collection of business model archetypes that can serve as sustainability plans for ITCR OSS development initiatives. The working group assembled models from the ITCR program, from other studies, and from the engagement of its extensive network of relationships with other organizations (eg, Chan Zuckerberg Initiative, Open Source Initiative, and Software Sustainability Institute) in support of this objective. METHODS: This paper reviews the existing sustainability models and describes 10 OSS use cases disseminated by the SIP-WG and others, including 3D Slicer, Bioconductor, Cytoscape, Globus, i2b2 (Informatics for Integrating Biology and the Bedside) and tranSMART, Insight Toolkit, Linux, Observational Health Data Sciences and Informatics tools, R, and REDCap (Research Electronic Data Capture), in 10 sustainability aspects: governance, documentation, code quality, support, ecosystem collaboration, security, legal, finance, marketing, and dependency hygiene. RESULTS: Information available to the public reveals that all 10 OSS have effective governance, comprehensive documentation, high code quality, reliable dependency hygiene, strong user and developer support, and active marketing. These OSS include a variety of licensing models (eg, general public license version 2, general public license version 3, Berkeley Software Distribution, and Apache 3) and financial models (eg, federal research funding, industry and membership support, and commercial support). However, detailed information on ecosystem collaboration and security is not publicly provided by most OSS. CONCLUSIONS: We recommend 6 essential attributes for research software: alignment with unmet scientific needs, a dedicated development team, a vibrant user community, a feasible licensing model, a sustainable financial model, and effective product management. We also stress important actions to be considered in future ITCR activities that involve the discussion of the sustainability and licensing models for ITCR OSS, the establishment of a central library, the allocation of consulting resources to code quality control, ecosystem collaboration, security, and dependency hygiene.
Schilling KG, Rheault F, Petit L, Hansen CB, Nath V, Yeh FC, Girard G, Barakovic M, Rafael-Patino J, Yu T, Fischi-Gomez E, Pizzolato M, Ocampo-Pineda M, Schiavi S, Canales-Rodríguez EJ, Daducci A, Granziera C, Innocenti G, Thiran JP, Mancini L, Wastling S, Cocozza S, Petracca M, Pontillo G, Mancini M, Vos SB, Vakharia VN, Duncan JS, Melero H, Manzanedo L, Sanz-Morales E, Peña-Melián Ángel, Calamante F, Attyé A, Cabeen RP, Korobova L, Toga AW, Vijayakumari AA, Parker D, Verma R, Radwan A, Sunaert S, Emsell L, De Luca A, Leemans A, Bajada CJ, Haroon H, Azadbakht H, Chamberland M, Genc S, Tax CMW, Yeh PH, Srikanchana R, Mcknight CD, Yang JYM, Chen J, Kelly CE, Yeh CH, Cochereau J, Maller JJ, Welton T, Almairac F, Seunarine KK, Clark CA, Zhang F, Makris N, Golby A, Rathi Y, Donnell LJO, Xia Y, Aydogan DB, Shi Y, Fernandes FG, Raemaekers M, Warrington S, Michielse S, Ramírez-Manzanares A, Concha L, Aranda R, Meraz MR, Lerma-Usabiaga G, Roitman L, Fekonja LS, Calarco N, Joseph M, Nakua H, Voineskos AN, Karan P, Grenier G, Legarreta JH, Adluru N, Nair VA, Prabhakaran V, Alexander AL, Kamagata K, Saito Y, Uchida W, Andica C, Abe M, Bayrak RG, Wheeler-Kingshott CAMG, Angelo ED, Palesi F, Savini G, Rolandi N, Guevara P, Houenou J, López-López N, Mangin JF, Poupon C, Román C, Vázquez A, Maffei C, Arantes M, Andrade JP, Silva SM, Calhoun VD, Caverzasi E, Sacco S, Lauricella M, Pestilli F, Bullock D, Zhan Y, Brignoni-Perez E, Lebel C, Reynolds JE, Nestrasil I, Labounek R, Lenglet C, Paulson A, Aulicka S, Heilbronner SR, Heuer K, Chandio BQ, Guaje J, Tang W, Garyfallidis E, Raja R, Anderson AW, Landman BA, Descoteaux M. Tractography Dissection Variability: What Happens When 42 Groups Dissect 14 White Matter Bundles on the Same Dataset?. Neuroimage. 2021;243:118502.
White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.
Novello L, Agarwal N, Vennarini S, Lorentini S, Zacà D, Mussano A, Pasternak O, Jovicich J. Longitudinal Changes in Brain Diffusion MRI Indices during and after Proton Beam Therapy in a Child with Pilocytic Astrocytoma: A Case Report. Diagnostics (Basel). 2021;12(1):26.
Proton beam therapy (PBT) is an effective pediatric brain tumor treatment. However, the resulting microstructural changes within and around irradiated tumors are unknown. We retrospectively applied diffusion tensor imaging (DTI) and free-water imaging (FWI) on diffusion-weighted magnetic resonance imaging (dMRI) data to monitor microstructural changes during the PBT and after 8 months in a pilocytic astrocytoma (PA) and normal-appearing white matter (NAWM). We evaluated the conventional MRI- and dMRI-derived indices from six MRI sessions (t0-t5) in a Caucasian child with a hypothalamic PA: at baseline (t0), during the PBT (t1-t4) and after 8 months (t5). The tumor voxels were classified as "solid" or "fluid" based on the FWI. While the tumor volume remained stable during the PBT, the dMRI analyses identified two different response patterns: (i) an increase in fluid content and diffusivity with anisotropy reductions in the solid voxels at t1, followed by (ii) smaller variations in fluid content but higher anisotropy in the solid voxels at t2-t4. At follow-up (t5), the tumor volume, fluid content, and diffusivity in the solid voxels increased. The NAWM showed dose-dependent microstructural changes. The use of the dMRI and FWI showed complex dynamic microstructural changes in the irradiated mass during the PBT and at follow-up, opening new avenues in our understanding of radiation-induced pathophysiologic mechanisms in tumors and the surrounding tissues.
Herberthson M, Boito D, Haije TD, Feragen A, Westin CF, Özarslan E. Q-Space Trajectory Imaging With Positivity Constraints (QTI+). Neuroimage. 2021;238:118198.
Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method’s robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique’s performance on shorter acquisitions could make it feasible in routine clinical practice.
Hsu TMH, Schawkat K, Berkowitz SJ, Wei JL, Makoyeva A, Legare K, DeCicco C, Paez N, Wu JSH, Szolovits P, Kikinis R, Moser AJ, Goehler A. Artificial Intelligence to Assess Body Composition on Routine Abdominal CT Scans and Predict Mortality in Pancreatic Cancer- a Recipe for Your Local Application. Eur J Radiol. 2021;142:109834.
BACKGROUND: Body composition is associated with mortality; however its routine assessment is too time-consuming. PURPOSE: To demonstrate the value of artificial intelligence (AI) to extract body composition measures from routine studies, we aimed to develop a fully automated AI approach to measure fat and muscles masses, to validate its clinical discriminatory value, and to provide the code, training data and workflow solutions to facilitate its integration into local practice. METHODS: We developed a neural network that quantified the tissue components at the L3 vertebral body level using data from the Liver Tumor Challenge (LiTS) and a pancreatic cancer cohort. We classified sarcopenia using accepted skeletal muscle index cut-offs and visceral fat based its median value. We used Kaplan Meier curves and Cox regression analysis to assess the association between these measures and mortality. RESULTS: Applying the algorithm trained on LiTS data to the local cohort yielded good agreement [>0.8 intraclass correlation (ICC)]; when trained on both datasets, it had excellent agreement (>0.9 ICC). The pancreatic cancer cohort had 136 patients (mean age: 67 ± 11 years; 54% women); 15% had sarcopenia; mean visceral fat was 142 cm2. Concurrent with prior research, we found a significant association between sarcopenia and mortality [mean survival of 15 ± 12 vs. 22 ± 12 (p < 0.05), adjusted HR of 1.58 (95% CI: 1.03-3.33)] but no association between visceral fat and mortality. The detector analysis took 1 ± 0.5 s. CONCLUSIONS: AI body composition analysis can provide meaningful imaging biomarkers from routine exams demonstrating AI’s ability to further enhance the clinical value of radiology reports.
Nilsson M, Eklund G, Szczepankiewicz F, Skorpil M, Bryskhe K, Westin CF, Lindh C, Blomqvist L, Jäderling F. Mapping Prostatic Microscopic Anisotropy Using Linear and Spherical B-Tensor Encoding: A Preliminary Study. Magn Reson Med. 2021;86(4):2025–33.
PURPOSE: Tensor-valued diffusion encoding provides more specific information than conventional diffusion-weighted imaging (DWI), but has mainly been applied in neuroimaging studies. This study aimed to assess its potential for the imaging of prostate cancer (PCa). METHODS: Seventeen patients with histologically proven PCa were enrolled. DWI of the prostate was performed with linear and spherical tensor encoding using a maximal b-value of 1.5 ms/µm2 and a voxel size of 3 × 3 × 4 mm3 . The gamma-distribution model was used to estimate the mean diffusivity (MD), the isotropic kurtosis (MKI ), and the anisotropic kurtosis (MKA ). Regions of interest were placed in MR-defined cancerous tissues, as well as in apparently healthy tissues in the peripheral and transitional zones (PZs and TZs). RESULTS: DWI with linear and spherical encoding yielded different image contrasts at high b-values, which enabled the estimation of MKA and MKI . Compared with healthy tissue (PZs and TZs combined) the cancers displayed a significantly lower MD (P < .05), higher MKI (P < 10-5 ), and lower MKA (P < .05). Compared with the TZ, tissue in the PZ showed lower MD (P < 10-3 ) and higher MKA (P < 10-3 ). No significant differences were found between cancers of different Gleason scores, possibly because of the limited sample size. CONCLUSION: Tensor-valued diffusion encoding enabled mapping of MKA and MKI in the prostate. The elevated MKI in PCa compared with normal tissues suggests an elevated heterogeneity in the cancers. Increased in-plane resolution could improve tumor delineation in future studies.